Orbiting Carbon Observatory-2 Observations Research Articles

Ocean-driven interannual variability in atmospheric CO2 quantified using OCO-2 observations and atmospheric transport simulations

Interannual variability (IAV) in the atmospheric CO2 growth rate is caused by variation in the balance between uptake by land and ocean and accumulation of anthropogenic emissions in the atmosphere. While variations in terrestrial fluxes are thought to drive most of the observed atmospheric CO2 IAV, the ability to characterize ocean impacts has been limited by the fact that most sites in the surface CO2 monitoring network are located on coasts or islands or within the continental interior. NASA’s Orbiting Carbon-Observatory 2 (OCO-2) mission has observed the atmospheric total column carbon dioxide mole fraction (XCO2) from space since September 2014. With a near-global coverage, this dataset provides a first opportunity to directly observe IAV in atmospheric CO2 over remote ocean regions. We assess the impact of ocean flux IAV on the OCO-2 record using atmospheric transport simulations with underlying gridded air-sea CO2 fluxes from observation-based products. We use three observation-based products to bracket the likely range of ocean air-sea flux contributions to XCO2 variability (over both land and ocean) within the GEOS-Chem atmospheric transport model. We find that the magnitude of XCO2 IAV generated by the whole ocean is between 0.08-0.12 ppm throughout the world. Depending on location and flux product, between 20-80% of the IAV in the simulations is caused by IAV in air-sea CO2 fluxes, with the remainder due to IAV in atmospheric winds, which modulate the atmospheric gradients that arise from climatological ocean fluxes. The Southern Hemisphere mid-latitudes and low-latitudes are the dominant ocean regions in generating the XCO2 IAV globally. The simulation results based on all three flux products show that even within the Northern Hemisphere atmosphere, Southern Hemisphere ocean fluxes are the dominant source of variability in XCO2. Nevertheless, the small magnitude of the air-sea flux impacts on XCO2 presents a substantial challenge for detection of ocean-driven IAV from OCO-2. Although the IAV amplitude arising from ocean fluxes and transport is 20 to 50% of the total observed XCO2 IAV amplitude of 0.4 to 1.6 ppm in the Southern Hemisphere and the tropics, ocean-driven IAV represents only 10% of the observed amplitude in the Northern Hemisphere. We find that for all three products, the simulated ocean-driven XCO2 IAV is weakly anti-correlated with OCO-2 observations, although these correlations are not statistically significant (p>0.05), suggesting that even over ocean basins, terrestrial IAV obscures the ocean signal.

Exploring CO2 anomalies in Brazilian biomes combining OCO-2 & 3 data: Linkages to wildfires patterns

Climate change is a challenge to the global community and one of the causes is the increase of greenhouse gases (GHG) in the atmosphere, especially carbon dioxide (CO2). The major emission source of this gas into the atmosphere comes from the burning of fossil fuels and biomass burning, on the other hand, the main sink comes from biochemical processes such as photosynthesis. Thus, the observation of CO2 is a key point to understanding sources and sinks. In this context, The Orbiting Carbon Observatory 2 (OCO-2) and 3 (OCO-3), are a NASA dedicated mission to monitor the column-averaged dry-air mole fraction of carbon dioxide (XCO2) on a global scale. We combined the OCO-2 and OCO-3 observations to study the spatial distribution of XCO2 anomalies and how some of these anomalies are related to fire occurrence in the Brazilian Biomes during 2020 and 2021 considering two different seasons, Dry and Wet. The fire occurrence was obtained from Fire Information for Resource Management System (FIRMS) that provides the data from the Moderate-Resolution Imaging Spectroradiometer (MODIS) product of active fires and thermal anomalies at Near-Real Time (MCD14DL, collection 6). The OCO-2/3 observations are affected by cloud formations in wet seasons, we observe that the dry period has more observations. The XCO2 anomaly values range from ∼ 7.0 ppm to −7.0 ppm and mostly positive anomalies occur in Amazon Biome, and this ecosystem has higher average values for all periods (∼0.9 ppm), compared to the other biomes. The fire occurrence was higher in dry periods, especially in 2020 when unprecedented fire outbreaks were registered in Brazil. The most affected biomes were Pantanal, Cerrado, and Amazon. XCO2 positive anomalies spatially agree with fire foci over some areas, and the correlation values between them ranged from 0.2 to 0.5 depending on the biome and season, and when considering observations with clouds the correlation is slightly higher. We point out for the first time the possibility of using OCO-2 and 3 combined, also, how positive XCO2 anomalies are related to fire occurrence in different ecosystems and periods, and the role of cloud detection in this relationship.

