Changes in surface hydrology, soil moisture and gross primary production in the Amazon during the 2015/2016 El Niño

The 2015/2016 El Niño event caused severe changes in precipitation across the tropics. This impacted surface hydrology, such as river run-off and soil moisture availability, thereby triggering reductions in gross primary production (GPP). Many biosphere models lack the detailed hydrological component required to accurately quantify anomalies in surface hydrology and GPP during droughts in tropical regions. Here, we take the novel approach of coupling the biosphere model SiBCASA with the advanced hydrological model PCR-GLOBWB to attempt such a quantification across the Amazon basin during the drought in 2015/2016. We calculate 30-40% reduced river discharge in the Amazon starting in October 2015, lagging behind the precipitation anomaly by approximately one month and in good agreement with river gauge observations. Soil moisture shows distinctly asymmetrical spatial anomalies with large reductions across the north-eastern part of the basin, which persisted into the following dry season. This added to drought stress in vegetation, already present owing to vapour pressure deficits at the leaf, resulting in a loss of GPP of 0.95 (0.69 to 1.20) PgC between October 2015 and March 2016 compared with the 2007-2014 average. Only 11% (10-12%) of the reduction in GPP was found in the (wetter) north-western part of the basin, whereas the north-eastern and southern regions were affected more strongly, with 56% (54-56%) and 33% (31-33%) of the total, respectively. Uncertainty on this anomaly mostly reflects the unknown rooting depths of vegetation.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.</p

Changes in surface hydrology, soil moisture and gross primary production in the Amazon during the 2015/2016 El Niño

The 2015/2016 El Niño event caused severe changes in precipitation across the tropics. This impacted surface hydrology, such as river run-off and soil moisture availability, thereby triggering reductions in gross primary production (GPP). Many biosphere models lack the detailed hydrological component required to accurately quantify anomalies in surface hydrology and GPP during droughts in tropical regions. Here, we take the novel approach of coupling the biosphere model SiBCASA with the advanced hydrological model PCR-GLOBWB to attempt such a quantification across the Amazon basin during the drought in 2015/2016. We calculate 30-40% reduced river discharge in the Amazon starting in October 2015, lagging behind the precipitation anomaly by approximately one month and in good agreement with river gauge observations. Soil moisture shows distinctly asymmetrical spatial anomalies with large reductions across the north-eastern part of the basin, which persisted into the following dry season. This added to drought stress in vegetation, already present owing to vapour pressure deficits at the leaf, resulting in a loss of GPP of 0.95 (0.69 to 1.20) PgC between October 2015 and March 2016 compared with the 2007-2014 average. Only 11% (10-12%) of the reduction in GPP was found in the (wetter) north-western part of the basin, whereas the north-eastern and southern regions were affected more strongly, with 56% (54-56%) and 33% (31-33%) of the total, respectively. Uncertainty on this anomaly mostly reflects the unknown rooting depths of vegetation.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.</p

Widespread reduction in sun-induced fluorescence from the Amazon during the 2015/2016 El Nino

The tropical carbon balance dominates year-to-year variations in the CO2 exchange with the atmosphere through photosynthesis, respiration and fires. Because of its high correlation with gross primary productivity (GPP), observations of sun-induced fluorescence (SIF) are of great interest. We developed a new remotely sensed SIF product with improved signal-to-noise in the tropics, and use it here to quantify the impact of the 2015/2016 El Nino Amazon drought. We find that SIF was strongly suppressed over areas with anomalously high temperatures and decreased levels of water in the soil. SIF went below its climatological range starting from the end of the 2015 dry season (October) and returned to normal levels by February 2016 when atmospheric conditions returned to normal, but well before the end of anomalously low precipitation that persisted through June 2016. Impacts were not uniform across the Amazon basin, with the eastern part experiencing much larger (10-15%) SIF reductions than the western part of the basin (2-5%). We estimate the integrated loss of GPP relative to eight previous years to be 0.34-0.48 PgC in the three-month period October-November-December 2015. This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Nino on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'

Evaluating atmospheric methane inversion model results for Pallas, northern Finland

A state-of-the-art inverse model, CarbonTracker Data Assimilation Shell (CTDAS), was used to optimize estimates of methane (CH4) surface fluxes using atmospheric observations of CH4 as a constraint. The model consists of the latest version of the TM5 atmospheric chemistry-transport model and an ensemble Kalman filter based data assimilation system. The model was constrained by atmospheric methane surface concentrations, obtained from the World Data Centre for Greenhouse Gases (WDCGG). Prior methane emissions were specified for five sources: biosphere, anthropogenic, fire, termites and ocean, of which bio-sphere and anthropogenic emissions were optimized. Atmospheric CH 4 mole fractions for 2007 from northern Finland calculated from prior and optimized emissions were compared with observations. It was found that the root mean squared errors of the posterior esti - mates were more than halved. Furthermore, inclusion of NOAA observations of CH 4 from weekly discrete air samples collected at Pallas improved agreement between posterior CH 4 mole fraction estimates and continuous observations, and resulted in reducing optimized biosphere emissions and their uncertainties in northern Finland

