Interannual variability in Australia's terrestrial carbon cycle constrained by multiple observation types

Abstract. Recent studies have shown that semi-arid ecosystems in Australia may be responsible for a significant part of the interannual variability in the global concentration of atmospheric carbon dioxide. Here we use a multiple constraints approach to calibrate a land surface model of Australian terrestrial carbon and water cycles, with a focus on interannual variability. We use observations of carbon and water fluxes at 14 OzFlux sites, as well as data on carbon pools, litterfall and streamflow. We include calibration of the function describing the response of heterotrophic respiration to soil moisture. We also explore the effect on modelled interannual variability of parameter equifinality, whereby multiple combinations of parameters can give an equally acceptable fit to the calibration data. We estimate interannual variability of Australian net ecosystem production (NEP) of 0.12–0.21 PgC yr−1 (1σ) over 1982–2013, with a high anomaly of 0.43–0.67 PgC yr−1 in 2011 relative to this period associated with exceptionally wet conditions following a prolonged drought. The ranges are due to the effect on calculated NEP anomaly of parameter equifinality, with similar contributions from equifinality in parameters associated with net primary production (NPP) and heterotrophic respiration. Our range of results due to parameter equifinality demonstrates how errors can be underestimated when a single parameter set is used.

[1]  Shaun Quegan,et al.  Model–data synthesis in terrestrial carbon observation: methods, data requirements and data uncertainty specifications , 2005 .

[2]  John Doherty,et al.  Two statistics for evaluating parameter identifiability and error reduction , 2009 .

[3]  P. Cook,et al.  Dynamics of component carbon fluxes in a semi‐arid Acacia woodland, central Australia , 2013 .

[4]  David J. Dahlstrom.  Calibration and Uncertainty Analysis for Complex Environmental Models , 2015 .

[5]  J. Doherty,et al.  Calibration‐constrained Monte Carlo analysis of highly parameterized models using subspace techniques , 2009 .

[6]  Ray Leuning,et al.  Diagnosing errors in a land surface model (CABLE) in the time and frequency domains , 2011 .

[7]  W. Meyer,et al.  Evaporation of perennial semi‐arid woodland in southeastern Australia is adapted for irregular but common dry periods , 2015 .

[8]  Jeffrey P. Walker,et al.  An introduction to the Australian and New Zealand flux tower network - OzFlux , 2016 .

[9]  D. Schimel,et al.  Remembrance of Weather Past: Ecosystem Responses to Climate Variability , 2005 .

[10]  Yaxing Wei,et al.  Toward more realistic projections of soil carbon dynamics by Earth system models , 2016 .

[11]  Q. Zhuang,et al.  Equifinality in parameterization of process‐based biogeochemistry models: A significant uncertainty source to the estimation of regional carbon dynamics , 2008 .

[12]  Randal D. Koster,et al.  On the Nature of Soil Moisture in Land Surface Models , 2009 .

[13]  D. Barrett,et al.  Spatial Estimates of Biomass in ‘Mature’ Native Vegetation , 2003 .

[14]  W. Parton,et al.  Intra‐annual and interannual variability of ecosystem processes in shortgrass steppe , 2000 .

[15]  Luis Guanter,et al.  Anomalous carbon uptake in Australia as seen by GOSAT , 2015 .

[16]  Li Zhang,et al.  Parameter identifiability, constraint, and equifinality in data assimilation with ecosystem models. , 2009, Ecological applications : a publication of the Ecological Society of America.

[17]  V. Haverd,et al.  Soil-Litter-Iso: A one-dimensional model for coupled transport of heat, water and stable isotopes in soil with a litter layer and root extraction , 2010 .

[18]  Edward A. King,et al.  Operational Delivery of Hydro-Meteorological Monitoring and Modeling Over the Australian Continent , 2009, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[19]  R. Leuning,et al.  Carbon and water fluxes over a temperate Eucalyptus forest and a tropical wet/dry savanna in Australia: measurements and comparison with MODIS remote sensing estimates , 2005 .

[20]  Ian G. Enting,et al.  A review of applications of model–data fusion to studies of terrestrial carbon fluxes at different scales , 2009 .