Comparison of Atmospheric Carbon Dioxide Concentrations Based on GOSAT, OCO-2 Observations and Ground-Based TCCON Data

Carbon dioxide (CO2) is one of the most significant greenhouse gases, and its concentration and distribution in the atmosphere have always been a research hotspot. To study the temporal and spatial characteristics of atmospheric CO2 globally, it is crucial to evaluate the consistency of observation data from different carbon observation satellites. This study utilizes data from the Total Carbon Column Observing Network (TCCON) to verify the column-averaged dry air mole fractions of atmospheric CO2 (XCO2) retrieved by satellites from October 2014 to May 2016, specifically comparing the XCO2 distributions from the Greenhouse Gases Observing Satellite (GOSAT) and Orbiting Carbon Observatory 2 (OCO-2). Our analysis indicates a strong correlation between the TCCON and both the GOSAT (correlation coefficient of 0.85) and OCO-2 (correlation coefficient of 0.91). Cross-validation further reveals that the measurements of the GOSAT and OCO-2 are highly consistent, with an average deviation and standard deviation of 0.92 ± 1.16 ppm and a correlation coefficient of 0.92. These differences remain stable over time, indicating that the calibration in the data set is reliable. Moreover, monthly averaged time-series and seasonal climatology comparisons were also performed separately over the six continents, i.e., Asia, North America, Europe, Africa, South America, and Oceania. The investigation of monthly XCO2 values across continents highlights greater consistency in Asia, North America, and Oceania (standard deviation from 0.15 to 0.27 ppm) as compared to Europe, South America, and Africa (standard deviation from 0.45 to 0.84 ppm). A seasonal analysis exhibited a high level of consistency in spring (correlation coefficient of 0.97), but lower agreement in summer (correlation coefficient of 0.78), potentially due to cloud cover and aerosol interference. Although some differences exist among the datasets, the overall findings demonstrate a strong correlation between the satellite measurements of XCO2. These results emphasize the importance of continued monitoring and calibration efforts to ensure the accurate assessment and understanding of atmospheric CO2 levels.

Spatial Statistical Prediction of Solar-Induced Chlorophyll Fluorescence (SIF) from Multivariate OCO-2 Data

Solar-induced chlorophyll fluorescence, or SIF, is a part of the natural process of photosynthesis. SIF can be measured from space by instruments such as the Orbiting Carbon Observatory-2 (OCO-2), making it a useful proxy for monitoring gross primary production (GPP), which is a critical component of Earth’s carbon cycle. The complex physical relationship between SIF and GPP is frequently studied using OCO-2 observations of SIF since they offer the finest spatial resolution available. However, measurement error (noise) and large gaps in spatial coverage limit the use of OCO-2 SIF to highly aggregated scales. To study the relationship between SIF and GPP across varying spatial scales, de-noised and gap-filled (i.e., Level 3) SIF data products are needed. Using a geostatistical methodology called cokriging, which includes kriging as a special case, we develop coSIF: a Level 3 SIF data product at a 0.05-degree resolution. As a natural secondary variable for cokriging, OCO-2 observes column-averaged atmospheric carbon dioxide concentrations (XCO2) simultaneously with SIF. There is a suggested lagged spatio-temporal dependence between SIF and XCO2, which we characterize through spatial covariance and cross-covariance functions. Our approach is highly parallelizable and accounts for non-stationary measurement errors in the observations. Importantly, each datum in the resulting coSIF data product is accompanied by a measure of uncertainty. Extant approaches do not provide formal uncertainty quantification, nor do they leverage the cross-dependence with XCO2.

Quantifying daily NOx and CO2 emissions from Wuhan using satellite observations from TROPOMI and OCO-2

Abstract. Quantification and control of NOx and CO2 emissions are essential across the world to limit adverse climate change and improve air quality. We present a new top-down method, an improved superposition column model to estimate day-to-day NOx and CO2 emissions from the large city of Wuhan, China, located in a polluted background. The latest released version 2.3.1 TROPOMI (TROPOspheric Monitoring Instrument) NO2 columns and version 10r of the Orbiting Carbon Observatory-2 (OCO-2)-observed CO2 mixing ratio are employed. We quantified daily NOx and CO2 emissions from Wuhan between September 2019 and October 2020 with an uncertainty of 31 % and 43 %, compared to 39 % and 49 % with the earlier v1.3 TROPOMI data, respectively. Our estimated NOx and CO2 emissions are verified against bottom-up inventories with minor deviations (<3 % for the 2019 mean, ranging from −20 % to 48 % on a daily basis). Based on the estimated CO2 emissions, we also predicted daily CO2 column mixing ratio enhancements, which match well with OCO-2 observations (<5 % bias, within ±0.3 ppm). We capture the day-to-day variation of NOx and CO2 emissions from Wuhan in 2019–2020, which does not reveal a substantial “weekend reduction” but does show a clear “holiday reduction” in the NOx and CO2 emissions. Our method also quantifies the abrupt decrease and slow NOx and CO2 emissions rebound due to the Wuhan lockdown in early 2020. This work demonstrates the improved superposition model to be a promising new tool for the quantification of city NOx and CO2 emissions, allowing policymakers to gain real-time information on spatial–temporal emission patterns and the effectiveness of carbon and nitrogen regulation in urban environments.

