key: cord-0864700-q0yv1o60 authors: Yadav, Vineet; Ghosh, Subhomoy; Mueller, Kimberly; Karion, Anna; Roest, Geoffrey; Gourdji, Sharon M.; Lopez‐Coto, Israel; Gurney, Kevin R.; Parazoo, Nicholas; Verhulst, Kristal R.; Kim, Jooil; Prinzivalli, Steve; Fain, Clayton; Nehrkorn, Thomas; Mountain, Marikate; Keeling, Ralph F.; Weiss, Ray F.; Duren, Riley; Miller, Charles E.; Whetstone, James title: The Impact of COVID‐19 on CO(2) Emissions in the Los Angeles and Washington DC/Baltimore Metropolitan Areas date: 2021-06-07 journal: Geophys Res Lett DOI: 10.1029/2021gl092744 sha: e5fbd3e0ff503899e361e03b33e2c7446f2e659c doc_id: 864700 cord_uid: q0yv1o60 Responses to COVID‐19 have resulted in unintended reductions of city‐scale carbon dioxide (CO(2)) emissions. Here, we detect and estimate decreases in CO(2) emissions in Los Angeles and Washington DC/Baltimore during March and April 2020. We present three lines of evidence using methods that have increasing model dependency, including an inverse model to estimate relative emissions changes in 2020 compared to 2018 and 2019. The March decrease (25%) in Washington DC/Baltimore is largely supported by a drop in natural gas consumption associated with a warm spring whereas the decrease in April (33%) correlates with changes in gasoline fuel sales. In contrast, only a fraction of the March (17%) and April (34%) reduction in Los Angeles is explained by traffic declines. Methods and measurements used herein highlight the advantages of atmospheric CO(2) observations for providing timely insights into rapidly changing emissions patterns that can empower cities to course‐correct CO(2) reduction activities efficiently. roles, yet each region has important differentiating characteristics related to its economy, infrastructure, climate, and demographics. As such, the impact of COVID-19 on CO 2 emissions is layered on top of other forces, for example, energy shifts or specific weather conditions, that drive emissions within these cities. We take a comprehensive approach by exploring the full value of CO 2 observations for understanding events like COVID-19 that drive urban-scale emissions. The CO 2 observations are used in three distinct analyses with increasing model dependency. The first detects the timing of emissions declines using atmospheric CO 2 observations, the second evaluates CO 2 enhancements to assess the persistence of reduced emissions, and the third employs an inverse model to estimate the relative change in emissions in 2020 compared to 2018 and 2019. In addition to detecting the timing and estimating the magnitude of emissions declines using atmospheric CO 2 observations, we also examine proxy activity data to tease apart the influence of COVID-19 from other factors that impact emissions. For this study, we consider Washington DC and Baltimore as a single metropolitan area (DC-Balt), given these cities' proximity and intersecting suburbs (US Census, 2010) . We also focus on the wider metropolitan extent of the Los Angeles area (LA). While topographically very different, both regions are similar in size (∼18,000 km 2 ) and are densely populated ( Figure S3 ). The two regions experience very different climates: LA's is characterized by seasonal changes in rainfall with a dry summer and a rainy winter, while DC-Balt experiences hot, humid summers and cold, wet winters. These climatological differences impact emissions from heating and cooling as well as urban and surrounding vegetation, which controls CO 2 uptake and release from the biosphere. The emissions magnitude and relative contribution of different sectors to emissions in each city reflect the economies in these metropolitan centers. Two existing 1 km 2 emissions data products give insight into sectoral distributions. Here we refer to the Hestia emissions product, tailored for LA (Gurney et al., 2019) , and Vulcan 3.0 , a national emissions product, in DC-Balt. The two emission products are constructed using similar methods and are consistent, such that the domain total for LA in the Vulcan national product is the same as in Hestia. According to these products, the total annual magnitude of in-boundary CO 2 emissions in LA is 43.7 megatons of carbon (MtC)/year, which is over twice that of DC-Balt, 20.5 MtC/year ( Figure S1 ). Although each city has its own sectoral emission mix, the on-road sector dominates in both LA and DC-Balt, with on-road sector emissions accounting for ∼45% of annual totals in both cities. The industrial sector is the second largest component in LA, accounting for an additional 23% of the yearly total. The electricity production sector in DC-Balt is the second most important contributor, accounting for 19% of annual emissions. The sectoral breakdown of fossil-fuel CO 2 (FFCO 2 ) emissions in DC-Balt exhibit much larger seasonal and annual variability compared to LA due to regional climate and energy demand ( Figure S1 ). The changes in emissions are evident during the transition from winter to spring, corresponding to our study period of January through May; DC-Balt emissions show wide ranges for the commercial (∼23%-4%) and residential (∼20%-4%) sector fractions due to reduced demand for heating in the warmer months compared with January and February ( Figure S2 ). Both the LA and DC-Balt metropolitan areas are being monitored for ambient CO 2 levels using networks of high-accuracy measurement stations. Each network has ∼9-13 sites with sensitivity to urban anthropogenic emissions ( Figure S3 ; Tables S1 and S2), although the number of sites and measurements available in any given month varies, depending on how each network was developed and maintained (Karion et al., 2020; Kort et al., 2013; Lopez-Coto et al., 2017; Mueller et al., 2018; Verhulst et al., 2017) . We use 2018 and 2019 as baseline years with which to compare 2020 observations and inferred emissions because we have a consistent number of observational sites for each domain during these times. The network sensitivity, or footprint, of the urban CO 2 measurements is also generally consistent in these 3 years ( Figure S4 ). Although we do not have a long-term record as a basis for comparison, using the average of two baseline years allows us to account for some year-to-year variability in emissions. Atmospheric CO 2 observations from the two urban networks contain signatures from biospheric and anthropogenic CO 2 sources and sinks within the footprint, which can include fluxes from both inside and outside the urban domain. We define FFCO 2 as all fluxes that result from combustion. As such, isolating the FFCO 2 contribution to CO 2 observations from each urban domain requires an understanding of biospheric fluxes and incoming CO 2 (i.e., background). Discerning FFCO 2 -driven enhancements in LA is more straightforward than in DC-Balt because, in part, the biosphere has less of an impact (Miller et al., 2020) . We use a non-year specific biospheric model, the Vegetation Photosynthesis Respiration Model (VPRM for July 2017 to June 2018 as described in Section S5.1.1; Figure S5a ), to model the biological fluxes. Modeled biospheric enhancements are consistent at the city scale with those presented in Miller et al. (2020) from January to May. We use measurements from the site on San Clemente Island (SCI) to represent incoming airflow (Verhulst et al., 2017) . Median monthly FFCO 2 enhancements estimated from atmospheric CO 2 observations within LA range between 8 µmol mol −1 and 30 µmol mol −1 (Section S6; Figure S6a ). Representing CO 2 concentration of incoming airflow for DC-Balt is more difficult because the domain is downwind of large agricultural areas, deciduous forests, and emission sources, including other metropolitan areas and power plants. We represent the incoming CO 2 based on observations from one of three near-boundary towers (BUC, SFD, or TMD, Figure S3 ) selected using modeled back-trajectories (described in Section S3.1) similar to a recent study in Boston (Karion et al., 2021; Sargent et al., 2018) , and described further in Section S4. In DC-Balt, we represent biospheric fluxes using year-specific (2018-2020) VPRM to ensure properly characterized variability in the modeled fluxes (Section S5.1.2). Average daily observed drawdown due to biological activity can be as high as 4.7 µmol mol −1 (as predicted by VPRM; Figure S5b ) during our study period. After accounting for inflow and the contribution of biological fluxes to observed CO 2 , the estimated median monthly FFCO 2 enhancements for DC-Balt are ∼1.5 µmol mol −1 -7 µmol mol −1 in January through May between 2018 and 2020, 3-4 times lower than those in LA ( Figure S6 ). We use three different types of analysis to detect the onset of the COVID-19 effect on CO 2 emissions and estimate the relative change of emissions. To simplify the interpretation of results, we keep most of the modeling and analysis components consistent between LA and DC-Balt. However, there are a few differences that cannot be avoided due to city-specific meteorology and availability of data. Our period of assessment includes the months of January-May for 2018, 2019, and 2020, with 2018 and 2019 observations and inferred emissions providing baselines for comparison with 2020. For the first analysis, we focus on the temporal variability of measured atmospheric CO 2 rather than the concentrations themselves to detect the timing of emissions changes. CO 2 enhancements contain large variability associated with weather conditions, making it difficult to isolate changes, while high-frequency (sub-hourly scale) variability is expected to directly correlate with emissions near the observation locations and are less affected by synoptic conditions (Umezawa et al., 2020) . To identify changes in variability in 2020 compared to previous years, we use the ratio of the standard deviations of the 1 min averages in 2020 to the averaged standard deviations from 2018 and 2019 in a change point analysis (Aminikhanghahi & Cook, 2017 ) (Section S6.1) and estimate dates associated with sudden emissions decline and recovery. In the second analysis, we examine the atmospheric CO 2 observations by assessing the mean daily cumulative afternoon (i.e., the time period when the atmosphere is well-mixed) FFCO 2 enhancement, comparing 2020 to other years. We calculate FFCO 2 enhancements ( FFCO 2 y ) as follows: where obs y are the CO 2 mole fractions in µmol CO 2 per mol of dry air, or ppm, as observed at the urban measurement sites; back y is the background CO 2 (Section S4); H, expressed in (µmol mol −1 )/(µmol m −2 s −1 ), represents sensitivities of the observations to surface fluxes as obtained from the atmospheric transport and dispersion model (Section S3.1); and bio y are the biological contributions to obs y as modeled through the convolution of H and the modeled VPRM fluxes ( VPRM s ; Section S5) with units µmol m −2 s −1 . Additional detail on the cumulative summation method is given in Section S6.2. In the final analysis, we used a geostatistical-Bayesian atmospheric inverse modeling approach with FFCO 2 y as observations to estimate posterior FFCO 2 emissions. We use Hestia (LA) and a temporally smoothed Vulcan 3.0 (DC-Balt) emission product for 2015, the most recent year available for both products, as covariates (priors) in the two inversions. Both priors are re-gridded to 0.02° × 0.02° for the domains in Figure S3 . To estimate the gridded posterior emissions, we adopt a moving window approach to estimate FFCO 2 emissions at daily temporal intervals by maintaining a one-day overlap between two 2-day inversions (for details see Section S7). The posterior FFCO 2 emissions are used to detect a relative change in emissions in March and April of 2020 in comparison to the mean of the baseline years, 2018 and 2019. The posteriors are further compared against adjusted domain total FFCO 2 emissions from Hestia and Vulcan 3.0 computed using information on statewide fuel sales, continuous emissions monitoring systems (CEMS) data from powerplants, natural gas consumption, and vehicle miles traveled (VMT) data (Section S8). The influence of the pandemic response on FFCO 2 emissions is evident in both the variability of concentrations (as indicated by their standard deviation) and the cumulative FFCO 2 enhancements in LA and DC-Balt ( Figure 1 ). In the observation analysis, we observe a reduction in the standard deviation of CO 2 within each hour, averaged daily over all sites, in mid-March 2020 relative to 2018 and 2019 in both metropolitan areas suggesting that emissions declined substantially at that time (Figures 1a and 1b) . The reduction in CO 2 variability is also apparent in the statistically significant smaller interquartile range of hourly averaged standard deviations after these dates, matching the timing of behavioral shifts indicated in Apple mobility data (Apple, 2020) . The corresponding uncertainties as indicated by the 95% credible intervals (Section S6.1) are 6 days (LA) and 21 days (DC-Balt). We also detect a second change point for LA on April 29 with a 95% credible interval of 9 days, associated with an increase in traffic (Figures 1a and S11 ) (California Dept. of Transportation, 2020) . Uncertainty bounds in LA are narrower due to the sharper relative drop in the standard deviation. We did not detect a second change point in the DC-Balt data, perhaps due to the ramp-up of the biosphere during this time period. The transition in LA is also evident when comparing the FFCO 2 cumulative enhancements before and after the decrease in mobility (Figure 1c) . Large interannual variability of the cumulative enhancements in LA is likely due to climatological or other extreme atmospheric events like a fire in early 2018. However, the reduced FFCO 2 cumulative enhancements after March 13, 2020 relative to 2018 and 2019 are apparent. Evaluating LA cumulative enhancements over additional years (2015-2020; Figure S7 ), we find that 2020 enhancements prior to March 13 are well within interannual variability, further indicating that the low post-March 13 enhancements in 2020 are likely caused by an emission change rather than year-specific meteorological variability. In DC-Balt, the 2020 cumulative enhancement is slightly higher in the beginning of the year than in 2018 and 2019 (Figure 1d ). After March 13, 2020 cumulative enhancements are lower than the previous 2 years, again indicating a change in emissions although less apparent than for LA. There may be several reasons why the FFCO 2 y decrease in 2020 is not larger, including lower overall enhancements in DC-Balt compared to LA (leading to smaller signals) or interannual differences in transport, which affect the relationship between enhancements and emissions. However, despite the difference in magnitude between the cities, at the onset of behavioral shifts (∼March 13, dashed vertical line in Figures 1c and 1d) , the enhancement growth for 2020 is lower suggesting a decrease in emissions-a conclusion consistent with the results of the LA FFCO 2 enhancement analysis. Posterior FFCO 2 emissions estimates from the inversion analysis also reflect the emissions decrease associated with the COVID-19 pandemic in both LA and DC-Balt (blue bars, Figure 2 ). In both March and April 2020, the posterior reflects significantly lower emissions than the inversion prior, and the April 2020 posterior is significantly lower than the posteriors for the pre-pandemic months in 2020 at 95% and 68% confidence intervals (CI) in LA ( Figure 2c ) and DC-Balt (Figure 2f ), respectively. April 2020 posteriors are also lower than any other January through March period in the previous 2 years in both cities, suggesting that the decrease in this month is outside of seasonal and interannual variability. The posterior in March 2020 does not reflect as large a decrease as April, given that the reduction in activity due to the pandemic response occurred in mid-March in both cities. In May 2020, the inversion FFCO 2 estimates are larger than those of April after accounting for uncertainties (Figures 2c and 2f) . In LA, the May increase in the posterior emissions is likely due to the rebound in traffic given the second change point identified in Figure 1a . Alternatively, in DC-Balt, the standard deviation analysis did not indicate any clear change point associated with a recovery period possibly due to the large influence of the biosphere on atmospheric CO 2 in May (Figure 1d ). In May, the biological contribution is nearly equal the For scaling, we use publicly available proxy data, such as reported emissions from electricity generation, natural gas usage, and fuel sales or VMT or vehicle hours traveled (VHT) (Section S8; Table S3 ). In DC-Balt and LA, all sectors contributing more than 10% of emissions in January through May were adjusted if public data was available. These year-specific estimates are not meant to provide accurate new emission products for these later years; rather, they provide a basis with which to compare our result for the whole-city domain (Figures 2a-2f , dark gray bars). Overall, we find that the estimated FFCO 2 posteriors are often (but not always) more consistent with year-specific activity-based BU totals than with the 2015 prior. The general good comparison gives confidence in our posterior estimates and helps explain variability in the emission estimates between years. This improved consistency indicates that the CO 2 observations within the inversion are driving the estimated emissions to better represent the true emissions for a particular year or month. For example, we expect more variability in posterior FFCO 2 emissions in DC-Balt than LA, largely due to climate. Indeed, in LA, the 2020 activity-based estimates are largely consistent with the 2015 prior for all months and years except for March, April, and May of 2020 (Figure 2, dark gray bars) . The monthly variation in DC-Balt activity-based BU estimates is large from year to year because natural gas consumption to heat residential and commercial buildings largely depends on weather conditions in any given winter or spring. In addition, emissions from electricity production vary significantly between months and years in this region (Section S8). For example, the February 2019 activity-based BU emissions in DC-Balt are 27% lower than the 2015 prior and agree well with the posterior for that month (Figure 2e ). We note that the posterior LA emissions in January and February 2018 do not agree well with the activity-based BU emissions for these months (Figure 2a) . The large posterior emissions in January 2018 may be caused by an active fire season in the region that stretched into the early days of 2018 (fire emissions are not included in the prior or activity-based BU estimates but would contribute to the posterior emissions as their signatures would be in obs y ). In February 2018, we suspect that poorly modeled meteorological YADAV ET AL. 10.1029/2021GL092744 6 of 10 conditions resulted in elevated emissions estimates, given the prevalence of low wind speeds in that month (Section S3.2). We compare March and April 2020 emissions to their respective means for 2018 and 2019 ( Figure 3 ) to estimate the relative decrease in emissions in these months. We also compare the March and April activity-based BU totals for 2020 to those of previous years to discern whether the observed reductions are consistent with changes in activity data. In LA, in addition to VMT, we also calculate BU reductions using VHT data as a second point of comparison for the estimated relative reduction (Section S8.4). VHT decreased more than VMT during the pandemic because of reduced traffic congestion resulting in higher average vehicular speeds ( Figure S11 ). We note that gasoline fuel sales, which is the best proxy data for on-road emissions, is only available at the state level, which may not reflect on-road activity in LA. We use both VMT and VHT in the analysis because they reflect different but related contributions to on-road emissions (distance traveled as well as time spent on the road). In LA, the posterior reductions are 17% ± 9% (0.57 MtC ± 0.30 MtC, 95% CI) in March and 34% ± 6% (1.09 MtC ± 0.21 MtC, 95% CI) in April (blue bars, Figure 3) . We also estimate a 28% ± 4% CO 2 emissions reduction for May in LA (not shown). The LA activity-based BU total for March and April suggests CO 2 emission reductions of 8% and 17%, respectively, all attributed to changes in VMT (light gray, Figure 3a ). Using VHT, the reductions are slightly larger at 11% and 21%, which are smaller than the posterior estimates; these are within our estimated 95% CI in March but not in April (medium gray, Figure 3a) . Overall, VMT/VHT in the LA area decreased by 17%/23% in March and 33%/42% in April 2020 compared to the mean of 2018 and 2019. Although VHT indicates a larger decrease than VMT ( Figure S11 ), VMT has been previously used to estimate emissions from on-road vehicles (Gately et al., 2015; . The industrial sector is the only other significant portion of LA's total emissions, at 23%. There are some indications that this sector was impacted by COVID-19 or other external forces (e.g., significant labor drops in employment for the manufacturing sector) (Bureau of Labor Statistics, 2020; California Energy Commission, 2020; Hydrocarbon Processing, 2020). However, without a reliable proxy for the industrial sector, we are unable to account for any quantitative changes. In DC-Balt, we estimate a slightly larger relative drop in the posterior emissions than in LA in March, that is, 25% ± 14% (0.45 MtC ± 0.25 MtC, 95% CI), while the April decrease is similar to LA, at 33% ± 11% (0.43 MtC ± 0.15 MtC, 95% CI) (blue bars, Figure 3b ). The estimated decrease for both months is broadly YADAV ET AL. 10.1029/2021GL092744 7 of 10 consistent with (although slightly larger than) the reductions predicted using the activity-based BU totals, which are 21% and 27% in March and April, respectively. However, only a 6% decrease in March is attributable to a reduction in the on-road sector relative to previous years, since only half of this month was affected by the slowdown in activity (medium gray bars, Figure 3b ). The remainder of the 21% March BU decrease is largely attributed to reduced natural gas use in the residential and commercial sectors due to decreased demand for heating, as March 2020 was warmer than March 2018 and 2019. Natural gas use in DC-Balt was not anomalous in 2020 after accounting for regional temperatures and follows the trend with heating degree days from previous years ( Figure S10 ), suggesting that behavioral changes due to COVID-19 were not responsible for the decline. In April 2020, the activity data suggest that most of the observed emissions reduction is indeed from on-road emissions, which are responsible for an estimated 23% decrease. The response to COVID-19 influenced CO 2 emissions in both Los Angeles and Washington DC/Baltimore, resulting in declines of over 30% in April relative to previous years. These reductions are consistent with a 30% emissions decline found in the San Francisco Bay area between the 6 weeks prior and the 6 weeks after that city's lockdown (Turner et al., 2020) . Our study demonstrates the detection of changes in urban-scale CO 2 emissions as small as 15%-30% during this time period. Overall, the atmospheric inversion-derived CO 2 emissions displayed significant drops in both cities due to factors (e.g., decreased traffic) impacted by COV-ID-19. Our three observation-based analyses allowed us to detect the timing of emissions declines, assess the continuation of reduced emissions, and estimate the relative change in emissions in March and April 2020 in comparison to 2018 and 2019. Beyond the reductions of CO 2 emissions due to COVID-19 impacts, we also found other factors that significantly influence monthly emissions using trends in publicly available activity data. In DC-Balt, both COVID-related traffic reduction and warmer weather caused the observed emissions decline; in LA, only a portion of the decline is attributable to a reduction in traffic, while some of the emissions declines remains unexplained. The methods used here show that emission changes (whether due to pandemic effects, as shown here, or any other cause) that may be missed with solely an activity-based approach can be detected using atmospheric measurements. Such timely emissions information allows cities to assess the effectiveness of policies designed to reduce emissions. Furthermore, accurate CO 2 emissions information provided at shorter timescales could be crucial to policy makers to course-correct their mitigation actions if necessary and propose new measures (e.g., redesigning transportation infrastructure) that will ultimately drive down CO 2 emissions. All data used in this analysis are available at data.nist.gov (https://doi.org/10.18434/mds2-2343). was partially funded by NIST grants 70NANB15H344, 70NANB14H322, 70NANB19H129, 70NANB19H130, 70NANB19H132, along with a commercial contract with EN (1333ND-19PNB600853). A portion of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). A survey of methods for time series change point detection Occupational employment and Wages in Los Angeles-Long Beach-Anaheim, state and area employment, hours, and earnings Caltrans performance measurement system (PeMS) Cities, traffic, and CO 2 : A multidecadal assessment of trends, drivers, and scaling relationships The Vulcan version 3.0 high-resolution fossil fuel CO 2 emissions for the United States Marathon Los Angeles Refinery cuts production due to demand drop Greenhouse gas observations from the Northeast Corridor tower network Background conditions for an urban greenhouse gas network in the Washington, D.C. and Baltimore metropolitan region Surface observations for monitoring urban fossil fuel CO 2 emissions: Minimum site location requirements for the Los Angeles megacity Time series analysis and its applications Temporary reduction in daily global CO 2 emissions during the COVID-19 forced confinement Near-real-time monitoring of global CO 2 emissions reveals the effects of the COVID-19 pandemic Tower-based greenhouse gas measurement network design-The national institute of standards and technology north east Corridor testbed Large and seasonally varying biospheric CO 2 fluxes in the Los Angeles megacity revealed by atmospheric radiocarbon Siting background towers to characterize incoming air for Uuban Greenhouse gas estimation A case study in the Anthropogenic and biogenic CO 2 fluxes in the Boston urban region Observed impacts of COVID-19 on urban CO 2 Emissions Statistical characterization of urban CO 2 emission signals observed by commercial airliner measurements Carbon dioxide and methane measurements from the Los Angeles Megacity Carbon Project -Part 1: Calibration, urban enhancements, and uncertainty estimates Satellite-based estimates of decline and rebound in China's CO 2 emissions during COVID-19 pandemic References From the Supporting Information Meteorological model evaluation for CalNex Assessment of an atmospheric transport model for annual inverse estimates of California greenhouse gas emissions A north American hourly assimilation and model forecast cycle: The rapid refresh Error statistics of Bayesian CO 2 flux inversion schemes as seen from GOSAT Climate warming impact on degree-days and building energy demand in Switzerland. Energy Conversion and Management A simplified framework for high-resolution urban vegetation classification with optical imagery in the Building energy demand assessment through heating degree days: The importance of a climatic dataset Global monthly CO 2 flux inversion with a focus over Heating and cooling building energy demand evaluation; a simplified model and a modified degree days approach Los Angeles megacity: A high-resolution land-atmosphere modeling system for urban CO 2 emissions Vulcan: High-resolution annual fossil fuel CO 2 emissions in USA The Hestia fossil fuel CO 2 emissions data product for the Los Angeles megacity (Hestia-LA) Hestia fossil fuel carbon dioxide emissions inventory for Los Angeles basin Intercomparison of atmospheric trace gas dispersion models: Barnett Shale case study Quasi-Linear geostatistical theory for inversing A city scale degree-day method to assess building space heating energy demands in Strasbourg Eurometropolis (France) Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model mcp: An R package for regression with multiple change points A satellite-based biosphere parameterization for net ecosystem CO 2 exchange: Vegetation Photosynthesis and Respiration Model (VPRM) Maximum likelihood estimation of covariance parameters for Bayesian atmospheric trace gas surface flux inversions CO 2 transport, variability, and budget over the Southern California air basin using the high-resolution WRF-VPRM model during the CalNex 2010 campaign Seasonal climatology of CO 2 across north America from aircraft measurements in the NOAA/ESRL Global Greenhouse Gas Reference Network Inverse Problem Theory and Methods for Model Parameter Estimation Monthly, state-level "Natural gas consumption by end use" for DC Monthly, state-level "Natural gas consumption by end use" for Maryland Monthly, state-level "Natural gas consumption by end use" for Virginia Monthly, state-level "Prime supplier sales volumes" for "motor gasoline" for DC Monthly, state-level "Prime supplier sales volumes" for "motor gasoline" for Maryland Monthly, state-level "Prime supplier sales volumes" for "motor gasoline" for Virginia Detecting urban emissions changes and events with a near-real-time-capable inversion system Spatio-temporally resolved methane fluxes from the Los Angeles Megacity On the relation between urban climate and energy performance of buildings. A three-years experience in Rome The authors thank Hratch Semerjian and Eleanor Waxman (NIST), Michael Stock, Elizabeth DiGangi, Bryan Biggs, and Charlie Draper (EN) for their contributions. This work