Search

The global positioning system (GPS) radio occultation (RO) technique has demonstrated the ability to precisely probe earth’s atmosphere globally with high vertical resolution. However, the lowermost troposphere still presents some challenges for the technique. Over moist marine areas, especially in subtropical regions, a very large negative moisture gradient often exists across the thermal inversion capping the marine boundary layer (MBL), which frequently causes superrefraction (SR), or ducting. In the presence of SR, the reconstruction of refractivity from RO data becomes an ill-posed inverse problem. This study shows that one given RO bending angle profile is consistent with a continuum (an infinite number) of refractivity profiles. The standard Abel retrieval gives the minimum refractivity solution of the continuum and thus produces the largest negative bias, consistent with a negative bias often present in the retrieved refractivity profiles in the moist lower troposphere. By applying a simple linear parameterization of the refractivity structure within and just below the SR layer, an analytical relation between the Abel-retrieved refractivity and a continuum of solutions is derived. Combining the Abel retrieval and the analytical relation with some physical constraints, a novel approach is developed to reconstruct the vertical refractivity structure within and below the SR layer. Numerical simulation studies in this paper have demonstrated the great potential of the reconstruction method to provide a much-improved retrieval in the presence of SR, and the method should greatly enhance the ability to measure the MBL structure globally using the GPS RO technique. [ABSTRACT FROM AUTHOR]

This paper describes the details of a fast, linear, forward-inverse refractive index mapping operator that can be used for assimilation of occultation data of various kinds into NWP models. Basically, the mapping consists of the integration of the refractive index along finite straight lines, mimicking the observational geometry as well as the subsequent retrieval of a refractive index profile, assuming spherical symmetry. Line integrals are discretized such that the refractivity is evaluated along the horizontal at fixed levels that can be chosen to coincide with the pressure levels of an NWP model. Integration of the hydrostatic equation at a large number of locations is thereby avoided. The mapping operator is tested using an idealized model of a weather front with large horizontal gradients. Mapped refractivity profiles are compared with retrieved refractivity profiles obtained via accurate 3D ray tracing simulations of GPS radio occultation events with ray path tangent points near the weather front. The simulations indicate that the mapping is a good representation of occultation measurements, including the influence large horizontal gradients have on retrieved refractivity profiles. To further the results, a simple ad hoc modification is introduced to approximately account for the ray path bending near the tangent points. The forward-inverse mapping allows for the near cancellation of otherwise crude approximations—for example, straight-line propagation—and the general concept could perhaps be adapted for the development of fast and accurate observation operators for the assimilation of other types of remote sensing data. [ABSTRACT FROM AUTHOR]