“…For the present purposes, a simple and approximate working model, guided by the results from the above references, has been alternatively implemented and used: (i) Equation is accepted for S 4 ≤ 0.6; (ii) Equation is used to determine δNi 0.6 and δNi 1.0 , at which S 4 = 0.6 and S 4 = 1.0, respectively; (iii) it is assumed that S 4 = 1.0 at [ δNi 1.0 + ( δNi 1.0 − δNi 0.6 )] = ( 2δNi 1.0 − δNi 0.6 ); (iv) the derivatives d S 4 / d (δNi) at δNi 0.6 and ( 2δNi 1.0 − δNi 0.6 ) are respectively equal to 0.6/ δNi 0.6 and 0; (v) a third‐degree polynomial is adjusted to the two points [ δNi 0.6 , S 4 = 0.6] and [(2 δNi 1.0 − δNi 0.6 ), S 4 = 1.0], considering the respective derivatives; and (vi) S 4 = 1.0 for δNi ≥ (2 δNi 1.0 − δNi 0.6 ). As an initial test, the approximate working model has been applied to the input data adopted by Bhattacharyya et al (2017) and Bhattacharyya et al (2019) (namely, L = 50 km, L o = 10 km , ν = 2.0, z S = 35,786 km , z = 350 km, and f = 251 MHz) and the corresponding results displayed a reasonable agreement, for small zenith angles. The next section will display the curves S 4 (δNi) for the parameters, frequencies, and zenith angles used in the present contribution, which exhibit smooth transitions from the weak to the strong scattering and saturation regimes similar to those observed in Figures 7a, 7c, and 7e of the paper by Carrano and Rino (2016), for spectral indices less than 2.5.…”