The Stochastic Time-Inverted Lagrangian Transport (STILT) model is comprised of a compiled Fortran executable that carries out advection and dispersion calculations as well as a higher-level code layer for simulation control and user interaction, written in the open-source data analysis language R. We introduce modifications to the STILT-R code base with the aim to improve the model's applicability to fine-scale (< 1 km) trace gas measurement studies. The changes facilitate placement of spatially distributed receptors and provide high-level methods for single- and multi-node parallelism. We present a kernel density estimator to calculate influence footprints and demonstrate improvements over prior methods. Vertical dilution in the hyper near field is calculated using the Lagrangian decorrelation timescale and vertical turbulence to approximate the effective mixing depth. This framework provides a central source repository to reduce code fragmentation among STILT user groups as well as a systematic, well-documented workflow for users. We apply the modified STILT-R to light-rail measurements in Salt Lake City, Utah, United States, and discuss how results from our analyses can inform future fine-scale measurement approaches and modeling efforts.
4 December 2018 in Geoscientific Model Development
Urban regions are responsible for emitting significant amounts of fossil fuel carbon dioxide (FFCO2), and emissions at the finer, city scales are more uncertain than those aggregated at the global scale. Carbon-observing satellites may provide independent top-down emission evaluations and compensate for the sparseness of surface CO2 observing networks in urban areas. Although some previous studies have attempted to derive urban CO2 signals from satellite column-averaged CO2 data (XCO2) using simple statistical measures, less work has been carried out to link upwind emission sources to downwind atmospheric columns using atmospheric models. In addition to Eulerian atmospheric models that have been customized for emission estimates over specific cities, the Lagrangian modeling approach – in particular, the Lagrangian particle dispersion model (LPDM) approach – has the potential to efficiently determine the sensitivity of downwind concentration changes to upwind sources. However, when applying LPDMs to interpret satellite XCO2, several issues have yet to be addressed, including quantifying uncertainties in urban XCO2 signals due to receptor configurations and errors in atmospheric transport and background XCO2. In this study, we present a modified version of the Stochastic Time-Inverted Lagrangian Transport (STILT) model, “X-STILT”, for extracting urban XCO2 signals from NASA's Orbiting Carbon Observatory 2 (OCO-2) XCO2 data. X-STILT incorporates satellite profiles and provides comprehensive uncertainty estimates of urban XCO2 enhancements on a per sounding basis. Several methods to initialize receptor/particle setups and determine background XCO2 are presented and discussed via sensitivity analyses and comparisons. To illustrate X-STILT's utilities and applications, we examined five OCO-2 overpasses over Riyadh, Saudi Arabia, during a 2-year time period and performed a simple scaling factor-based inverse analysis. As a result, the model is able to reproduce most observed XCO2 enhancements. Error estimates show that the 68 % confidence limit of XCO2 uncertainties due to transport (horizontal wind plus vertical mixing) and emission uncertainties contribute to ∼33 % and ∼20 % of the mean latitudinally integrated urban signals, respectively, over the five overpasses, using meteorological fields from the Global Data Assimilation System (GDAS). In addition, a sizeable mean difference of −0.55 ppm in background derived from a previous study employing simple statistics (regional daily median) leads to a ∼39 % higher mean observed urban signal and a larger posterior scaling factor. Based on our signal estimates and associated error impacts, we foresee X-STILT serving as a tool for interpreting column measurements, estimating urban enhancement signals, and carrying out inverse modeling to improve quantification of urban emissions.
A growing constellation of satellites is providing near‐global coverage of column‐averaged CO2 observations. Launched in 2019, NASA’s OCO‐3 instrument is set to provide XCO2 observations at a high spatial and temporal resolution for regional domains (100 × 100 km). The atmospheric column version of the Stochastic Time‐Inverted Lagrangian Transport (X‐STILT) model is an established method of determining the influence of upwind sources on column measurements of the atmosphere, providing a means of analysis for current OCO‐3 observations and future space‐based column‐observing missions. However, OCO‐3 is expected to provide hundreds of soundings per targeted observation, straining this already computationally intensive technique. This work proposes a novel scheme to be used with the X‐STILT model to generate upwind influence footprints with less computational expense. The method uses X‐STILT generated influence footprints from a key subset of OCO‐3 soundings. A nonlinear weighted averaging is applied to these footprints to construct additional footprints for the remaining soundings. The effects of subset selection, meteorological data, and topography are investigated for two test sites: Los Angeles, California, and Salt Lake City, Utah. The computational time required to model the source sensitivities for OCO‐3 interpretation was reduced by 62% and 78% with errors smaller than other previously acknowledged uncertainties in the modeling system (OCO‐3 retrieval error, atmospheric transport error, prior emissions error, etc.). Limitations and future applications for future CO2 missions are also discussed.
Observing air quality from sensors onboard light rail cars in Salt Lake County, Utah began as a pilot study in 2014 and has now evolved into a five-year, state-funded program. This metropolitan region suffers from both elevated ozone levels during summer and high PM2.5 events during winter. Pollution episodes result predominantly from local anthropogenic emissions but are also impacted by regional transport of dust, chemical precursors to ozone, and wildfire smoke, as well as being exacerbated by the topographical features surrounding the city. Two electric light-rail train cars from the Utah Transit Authority light-rail Transit Express (“TRAX”) system were outfitted with PM2.5 and ozone sensors to measure air quality at high spatial and temporal resolutions in this region. Pollutant concentration data underwent quality control procedures to determine whether the train motion affected the readings and how the sensors compared against regulatory sensors. Quality assurance results from data obtained over the past year show that TRAX Observation Project sensors are reliable, which corroborates earlier preliminary validation work. Three case studies from summer 2019 are presented to illustrate the strength of the finely-resolved air quality observations: (1) an elevated ozone event, (2) elevated particulate pollution resulting from 4th of July fireworks, and (3) elevated particle pollution during a winter time inversion event. The mobile observations were able to capture spatial gradients, as well as pollutant hotspots, during both of these episodes. Sensors have been recently added to a third light rail train car, which travels on a north–south oriented rail line, where air quality was unable to be monitored previously. The TRAX Observation Project is currently being used to provide reliable pollutant data for health studies and inform urban planning efforts. Links to real-time data displays and updated information on the quality-controlled data from this study are available on the webpage for the Department of Atmospheric Sciences at the University of Utah.
4 December 2017 in Journal of Geophysical Research: Atmospheres