Estimating CO2 Emissions from Large Scale Coal-Fired Power Plants Using OCO-2 Observations and Emission Inventories

The concentration of atmospheric carbon dioxide (CO2) has increased rapidly worldwide, aggravating the global greenhouse effect, and coal-fired power plants are one of the biggest contributors of greenhouse gas emissions in China. However, efficient methods that can quantify CO2 emissions from individual coal-fired power plants with high accuracy are needed. In this study, we estimated the CO2 emissions of large-scale coal-fired power plants using Orbiting Carbon Observatory-2 (OCO-2) satellite data based on remote sensing inversions and bottom-up methods. First, we mapped the distribution of coal-fired power plants, displaying the total installed capacity, and identified two appropriate targets, the Waigaoqiao and Qinbei power plants in Shanghai and Henan, respectively. Then, an improved Gaussian plume model method was applied for CO2 emission estimations, with input parameters including the geographic coordinates of point sources, wind vectors from the atmospheric reanalysis of the global climate, and OCO-2 observations. The application of the Gaussian model was improved by using wind data with higher temporal and spatial resolutions, employing the physically based unit conversion method, and interpolating OCO-2 observations into different resolutions. Consequently, CO2 emissions were estimated to be 23.06 ± 2.82 (95% CI) Mt/yr using the Gaussian model and 16.28 Mt/yr using the bottom-up method for the Waigaoqiao Power Plant, and 14.58 ± 3.37 (95% CI) and 14.08 Mt/yr for the Qinbei Power Plant, respectively. These estimates were compared with three standard databases for validation: the Carbon Monitoring for Action database, the China coal-fired Power Plant Emissions Database, and the Carbon Brief database. The comparison found that previous emission inventories spanning different time frames might have overestimated the CO2 emissions of one of two Chinese power plants on the two days that the measurements were made. Our study contributes to quantifying CO2 emissions from point sources and helps in advancing satellite-based monitoring techniques of emission sources in the future; this helps in reducing errors due to human intervention in bottom-up statistical methods.

Linking global terrestrial CO<sub>2</sub> fluxes and environmental drivers: inferences from the Orbiting Carbon Observatory 2 satellite and terrestrial biospheric models

Abstract. Observations from the Orbiting Carbon Observatory 2 (OCO-2) satellite have been used to estimate CO2 fluxes in many regions of the globe and provide new insight into the global carbon cycle. The objective of this study is to infer the relationships between patterns in OCO-2 observations and environmental drivers (e.g., temperature, precipitation) and therefore inform a process understanding of carbon fluxes using OCO-2. We use a multiple regression and inverse model, and the regression coefficients quantify the relationships between observations from OCO-2 and environmental driver datasets within individual years for 2015–2018 and within seven global biomes. We subsequently compare these inferences to the relationships estimated from 15 terrestrial biosphere models (TBMs) that participated in the TRENDY model inter-comparison. Using OCO-2, we are able to quantify only a limited number of relationships between patterns in atmospheric CO2 observations and patterns in environmental driver datasets (i.e., 10 out of the 42 relationships examined). We further find that the ensemble of TBMs exhibits a large spread in the relationships with these key environmental driver datasets. The largest uncertainty in the models is in the relationship with precipitation, particularly in the tropics, with smaller uncertainties for temperature and photosynthetically active radiation (PAR). Using observations from OCO-2, we find that precipitation is associated with increased CO2 uptake in all tropical biomes, a result that agrees with half of the TBMs. By contrast, the relationships that we infer from OCO-2 for temperature and PAR are similar to the ensemble mean of the TBMs, though the results differ from many individual TBMs. These results point to the limitations of current space-based observations for inferring environmental relationships but also indicate the potential to help inform key relationships that are very uncertain in state-of-the-art TBMs.

Five years of variability in the global carbon cycle: comparing an estimate from the Orbiting Carbon Observatory-2 and process-based models