[21]  John Doherty,et al.  Reply to discussion by Carl W. Chen "Methodologies for calibration and predictive analysis of a watershed model" , 2003 .

[22]  Eva van Gorsel,et al.  OzFlux data: network integration from collection to curation , 2016 .

[23]  Tim R. McVicar,et al.  Climate‐related trends in Australian vegetation cover as inferred from satellite observations, 1981–2006 , 2009 .

[24]  Edwin W. Pak,et al.  An extended AVHRR 8‐km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data , 2005 .

[25]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[26]  A. Pitman,et al.  Evaluating the Performance of Land Surface Models , 2008 .

[27]  Jason Beringer,et al.  Patterns and processes of carbon, water and energy cycles across northern Australian landscapes: From point to region , 2011 .

[28]  John Doherty,et al.  Efficient accommodation of local minima in watershed model calibration , 2006 .

[29]  E. Davidson,et al.  Estimating parameters of a forest ecosystem C model with measurements of stocks and fluxes as joint constraints , 2010, Oecologia.

[30]  J. Canadell,et al.  The Australian terrestrial carbon budget , 2012 .

[31]  Atul K. Jain,et al.  Global Carbon Budget 2018 , 2014, Earth System Science Data.

[32]  P. Rayner,et al.  Interannual variability of the global carbon cycle (1992–2005) inferred by inversion of atmospheric CO2 and δ13CO2 measurements , 2008 .

[33]  G. Abramowitz,et al.  Sensitivity of net ecosystem exchange and heterotrophic respiration to parameterization uncertainty , 2013 .

[34]  C. M Holgatea,et al.  Comparison of remotely sensed and modelled soil moisture data sets across Australia , 2016 .

[35]  E. Kowalczyk,et al.  The CSIRO Atmosphere Biosphere Land Exchange (CABLE) model for use in climate models and as an offline model , 2006 .

[36]  M. R. R A U Pa C H,et al.  Model – data synthesis in terrestrial carbon observation : methods , data requirements and data uncertainty specifications , 2005 .

[37]  Tim R. McVicar,et al.  Evaluation of the remote-sensing-based DIFFUSE model for estimating photosynthesis of vegetation , 2014 .

[38]  Benjamin Smith,et al.  Process contributions of Australian ecosystems to interannual variations in the carbon cycle , 2016 .

[39]  Benjamin Smith,et al.  Dryland vegetation response to wet episode, not inherent shift in sensitivity to rainfall, behind Australia's role in 2011 global carbon sink anomaly , 2016, Global change biology.

[40]  S. Wofsy,et al.  Factors controlling CO2 exchange on timescales from hourly to decadal at Harvard Forest , 2007 .

[41]  R. Isbell,et al.  A description of Australian soils. , 1975 .

[42]  D. Jones,et al.  High-quality spatial climate data-sets for Australia , 2009 .

[43]  L. Hutley,et al.  Dynamic INtegrated Gap-filling and partitioning for OzFlux (DINGO) , 2016 .

[44]  Huimin Wang,et al.  Complementarity of flux- and biometric-based data to constrain parameters in a terrestrial carbon model , 2015 .

[45]  Peter R. J. North,et al.  New Vegetation Albedo Parameters and Global Fields of Soil Background Albedo Derived from MODIS for Use in a Climate Model , 2009 .

[46]  J. Canadell,et al.  Multiple observation types reduce uncertainty in Australia's terrestrial carbon and water cycles , 2013 .

[47]  R.A.M. de Jeu,et al.  Comparison of remotely sensed and modelled soil moisture data sets across Australia , 2016 .

[48]  Alessandro Anav,et al.  Global Data Sets of Vegetation Leaf Area Index (LAI)3g and Fraction of Photosynthetically Active Radiation (FPAR)3g Derived from Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) for the Period 1981 to 2011 , 2013, Remote. Sens..

[49]  Benjamin Smith,et al.  Coupling carbon allocation with leaf and root phenology predicts tree-grass partitioning along a savanna rainfall gradient , 2016 .

[50]  Jason Beringer,et al.  Carbon cycling in a mountain ash forest: Analysis of below ground respiration , 2007 .

