This study revisits previous formulations of atmospheric refractivity at L-band frequencies, focusing on signals from Global Navigation Satellite Systems (GNSS). A refined model expression is proposed as a function of air density, temperature, and composition, evaluated using a comprehensive set of existing laboratory and atmospheric measurements. The key measurements that most affect the final accuracy are identified, establishing traceable error bounds and indicating where further experimental work could confirm or improve the model.
Recent studies on the use of large volumes of GNSS radio occultation (GNSSRO) observations in Numerical Weather Prediction (NWP) show that the precise formulation of refractivity becomes increasingly critical as data volumes grow. Although the revision is modest, its impact lies within the range where NWP sensitivity becomes non-negligible.
Compared to earlier work, this study (1) incorporates updated fundamental measurements, (2) accounts for the small but measurable variability in atmospheric composition, mainly increasing CO2 and decreasing O2, emphasizing that refractivity traceability is composition-dependent, and (3) extends the model to include hydrometeors. A simplified formulation based on hydrometeor oblateness is proposed, suitable for NWP applications where only limited hydrometeor information is available. Nonspherical hydrometeors tend to align during fall, introducing weak birefringence that can be detected during GNSS occultations with dual-polarization receivers.
The resulting refractivity expression is presented as a function of air density, temperature, moisture, and composition, and (using a simplified model of atmospheric evolution) also as a function of density, temperature, moisture, and time. It is concluded that the result meets the target accuracy of 0.01%, with all remaining uncertainty sources below this threshold. With respect to earlier work, composition and hydrometeor additions do not modify the largest part of the bulk refractivity above this threshold.
Global Navigation Satellite System (GNSS) Radio Occultation (GNSSRO) is a technique for remote sensing of the Earth's atmosphere . It has proven highly useful for operational numerical weather prediction, (e.g. ), and climate research, (e.g. ), providing high-resolution vertical profiles of temperature and humidity.
A major reason for this success is the excellent absolute calibration of GNSSRO. The data are derived from highly accurate and traceable frequency observations, a procedure that is intrinsically robust against instrumental bias. For thermodynamic measurements of the atmosphere, GNSSRO has established a qualitatively new traceability link to the International System of Units (SI), complementing in-situ instruments such as surface thermometers, radiosondes, and aircraft sensors.
Although no measurement is entirely free from error or bias, radio occultation provides observations of exceptionally high raw precision. How much of this precision can be translated into accuracy in the retrieved atmospheric variables remains an open question. The issue is linked, on one hand, to the interpretation of raw measurements in terms of specific propagation paths (e.g. ), and on the other, to the connection between the radio-optical and thermodynamic properties of air, expressed through the refractive index as a constitutive relation). This latter connection is the focus of this study.
A previous investigation , hereafter referred to as AL11, concluded that for this constitutive relation, an average accuracy of 0.1% can be routinely achieved in linking the radio-optical and thermodynamic properties of air (i.e. refractivity, pressure, temperature and moisture). It also concluded that current understanding of the relevant physics should allow an accuracy of 0.01%. Accurate laboratory measurements of air's constitutive properties already exist, but to exploit these and achieve accuracy better than 0.1%, proper attention must be given to the precise definitions of thermodynamic quantities, to the equation of state of moist air and its deviations from ideal-gas behaviour, to the various standard formulations of atmospheric humidity (specific and relative humidity, vapour pressure, dew-point depression, etc.), and to the strongly non-ideal behaviour of water vapour. Such considerations are essential not only when formulating constitutive relations and their parameters but also, and most importantly, during their routine use in NWP.
For these reasons, the simplest approach for a user may be to accept that an accuracy better than 0.1% is unlikely to be achieved. In this case, the thermodynamic details of the air equation of state or the representation of water vapour concentration are not especially demanding. However, achieving an additional order of magnitude in accuracy is realistic and can be valuable for many scientific and operational applications. Such a goal has been expressed, for example, by , to improve the traceability of radio occultation data for both weather forecasting and climate research. The motivation is twofold: first, to ensure that the use of GNSSRO data is not limited by uncertainties in the constitutive relations, and second, to provide SI-traceable estimates of the accuracy of these relations.
It is therefore assumed hereafter that exceeding the 0.1% accuracy level is of interest. To ensure traceability, it is not only desirable that the resulting expressions are as accurate as possible, but also that their uncertainties be quantifiable and attributable to specific sources. On the other hand, achieving accuracy better than 0.01% does not appear realistic at present. Accordingly, the analysis presented here focuses on this intermediate range of accuracy.
Recently, in connection with the Radio Occultation Meteorology Experiment (ROMEX, ), several studies have addressed the expansion of radio-optical GNSS measurements to volumes significantly higher than those available over the past two decades. The number of global GNSSRO profiles has increased from about 10 000 profiles per day, typical in major NWP centers since the launch of the COSMIC-2 mission , to between 20 000 and 40 000 global profiles per day. The impact of GNSSRO data, already substantial, generally increases with the higher volume, but also becomes more sensitive to details of data processing and assimilation.
NWP systems are highly sensitive to biases between expected and observed data. While instrumental biases may sometimes need correction, it is common practice to apply bias correction simply to ensure statistical consistency between forecasts and observations (see ). GNSSRO measurements, however, are known to exhibit very low intrinsic bias. Therefore, they should not be bias-corrected solely for statistical alignment; rather, it is preferable to ensure that the traceability link between the radio-optical and thermodynamic properties of air remains as bias-free as possible.
In the present context, “sufficiently free of bias” is taken to mean “accurate enough that NWP systems do not perceive a bias originating from GNSSRO data”, even when used at the highest practical data volumes, currently about 40 000 profiles per day. Both the earlier AL11 formulation and the expression presented here have been tested with these enhanced volume experiments in the global NWP system at Environment and Climate Change Canada (ECCC), as part of the ROMEX studies. In both cases, the results were found to be sufficiently free of bias in the sense defined above, albeit with some differences discussed below. Since the GNSSRO perceived bias against a model is a complex combination of many elements, and the refractivity model is only one of them, a summary of the prime among those in these tests is given in the discussion section. Within ECCC, refractivity, either AL11 or this work, coupled with a suitable equation of state, is not currently viewed as a source of concern towards GNSSRO bias, while it had been in the past .
An ansatz expression is developed here as an extension of AL11, together with an associated error budget. AL11 has already proven valuable in operational forecasting environments, improving the consistency among different SI-traceability links of atmospheric measurements, notably against radiosondes and aircraft data , and leading to enhanced forecast performance. Both AL11 and the updated ansatz expression proposed here are similar to existing ones in the literature but differ in recommending the use of density as an independent variable rather than pressure. This choice allows a refractivity expression independent of the equation of state, since air shows measurable deviations from ideal-gas behaviour, already noticeable in dry air and increasingly relevant when water vapour is present.
The motivation for this expansion in scope beyond AL11 is twofold. First, some properties of air are imperfectly known; second, certain properties are known not to be static. Both must be included in the error budget. Regarding imperfect knowledge, the dielectric response of air's constituents (mainly N2, O2, etc.) to an electromagnetic wave is known with finite accuracy. Regarding variability, even if the molecular properties were perfectly known, the composition of air is not fixed in space and time. The most obvious source of variability is the moisture content, already explicit in standard refractivity formulations. However, the composition of dry air is also variable, albeit to a lesser extent. Records of CO2 and O2 concentration indicate that dry air is not an immutable mixture. Several modes of variability can be identified in these records: global drift, seasonal, geographical, and vertical. Most are small in amplitude for the purpose of this study, and their effect on refractivity is below the target accuracy defined above. Of all the variability modes in the composition of dry air, this work focuses on the largest one, the long-term global drift associated with the measurable accumulation of CO2 in the atmosphere and the associated decrease in O2.
The main microphysical model describing the molecular response of gas-phase air to an electromagnetic wave remains unchanged from . However, the sources of information have been revised to provide error bounds traceable to specific laboratory measurements. Some microphysical elements have been incorporated, notably allowing air to be treated as a multiphase heterogeneous mixture, described as an effective medium: a gaseous matrix that may contain small inclusions of liquid (drops), or solid water (ice, snow). As long as these heterogeneities are small, the medium is effectively homogeneous at the scale of the test electromagnetic field. With GNSS wavelengths in the range λGNSS=0.18–0.26 m, this approximation applies to rain drops, ice particles, and snowflakes in the atmosphere, all small with respect to λGNSS.
Compared to AL11, the accuracy of the various input parameters has been reassessed, and the most critical ones identified, to determine which limit the overall accuracy of the final model and to apportion the net uncertainty among contributing sources. An updated refractivity expression is then presented. Although the differences relative to the earlier AL11 expression are small, they may be significant for certain applications. Most of the dry-air difference arises from a slightly revised composition of dry air. The new expression and its accompanying error analysis enable several additional applications.
First, specific laboratory measurements that are critical to the model's accuracy are identified, linking the reliability of the refractivity expression directly to the accuracy of those experimental data. Second, the dependence of refractivity on atmospheric composition is made explicit, which may be relevant to decadal-scale climate studies. Since composition is slightly time-dependent, the thermodynamic relationship between refractivity, pressure, moisture, and temperature effectively becomes time-dependent as well. Third, hydrometeors are included. Expressions for their contribution have been presented previously (e.g. ), and the main results obtained here are consistent with those studies. However, the dependence on hydrometeor shape is explicitly explored, as nonspherical particles tend to reorient during fall, leading to anisotropic shape distributions that interact differently with each polarization component of the signal. Whereas specific microphysical studies of refractivity and polarization effects for individual ice-crystal habits are ongoing (, and references therein), this work provides an estimate of refractivity dependence on general hydrometeor shape properties. Because NWP systems typically have limited information about hydrometeors, this analysis focuses on two key parameters: the density of water in hydrometeor form (as a fraction of air) and the hydrometeor oblateness as perceived by the electromagnetic signal. Including this anisotropy enables estimates of polarimetric effects. Users should evaluate whether these two parameters suffice for their purposes or whether more detailed scattering models for specific ice or snow shapes are required.
The various aspects of the microphysical model, relating molecular properties, atmospheric composition, hydrometeor presence, and refractivity are developed in Sect. . This model is then applied under a range of conditions to produce a fit suitable for practical use, in Sect. . This fit contain several adjusted parameters whose precise value depend on input measurements. An uncertainty analysis of the fit attributes the main sources of parameter uncertainty and identifies which measurements are most critical, indicating where updated experimental data could improve model performance. This analysis is discussed in Sect. , along with the differences relative to the earlier AL11 model and to other existing expressions. Finally, in Sect. , user expressions are recommended.
Microphysical model of the air
The expression of atmospheric refractivity combines two largely independent aspects: (1) the relationship between the electromagnetic response and the number of molecules in an air parcel, and (2) the relationship between the number of molecules and the parcel's mass or density. The first depends on molecular polarizability and dipole, whereas the second is governed by molecular mass.
A microphysical model linking molecular properties (polarizability, electric and magnetic dipoles, and molecular mass) to macroscopic refractivity was described in AL11. The present study follows a similar approach, updating the selection of experimental data. AL11 described a single-phase gaseous medium, which with minor modifications is retained here as the gas phase of atmospheric air. In the present work, the model is extended to a potentially multiphase mixture, allowing for small liquid or solid inclusions (droplets or snow, not necessarily spherical) may be present. Inclusions (solid or liquid) are assumed to occupy small fractions of the total volume and mass, and to be much smaller than the electromagnetic wavelength. Under these conditions the mixture can be treated as a homogeneous effective medium .
The microphysical model is here reviewed briefly. The refractive index of the effective medium will follow Maxwell's equation (, Sect. 1.2):
n=ϵrμr
where ϵr and μr are the relative electric permittivity and magnetic permeability. Although ϵr and μr are generally frequency-dependent, at GNSS frequencies (1.1–1.6 GHz), which lie well below the lowest vibrational and rotational resonances of air and water vapour , they can be accurately described by their static limits. Most of this response is electric, and only a small fraction is magnetic.
Air, of relative permittivity and permeability (ϵr,μr), is here considered as a mixed phase medium containing:
A matrix gas phase, with (ϵrg,μrg), containing dry air and water vapour.
A liquid phase (water drops), with ϵrl, while assuming μrl=1.
A solid phase (ice or snow), with ϵri, while assuming μri=1.
The magnetic response of water and ice hydrometeors is thus neglected. The single-phase relative permittivities (ϵrg, ϵrl, ϵri) will be evaluated, and from them the effective relative permittivity of the composite, ϵr. Similarly, for the permeabilities, but since the magnetic response of water and ice is neglected, the permeability of the composite will be μr=μrg.
Dry air is assumed to be well mixed at the spatial scales relevant to GNSSRO. However, the model allows for spatial and temporal variability at much larger scales (Earth's hemispheres, seasons, decades). This variability is primarily driven by the increasing concentration of CO2, and the generally anticorrelated concentration of O2.
A list of values for the constituents' concentrations was presented in AL11, which were assumed to be static. The four leading substances in dry air (N2, O2, Ar and CO2) are revisited here, and based on additional observations. Firstly, their concentrations are not considered here constant parameters. The dielectric and magnetic properties of these substances have been slightly updated from AL11, and we follow the expressions offered by some laboratory measurements that have been here selected as particularly accurate. The use of the provided expressions allows better linking between refractivity and measurement, for error traceability. All trace constituents (those found in air at concentration below that of CO2, such as Ne, He, etc.) are retained, as in AL11, unmodified. Since concentrations are not assumed to be fixed, the description of the state of a parcel of air will therefore include density, temperature, and moisture, as usual, but also some composition details (concentration of O2 and CO2).
Each chemical substance is assumed to have a fixed isotopic composition, and its properties assumed to be those of the average isotopic mixture found in nature. The isotopic composition affects nearly only molar mass, as molecular polarizability or magnetizability show very small isotopic dependence for dry air substances. Only water is split into a fixed isotopic mix, since besides mass, molecular electric dipole can be affected, which is one of the major sources of air refractivity. In the case of water, the isotopic composition is kept constant across the vapour, liquid, and solid phases.
The gas phase
In the gas phase, intermolecular interactions are sufficiently weak that the standard Lorentz-Lorenz formulation applies :
ϵrg=1+4π2+ϵr3∑iniαi+μi23kBTf(ϵr)
where ϵrg is the relative permittivity of the gas phase, and the sum is over all pure constituents i of the mixture. This expression relates the macroscopic dielectric response to microscopic properties of each species, such as molecular polarizability αi, permanent electric dipole moment μi, and particle number density ni.
The structure of Eq. () shows that the refractivity associated with molecular polarizability scales with ni and therefore with the constituent mass density ρi. The refractivity arising from the permanent dipoles μi scales as ρi/T. We emphasize that all microscopic contributions depend on the effective medium in which the molecules are embedded, through its effective permittivity ϵr. The factors 2+ϵr3 and f(ϵr) represent these medium interactions. Because the environment includes not only the gas phase but also any coexisting condensed phases, Eq. () is implicit in ϵr. The weak interactions in air allow this dependence to be resolved iteratively, starting from ϵr=1, with very fast convergence.
The function f(ϵr) represents the interaction between a molecular dipole and its polarizable surroundings and can be written :
f(ϵr)=g9ϵr(2ϵr+1)(ϵr+2)
where the correlation factor g quantifies the difference between the local environment of a polar molecule, and that of an average air parcel. When molecular interactions are weak, both f and g are close to unity. Since water vapour is the only major atmospheric constituent with a significant permanent electric dipole moment, the dipole term in Eq. () applies almost entirely to H2O. Thus, the functions f and g of interest are essentially water-specific. An explicit expression for g is given in , which also enables evaluation of f. Although g can become large in condensed water phases and is easy measured at high molecular densities , here it is required only at the very low densities of atmospheric water vapour, where g-1 is small but not negligible. Existing parameterizations interpolate between the zero-density limit (g=1) and condensed-phase values, leaving uncertainty in g at the low densities relevant for atmospheric moist refractivity.
The magnetic response is always very weak and follows a Lorentz-Lorenz relation analogous to the electric case. The molecular dipole contribution arises almost entirely from O2. Because magnetic dipole interactions are much weaker than those of water electric dipoles, the associated correlation factors are negligible: f(μr)≈1, and we take (2+μr)/3≈1. The magnetic permeability of the gas phase then reduces to
μrg=1+4π∑iniχi+νi23kBT
where χi and νi are the molecular magnetizability and permanent magnetic dipole moment of species i. Unlike the electric case, Eq. () does not require iterative evaluation. Furthermore, water – vapour, liquid, or solid – has a very weak magnetic response, so the relative permeability of air may be taken to be that of its gas phase alone: μr=μrg.
Concentration of dry substances
Excluding water in any phase, and assuming that other particulates such as dust contribute negligibly, the remaining atmospheric constituents are referred to as dry air. A representative list of dry-air components and their approximate properties was provided in AL11. Here we examine in more detail the four dominant species (N2, O2, Ar and CO2). For all other constituents, we assume that their concentrations and electromagnetic properties are known with sufficient accuracy for the present purpose.
The concentration of CO2 varies , most prominently as a long-term drift, and also following seasonal and latitudinal patterns. To accommodate this variability, the present model explicitly treats the molar fraction of carbon dioxide, xCO2, as a free variable. Measurements and reconstructions over the industrial period indicate 275<xCO2/10-6<430, with the upper end representative of recent values.
It is also well established that the O2 concentration exhibits a small but measurable variability. Its value is approximately xO2≈0.2095, with a decreasing trend that is, to a large extent, negatively correlated with xCO2. For this reason, refractivity calculations in this study do not assume fixed values of xO2 or xCO2. Instead, both are treated as independent variables, analogous to density, temperature, and water vapour content.
The molar fraction of Ar, and its uncertainty, has been determined (see , and references therein) as xAr=0.009332(3). It will be assumed in this study that this value and error estimate are accurate, and that this concentration is stable. Similarly, the molar fractional concentration of all traces (Ne and other minor species) are taken as fixed at their AL11 values. Their abundances are low enough that any natural variability would have a negligible effect on L-band refractivity.
Finally, because N2 behave as a nearly passive background gas, its molar fraction is obtained by normalization within the dry-air mixture:
xN2=1-∑i≠N2xi
Electromagnetic properties of dry substances
Accurate measurements of the electromagnetic response of the four dominant dry-air constituents, N2, O2, Ar and CO2 have been obtained using high-sensitivity cavity resonators . These measurements do not directly provide the microscopic parameters required in Eqs. () and (). Instead, they yield closely related macroscopic quantities, namely the virial coefficients describing the dielectric or paramagnetic response of the pure substances.
The polarizability α0 of a molecule is related to the dielectric virial Bϵ of a pure substance through:
α0=14π3×10-24⋅NABϵ
where NA is Avogadro's number. The measurements of show that Bϵ is not strictly constant but exhibits small dependencies on temperature and molar density. Their empirical fit is:
Bϵ(dm,T)=Aϵ1+b⋅dm+c⋅dm2+AτT273.16K-1+q⋅dm273.16KT-1
with dm the molar density. The coefficients for the relevant substances are listed in Table ; the values were provided with very high nominal accuracy.
It is noteworthy that the experimentally determined dielectric virial Bϵ depends on the thermodynamic state, despite being proportional in principle to the molecular property α0, which is often assumed constant. In AL11, the α0 were indeed treated as constant. Part of the observed variation may reflect macroscopic effects such as uncertainties in the equation state used during data reduction in . However, small physical variations of molecular polarizability with density or temperature are also plausible, due to intermolecular forces or molecular deformation, vibration or stretching.
Because measurements of are representative of highest current or foreseeable achievable accuracy for this type of experiment, we retain Eqs. () and () in the present study. This choice accepts that the molecular polarizabilities may not be strictly constant and that Table captures these subtle effects. Including the tabulated coefficients preserves direct traceable connection between laboratory measurements and the refractivity model.
Virial expression of the dielectric properties of the four leading constituents of dry air. Units of Aϵ, b, and Aτ in cm3 mol−1. Units of c and q in cm6 mol−2.
Most dry-air constituents exhibit only very weak magnetic response. To verify the order of magnitude of the contributions, simple estimates of molecular magnetizability were used . The sole non-negligible contribution is from O2, which is paramagnetic and possesses a comparatively large magnetic dipole moment νi. Following AL11, a simplified model is adopted for all species except O2, for which a more accurate expression from is used to represent the magnetic dipole νO2. The contribution of O2 to Eq. () is:
4πχO2+νO223kBT=-0.13×10-9+NA⋅μ0⋅μm23kBT
where μ0≡4π×10-7 N A−2 is the vacuum permeability, and μm=2⋅2.00386⋅μB is the molecular magnetic moment, with μB=9.274×10-24 J T−1 the Bohr magneton. The comparatively large bulk permeability arises from the substantial magnetic dipole moment νO2, rather than a large molecular susceptibility χO2.
Water fraction
Water contributes to atmospheric refractivity primarily through its permanent electric dipole μH2O, which is comparatively large among small molecules; by contrast most other air constituents are non-polar. The static polarizability αH2O, however, is similar in magnitude to that of the other major atmospheric species.
Because water is a major contributor to refractivity, and for consistency with AL11, several isotopologues are considered. These exhibit slightly different electric responses, mainly through variations in their molecular dipole moments. The overall effect on refractivity is minor, since the most abundant isotopologue (light water) comprises approximately 99.8 % of naturally occurring water. Although various processes in the hydrological cycle fractionate the isotopologues away from the global mean composition, this effect is neglected here, because it influences only a small fraction of water molecules. The isotopic fractions and microscopic properties used in this study follow AL11 and are based on and .
The average molar mass of water is determined from the assumed isotopic composition and is taken to be constant across vapour, liquid and solid phases. The adopted value is mw=18.01525.
Precipitated inclusions
Besides water vapour, atmospheric water also occurs in liquid (drops) and solid form (ice and snow). Under such conditions, condensed water contributes to the refractivity of the air matrix . This contribution may be represented by treating condensed water as small inclusions embedded in the air matrix (see ), or, more symmetrically between phases, by an effective medium formulation that mixes multiple phases . The resulting refractivity depends not only on the fraction of water present as hydrometeors but also on their shape and orientation.
Here, it is assumed that drops, ice crystals, and snowflakes may deviate from perfect sphericity, but only the simplest mode of asphericity is considered: rotational ellipsoids, with symmetry axes aligned with gravity.
Real hydrometeors are neither perfectly ellipsoidal, nor perfectly aligned as they fall. Random tilt of the principal axes is to be expected, and aerodynamic forcing (e.g. wind shear) may induce systematic tilt of the hydrometeors (known as canting). Moreover, drops, flakes and ice particles exhibit wide variability in size, shape, and orientation, with substantial spatial and temporal variability. Therefore, the present exploration of non-sphericity is not intended to be exhaustive; its purpose is to estimate the approximate effect on refractivity of the dominant main modes of asphericity. For the same reason, this model is not designed to retain high accuracy when hydrometeor concentrations are large, and the traceability link is not expected to remain robust in such cases. Although the accuracy of the nominal hydrometeor refractivity for a perfectly known distribution of hydrometeors should be better than 0.1% of the hydrometeor contribution, the large variability in the atmosphere of hydrometeor concentration, size, shape, and orientation, suggests that use of this expression should allow a much larger uncertainty deriving from this variability. A safe value for this uncertainty in the presence of hydrometeors may be as large as 50% of the estimated hydrometeor contribution.
An ellipsoidal dielectric inclusion embedded in a homogeneous matrix is polarizable and responds to an external test field (e.g. ). The response includes a depolarizing field within the inclusion, which partially screens the interior from the applied field. The magnitude of this depolarization depends strongly on the object's shape and orientation relative to the field, and conversely, for a fixed shape and orientation of the object, the response depends on the direction of the test field with respect to the object.
All phases considered here – the matrix gas, liquid water, and to a very good approximation ice – are optically isotropic, responding identically to all electromagnetic polarizations. However, the geometric orientation of the hydrometeors may introduce anisotropy if they assume preferential orientations while falling. This constitutes a case of form anisotropy (, Sect. 15.5.2). Falling drops deform into oblate shapes , with the minor axis roughly vertical and two equal major axes lying horizontally. Ice and snow particles also tend to reorient.
Because depolarization depends on the relative alignment between the test field and the oblate hydrometeor, the effective medium consisting of the air matrix and the inclusions becomes anisotropic. Two principal configurations arise: the electric vector is aligned with the major axis of the hydrometeor, or with the minor axis (see Fig. ). In the atmosphere, for horizontally propagating electromagnetic signals and vanishing canting, the electric vector aligns with the major axis for linear horizontal polarization (H, in Fig. ), and with the minor axis for linear vertical polarization (V). All other combinations of drop orientation and polarization can be expressed as mixtures of these two cases.
Scheme for a signal propagating along P, through the Earth's atmosphere. The signal finds flattened drops, of axes a=b>c in its path. The electric field oscillates within the polarization plane P. Due to the symmetry of the drops, there is always in P a direction H that lies inside the symmetry plane (a, b). Oscillations along this direction H meet the drops along their major axis. Oscillations in the perpendicular direction V meet the drops at some other angle, that depends on the elevation angle of the ray and the canting angle of the drops. In the case depicted, a horizontal ray and zero canting, V is oriented along the minor axis of the drops.
Interestingly, the dielectric response of the inclusions depends on the amount of condensed matter that they contain, as well as on their shape, but not directly on their size. Size and shape are correlated due to the fall , and also through their association with meteorological variables such as rain rate. Because the aim is to characterise the optical properties of air, the independent variables are chosen to be those directly entering the refractivity expression – namely the amount of condensed water (in drops, ice or snow) and the particle shape –, rather than meteorological quantities that correlate with them.
Shape will be expressed as the axis ratio between the vertical (V) and horizontal (H) axes, denoted al for liquid drops and ai for frozen particles. This ratio is 1 if they are spherical. Falling drops tend to become oblate, so typically al≤1.
A volume fraction ηl is occupied by liquid water, whose bulk electric permittivity many be represented by :
log10ϵl=1.94315-0.0019720⋅T-273.15K
These authors report an uncertainty of σϵl=0.050. This permittivity is very large, and produces large depolarization fields in the drops, strongly screening the molecules inside from the test field. Since the overall contribution of liquid water to air's refractivity is small, this expression is here considered to be sufficiently accurate.
Similarly, a volume fraction ηi is occupied by ice or snow, with axis ratio ai, which may differ from unity. For ice, the electric permittivity from is chosen here:
ϵi=3.1884+9.1×10-4⋅T-273.15K
where they estimate the error to be σϵi=0.0050. This permittivity is also large, though less than that of liquid water. This expression is considered adequate, since even under extreme conditions (e.g. intense hail) ice contributes only a small fraction air's refractivity.
In summary, the independent variables introduced here are the precipitated mass of liquid and solid water per unit volume of air, ρl and ρi. These correspond to the liquid water content (LWC) and ice water content (IWC), typically ranging 0≤ρl<10-2 and 0≤ρi≤4×10-3, in kg m−3. Axis ratios (vertical/horizontal) for liquid and solid ellipsoids al and ai have been considered between 0.5 and 1.25.
Multiphase model
Air is represented here as an effective medium composed of three phases: a gas matrix occupying most of the volume, and possibly liquid and ice phases present as small inclusions. Their respective volume fractions are ηg, ηl, and ηi, which are subject to ηg+ηl+ηi=1, and typically ηg≫ηl+ηi.
The phases have permittivities ϵg, ϵl, and ϵi. These are assumed to be known: ϵg from the microscopic gas-phase model described above, and ϵl and ϵi from published measurements . At scales large enough for the mixture to appear homogeneous, the ensemble of the gas matrix and embedded inclusions exhibits an effective permittivity ϵr, which we aim to determine.
Following (and references therein), the effective medium corresponds to a certain average of the different phases k. The average can be obtained by solving:
∑kηkϵk-ϵrϵk+pki⋅ϵr=0
which is approximately a volume-weighted average. Here pki=1/Pki-1 is a shape factor that depends on the depolarizing factor Pki. The depolarizing factors describe the geometric exposure of the material inside the inclusions of phase k, to a test field oriented along direction i direction, and satisfy, for each phase k: ∑i∈(x,y,z)Pki=1. Because the depolarization depends on both the shape of the inclusions and their relative orientation relative to the test field , the effective dielectric response ϵr is direction-dependent. Thus, instead of a single scalar permittivity, the model yields directional components ϵri, corresponding for example to different polarization states of the test field.
A rearrangement of Eq. (), using the dominance of the gas phase and allowing for direction-dependent effective permittivities ϵri, gives:
ϵri=ϵg+ϵg+pgi⋅ϵriηg∑k≠gηk⋅ϵk-ϵriϵk+pki⋅ϵri
As noted in Sect. (), ϵg is implicitly defined and depends weakly on the dielectric response of the surrounding medium, previously denoted ϵr, but which becomes direction-dependent ϵri in the presence of hydrometeors. The nonlinearities in Eq. (), for ϵg, and in Eq. (), for ϵri, suggest an iterative solution as the most appropriate approach. Starting with an initial guess ϵri=1 leads to rapid convergence.
The simplest inclusions are small spheres, representative for example of small droplets. Their symmetry implies Pi=1/3 for all i. Here, the possible inclusions are generalized to rotational ellipsoids whose axis of rotation is nearly vertical, possibly with some canting. Because particles fall under gravity, the two quasi-horizontal directions are expected to be equivalent. The degree of asphericity is represented by the axial ratio a, defined as the ratio of the quasi-vertical to quasi-horizontal axes. Typical atmospheric values lie between 0.5–1 . Exact closed-form expressions exist for the depolarization factors Pi of arbitrary ellipsoids (e.g. ), but given the variability of real hydrometeors in shape and orientation, a simpler but accurate approximation is sufficient.
We therefore use two depolarization factors for each inclusion type: Pv for test fields along the quasi-vertical axis (the axis of symmetry), and Ph for fields in the quasi-horizontal plane. These must satisfy Pv+2⋅Ph=1. The following approximate expressions, proposed by , satisfy the closure relations and reproduce the limiting behaviours:
Pv=11+1.6⋅a+0.4⋅a2
and Ph=0.5⋅1-Pv. For spherical hydrometeors (a=1), these reduce to Pv=Ph=1/3. For aligned oblate spheroids, as drops, where a<1, they satisfy Pv>1/3>Ph.
Axial ratios
Liquid and frozen hydrometeors (drops, ice particles, and snowflakes) have been treated as inclusions within an effective-medium, the matrix being the gas phase of air, which includes dry air and water vapour. The effective index of refraction of the mixture depends primarily on mass fraction of each phase (gas, liquid, and ice), and the dielectric properties of those phases, with an additional, weaker dependence on the shapes of the inclusions. Even if each individual phase is dielectrically isotropic, preferential hydrometeor shapes introduce geometric anisotropy into the mixture, resulting in a macroscopically anisotropic optical response.
When liquid water and snow particles are preferentially aligned around a quasi-vertical direction, two refraction indices must be evaluated, one for electric test fields along the alignment direction (“vertical”), and another for fields perpendicular to it (“horizontal”). While the alignment need not be perfectly vertical or horizontal because of canting, these terms remain convenient. These resulting refraction indices depend on the amount of precipitated water, its thermodynamic phase, and the mean axial ratios of the inclusions, but not directly on particle size. Indirectly, however, size does influence the axial ratios: larger raindrops and higher rain rates are associated with more pronounced flattening (see , and references therein).
Numerous studies of raindrop shape indicate that liquid hydrometeors can be represented as oblate spheroids with axial ratios between 0.5 and 1. Several empirical relationships have been proposed linking axial ratio, drop size, and rain rate. A simple, order of magnitude example, suggests:
al=min(1.030-0.124r0,1)
where r0 is the drop-equivalent (volume-based) radius. In this study, however, we focus solely on the optical coupling between drops and the electromagnetic field, which does not depend explicitly on drop size.
Ice particles are not flexible to adjust their shapes as they fall, but they do exhibit preferential orientation. Depending on the specific hydrometeor type (pellets, hail, crystals, etc.), both oblate and prolate forms are commonly observed. For the purposes of this work, limited to only the axial ratio, we explore a representative range of 0.5<ai<1.25.
Approximation to the evolution of composition
Long-term atmospheric composition records show that, beyond seasonal patterns, the atmospheric fraction of CO2 has been steadily increasing, while the fraction of O2 has been decreasing (e.g. ). Both trends reflect carbon-oxidation processes. The quantitative relationship between the two is not fixed, and depends on the stoichiometry of oxidation processes, as well as on additional sources and sinks of carbon and oxygen.
The predominant forms of carbon oxidation would suggest that substitution of O2 to CO2 in air approximately a 1:1 molar ratio, implying that the molar sum of O2 and CO2 should remain nearly constant. This behaviour is often observed in short-term laboratory studies (e.g. ), and in urban environments (e.g. ). However, long-term background observations of molar fraction in well-mixed non-urban areas, show that the decline in O2 is significantly larger – by roughly a factor of two, () – than the increase of CO2. This disparity is interpreted as evidence that a substantial fraction of the newly produced atmospheric CO2 is transferred to land and ocean reservoirs, where it accumulates – particularly in dissolved form within the oceans – while atmospheric O2 experiences comparatively limited exchange .
Given the availability of multi-decadal background observations from numerous stations (see ), we may produce fits to the composition over recent decades. A quadratic function of time provides a good fit to the deseasonalized record from 1958 to the present. Because there is a significant latitudinal gradient, a hemispheric asymmetry term is included, proportional to the sine of latitude λ. Since the period of interest for the present study corresponds to the availability of GNSS observations (approximately 1994–present), time is expressed relative to the year 2000, as y2000=y-2000. The residuals of the fit are small, suggesting that expression can reasonably be extrapolated into the near future.
xCO2(y,λ)=10-6⋅[368.625+1.798⋅y2000+0.0118⋅y20002+2.224⋅sin(λ)]
which captures both a linear trend and a measurable acceleration.
Oxygen is not measured by in absolute units, but relative to a reference sample. The deseasonalized trends from this relative dataset constrain the rate of decline, its acceleration, and the latitudinal dependence. The absolute level is set using measurements by , who report xO2=0.209392±3, at the year 2000 and at latitude 24N. The constant term in Eq. () is chosen to reproduce that value, with the remaining terms reproduce the relative measurements:
xO2(y,λ)=10-6⋅[209393-3.953⋅y2000-0.0363⋅y20002-3.064⋅sin(λ)]
Both fits are intended to capture the long-term background evolution. These measurements reflect well-mixed tropospheric air, which becomes homogenized over a time scale of a few months. Stratospheric air exchanges more slowly with the troposphere, leading to detectable differences in CO2 concentration (e.g. ) and to a typical lag of about 4 years in stratospheric CO2 relative to the deseasonalized tropospheric mean. This lag varies with latitude and altitude, typically between 1–6 years . A further refinement could incorporate the seasonal cycle in the troposphere – especially pronounced in the northern hemisphere and reaching a few ppm – but for the purposes of this study the above expressions are sufficiently accurate. In particular, the absolute abundance of O2 is not known to within better than about 3 ppm, comparable to the amplitude of the seasonal cycle.
Proposed expression
An ansatz analogous to the one proposed in is adopted in this study. Relative to that earlier formulation, two principal modifications are introduced. First, the main dry-air term is no longer treated as constant, because the atmospheric composition is not assumed to be static. Although the resulting dependence is small for compositions observed in recent decades, it is included here to assess the long-term signature of composition changes in refractivity. Second, additional terms are incorporated to represent the contributions of precipitated water. These are known to be moderately small , but particular attention is paid in order to explore the optical anisotropy that they may develop. The proposed ansatz is:
17N≡(n-1)⋅106=N0⋅1+10-66N018N0=q1xO2,xCO2+q2⋅τρd+19q3+q4⋅τρw+20q5⋅fl(al;p)⋅ρl+21q6⋅fi(ai;p)⋅ρi
where the refractivity N depends on the amount of matter per unit volume in each component: dry air ρd, water vapour ρw, liquid water ρl, and frozen water ρi. Refractivity also depends on the absolute temperature T through τ=273.15K/T-1. The function q1, which is nearly constant, is assumed to take the form:
q1xO2,xCO2=q10+q11⋅xO2-0.2095+q12⋅xCO2
where the first term represents the average composition of dry air, in the limit of vanishing CO2. The second accounts for small deviations of the O2 molar fraction relative to the reference value of xO2=0.2095, close to present-day conditions. The third term represents the contribution of CO2. A linear dependence is considered sufficient for both the second and third terms, because variations of O2 around the reference value are small, and CO2 remains a trace constituent.
Since Eq. () is expressed in terms of density, and its conversion to pressure is frequently required, we provide the corresponding expression for the mean molar mass of dry air that is consistent with the assumed composition.
mdxO2,xCO2=md0+md1⋅xO2-0.2095+md2⋅xCO2
The polarization functions are modelled as quadratic functions of the axial ratio and depend on the linear polarization state p, which may be aligned with the rotation axis of the hydrometeor ellipsoids (denoted vertical or V), or perpendicular to it (horizontal, H), see Fig. ). For liquid inclusions:
24fl(al;V)=1+c1V⋅(al-1)+c2V⋅(al-1)225fl(al;H)=1+c1H⋅(al-1)+c2H⋅(al-1)2
and solid:
26fi(ai;V)=1+d1V⋅(ai-1)+d2V⋅(ai-1)227fi(ai;H)=1+d1H⋅(ai-1)+d2H⋅(ai-1)2
The closure relation for the depolarization coefficients imposes constraints between the two polarizations, namely: c1V+2⋅c1H=0, and d1V+2⋅d1H=0.
Fit to a user expression
The coefficients introduced above are determined by constructing a large ensemble of physical conditions – spanning pressure, temperature, moisture, dry air composition, and hydrometeor shape – that reflects a broad range of realistic atmospheric states.
Pressures are selected to represent conditions from the surface to the mid-stratosphere. Surface temperatures range from -60 to 40 °C, and are extrapolated to upper-air temperatures using typical lapse rates . Water vapour is varied from zero up to 25 % above saturation at the local conditions. To avoid unrealistically moist stratospheric states, a constraint is imposed to prevent the water fraction from increasing with altitude. In supersaturated cases, water vapour is capped at the saturation density, with any excess assigned to hydrometeors. Hydrometeor mass therefore reaches up to 25 % of the saturation density of the local saturation vapour density, partitioned liquid or solid according to temperature.
Carbon dioxide concentrations range from 300 and 450 ppm, and the O2 molar fraction from 0.209 and 0.210. Axial ratios (vertical/horizontal) of hydrometeor ellipsoids span the interval 0.5–1.25.
The full ensemble comprises more than 150 000 atmospheric states. For each state, the microphysical model is used to compute the refractivity for both linear polarizations, and the molar mass of dry air is evaluated.
To ensure numerical stability, the fitting procedure is carried out in two stages. In the first stage, all coefficients not related to hydrometeor asphericity are adjusted. This stage involves the expressions between Eqs. () and (), and determines coefficients q1 to q6, including q10 to q12, and md0 to md2. Only the subset of conditions that contain either no hydrometeors, or hydrometeors of spherical shape is included. This first fit yields the following expression:
28N0=q1(xO2,xCO2)+0.097⋅τρd+296703.497+6393.484⋅τρw+301447.827⋅fl(al;p)⋅ρl+31686.944⋅fi(ai;p)⋅ρi
with refractivity expressed in N-units, densities in kg m−3, and τ=(273.15K/T-1). The function q1 is determined from the first-stage fit to be:
q1(xO2,xCO2)=222.637-51.817⋅xO2-0.2095+30.266⋅xCO2
To preserve the intended accuracy of the formulation, the expression must be used together with a consistent relationship between mass and molar density, and therefore with a consistent value of the mean molar mass of dry air, which is likewise obtained from the fit:
md(xO2,xCO2)=28.95949+3.985⋅xO2-0.2095+15.996⋅xCO2
As noted in Sect. , the molar mass of water is taken as constant in this model, with value:
mw=18.01525
With all these coefficients fixed, a second fit is performed to determine the shape-dependent and polarization-dependent functions, now extending the physical conditions to include aspherical hydrometeors. For liquid hydrometeors, the functions for the two linear polarizations are:
35fl(al;V)=1+0.743⋅(al-1)+0.043⋅(al-1)236fl(al;H)=1-0.371⋅(al-1)+0.753⋅(al-1)2
and for ice:
37fi(ai;V)=1+0.330⋅(ai-1)-0.125⋅(ai-1)238fi(ai;H)=1-0.165⋅(ai-1)+0.215⋅(ai-1)2
By construction, this second fit satisfies the closure condition on the coefficients that are linear in (a-1).
Across the full range of atmospheric conditions considered – spanning and exceeding those expected in practice – the fitted expressions reproduce the microphysical model with an average absolute residual of 0.0016%, and a maximum residual of 0.009%. This performance exceeds both the estimated accuracy of the microphysical model and the intended accuracy of this work. The overall accuracy is therefore limited solely by the microphysical model itself, not by the fitting procedure.
A time-dependent form of the proposed expression
The fit to the microphysical model spans a range of physical conditions and dry-air compositions representative of recent decades and the foreseeable future. Two composition parameters are retained as independent variables, appearing in the functions q1xO2,xCO2 and mdxO2,xCO2. Because measurements exist for both quantities, and can be described, for example, by the functions introduced in Sect. , it is possible to further reduce q1 and md to explicit functions of time. In this section, we illustrate this procedure using Eqs. () and (). The user must judge whether adopting only the long-term trends of O2 and CO2 is sufficient, or whether latitudinal, seasonal, or altitude variations should also be incorporated.
Using the temporal fits of molar fractions yields:
q1xO2,xCO2≃q1(y)=222.654+0.000259⋅y2000+2.24×10-6⋅y20002
and
mdxO2,xCO2≃md(y)=28.96496+1.30×10-5⋅y2000+4.41×10-8⋅y20002
A time-dependent user expression is therefore obtained by combining Eqs. (), () and ().
A particularly relevant limit is that of pure dry air, often invoked in a variety of atmospheric applications:
Ndry=q1(xO2,xCO2)+0.097⋅τρd
It is useful to compare this expression with others commonly found in the literature, such as the classical form in :
Ndrya=k1PT
in which the parameter here labeled k1 is not exactly a constant, but is modulated by the compressibility factor Z, for instance in :
NdryTh74=77.60ZPT
An accurate expression of the compressibility factor from is used in the microphysical model. For the purposes of this section, however, it is sufficient to adopt the approximation:
Zapp≈1-P⋅5.789×10-9-3.512×10-8⋅τ
with P in pascals. In the dry-air limit – appropriate for instance to the stratosphere, and to good approximation, in dry tropospheric air where deviations of Z from unity remain small –, the following holds:
k1=q1⋅md10⋅RZ≃77.5655Zapp+1.25⋅10-4⋅y2000+8.98⋅10-7⋅y20002
If we compare this with the formulation of AL11, which used a composition approximately representative of 2010, we would obtain here k1≃77.5668Zapp, against an asymptotic low-density value of AL11 is k1≃77.575Zapp. For the ROMEX experiment, which is performed with measurements from 2022, the asymptotic low-density behaviour is k1≃77.5687Zapp. The difference between AL11 and the current model arises primarily from the use of a slightly different fit in to the same underlying laboratory data for dielectric gas properties, and also to a different choice of atmospheric conditions to reduce the microphysical model to a practical expression. Both formulations remain compatible within their intended target of accuracy. Among these available dielectric fits, the expression adopted here is preferred because it provides uncertainty estimates for all coefficients, thereby enabling traceability of the error budget back to the laboratory measurements, which is one of the purposes of this study.
Double refraction
The refractivity model in Eq. () describes an isotropic medium in the absence of hydrometeors, or when these are spherical (alρl=0 and aiρi=0). In contrast, aligned aspherical hydrometeors introduce anisotropy, yielding different refractivities for test fields (i.e. linear polarizations) aligned with the major and minor axes of the hydrometeors. In this section we estimate, at an order of magnitude level, the resulting impact on wave propagation in the two polarizations, focusing on the geometry of a radio occultation.
Because GNSS applications are a primary motivation for the present formulation, it is useful to recall that the carrier of a standard GPS signal is, upon emission, right-hand circularly polarized. This can be viewed as the superposition of two orthogonal linearly polarized components, one of them delayed by a quarter wavelength with respect to the other.
The electric-field polarization is defined in the plane perpendicular to the local Poynting vector P. In an occultation (a limb observation) this vector is approximately horizontal, and the normal plane can be described with a basis of a horizontal and a vertical orthogonal unit vectors, which align to the major and minor axes of the hydrometeors; see Fig. . In contrast, in a quasi-zenithal observation, as with Ground-Based GNSS receivers, P is aligned to the hydrometeors' rotation axis. The plane perpendicular to P is then symmetric concerning Eq. (), and all polarizations behave as the H component in that equation, with identical refractivity.
Polarization components in a GNSS observation may therefore experiance different refractive indices. In the limb-viewing geometry, the plane of the electric field contains a subspace aligned with a horizontal axis of the hydrometeors. Electic activity in this subspace senses a refraction index nH, a component referred to as ordinary (, Sect. 15.3.2). Activity in the perpendicular subspace of the same plane, which is aligned with the vertical minor axis V, senses instead nV, a component referred as extraordinary.
If both polarizations encounter the same refractive index, they follow identical trajectories and retain their relative phase delay. If instead one polarization experiences a higher refractive index, two effects arise:
a modification of the relative delay between the two components, and
a differential geometric propagation path (i.e. double refraction) between transmitter and receiver.
Although both effects vanish in many atmospheris situations (no hydrometeors, or hydrometeors of spherical average shape), we estimate here the magnitude they may reach under extreme but still realistic conditions during an occultation in the low troposphere:
a vertical/horizontal axis ratio al=0.5 for liquid or ai=1.25 for ice, representative of large raindrops and hail ;
maximum precipitated densities of ρl=10-2 kg m−3 for rain and ρi=4×10-3 kg m−3 for frozen hydrometeors;
a horizontal propagation distance of 50km under these conditions (, and references therein).
The rain case suggests a liquid contribution to the total refractivity (in N-units) of about Nl≈15, with a large asymmetry between polarizations: NlH≈20 and NlV≈10. Here, NlH>NlV, because the drops are substantially oblate . In the hail case, contribution to refractivity is smaller <3, with moderate asymmetry, NiH≈2.7 and NiV≈3.0, i.e. NiH<NiV, for the assumed prolateness.
Accumulated over the 50km of horizontal propagation, the refractivity asymmetry yields an optical path difference (at equal geometric paths) of ΔLH-V≈0.5m (H minus V polarizations) for the rain case, and ΔLH-V≈-0.015m for the hail case. Rain is both more refractive and can reach larger extreme mass densities than ice, but even in these extreme cases the resulting differences remain moderate. These path differences are detectable by receivers capable of discriminating both polarizations, see for instance .
In addition to producing a differential optical path length, a difference in refractivity between the two polarization components also causes their propagation paths to split (birefringence). Under typical occultation geometry, and given an assumed typical background atmospheric gradient, the rain example implies that the horizontal polarization component propagates above the vertical polarization, as the horizontal component encounters a larger refractivity at the same altitude, in this example by about 200m. The example is chosen as an extreme case, derived from large but still realistic rain density, sustained over a large geographic area (infrequent under intense rain). Smaller values of the tangent point split, of 10s of meters, should be more common but still detectable. At fixed reception time, an occulting signal containing both polarizations (e.g. a circularly polarized GPS signal) would therefore exhibit different tangent point heights in each polarization, separated by roughly this amount. This is illustrated in Fig. .
Schematic representation of the tangent point (TP) separation between polarizations in the case of oblate hydrometeors, typical of rain. In the path between two Earth-orbiting satellites, the horizontal polarization component (H) senses a higher refractive index, and follows a longer (higher) path than the vertical polarization (V). The scale is exaggerated for clarity.
During the occultation scan, the tangent point reaches a given altitude with a time delay between the two components; with a typical vertical scan rate of the order of 1kms-1, the delay reaches about 0.2s in this extreme scenario – still smaller than the Fresnel diameter. The reduction of scan rate in the low atmosphere would favor the detectability. In the hail example case, the path split is correspondingly smaller, about 10m, and the order of the polarizations is reversed: the vertical component lies above the horizontal, consistent with the assumed prolateness. The resulting time difference between measurements at fixed tangent altitude should also be detectable, in addition to the path delay, by receivers capable of discriminating both linear polarizations.
A dual linearly polarized antenna may thus detect certain sharp features such as the Planetary Boundary Layer (PBL) with a small time difference if nonspherical hydrometeors are present. With a standard circularly polarized antenna, the time difference may instead blur the feature.
List of microphysical input quantities
We summarize here the quantities – and their sources – that play a significant role in the microphysical model and thus in the determination of the final refractivity expression. These quantities are obtained from external experiments and come with associated uncertainties. They are considered sufficiently accurate for deriving a useful refractivity model and for assessing the influence of individual physical processes on the observable refractivity. Nonetheless, the expression retains a non-negligible residual uncertainty and, if SI-traceability is desired, it becomes important to associate this final uncertainty with the underlying experimental inputs.
The microphysical formulation contains a set of parameters pj (molecular polarizabilities, electric dipoles, standard concentrations, magnetic susceptibilities, etc). A final user-oriented expression was obtained by numerical optimization, yielding a set of coefficients qi, fitted to the microphysical model over a wide range of realistic atmospheric conditions. Realistic deviations of pj within their experimental uncertainties lead, over the same set of physical conditions, to modified best-fit coefficients qi. To quantify this dependence, derivatives
dij=dqidpj
were evaluated.
Because of the chosen ansatz for the macroscopic expression, most of the derivatives are negligible, and most experimental inputs influence only a single coefficient. The input parameters whose uncertainties can produce significant changes in the fitted qi are collected in Table .
Collected list of input quantities that have a sizeable impact in the final user expression
Table collects the list of primary inputs entering the microphysical model, and correspond to quantities that are either mdirectly measurable in laboratory conditions or determined through dedicated measurement campaigns. In total, the list contains 19 real-valued primary inputs.
The user expressions ()–() contain a number of coefficients obtained by fitting to the microphysical model constructed upon those primary inputs. Consequently, each user-level coefficient inherits the dependence on the numerical values, and uncertainties of the primary inputs from which the microphysical model is built.
The user expression contains 18 coefficients, grouped as follows:
4 coefficients to describe the reference gas phase –namely refractivity of dry air and water vapour, and associated temperature dependences– following the structure of .
2 coefficients describing the refractivity sensitivity dependence to atmospheric composition, specifically linked to variations in CO2 and O2.
4 coefficients related to average molar masses, three for the dry-air fraction (allowing variable CO2 and O2 composition), and one for water vapour.
8 coefficients describing contributions of hydrometeors – four for each condensed phase (liquid and solid) – including their dependence on hydrometeor shape under vertical and horizontal field polarizations.
In principle, propagating the dependence of these fitted coefficients on the 19 primary inputs would require the full matrix of partial derivatives, of dimension (18×19). However, the structure of the ansatz was chosen explicitly to minimize cross-dependencies between primary inputs and fitted parameters. As a result, most primary measurements influence only a single coefficient, and their impact on the remaining is either very small or entirely negligible.
The practical traceability table is therefore far simpler than the full derivative matrix. The primary inputs that significantly affect each fitted coefficient are identified in the last column of Table . The reduced set of non-negligible dependencies is presented in Tables and .
List of significant dependencies of the dry user coefficients q1 and q2, with respect to primary inputs. User coefficients are in SI units (N-units ⋅ kg−1 m3). Units of the inputs and derivatives, follow the source units from their respective references, and are diverse: Aϵ and Aτ in cm3 mol, molecular mass in amu, gO2 adimensional, polarizability α in 10-24cm3, electric dipole in Debye units. Uncertainty statistics σ are those that propagate to the coefficients, arising from the estimated uncertainties of the inputs.
As Table , for moist user coefficients q3 and q4, precipitated terms q5 and q6, and main molar mass coefficients md0 and mw.
QuantityCoef valueDependenciesDerivativeDeriv. valueσCoefMain moist constant q36703.497Polarizability H2O∂q3/∂αH2O312.960.03––Elec. Dipole H2O∂q3/∂μH2O6869.420.62––Harris Alder corr.0.03––Composed0.62Moist constant T-dep q46393.484Elec. Dipole H2O∂q4/∂μH2O6868.560.62––Harris Alder corr.0.8––Composed1.0Main liquid term q51447.827Permit. liq. H2O∂q5/∂ϵl,H2O0.690.034Main solid term q6686.943Permit. sol. H2O∂q6/∂ϵi,H2O182.670.91Main dry molar mass md028.95949Mass N2∂md0/∂mN20.780.00031–—Mass O2∂md0/∂mN20.210.00015––Composed0.00035Moist molar mass mw18.01525Fraction 18O/O∂mw/∂x18O2.000.00028Discussion and analysis
The microphysical model and the corresponding fitted user expression have been evaluated using best estimates of the primary laboratory or field-measured parameters listed in Table . These quantities carry uncertainties, that propagate into the user-level coefficients. This propagation serves two purposes: (i) assigning a total uncertainty to each fitted parameter, and (ii) identifying the dominant experimental sources of these uncertainties. The detailed propagation is summarized in Tables (dry-air terms) and (all other terms). As customary, uncertainties are expressed as std errors affecting the last digits. We discuss the fitted coefficients in turn:
Main dry parameter q10: The uncertainty is dominated by the dielectric response of N2. The next most relevant contribution originates from the molar mass of N2 (via the 15N/14N ratio), owing to the explicit dependence of the ansatz on density. The dielectric response of O2, assuming the quoted uncertainty is reliable, contributes very little, whereas the mass uncertainty of O2 (through the 18O/16O ratio) contributes more strongly than its dielectric term. The resulting coefficient, q10=222.637(7), is sufficiently accurate for the intended applications (fractional uncertainty 3×10-5). Improved dielectric measurements of N2 could modestly reduce the uncertainty (down to 1.3×10-5), though at that level limitations of the ansatz (e.g. missing nonlinearities) would likely dominate. This term contributes most to the refractivity.
O2 dry parameter q11: For O2, which is close in mass and polarizability to the mean air molecule, the uncertainty is shared among the dielectric responses and molar masses of both N2 and O2. The fitted value q11=-51.817(12) is sufficiently accurate given the relatively small overall contribution of this term.
CO2 dry parameter q12: The uncertainty is governed by the polarizability and mass of the CO2 molecule. The coefficient q12=30.266(33) has adequate accuracy, as its contribution to refractivity is small.
Temperature-dependent dry term q2: This coefficient is small and would vanish for an ideal gas of nonpolar molecules. It arises from temperature-dependent polarizability (intermolecular interactions), and the magnetic susceptibility of O2. The latter is known very precisely, hence the quoted uncertainty is dominated by the temperature-dependent polarizability of N2, attributed by to anharmonic vibrational effects. The coefficient q2=0.097(6) is sufficiently accurate given its modest contribution.
Main moist term q3: The term includes contributions from the electric polarizability and permanent dipole moment of water, toguether with a dependence on molar mass. The uncertainty is dominated by the dipole moment of water, even though its measured value (Table ) is already precise. The coefficient is q3=6703.5(6), and this term is the second largest contributor to refractivity, though only in moist conditions.
Temperature-dependent moist term q4: The uncertainty is again dominated by the electric dipole moment of water, but mostly through the correlation factor g. The dipole is enhanced by molecular interactions described by a Buckingham factor and, critically, by a density and temperature-dependent correlation factor g(ρ,T) which is very specific to the interaction of water with other near molecules. Water forms transient aggregates such as dimers , which increase the effective dipole moment beyond the single-molecule value. While the high-density behaviour of g is well constrained , the low-density regime relevant here relies on the approximation g(ρ,T)≈1+h(T)ρw, which carries significant uncertainty, and dominates the error estimate. The resulting coefficient is q4=6393.5(10). Although moderately important in moist conditions, the current uncertainty is acceptable.
Liquid water term q5: Somewhat surprisingly, this term depends only weakly on the permittivity of liquid water ϵr≈90. With this strong response, water screens itself from the test field, so the uncertainty of ϵr has a modest effect. Most of the remaining uncertainty is aliased from correlations with water vapour in the fitting dataset. The fitted value q5=1447.83(13) is satisfactory given the typically small amounts of liquid water.
Solid water term q6: Ice and snow have much lower permittivity ϵr≈3 and therefore exhibit weaker self-screening. The refractivity contribution is thus more sensitive to the dielectric uncertainty, which dominates the error budget of this term. The coefficient q6=686.94(91) remains adequate given the small contribution of frozen water to refractivity.
A key point is that the dominant uncertainty in the main dry term, q10 derives from the physical properties of the nitrogen molecule. This uncertainty is effectively constant in time to a precision much better than the quoted error. As discussed in Sect. , the dry-air term exhibits a slow temporal drift due to evolving atmospheric composition. However, this drift does not involve N2, thus we can conclude that q10 is constant to high accuracy.
The user model also includes several parameters describing the asphericity of hydrometeors. The model assumes uniform asphericity and a single canting angle, which are not expected to hold with high accuracy. Realistic hydrometeor populations involve heterogeneous and poorly known shape and orientation distributions, far more complex than an ensemble of aligned rotational ellipsoids. Consequently, SI-traceability cannot be guaranteed when aspherical hydrometeors contribute appreciably. For this reason, the propagation of the microphysical error has not been carried to the functions fl(al;p) and fi(ai;p), nor to the eight parameters that govern them. Their purpose is not to yield precise refractivities for hydrometeors of known shape, but rather to bound the plausible range of refractivities when shapes and orientations are inherently uncertain, and to estimate the order of magnitude of hydrometeor-related effects such as polarization-dependent differential delay or the splitting of propagation paths.
We may compare the current proposed formulation with several others in the literature. To illustrate the basic behaviour of dry-air refractivity, it is useful to examine N⋅T/P, which is very close to 77.6 N-units K hPa−1, and for some formulations exactly constant and denoted k1. This comparison is shown in Fig. . Representative atmospheric states reveal a temperature and density dependence: at lower density (lower pressure, higher altitude), the asymptotic value of k1 decreases; at lower altitude, and particularly in colder conditions, k1 increases. Since GNSSRO impact in NWP peaks at around 200–300 hPa (approximately 8–10 km of altitude), the mean effective k1 is likely closer to its lower-density value.
Quantity Ndry⋅T/P evaluated with various available formulations, and including the present one. For realistic examples, the atmospheric states follow the lapse rates of a US Standard atmosphere profile, illustrating typical behavior at different altitudes and latitudes. Warm and cold examples are shown, with surface temperatures of Ts=-30 °C, and Ts=+30 °C. The compared expressions are SW53, TH74, RU02, AL11, and the present work. The current formulation depends on the composition; for consistency with the ROMEX experiment, 2022 values have been adopted.
This work also confirms a refractivity generally in the lower band within the literature formulations, except in colder denser conditions, where suggests results in the upper band. In general, these results suggest that the formulation by appears to be too high by between 0.05 %–0.1 %, which agrees with the appearence of NWP bias if this formulation is used . The formulation reproduces significantly the qualitative dry-air behaviour indicated in and here, although also appears systematically high. The excess, though, is smaller, about 0.03% in dry-air refractivity. However, although the dry-air behaviour of is good, it does not clarify how moist refractivity should be calculated, given the strongly nonideal behaviour of water vapor, which renders the expression ambiguous . The formulation in does not reproduce the temperature or density dependence, but appears to be an average compromise except at the densest conditions (cold conditions in the low troposphere). Finally, compared to AL11, the ensemble of modifications proposed here reduces slightly the average refractivity, by about 0.01%, but maintains most of the density and temperature dependence.
In Fig. it is also interesting to note that, whereas TH74 converges to a constant with altitude (its compressibility factor Z converging to 1), both AL11 and this work continue to show some descent. The compressibility factor in these latter also converges to 1, but these expressions contain some small but non-zero temperature dependence of NT/P, stemming from the magnetic dipole of the O2 molecule, and expressed as term q2 in the proposed expression. Figure also indicates that, although overestimates dry refractivity, it is better at low altitude. In retrospect, this is consistent with the primary application of that paper, which was for electronic distance measurements at surface level. Similarly, underestimates refractivity in most of the troposphere, yet is better in the upper troposphere and low stratosphere, leading to a good average behaviour for GNSSRO.
Within ECCC, neither AL11 nor this formulation introduce a concerning bias. Upon significant increase of the GNSSRO volume, as for instance ROMEX data, thermodynamic indicators such as the bias profiles of temperature or geopotential against radiosondes are not dragged by more than 1 part in 10 000, consistent with the estimation of absolute accuracy of 0.01%. However, initial GNSSRO use at ECCC (e.g. ), while positive, did drag thermodynamic profiles at the level of 0.1%. Since GNSSRO bias is not only a function of the refractivity, it will be summarily mentioned that at ECCC, assimilation of GNSSRO also
uses GNSSRO data in the form of refractivity profiles
accounts for the deviation of air from ideal gas , both in the refractivity/pressure relationship, and in the hypsometric relationship
follows the horizontal drift of the tangent point over a profile
the hypsometric equation considers the variability of the strength of gravity with latitude and altitude, and through a wind-dependent vertical component of the Coriolis force
all of which have some impact in the GNSSRO observation minus background statistics.
Of all these elements, and concerning the refractivity expression and the bias between GNSSRO and an NWP system, the dry component (here, q1 with its subcomponents, and q2) is the most important, as it detemines the behaviour in the core region of GNSSRO, between 300 and 50 hPa, where these data are more accurate and have higher impact). Among the other elements besides refractivity, the non-ideal gas effects are judged as next in importance, with the Coriolis force being small but systematic, thus affecting bias.
Conclusions
A previous study has been revisited, with particular attention to re-evaluating the error budget and to identifying the main limitations to improved accuracy. Beyond several minor refinements, the atmospheric composition is identified as one of the elements that constrain accuracy within the desired performance range. In view of the secular evolution of O2 and CO2 abundances, their fractions are made explicit independent variables.
47N=q1+0.097(6)⋅τρd+486703.5(6)+6393.5(10)⋅τρw+491447.83(13)⋅fl(al;p)⋅ρl+50686.94(91)⋅fi(ai;p)⋅ρi
where q1 may be expressed as a function of O2 and CO2, if these are available:
q1(xO2,xCO2)=222.637(7)-51.817(12)⋅xO2-0.2095+30.266(33)⋅xCO2
In all cases, the parentheses express the uncertainty in the coefficients as stemming from the uncertainty in the experimental sources collected in the microphysical model.
Alternatively, using fits to recent decades of measured atmospheric composition, the expression for q1 may be written instead as a function of time. Such a fit to the composition is suggested in Eqs. () and () and will be valid as long as it continues to represent the main trend of the observed evolution, presumably at least several decades. The time-dependent q1 may then be expressed as:
q1(y)=222.654(7)+0.000259⋅y2000+2.24×10-6⋅y20002
with y2000 time in years since 2000. This is the only coefficient in Eq. () that shows a significant dependence on composition. At Environment and Climate Change Canada, Eqs. () and () have been implemented operationally as an update to AL11, which was in use since 2011. For practical operational use, it was selected to use the time-dependent form of q1, which implicitly describe composition. As noted above, the uncertainty associated with seasonal and latitudinal variability in composition is not dominant relative to the uncertainty in absolute concentrations. Under these circumstances, the expression q1(y), which has been simplified, is not inferior to a more complex q1(xO2,xCO2).
The refractivity expression exhibits a small dependence on the polarization state when hydrometeors (rain, snow, or ice) are non-spherical. This dependence is captured through the shape factors fl(al;p) and fi(ai;p). The approximate magnitude of their impact has been evaluated: it is modest, but is measurable as differential path delay, and should also be measurable as a form of double refraction for rays traversing extended regions of dense precipitation during an occultation. Measuring this effect requires a dual-polarization receiver. Such receivers, sensitive to both linear polarizations independently, have been launched in recent years.
For many applications of the refractivity expression – such as its routine use in the assimilation of refractivity in Numerical Weather Predicition –, this update introduces only small differences, and the changes with respect to are minor. As observed at other operational centers, ECCC identified that the absence of bias between GNSSRO and model fields becomes particularly critical at higher observational volumes. This motivated the implementation of the present update, even though AL11 had remained generally adequate.
The comparison between the different existing formulations suggests that the threshold of 0.1% uncertainty (that is, fractional uncertainty of 10-3) is exceeded, provided a non-ideal gas formulation is adopted. The difference between existing formulations, with some differences between the implementation details, in the assumptions about the nature of atmospheric air, and the link to fundamental measurements, indicates rather a difference of 0.01% (thus 10-4) between studies, provided they account for non-ideal gas effects. Finally, the formal uncertainty as propagated in this study from fundamental measurements indicates a fractional dry-air uncertainty of about 3×10-5 of dielectric origin, thus if the variabilities in composition are accounted.
This study suggests that, whereas a confirmation of the low dielectric uncertainty would be welcome, the remaining sources of refractive uncertainty are, first, of thermodynamic origin, in the non-ideal gas behaviour of air, and particularly water vapour, and second, in details of composition, from the evolution of CO2 and O2 to uncertainty in the molar mass, by uneven isotopic composition, where notably water vapour and N2 should be mentioned. These diverse residual sources of uncertainty – of the order of 10-5 –, limit our knowledge of refractivity to not be better than 10-4 in fractional uncertainty unless a model of considerably more complex nature was prepared. However, these residual sources of uncertainty do not prevent reaching that stated level of 10-4, or 0.01%, which we judge as sufficient for NWP purposes. Concerning hydrometeors, the model does maintain this stated accuracy even under the highest densities of hydrometeors that are likely in the atmosphere, for a given amount of liquid water, ice, their shape and orientation. It is however unlikely that the atmosphere contains amounts of water, ice, shapes and orientations that are sufficiently uniform for this to be meaningful in practice. The natural heterogeneity of hydrometeors introduces an uncertainty in the use of an electromagnetic model, above the physical uncertainty of that model. In the atmosphere, that may amount to perhaps as much as 50% of the hydrometeor contribution to refractivity. However, the hydrometeor contribution is generally less than 1% of the total refractivity, and in the upper atmosphere much less. Therefore this heterogeneity may dominate uncertainty only when the hydrometeor contribution is large, which is not the case in most of the upper atmosphere, and may happen only at hydrometeor densities above 0.02 (grams of hydrometeor per kg of air). It is thus recommended to allow for a hydrometeor-dependent uncertainty estimate, composed of the mentioned 10-4 fractional uncertainty of the electromagnetic model, plus a fraction of the estimated hydrometeor contribution, derived not from electromagnetic uncertainty but from their large variability.
Several other applications may benefit from the expressions presented here. These include:
Radio-optical measurements, such as GNSSRO, in decadal climate studies, where cumulative changes in atmospheric composition produce a non-negligible drift.
Situations where traceability between radio-optical and thermodynamic measurements is essential, since traceability requires a quantified uncertainty chain.
Investigations of the polarization dependence introduced by hydrometeors.
Estimation of atmospheric delay, such as GNSS zenith total delay (ZTD), or subtraction of this delay, such as radar altimeters
Code availability
The code expressing the microphysical model, and to evaluate the fit to user expressions, can be obtained from the author upon request.
Data availability
The data that were used as part of the microphysical model, fit to compositions, and meteorological states can be obtained from the author upon request. The composition data were obtained from the University of California San Diego databases 10.6075/J0542KSG and 10.6075/J0WS8RJR.
Author contributions
JMA prepared the conceptual plan, assembled the needed information, developed the required tracing software, and wrote the manuscript.
Competing interests
The author has declared that there are no competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Special issue statement
This article is part of the special issue “The Radio Occultation Modeling EXperiment (ROMEX): observational quality, processing, and numerical weather prediction (NWP) applications”. It is not associated with a conference.
Acknowledgements
The author acknowledges comments received from Maria Teresa Hernandez (Sinia Research) on uncertainty propagation, from Ulrich Foelshe and Gottfried Kirchengast (Wegener Center) on the evolution of atmospheric composition, as well as general comments from Louis Garand, Stéphane Laroche and Jean-François Cossette (Environment Canada), Anthony Mannucci (JPL), and Robert Kursinski (PlanetIQ).
Review statement
This paper was edited by Richard Anthes and reviewed by two anonymous referees.
ReferencesAndrews, A. E., Boering, K. A., Daube, B. C., Wofsy, S. C., Loewenstein, M., Jost, H., Podolske, J. R., Webster, C. R., Herman, R. L., Scott, D. C., Flesch, G. J., Moyer, E. J., Elkins, J. W., Dutton, G. S., Hurst, D. F., Moore, F. L., Ray, E. A., Romashkin, P. A., and Strahan, S. E.: Mean ages of stratospheric air derived from in situ observations of CO2, CH4, and N2O, J. Geophys. Res.-Atmos., 106, 32295–32314, 10.1029/2001JD000465, 2001. Anthes, R. A., Rocken, C., and Kuo, Y. H.: Applications of COSMIC to Meteorology and Climate, Terr. Atmos. Ocean. Sci., 11, 115–156, 2000.Anthes, R. A., Marquardt, C., Ruston, B., and Shao, H.: Radio Occultation Modeling Experiment (ROMEX): Determining the Impact of Radio Occultation Observations on Numerical Weather Prediction, B. Am. Meteoror. Soc., 105, E1552–E1568, 10.1175/BAMS-D-23-0326.1, 2024.Ao, C. O., Meehan, T. K., Hajj, G. A., Mannucci, A. J., and Beyerle, G.: Lower troposphere refractivity bias in GPS occultation retrievals, J. Geophys. Res.-Atmos., 108, 4577, 10.1029/2002JD003216, 2003.Aparicio, J. M. and Deblonde, G.: Impact of the Assimilation of CHAMP Refractivity Profiles on Environment Canada Global Forecasts, Mon. Weather Rev., 133, 257–275, 10.1175/2007MWR1951.1, 2008.Aparicio, J. M. and Laroche, S.: An evaluation of the expression of the atmospheric refractivity for GPS signals, J. Geophys. Res.-Atmos., 116, D11104, 10.1029/2010JD015214, 2011.Aparicio, J. M. and Laroche, S.: Estimation of the Added Value of the Absolute Calibration of GPS Radio Occultation Data for Numerical Weather Prediction, Mon. Weather Rev., 143, 1259–1274, 10.1175/MWR-D-14-00153.1, 2015.Aparicio, J. M., Deblonde, G., Garand, L., and Laroche, S.: Signature of the atmospheric compressibility factor in COSMIC, CHAMP, and GRACE radio occultation data, J. Geophys. Res.-Atmos., 114, D16114, 10.1029/2008JD011156, 2009.Bani-Shahabadi, M., Aparicio, J. M., and Garand, L.: Impact of Slant-Path Radiative Transfer in the Simulation and Assimilation of Satellite Radiances in Environment Canada's Weather Forecast System, Mon. Weather Rev., 146, 4357–4372, 10.1175/MWR-D-18-0126.1, 2018.Bani-Shahabadi, M., Buehner, M., Aparicio, J. M., and Garand, L.: Implementation of Slant-Path Radiative Transfer in Environment Canada's Global Deterministic Weather Prediction System, Mon. Weather Rev., 148, 4231–4245, 10.1175/MWR-D-20-0060.1, 2020.Beard, K. V. and Chuang, C.: A New Model for the Equilibrium Shape of Raindrops, J. Atmos. Sci., 44, 1509–1524, 10.1175/1520-0469(1987)044<1509:ANMFTE>2.0.CO;2, 1987.Beard, K. V., Bringi, V. N., and Thurai, M.: A new understanding of raindrop shape, Atmos. Res., 97, 396–415, 10.1016/j.atmosres.2010.02.001, 2010. Born, M. and Wolf, E.: Principles of Optics, Cambridge University Press, Cambridge, UK, 7th edn., ISBN 9780521642224, 1999.Bruggeman, D. A. G.: Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen. I. Dielektrizitätskonstanten und Leitfähigkeiten der Mischkörper aus isotropen Substanzen, Ann. Phys., 416, 665–679, 10.1002/andp.19354160802, 1935.Buckingham, A. D.: A Theory of the Dielectric Polarization of Polar Substances, Proc. R. Soc. Lond. A Mat., 238, 235–244, 10.1098/rspa.1956.0216, 1956.Cardellach, E., Tomás, S., Oliveras, S., Padullés, R., Rius, A., de la Torre-Juárez, M., Turk, F. J., Ao, C. O., Kursinski, E. R., Schreiner, B., Ector, D., and Cucurull, L.: Sensitivity of PAZ LEO polarimetric GNSS radio-occultation experiment to precipitation events, IEEE T. Geosci. Remote, 53, 190–206, 10.1109/TGRS.2014.2320309, 2015.CGMS: Report of the 40th meeting, Tech. rep., Coordination Group for Meteorological Satellites, Lugano, Switzerland, https://www.cgms-info.org/wp-content/uploads/2021/10/cgms-40-report.pdf (last access: 6 August 2026), 2012.Choy, T. C.: Effective Medium Theory: Principles and Applications, Oxford University Press, Oxford, UK, 10.1093/acprof:oso/9780198705093.001.0001, 1999.Cucurull, L., Derber, J. C., Treadon, R., and Purser, R. J.: Assimilation of Global Positioning System Radio Occultation Observations into NCEP's Global Data Assimilation System, Mon. Weather Rev., 135, 3174–3193, 10.1175/MWR3461.1, 2007.Eötvös, R.: Experimenteller Nachweis der Schwereänderung, die ein auf normal geformter Erdoberfläche in östlicher oder westlicher Richtung bewegter Körper durch diese Bewegung erleidet, Ann. Phys., 364, 743–752, 10.1002/andp.19193641604, 1919.Eyre, J.: A bias correction scheme for simulated TOVS brightness temperatures, ECMWF Technical Memoranda, 35, 10.21957/tmhrqv5cp, 1992.Eyre, J. R.: Observation bias correction schemes in data assimilation systems: a theoretical study of some of their properties, Q. J. Roy. Meteor. Soc., 142, 2284–2291, 10.1002/qj.2819, 2016.FAA: Title 14 Code of Federal Regulations: Part 33 – Airworthiness Standards: Aircraft Engines, Appendix B, https://www.ecfr.gov/current/title-14/chapter-I/subchapter-C/part-33 (last access: 6 August 2026), 2016.Fernández, D. P., Goodwin, A. R. H., Lemmon, E. W., Sengers, J. M., and Williams, R. C.: A formulation for the static permittivity of water and steam at temperatures from 238 K to 873 K at pressures up to 1200 MPa, including derivatives and Debye–Hückel coefficients, J. Phys. Chem. Ref. Data, 26, 1125–1166, 10.1063/1.555997, 1997.Garnett, J. C. M.: Colours in Metal Glasses and in Metallic Films, Proc. R. Soc. Lond. A Mat., 203, 385–420, 10.1098/rsta.1904.0024, 1904.Georgii, H. W. and Jost, D.: Concentration of CO2 in the upper troposphere and lower stratosphere, Nature, 221, 1040, 10.1038/2211040a0, 1969.Healy, S. B.: Refractivity coefficients used in the assimilation of GPS radio occultation measurements, J. Geophys. Res., 116, D01106, 10.1029/2010JD014013, 2011.Healy, S. B. and Thépaut, J. N.: Assimilation experiments with CHAMP GPS radio occultation measurements, Q. J. Roy. Meteor. Soc., 132, 605–623, 10.1256/qj.04.182, 2006.IAPWS: Release on the Static Dielectric Constant of Ordinary Water Substance for Temperatures from 238 K to 873 K and Pressures up to 1000 MPa, Tech. rep., International Association for the Properties of Water and Steam, Erlangen, Germany, https://iapws.org/documents/release/Dielec (last access: 6 August 2026), 1997.IAPWS: Guideline on the Use of Fundamental Physical Constants and Basic Constants of Water, Tech. rep., International Association for the Properties of Water and Steam, Gaithersburg, MD, USA, https://iapws.org/documents/release/fundam (last access: 6 August 2026), 2001.Jones, S. B. and Friedman, S. P.: Particle shape effects on the effective permittivity of anisotropic or isotropic media consisting of aligned or randomly oriented ellipsoidal particles, Water Resour. Res., 36, 2821–2833, 10.1029/2000WR900198, 2000.Keeling, C. D., Piper, S. C., Bacastow, R. B., Wahlen, M., Whorf, T. P., Heimann, M., and Meijer, H. A.: Exchanges of atmospheric CO2 and 13CO2 with the terrestrial biosphere and oceans from 1978 to 2000. I. Global aspects, Tech. Rep. 01–06, Scripps Institution of Oceanography, San Diego, CA, USA, https://escholarship.org/uc/item/09v319r9 (last access: 6 August 2026), 2001.Keeling, R. F.: Measuring correlations between atmospheric oxygen and carbon dioxide mole fractions: A preliminary study in urban air, J. Atmos. Chem., 7, 153–176, 10.1007/BF00048044, 1988.Keeling, R. F.: Atmospheric CO2 data, 10.6075/J0542KSG, 2015a.Keeling, R. F.: Atmospheric O2 data, 10.6075/J0WS8RJR, 2015b.Keeling, R. F. and Morgan, E. J.: Scripps O2 Program Data, UC San Diego Library Digital Collections [data set], 10.6075/J0WS8RJR, 2019.Keeling, R. F., Morgan, E. J., and Keeling, C. D.: Scripps CO2 Program Data, UC San Diego Library Digital Collections [data set], 10.6075/J0542KSG, 2017.Kessler, E.: Condensate Content in Relation to Sloping Updraft Parameters, J. Atmos. Sci., 32, 443–443, 10.1175/1520-0469(1975)032<0443:CCIRTS>2.0.CO;2, 1975.Knight, C. A. and Knight, N. C.: Very large hailstones from Aurora, Nebraska, B. Am. Meteor. Soc., 86, 1773–1781, 10.1175/BAMS-86-12-1773, 2005.Kursinski, E. R., Hajj, G. A., Schofield, J. T., Linfield, R. P., and Hardy, K. R.: Observing Earth's atmosphere with radio occultation measurements using the Global Positioning System, J. Geophys. Res.-Atmos., 102, 23429–23465, 10.1029/97JD01569, 1997.Landau, L. D., Lifshitz, E. M., and Pitaevskii, L. P.: Electrodynamics of Continuous Media, Pergamon Press, Oxford, UK, 2nd edn., 10.1016/B978-0-08-030275-1.50007-2, 1984. Lide, D. R.: CRC Handbook of Chemistry and Physics, CRC Press, New York, USA, 82nd edn., ISBN: 978-0-8493-0482-8, 2001.Liebe, H. J., Hufford, G. A., and Cotton, M. G.: Propagation modeling of moist air and suspended water/ice particles at frequencies below 1000 GHz, in: AGARD Conference Proceedings 542: Atmospheric Propagation Effects Through Natural and Man-Made Obscurants for Visible to mm-Wave Radiation, NATO, https://its.ntia.gov/software/millimeter-wave-propagation-model-mpm (last access: 6 August 2026), 1993.Malmberg, C. G. and Maryott, A. A.: Dielectric Constant of Water from 0 to 1000 C, J. Res. Nat. Bur. Stand., 56, 1–8, 10.6028/jres.056.001, 1956.Manning, A. C. and Keeling, R. F.: Global oceanic and land biotic carbon sinks from the Scripps atmospheric oxygen flask sampling network, Tellus B, 58, 95–116, 10.1111/j.1600-0889.2006.00175.x, 2006.Mätzler, C. and Wegmüller, U.: Dielectric properties of freshwater ice at microwave frequencies, J. Phys. D: Appl. Phys., 20, 1623, 10.1088/0022-3727/20/12/013, 1987.May, E. F., Moldover, M. R., and Schmidt, J. W.: Electric and magnetic susceptibilities of gaseous oxygen: Present data and modern theory compared, Phys. Rev. A, 78, 32522–32536, 10.1103/PhysRevA.78.032522, 2008.Melbourne, W. G., Davis, E. S., Duncan, C. B., Hajj, G. A., Hardy, K., Kursinski, E. R., Meehan, T. K., Young, L. E., and Yunck, T. P.: The Application of Spaceborne GPS to Atmospheric Limb Sounding and Global Change Monitoring, 94, JPL, Pasadena, CA, https://ntrs.nasa.gov/api/citations/19960008694/downloads/19960008694.pdf (last access: 6 August 2026), 1994.NOAA, NASA, and USAF: U.S. Standard Atmosphere 1976, U.S. Government Printing Office, Washington, DC, USA, https://ntrs.nasa.gov/citations/19770009539 (last access: 6 August 2026), 1976.Padullés, R., Ao, C. O., Turk, F. J., de la Torre Juárez, M., Iijima, B., Wang, K. N., and Cardellach, E.: Calibration and validation of the Polarimetric Radio Occultation and Heavy Precipitation experiment aboard the PAZ satellite, Atmos. Meas. Tech., 13, 1299–1313, 10.5194/amt-13-1299-2020, 2020.Padullés, R., Cardellach, E., and Turk, F. J.: On the global relationship between polarimetric radio occultation differential phase shift and ice water content, Atmos. Chem. Phys., 23, 2199–2214, 10.5194/acp-23-2199-2023, 2023.Padullés, R., Cardellach, E., Paz, A., and Burger, T.: Initial Polarimetric Radio Occultation Results from Spire's Nanosatellite Constellation: Independent Assessment and Potential Applications, B. Am. Meteor. Soc., 106, E735–E751, 10.1175/BAMS-D-23-0322.1, 2025.Park, S. Y., Kim, J. S., Lee, J. B., Esler, M. B., Davis, R. S., and Wielgosz, R. I.: A redetermination of the argon content of air for buoyancy corrections in mass standard comparisons, Metrologia, 41, 387, 10.1088/0026-1394/41/6/005, 2004.Picard, A., Davis, R. S., Gläser, M., and Fujii, K.: Revised formula for the density of moist air (CIPM-2007), Metrologia, 45, 149–155, 10.1088/0026-1394/45/2/004, 2008.Poli, P., Moll, P., Puech, D., Rabier, F., and Healy, S. B.: Quality Control, Error Analysis, and Impact Assessment of FORMOSAT-3/COSMIC in Numerical Weather Prediction, Terr. Atmos. Ocean. Sci., 20, 101–113, 10.3319/TAO.2008.01.21.02(F3C), 2009.Pruppacher, H. R. and Beard, K. V.: A wind tunnel investigation of the internal circulation and shape of water drops falling at terminal velocity in air, Q. J. Roy. Meteor. Soc., 96, 247–256, 10.1002/qj.49709640807, 1970.Pruppacher, H. R. and Pitter, R. L.: A Semi-Empirical Determination of the Shape of Cloud and Rain Drops, J. Atmos. Sci., 28, 86–94, 10.1175/1520-0469(1971)028<0086:ASEDOT>2.0.CO;2, 1971.Rennie, M.: The assimilation of GPS radio occultation measurements at the Met Office, in: GRAS SAF Workshop on Applications of GPSRO measurements, 84–94, European Centre for Medium-Range Weather Forecast, Reading, UK, https://www.ecmwf.int/sites/default/files/ (last access: 6 August 2026), 2008.Rüeger, J. M.: Refractive Index Formulae for Electronic Distance Measurement with Radio and Millimetre Waves, Tech. Rep. S-68, School of Surveying and Spatial Information Systems, University of New South Wales, Sydney, Australia, https://www.unsw.edu.au/content/dam/pdfs/engineering/ (last access: 6 August 2026), 2002.Schmidt, J. W. and Moldover, M. R.: Dielectric permittivity of eight gases measured with cross capacitors, Int. J. Thermophys., 24, 375–403, 10.1023/A:1022963720063, 2003.Schreiner, W., Weiss, J., Anthes, R., Braun, J., Chu, V., Fong, J., Hunt, D., Kuo, Y.-H., Meehan, T., Serafino, W., Sjoberg, J., Sokolovskiy, S., Talaat, E., Wee, T., and Zeng, Z.: COSMIC-2 Radio Occultation Constellation: First Results, Geophys. Res. Lett., 47, e2019GL086841, 10.1029/2019GL086841, 2020.Scribano, Y., Goldman, N., Saykally, R. J., and Leforestier, C.: Water dimers in the atmosphere III: Equilibrium constant from a flexible potential, J. Phys. Chem. A, 110, 5411–5419, 10.1021/jp056759k, 2006. Shostak, S. L., Ebenstein, W. L., and Muenter, J. S.: The dipole moment of water. I. Dipole moments and hyperfine properties of H2O and HDO in the ground and excited vibrational states, J. Chem. Phys., 94, 5875–5882, 10.1063/1.460471, 1991.Smith, E. K. and Weintraub, S.: The constants in the equation for the atmospheric refractive index at radio frequencies, P. IRE, 41, 1035–1037, 10.1109/JRPROC.1953.274297, 1953.Sokolovskiy, S.: Effect of superrefraction on inversions of radio occultation signals in the lower troposphere, Radio Sci., 38, 1058, 10.1029/2002RS002728, 2003.Steiner, A. K., Lackner, B. C., Ladstädter, F., Scherllin-Pirscher, B., Foelsche, U., and Kirchengast, G.: GPS radio occultation for climate monitoring and change detection, Radio Sci., 46, RS0D24, 10.1029/2010RS004614, 2011.Stratton, J. A.: Electromagnetic Theory, McGraw-Hill, New York, USA, 10.1002/9781119134640, 1941.Talpe, M. J., Nguyen, V. A., and Tomás, S.: Initial Polarimetric Radio Occultation Results from Spire's Nanosatellite Constellation: Satellite Payload, Collection, and Calibration, B. Am. Meteor. Soc., 106, E714–E734, 10.1175/BAMS-D-23-0314.1, 2025.Thayer, G. D.: An improved equation for the radio refractive index of air, Radio Sci., 9, 803–807, 10.1029/RS009i010p00803, 1974.Tohjima, Y., Machida, T., Watai, T., Akama, I., Amari, T., and Moriwaki, Y.: Preparation of gravimetric standards for measurements of atmospheric oxygen and reevaluation of atmospheric oxygen concentration, J. Geophys. Res.-Atmos., 110, D11302, 10.1029/2004JD005595, 2005.Xie, F., Wu, D. L., Ao, C. O., Kursinski, E. R., Mannucci, A. J., and Syndergaard, S.: Super-refraction effects on GPS radio occultation refractivity in marine boundary layers, Geophys. Res. Lett., 37, L11805, 10.1029/2010GL043299, 2010.Zou, X., Yang, S., and Ray, P. S.: Impacts of Ice Clouds on GPS Radio Occultation Measurements, J. Atmos. Sci., 69, 3670–3682, 10.1175/JAS-D-11-0199.1, 2012.