Year-to-year variability in CO2 fluxes can yield insight into climate-carbon cycle relationships, a fundamental yet uncertain aspect of the terrestrial carbon cycle. In this study, we use global observations from NASA’s Orbiting Carbon Observatory-2 (OCO-2) satellite for years 2015–2019 and a geostatistical inverse model to evaluate 5 years of interannual variability (IAV) in CO2 fluxes and its relationships with environmental drivers. OCO-2 launched in late 2014, and we specifically evaluate IAV during the time period when OCO-2 observations are available. We then compare inferences from OCO-2 with state-of-the-art process-based models (terrestrial biosphere model, TBMs). Results from OCO-2 suggest that the tropical grasslands biome (including grasslands, savanna, and agricultural lands within the tropics) makes contributions to global IAV during the 5 year study period that are comparable to tropical forests, a result that differs from a majority of TBMs. Furthermore, existing studies disagree on the environmental variables that drive IAV during this time period, and the analysis using OCO-2 suggests that both temperature and precipitation make comparable contributions. TBMs, by contrast, tend to estimate larger IAV during this time and usually estimate larger relative contributions from the extra-tropics. With that said, TBMs show little consensus on both the magnitude and the contributions of different regions to IAV. We further find that TBMs show a wide range of responses on the relationships of CO2 fluxes with annual anomalies in temperature and precipitation, and these relationships across most of the TBMs have a larger magnitude than inferred from OCO-2. Overall, the findings of this study highlight large uncertainties in process-based estimates of IAV during recent years and provide an avenue for evaluating these processes against inferences from OCO-2.

Evaluation of earth system model and atmospheric inversion using total column CO2 observations from GOSAT and OCO-2

The measurements of one of the major greenhouse gases, carbon dioxide (CO2), are being made using dedicated satellite remote sensing since the launch of the greenhouse gases observing satellite (GOSAT) by a three-way partnership between the Japan Aerospace Exploration Agency (JAXA), the Ministry of Environment (MoE) and the National Institute for Environmental Studies (NIES), and the National Aeronautics and Space Administration (NASA) Orbiting Carbon Observatory-2 (OCO-2). In the past 10 years, estimation of CO2 fluxes from land and ocean using the earth system models (ESMs) and inverse modelling of in situ atmospheric CO2 data have also made significant progress. We attempt, for the first time, to evaluate the CO2 fluxes simulated by an earth system model (MIROC-ES2L) and the fluxes estimated by an inverse model (MIROC4-Inv) using in situ data by comparing with GOSAT and OCO-2 observations. Both MIROC-ES2L and MIROC4-Inv fluxes are used in the MIROC4-atmospheric chemistry transport model (referred to as ACTM_ES2LF and ACTM_InvF, respectively) for calculating total column CO2 mole fraction (XCO2) that are sampled at the time and location of the satellite measurements. Both the ACTM simulations agreed well with the GOSAT and OCO-2 satellite observations, within 2 ppm for the spatial maps and time evolutions of the zonal mean distributions. Our results suggest that the inverse model using in situ data is more consistent with the OCO-2 retrievals, compared with those of the GOSAT XCO2 data due to the higher accuracy of the former. This suggests that the MIROC4-Inv fluxes are of sufficient quality to evaluate MIROC-ES2L simulated fluxes. The ACTM_ES2LF simulation shows a slightly weaker seasonal cycle for the meridional profiles of CO2 fluxes, compared with that from the ACTM_InvF. This difference is revealed by greater XCO2 differences for ACTM_ES2LF vs GOSAT, compared with those of ACTM_InvF vs GOSAT. Using remote sensing–based global products of leaf area index (LAI) and gross primary productivity (GPP) over land, we show a weaker sensitivity of MIROC-ES2L biospheric activities to the weather and climate in the tropical regions. Our results clearly suggest the usefulness of XCO2 measurements by satellite remote sensing for evaluation of large-scale ESMs, which so far remained untested by the sparse in situ data.

A fast cloud geometrical thickness retrieval algorithm for single-layer marine liquid clouds using OCO-2 oxygen A-band measurements