[51]  Jin Teng,et al.  Rainfall-Runoff Modelling Across Southeast Australia: Datasets, Models and Results , 2011 .

[52]  David Cook,et al.  Quantification of terrestrial ecosystem carbon dynamics in the conterminous United States combining a process-based biogeochemical model and MODIS and AmeriFlux data , 2011 .

[53]  Philippe Ciais,et al.  Terrestrial biosphere model performance for inter‐annual variability of land‐atmosphere CO2 exchange , 2012 .

[54]  F. Zhao,et al.  Evaluation of AWRA-L: the Australian Water Resource Assessment model , 2015 .

[55]  Atul K. Jain,et al.  The dominant role of semi-arid ecosystems in the trend and variability of the land CO2 sink , 2015, Science.

[56]  R. Aalderink,et al.  Identification of the parameters describing primary production from continuous oxygen signals. , 1997 .

[57]  Yi Y. Liu,et al.  Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle , 2014, Nature.

[58]  Matthew F. McCabe,et al.  Recent reversal in loss of global terrestrial biomass , 2015 .

[59]  Ankur R. Desai,et al.  Climatic and phenological controls on coherent regional interannual variability of carbon dioxide flux in a heterogeneous landscape , 2010 .

[60]  Steven W. Running,et al.  The global NPP dependence on ENSO: La Niña and the extraordinary year of 2011 , 2013 .

[61]  Marko Scholze,et al.  On the capability of Monte Carlo and adjoint inversion techniques to derive posterior parameter uncertainties in terrestrial ecosystem models , 2012 .

[62]  M. Hutchinson,et al.  Integrating a global agro‐climatic classification with bioregional boundaries in Australia , 2005 .

[63]  Atul K. Jain,et al.  Global Carbon Budget 2018 , 2014, Earth System Science Data.

[64]  H. Keith,et al.  Dynamics of carbon exchange in a Eucalyptus forest in response to interacting disturbance factors , 2012 .

[65]  John Doherty,et al.  Approaches to Highly Parameterized Inversion: A Guide to Using PEST for Model-Parameter and Predictive-Uncertainty Analysis , 2014 .

[66]  Keith Beven,et al.  A manifesto for the equifinality thesis , 2006 .

[67]  E. Davidson,et al.  Using model‐data fusion to interpret past trends, and quantify uncertainties in future projections, of terrestrial ecosystem carbon cycling , 2012 .

[68]  D. Barrett,et al.  NPP Multi-Biome: VAST Calibration Data, 1965-1998, R1 , 2001 .

[69]  R. Leuning,et al.  Nocturnal carbon efflux: reconciliation of eddy covariance and chamber measurements using an alternative to the u.-threshold filtering technique , 2007 .

[70]  Yujie He,et al.  Toward more realistic projections of soil carbon dynamics by Earth system models , 2016 .

[71]  V. Haverd,et al.  CSIRO Marine and Atmospheric Research Component: Final Report for Phase 3 , 2008 .

[72]  J. Doherty,et al.  Uncertainty Analysis of a Groundwater Flow Model in East‐Central Florida , 2015, Ground water.

[73]  E. Davidson,et al.  Sensitivity of decomposition rates of soil organic matter with respect to simultaneous changes in temperature and moisture , 2015 .

[74]  L. Hutley,et al.  The utility of the eddy covariance techniques as a tool in carbon accounting: tropical savanna as a case study , 2005 .

[75]  L. Hutley,et al.  Savanna fires and their impact on net ecosystem productivity in North Australia , 2004 .

[76]  M. Raupach,et al.  Decomposition of vegetation cover into woody and herbaceous components using AVHRR NDVI time series , 2003 .

[77]  A. Pitman,et al.  Examining soil carbon uncertainty in a global model: response of microbial decomposition to temperature, moisture and nutrient limitation , 2013 .

[78]  Rachel M. Law,et al.  A global model of carbon, nitrogen and phosphorus cycles for the terrestrial biosphere , 2009 .

[79]  J. Gove,et al.  The REFLEX project: Comparing different algorithms and implementations for the inversion of a terrestrial ecosystem model against eddy covariance data , 2009 .