CTDAS-Lagrange v1.0:a high-resolution data assimilation system for regional carbon dioxide observations

We have implemented a regional carbon dioxide data assimilation system based on the CarbonTracker Data Assimilation Shell (CTDAS) and a high-resolution Lagrangian transport model, the Stochastic Time-Inverted Lagrangian Transport model driven by the Weather Forecast and Research meteorological fields (WRF-STILT). With this system, named CTDAS-Lagrange, we simultaneously optimize terrestrial biosphere fluxes and four parameters that adjust the lateral boundary conditions (BCs) against CO2 observations from the NOAA ESRL North America tall tower and aircraft programmable flask packages (PFPs) sampling program. Least-squares optimization is performed with a time-stepping ensemble Kalman smoother, over a time window of 10 days and assimilating sequentially a time series of observations. Because the WRF-STILT footprints are pre-computed, it is computationally efficient to run the CTDAS-Lagrange system.To estimate the uncertainties in the optimized fluxes from the system, we performed sensitivity tests with various a priori biosphere fluxes (SiBCASA, SiB3, CT2013B) and BCs (optimized mole fraction fields from CT2013B and CTE2014, and an empirical dataset derived from aircraft observations), as well as with a variety of choices on the ways that fluxes are adjusted (additive or multiplicative), covariance length scales, biosphere flux covariances, BC parameter uncertainties, and model-data mismatches. In pseudo-data experiments, we show that in our implementation the additive flux adjustment method is more flexible in optimizing net ecosystem exchange (NEE) than the multiplicative flux adjustment method, and our sensitivity tests with real observations show that the CTDAS-Lagrange system has the ability to correct for the potential biases in the lateral BCs and to resolve large biases in the prior biosphere fluxes.Using real observations, we have derived a range of estimates for the optimized carbon fluxes from a series of sensitivity tests, which places the North American carbon sink for the year 2010 in a range from -0.92 to -1.26 PgC yr( -1). This is comparable to the TM5-based estimates of CarbonTracker (version CT2016, -0.91 +/- 1.10 PgC yr (-1)) and CarbonTracker Europe (version CTE,2016, -0.91 +/- 0.31 PgC yr(-1)). We conclude that CTDAS-Lagrange can offer a versatile and computationally attractive alternative to these global systems for regional estimates of carbon fluxes, which can take advantage of high-resolution Lagrangian footprints that are increasingly easy to obtain

Development of CarbonTracker Europe-CH4 – part 1 : system set-up and sensitivity analyses

CarbonTracker Europe-CH4 (CTE-CH4) inverse model versions 1.0 and 1.1 are presented. The model optimizes global surface methane emissions from biosphere and anthropogenic sources using an ensemble Kalman filter (EnKF) based optimization method, using the TM5 chemistry transport model as an observation operator, and assimilating global in-situ atmospheric methane mole fraction observations. In this study, we examine sensitivity of our CH4 emission estimates on the ensemble size, covariance matrix, prior estimates, observations to be assimilated, assimilation window length, convection scheme in TM5, and model structure in the emission estimates by performing CTE-CH4 with several set-ups. The analyses show that the model is sensitive to most of the parameters and inputs that were examined. Firstly, using a large enough ensemble size stabilises the results. Secondly, using an informative covariance matrix reduces uncertainty estimates. Thirdly, agreement with discrete observations became better when assimilating continuous observations. Finally, the posterior emissions were found sensitive to the choice of prior estimates, convection scheme and model structure, particularly to their spatial distribution. The distribution of posterior mole fractions derived from posterior emissions is consistent with the observations to the extent prescribed in the various covariance estimates, indicating a satisfactory performance of our system.peer-reviewe

Global 3-D Simulations of the Triple Oxygen Isotope Signature Δ17O in Atmospheric CO2

The triple oxygen isotope signature Δ¹⁷O in atmospheric CO₂, also known as its “¹⁷O excess,” has been proposed as a tracer for gross primary production (the gross uptake of CO₂ by vegetation through photosynthesis). We present the first global 3-D model simulations for Δ¹⁷O in atmospheric CO₂ together with a detailed model description and sensitivity analyses. In our 3-D model framework we include the stratospheric source of Δ¹⁷O in CO₂ and the surface sinks from vegetation, soils, ocean, biomass burning, and fossil fuel combustion. The effect of oxidation of atmospheric CO on Δ¹⁷O in CO2 is also included in our model. We estimate that the global mean Δ¹⁷O (defined as Δ¹⁷O = ln( ¹⁷O + 1) − RL · ln( ¹⁸O + 1) with RL = 0.5229) of CO₂ in the lowest 500 m of the atmosphere is 39.6 per meg, which is ∼20 per meg lower than estimates from existing box models. We compare our model results with a measured stratospheric Δ¹⁷O in CO₂ profile from Sodankylä (Finland), which shows good agreement. In addition, we compare our model results with tropospheric measurements of Δ¹⁷O in CO₂ from Göttingen (Germany) and Taipei (Taiwan), which shows some agreement but we also find substantial discrepancies that are subsequently discussed. Finally, we show model results for Zotino (Russia), Mauna Loa (United States), Manaus (Brazil), and South Pole, which we propose as possible locations for future measurements of Δ¹⁷O in tropospheric CO₂ that can help to further increase our understanding of the global budget of Δ¹⁷O in atmospheric CO₂