Knowledge of cloud geometrical thickness (H) is of importance for the study of radiative balance and cloud microphysics. However, the retrieval of H remains challenging, especially for passive instruments. In this work, we derive a semi-analytical algorithm for retrieving H of single-layer liquid cloud based on oxygen A-band hyperspectral measurements from NASA's Orbiting Carbon Observatory-2 (OCO-2). In this algorithm, a high-order correction is introduced to the approximation formula to accurately calculate oxygen A-band hyperspectral cloud reflectance. The algorithm can retrieve H in real time as it does not require the use of the time-consuming radiative transfer model for radiation calculation during each retrieval. In addition, the algorithm currently requires cloud optical depth, cloud top height, and aerosol properties measured by other instruments as input. In idealized simulations using ten thousand A-band spectra spanning a range of cloud cases, the root mean squared error (RMSE) of the retrieval is approximately 2.0 hPa (for low clouds, 1 hPa is about 10 m). We also retrieve H based on millions of real OCO-2 observations and compare the retrieval results with the cloud product from CloudSat/CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations). After abnormal samples are removed, the correlation coefficient is 0.716, the average bias is −15.6 hPa, and the RMSE is 27.4 hPa. The statistical results show that the absolute bias increases systematically with the reference cloud geometrical thickness, which may be caused by the unrealistic vertical homogeneous cloud assumption. The similar phenomenon was also found in comparison with OCO2CLD-LIDAR-AUX, a joint retrieval product based on OCO-2, CALIPSO, and accurate radiative transfer model. The fast algorithm shows a similar distribution of retrieved cloud bottom pressure and a smaller bias on cloud geometrical thickness retrieval compared to the OCO2CLD-LIDAR-AUX products.

Optimization of the OCO-2 Cloud Screening Algorithm and Evaluation against MODIS and TCCON Measurements over Land Surfaces in Europe and Japan

A method to tighten the cloud screening thresholds based on local conditions is used to provide more stringent schemes for Orbiting Carbon Observatory-2 (OCO-2) cloud screening algorithms. Cloud screening strategies are essential to remove scenes with significant cloud and/or aerosol contamination from OCO-2 observations, which helps to save on the data processing cost and ensure high quality retrievals of the column-averaged CO2 dry air mole fraction (XCO2). Based on the radiance measurements in the 0.76 μm O2A band, 1.61 μm (weak), and 2.06 μm (strong) CO2 bands, the current combination of the A-Band Preprocessor (ABP) algorithm and Iterative Maximum A Posteriori (IMAP) Differential Optical Absorption Spectroscopy (DOAS) Preprocessor (IDP) algorithm passes around 20%–25% of all soundings, which means that some contaminated scenes also pass the screening process. In this work, three independent pairs of threshold parameters used in the ABP and IDP algorithms are sufficiently tuned until the overall pass rate is close to the monthly clear-sky fraction from the MODIS cloud mask. The tightened thresholds are applied to observations over land surfaces in Europe and Japan in 2016. The results show improvement of agreement and positive predictive value compared to the collocated MODIS cloud mask, especially in summer and fall. In addition, analysis indicates that XCO2 retrievals with more stringent thresholds are in closer agreement with measurements from collocated Total Carbon Column Observing Network (TCCON) sites.

The impact of improved satellite retrievals on estimates of biospheric carbon balance

Abstract. The Orbiting Carbon Observatory 2 (OCO-2) is NASA's first satellite dedicated to monitoring CO2 from space and could provide novel insight into CO2 fluxes across the globe. However, one continuing challenge is the development of a robust retrieval algorithm: an estimate of atmospheric CO2 from satellite observations of near-infrared radiation. The OCO-2 retrievals have undergone multiple updates since the satellite's launch, and the retrieval algorithm is now on its ninth version. Some of these retrieval updates, particularly version 8, led to marked changes in the CO2 observations, changes of 0.5 ppm or more. In this study, we evaluate the extent to which current OCO-2 observations can constrain monthly CO2 sources and sinks from the biosphere, and we particularly focus on how this constraint has evolved with improvements to the OCO-2 retrieval algorithm. We find that improvements in the CO2 retrieval are having a potentially transformative effect on satellite-based estimates of the global biospheric carbon balance. The version 7 OCO-2 retrievals formed the basis of early inverse modeling studies using OCO-2 data; these observations are best equipped to constrain the biospheric carbon balance across only continental or hemispheric regions. By contrast, newer versions of the retrieval algorithm yield a far more detailed constraint, and we are able to constrain CO2 budgets for seven global biome-based regions, particularly during the Northern Hemisphere summer when biospheric CO2 uptake is greatest. Improvements to the OCO-2 observations have had the largest impact on glint-mode observations, and we also find the largest improvements in the terrestrial CO2 flux constraint when we include both nadir and glint data.

Prior biosphere model impact on global terrestrial CO<sub>2</sub> fluxes estimated from OCO-2 retrievals

