Microwave Limb Sounder
Lagrangian Trajectory Diagnostics
>Download this as a PDF
IntroductionThe Aura Microwave Limb Sounder (MLS) Lagrangian Trajectory Diagnostics (LTDs) are support products for MLS, developed under the NASA Aura Science Team program. More details on MLS are given in [Waters et al., 2006]. Details of the MLS Level 2 data that the LTDs support are given in [Livesey et al., 2011]. The diagnostics describe the motion of the airmasses observed by MLS from 15-days prior to 15-days following each MLS observation. The LTDs are provided to aid scientific studies using MLS observations of atmospheric properties, and are generated for MLS observations from the upper troposphere to the lower mesosphere.
The main files (known as "combined" LTDs) contain trajectories running from 15 days prior to 15 days following a "launch day", being a single day of MLS observations. Note that these files represent a "merge" of separate forward and reverse Lagrangian trajectory computations. Accordingly, in these files, the airmasses start at some "arbitrary" location 15 days prior to the launch day, pass through a region observed by MLS at the time of that observation on the launch day then travel on a further 15 days. For user convenience, two smaller subsets describe only the launch day and the 5 days either following (5daysfwd files) or prior (5daysrev).
The LTDs describe diabatic parcel trajectories computed using NASA Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System (GEOS) model version 5.2 winds and temperature tendencies [Reinecker et al., 2008]. This limits the temporal span of the LTDs to the period from the July 2004 Aura launch to April 30, 2013. Planned updates to the LTD products will use the new GEOS 5.9 record, enabling extension to more recent dates. Parcel advection is computed using a fourth-order Runge Kutta method with a five-minute timestep, based on the algorithm described in Manney et al. [1994b].
The specific LTDs described herein are a user-friendly "reduced" subset of a set of "standard" LTD products generated at the MLS Science Computing Facility (SCF). Differences between the reduced and full sets (mainly decreased temporal resolution) are described below. Some of the quantities within the LTD files, although only relevant for the full files, are retained in the reduced files in order to maintain compatibility and interoperability for user software. A description of all the LTD variables, including these less-relevant ones, is included in this document for completeness.
LTD files and their contentsThe LTD file naming convention, is best illustrated by the following example:
MLS-Aura_L2LTD-G5D-RedV1p2a_v03-30-c01-ltdv01p0_2005d001-combined.ncIn compliance with the Aura-wide file naming convention [Craig et al., 2003] the name consists of four main fields -- instrument, product, version, date -- separated by underscores, with a file type (nc, indicating NetCDF) following the only period. "Words" within the individual fields are delineated by dashes. The first field, MLS-Aura, indicates that this is an Aura MLS product.
The second field, L2LTD-GSD-RedV1p2a denotes the product type. L2LTD denotes that these are Level 2 Lagrangian Trajectory Diagnostics. G5D indicates that the file contains diabatic trajectories (D), computed using the GEOS 5.2 wind fields (G5). The final word describes the particular trajectory run, with RedV1p2a being the current publicly-available run. (The Red indicates that these are "reduced" trajectories, as described above; the original "full" trajectories are called StdV1p2a).
The third, "version", field mirrors that of the MLS Level 2 data from which the trajectories were launched (e.g., v03-30-c01), with additional information on the version of the software used to generate the trajectories (ltdv01p0). The final, "date", field gives the Aura format date for the launch day (2005d001 indicating January 1, 2005), with -combined indicating that the file contains a combination of forward and reverse trajectories from 15 days prior to 15 days subsequent to the observations. Alternative smaller files contain only 5-day forward (5daysfwd) or reverse (5daysrev) trajectory information (though note that the 5daysfwd file will contain reverse trajectories back to midnight universal time on the launch day, and 5daysrev forward trajectories to midnight the following day.
The LTD files are stored in the widely used NetCDF format. Each file contains two main sets of information: details of the trajectory launch locations, and information on the motion of each air parcel prior to and following its launch. The table below details each quantity in the LTD file and its meaning. The launch information describes, for each trajectory, the latitude, longitude, time and height (actually - log10[pressure / hPa]) of its launch, along with arrays relating each launch to a specific MLS profile in the parent MLS Level 2 (L2GP) file.
The trajectory information is simply a set of arrays describing the location (latitude, longitude and potential temperature, stored as log10[θ / K]) of each air parcel at each timestep. In addition to location and local atmospheric temperature, additional information can be stored for each trajectory in fields prefixed with extra_. Currently only potential vorticity (SPV, scaled as in Manney et al. [1994a]), equivalent latitude [Butchart and Remsberg, 1986] derived from the SPV, and a diagnostic of trajectory dispersion are so stored; others may be added later. Units for all fields are described in the supporting metadata. All latitudes, longitudes and equivalent latitudes are in degrees, temperatures are in Kelvin, and distances are in meters. Times are given in the EOS-standard "TAI93" format, being the number of seconds since midnight universal time on January 1, 1993 (including the subsequent leap seconds). Other units are given in the table below.
The large arrays associated with the trajectory information are stored "packed" (i.e., as scaled and offsetted two-byte integers). The required unpacking is handled transparently to the user by some NetCDF reading software (e.g., Matlab) but not others (most notably IDL where the user has to apply the offset and scaling themselves). The quantization imposed by this storage equates to ~500 m horizontal and ~1.5 m vertical accuracy.
Obtaining the LTD filesThe LTD files are available from the MLS Science Computing Facility (SCF) at:
People are encouraged to register as an LTD user, in order to be kept informed of updates and other issues. The registration page is at http://mls.jpl.nasa.gov/forms/regltd.php. In addition, the MLS team would appreciate being informed (via an e-mail to email@example.com) of any publications using the LTD products.
Notes on "reduced" vs. "standard" trajectories and other issuesThese "reduced" LTD products described here require ~250 Gb of storage per year of MLS observations for the combined set, with ~50GB per year for each of the 5dayfwd and 5dayrev sets. These are all derived from a significantly more voluminous record of "standard" LTDs that require 2.1 Tb per year. Differences between standard and reduced LTDs are twofold. Firstly, standard LTDs record the parcel location every 20 minutes rather than every two hours. Secondly, standard LTDs actually describe two sets of trajectories. One set is launched from along the MLS track, as for the reduced LTDs, the other set is launched 20 km to the right of the MLS track (i.e., to the east for ascending parts of the orbit and to the west for descending). The divergence of each parcel pair is one measure of the reliability of the trajectory calculation. The reduced LTDs retain the locations of only the "along track" parcels, but also contain an additional field (extra_flankseparation) that records the distance (in meters) between the two "parallel" trajectories at each output time.
Although the launchTime field denotes the TAI93-format time of each launch (i.e., each MLS profile observation), the actual effective trajectory launch times are rounded to the nearest 20-minute timestep. Note that, unless that timestep coincides with one of the two-hour timesteps in the reduced files, the trajectory will not be "seen" to land directly on the MLS observation location (as recorded in launchLat and launchLon). In any case, even for those that are launched 10 minutes either side of a two-hour timestep, the quantization introduced by the two-byte integer storage means that lat and lon will likely not coincide exactly with launchLat and launchLon.
Contents of the LTD files
† Denotes arrays stored as packed (scaled and offset) 2-byte integers.
Work at the Jet Propulsion Laboratory, California Institute of Technology was performed under contract with the National Aeronautics and Space Administration.
Butchart, N., and E. E. Remsberg, The area of the stratospheric vortex as a diagnostic for tracer transport on an isentropic surface, J. Atmos. Sci., 43, 1319–1339, 1986.
Craig, C., K. Stone, D. Cuddy, S. Lewicki, P. Veefkind, P. Leonard, A. Fleig, and P. Wagner, HDF-EOS Aura file format guidelines, Tech. rep., National Center For Atmospheric Research, 2003.
Livesey, N. J., et al., EOS MLS version 3.3 Level 2 data quality and description document, Tech. rep., Jet Propul- sion Laboratory, available from http://mls.jpl.nasa.gov/, 2011.
Manney, G. L., R. W. Zurek, M. E. Gelman, A. J. Miller, and R. Nagatani, The anomalous Arctic lower strato- spheric polar vortex of 1992–1993, Geophys. Res. Lett., 21(22), 2405–2408, 1994a.
Manney, G. L., R. W. Zurek, A. O’Neill, and R. Swinbank, On the motion of air through the stratospheric polar vortex, J. Atmos. Sci., 51, 2973–2994, 1994b.
Reinecker, M. M., et al., The GEOS-5 data assimilation system: A documentation of GEOS-5.0, Tech. Rep. 104606 V27, NASA, 2008.
Waters, J. W., et al., The Earth Observing System Microwave Limb Sounder (EOS MLS) on the Aura satellite, IEEE Trans. Geosci. Remote Sens., 44(5), 1075–1092, 2006.