CTDAS-Lagrange v1.0 : A high-resolution data assimilation system for regional carbon dioxide observations

We have implemented a regional carbon dioxide data assimilation system based on the CarbonTracker Data Assimilation Shell (CTDAS) and a high-resolution Lagrangian transport model, the Stochastic Time-Inverted Lagrangian Transport model driven by the Weather Forecast and Research meteorological fields (WRF-STILT). With this system, named CTDAS-Lagrange, we simultaneously optimize terrestrial biosphere fluxes and four parameters that adjust the lateral boundary conditions (BCs) against CO2 observations from the NOAA ESRL North America tall tower and aircraft programmable flask packages (PFPs) sampling program. Least-squares optimization is performed with a time-stepping ensemble Kalman smoother, over a time window of 10 days and assimilating sequentially a time series of observations. Because the WRF-STILT footprints are pre-computed, it is computationally efficient to run the CTDAS-Lagrange system. To estimate the uncertainties in the optimized fluxes from the system, we performed sensitivity tests with various a priori biosphere fluxes (SiBCASA, SiB3, CT2013B) and BCs (optimized mole fraction fields from CT2013B and CTE2014, and an empirical dataset derived from aircraft observations), as well as with a variety of choices on the ways that fluxes are adjusted (additive or multiplicative), covariance length scales, biosphere flux covariances, BC parameter uncertainties, and model-data mismatches. In pseudo-data experiments, we show that in our implementation the additive flux adjustment method is more flexible in optimizing net ecosystem exchange (NEE) than the multiplicative flux adjustment method, and our sensitivity tests with real observations show that the CTDAS-Lagrange system has the ability to correct for the potential biases in the lateral BCs and to resolve large biases in the prior biosphere fluxes. Using real observations, we have derived a range of estimates for the optimized carbon fluxes from a series of sensitivity tests, which places the North American carbon sink for the year 2010 in a range from -0.92 to -1.26PgCyr-1. This is comparable to the TM5-based estimates of CarbonTracker (version CT2016, -0.91±1.10PgCyr-1) and CarbonTracker Europe (version CTE2016, -0.91±0.31PgCyr-1). We conclude that CTDAS-Lagrange can offer a versatile and computationally attractive alternative to these global systems for regional estimates of carbon fluxes, which can take advantage of high-resolution Lagrangian footprints that are increasingly easy to obtain.</p

Warm spring reduced carbon cycle impact of the 2012 US summer drought

The global terrestrial carbon sink offsets one-third of the world’s fossil fuel emissions, but the strength of this sink is highly sensitive to large-scale extreme events. In 2012, the contiguous United States experienced exceptionally warm temperatures and the most severe drought since the Dust Bowl era of the 1930s, resulting in substantial economic damage. It is crucial to understand the dynamics of such events because warmer temperatures and a higher prevalence of drought are projected in a changing climate. Here, we combine an extensive network of direct ecosystem flux measurements with satellite remote sensing and atmospheric inverse modeling to quantify the impact of the warmer spring and summer drought on biosphere-atmosphere carbon and water exchange in 2012. We consistently find that earlier vegetation activity increased spring carbon uptake and compensated for the reduced uptake during the summer drought, which mitigated the impact on net annual carbon uptake. The early phenological development in the Eastern Temperate Forests played a major role for the continental-scale carbon balance in 2012. The warm spring also depleted soil water resources earlier, and thus exacerbated water limitations during summer. Our results show that the detrimental effects of severe summer drought on ecosystem carbon storage can be mitigated by warming-induced increases in spring carbon uptake. However, the results also suggest that the positive carbon cycle effect of warm spring enhances water limitations and can increase summer heating through biosphere–atmosphere feedbacks.Data deposition: The eddy-covariance data are available in the AmeriFlux data archive at the Carbon Dioxide Information Analysis Center at the Oak Ridge National Laboratory (cdiac.ornl.gov/ftp/ameriflux/data). This is the publisher’s final pdf. The article is published by National Academy of Sciences and can be found at: http://www.pnas.org/Keywords: ecosystem fluxes, eddy covariance, carbon uptake, biosphere, seasonal climate anomalies, atmosphere feedback