Abstract. This study assesses the impact of different state of the art global biospheric CO2 flux models, when applied as prior information, on inverse model “top-down” estimates of terrestrial CO2 fluxes obtained when assimilating Orbiting Carbon Observatory 2 (OCO-2) observations. This is done with a series of observing system simulation experiments (OSSEs) using synthetic CO2 column-average dry air mole fraction (XCO2) retrievals sampled at the OCO-2 satellite spatiotemporal frequency. The OSSEs utilized a 4-D variational (4D-Var) assimilation system with the GEOS-Chem global chemical transport model (CTM) to estimate CO2 net ecosystem exchange (NEE) fluxes using synthetic OCO-2 observations. The impact of biosphere models in inverse model estimates of NEE is quantified by conducting OSSEs using the NASA-CASA, CASA-GFED, SiB-4, and LPJ models as prior estimates and using NEE from the multi-model ensemble mean of the Multiscale Synthesis and Terrestrial Model Intercomparison Project as the “truth”. Results show that the assimilation of simulated XCO2 retrievals at OCO-2 observing modes over land results in posterior NEE estimates which generally reproduce “true” NEE globally and over terrestrial TransCom-3 regions that are well-sampled. However, we find larger spread among posterior NEE estimates, when using different prior NEE fluxes, in regions and seasons that have limited OCO-2 observational coverage and a large range in “bottom-up” NEE fluxes. Seasonally averaged posterior NEE estimates had standard deviations (SD) of ∼10 % to ∼50 % of the multi-model-mean NEE for different TransCom-3 land regions with significant NEE fluxes (regions/seasons with a NEE flux ≥0.5 PgC yr−1). On a global average, the seasonally averaged residual impact of the prior model NEE assumption on the posterior NEE spread is ∼10 %–20 % of the posterior NEE mean. Additional OCO-2 OSSE simulations demonstrate that posterior NEE estimates are also sensitive to the assumed prior NEE flux uncertainty statistics, with spread in posterior NEE estimates similar to those when using variable prior model NEE fluxes. In fact, the sensitivity of posterior NEE estimates to prior error statistics was larger than prior flux values in some regions/times in the tropics and Southern Hemisphere where sufficient OCO-2 data were available and large differences between the prior and truth were evident. Overall, even with the availability of spatiotemporally dense OCO-2 data, noticeable residual differences (up to ∼20 %–30 % globally and 50 % regionally) in posterior NEE flux estimates remain that were caused by the choice of prior model flux values and the specification of prior flux uncertainties.

Estimating power plant CO2 emission using OCO-2 XCO2 and high resolution WRF-Chem simulations

Anthropogenic CO2 emission from fossil fuel combustion has major impacts on the global climate. The Orbiting Carbon Observatory 2 (OCO-2) observations have previously been used to estimate individual power plant emissions with a Gaussian plume model assuming constant wind fields. The present work assesses the feasibility of estimating power plant CO2 emission using high resolution chemistry transport model simulations with OCO-2 observation data. In the new framework, 1.33 km Weather Research and Forecasting-Chem (WRF)-Chem simulation results are used to calculate the Jacobian matrix, which is then used with the OCO-2 XCO2 data to obtain power plant daily mean emission rates, through a maximum likelihood estimation. We applied the framework to the seven OCO-2 observations of near mid-to-large coal burning power plants identified in Nassar et al (2017 Geophys. Res. Lett. 44, 10045–53). Our estimation results closely match the reported emission rates at the Westar power plant (Kansas, USA), with a reported value of 26.67 ktCO2/day, and our estimated value at 25.82–26.47 ktCO2/day using OCO-2 v8 data, and 22.09–22.80 ktCO2/day using v9 data. At Ghent, KY, USA, our estimations using three versions (v7, v8, and v9) range from 9.84–20.40 ktCO2/day, which are substantially lower than the reported value (29.17 ktCO2/day). We attribute this difference to diminished WRF-Chem wind field simulation accuracy. The results from the seven cases indicate that accurate estimation requires accurate meteorological simulations and high quality XCO2 data. In addition, the strength and orientation (relative to the OCO-2 ground track) of the XCO2 enhancement are important for accurate and reliable estimation. Compared with the Gaussian plume model based approach, the high resolution WRF-Chem simulation based approach provides a framework for addressing varying wind fields, and possible expansion to city level emission estimation.

Exploring SMAP and OCO-2 observations to monitor soil moisture control on photosynthetic activity of global drylands and croplands

Sun-induced chlorophyll fluorescence (SIF) has proven valuable as a proxy of terrestrial ecosystem photosynthesis that can be used to study the impact of environmental limiting factors on terrestrial ecosystem productivity. In global drylands and croplands, soil moisture is often a primary limiting factor to plant growth. Here we use the SIF retrievals from the Orbiting Carbon Observatory-2 (OCO-2) and soil moisture (SM) retrievals from the Soil Moisture Active Passive (SMAP) satellites to study the soil moisture control on photosynthetic activity of global drylands and croplands. We primarily analyze the predawn SM control on daily peak photosynthetic activity using the morning SM observations from SMAP and the noon SIF observations from OCO-2. On global drylands (i.e., dry sub-humid, semi-arid, and arid zones) and majority of croplands, soil water content is typically low, moisture in the topsoil is critical for plant growth, and satellite-based SIF and SM observations are expected to have high retrieval accuracy. The SMAP SM retrievals show positive daily relationships with OCO-2 SIF for drylands and croplands of the tropics and Australia, where SM is expected to limit plant growth and concurrent data records are sufficient to make statistical inferences. Negative relationships between SIF and SM were observed in forested areas of mid-latitude dry subhumid zones with high average annual SM for which SIF showed a positive relationship with air temperature. The relationships between SM and SIF follow the expected spatial patterns, positive and stronger in drier and warmer regions whereas weaker and often negative in colder and wetter regions. We find strong evidence that the OCO-2 SIF is accurately capturing monthly SMAP SM dynamics particularly for regions with distinct seasonality of rainfall such as the Sub-Saharan North Africa, Indian subcontinent, and South African countries. These results indicate new opportunities for monitoring dryland and cropland productivity solely using satellite observed SIF and soil moisture data.

Solar-induced chlorophyll fluorescence is strongly correlated with terrestrial photosynthesis for a wide variety of biomes: First global analysis based on OCO-2 and flux tower observations.

Solar-induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite-observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3km×2.25km) measured from the Orbiting Carbon Observatory-2 (OCO-2) provides the first opportunity to examine the SIF-GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO-2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO-2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757nm (R2 =0.72, p<0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2 =0.57-0.79, p<0.0001) except evergreen broadleaf forests (R2 =0.16, p<0.05) at the daily timescale. A higher slope was found for C4 grasslands and croplands than for C3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome-specific SIF-GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO-2 SIF generally had a better performance for predicting GPP than satellite-derived vegetation indices and a light use efficiency model. The universal SIF-GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO-2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

Characterizing biospheric carbon balance using CO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt; observations from the OCO-2 satellite

Abstract. NASA's Orbiting Carbon Observatory 2 (OCO-2) satellite launched in summer of 2014. Its observations could allow scientists to constrain CO2 fluxes across regions or continents that were previously difficult to monitor. This study explores an initial step toward that goal; we evaluate the extent to which current OCO-2 observations can detect patterns in biospheric CO2 fluxes and constrain monthly CO2 budgets. Our goal is to guide top-down, inverse modeling studies and identify areas for future improvement. We find that uncertainties and biases in the individual OCO-2 observations are comparable to the atmospheric signal from biospheric fluxes, particularly during Northern Hemisphere winter when biospheric fluxes are small. A series of top-down experiments indicate how these errors affect our ability to constrain monthly biospheric CO2 budgets. We are able to constrain budgets for between two and four global regions using OCO-2 observations, depending on the month, and we can constrain CO2 budgets at the regional level (i.e., smaller than seven global biomes) in only a handful of cases (16 % of all regions and months). The potential of the OCO-2 observations, however, is greater than these results might imply. A set of synthetic data experiments suggests that retrieval errors have a salient effect. Advances in retrieval algorithms and to a lesser extent atmospheric transport modeling will improve the results. In the interim, top-down studies that use current satellite observations are best-equipped to constrain the biospheric carbon balance across only continental or hemispheric regions.

Chlorophyll fluorescence observed by OCO-2 is strongly related to gross primary productivity estimated from flux towers in temperate forests

Solar-induced chlorophyll fluorescence (SIF) opens a new perspective on the monitoring of vegetation photosynthesis from space, and has been recently used to estimate gross primary productivity (GPP). However, previous studies on SIF were mainly based on satellite observations from the Greenhouse Gases Observing Satellite (GOSAT) and Global Ozone Monitoring Experiment-2 (GOME-2), and the evaluation of these coarse-resolution SIF measurements using GPP derived from eddy covariance (EC) flux towers has been hindered by the scale mismatch between satellite and tower footprints. We use new far-red SIF observations from the Orbiting Carbon Observatory-2 (OCO-2) with much finer spatial resolution and GPP data from EC flux towers from 2014 to 2016 to examine the relationship between GPP and SIF for temperate forests. The OCO-2 SIF tracked tower GPP well, and had strong correlation with tower GPP at both retrieval bands (757nm and 771nm) and both instantaneous (mid-day) and daily timescales, with the strongest relationship for the 757nm band (SIF757) at the daily timescale. Daily SIF757 exhibited much stronger correlation with tower GPP compared to the enhanced vegetation index (EVI) and normalized difference vegetation index (NDVI) based on the surface reflectance product derived from the moderate resolution imaging spectroradiometer (MODIS) onboard either Terra or Aqua satellite and had a similarly strong relationship as EVI based on the bidirectional reflectance distribution function (BRDF) corrected reflectance product (MCD43A4) based on data from both Terra and Aqua. Absorbed photosynthetically active radiation (APAR) explained 85% of the variance in SIF757, while the product of APAR and two environmental scalars - fTmin and fVPD (representing minimum temperature stress and water stress) explained higher variance (92%) in SIF757. This suggests that SIF mainly depends on APAR and is also affected by environmental stresses that determine photosynthetic light use efficiency. The BRDF-corrected EVI also has advantages over SIF in that EVI has spatially and temporally continuous coverage and has a much longer record (2000-present). The OCO-2 SIF757 estimated GPP well (R2=0.81, p<0.0001; RMSE=1.11gCm−2d−1), and its performance was comparable to or slightly better than that of the MODIS GPP algorithm, a GPP model based on the light use efficiency logic. With intensive OCO-2 observations acquired in the Target mode around Park Falls, a calibration/validation site for OCO-2, we also generated gridded SIF and then estimated GPP on a per-pixel basis at the landscape scale. Our findings demonstrate the strong ability of chlorophyll fluorescence observed by OCO-2 in estimating GPP for temperate forests and reveal its potential and limitations in future ecosystem functioning and carbon cycling studies.

Influence of El Niño on atmospheric CO2 over the tropical Pacific Ocean: Findings from NASA's OCO-2 mission.

Spaceborne observations of carbon dioxide (CO2) from the Orbiting Carbon Observatory-2 are used to characterize the response of tropical atmospheric CO2 concentrations to the strong El Niño event of 2015-2016. Although correlations between the growth rate of atmospheric CO2 concentrations and the El Niño-Southern Oscillation are well known, the magnitude of the correlation and the timing of the responses of oceanic and terrestrial carbon cycle remain poorly constrained in space and time. We used space-based CO2 observations to confirm that the tropical Pacific Ocean does play an early and important role in modulating the changes in atmospheric CO2 concentrations during El Niño events-a phenomenon inferred but not previously observed because of insufficient high-density, broad-scale CO2 observations over the tropics.

Orbiting Carbon Observatory-2 (OCO-2) cloud screening algorithms: validation against collocated MODIS and CALIOP data

Abstract. The objective of the National Aeronautics and Space Administration's (NASA) Orbiting Carbon Observatory-2 (OCO-2) mission is to retrieve the column-averaged carbon dioxide (CO2) dry air mole fraction (XCO2) from satellite measurements of reflected sunlight in the near-infrared. These estimates can be biased by clouds and aerosols, i.e., contamination, within the instrument's field of view. Screening of the most contaminated soundings minimizes unnecessary calls to the computationally expensive Level 2 (L2) XCO2 retrieval algorithm. Hence, robust cloud screening methods have been an important focus of the OCO-2 algorithm development team. Two distinct, computationally inexpensive cloud screening algorithms have been developed for this application. The A-Band Preprocessor (ABP) retrieves the surface pressure using measurements in the 0.76 µm O2 A band, neglecting scattering by clouds and aerosols, which introduce photon path-length differences that can cause large deviations between the expected and retrieved surface pressure. The Iterative Maximum A Posteriori (IMAP) Differential Optical Absorption Spectroscopy (DOAS) Preprocessor (IDP) retrieves independent estimates of the CO2 and H2O column abundances using observations taken at 1.61 µm (weak CO2 band) and 2.06 µm (strong CO2 band), while neglecting atmospheric scattering. The CO2 and H2O column abundances retrieved in these two spectral regions differ significantly in the presence of cloud and scattering aerosols. The combination of these two algorithms, which are sensitive to different features in the spectra, provides the basis for cloud screening of the OCO-2 data set.To validate the OCO-2 cloud screening approach, collocated measurements from NASA's Moderate Resolution Imaging Spectrometer (MODIS), aboard the Aqua platform, were compared to results from the two OCO-2 cloud screening algorithms. With tuning of algorithmic threshold parameters that allows for processing of ≃ 20–25 % of all OCO-2 soundings, agreement between the OCO-2 and MODIS cloud screening methods is found to be ≃ 85 % over four 16-day orbit repeat cycles in both the winter (December) and spring (April–May) for OCO-2 nadir-land, glint-land and glint-water observations.No major, systematic, spatial or temporal dependencies were found, although slight differences in the seasonal data sets do exist and validation is more problematic with increasing solar zenith angle and when surfaces are covered in snow and ice and have complex topography. To further analyze the performance of the cloud screening algorithms, an initial comparison of OCO-2 observations was made to collocated measurements from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). These comparisons highlight the strength of the OCO-2 cloud screening algorithms in identifying high, thin clouds but suggest some difficulty in identifying some clouds near the surface, even when the optical thicknesses are greater than 1.