ALBERT

Description: The sensitivity of the process parameters of the biosphere model BETHY (Biosphere Energy Transfer HYdrology) to choices of atmospheric concentration network, high frequency terrestrial fluxes, and the choice of flux measurement network is investigated by using a carbon cycle data assimilation system. Results show that monthly mean or low-frequency observations of CO2 concentration provide strong constraints on parameters relevant for net flux (NEP) but only weak constraints for parameters controlling gross fluxes. The use of high-frequency CO2 concentration observations, which has allowed a great refinement of spatial scales in direct inversions, adds little to the observing system in this case. This unexpected result is explained by the fact that the stations of the CO2 concentration network we are using are not well placed to measure such high frequency signals. Indeed, CO2 concentration sensitivities relevant for such high frequency fluxes are found to be largely confined in the vicinity of the corresponding fluxes, and are therefore not well observed by background monitoring stations. In contrast, our results clearly show the potential of flux measurements to better constrain the model parameters relevant for gross primary productivity (GPP) and net primary productivity (NPP). Given uncertainties in the spatial description of ecosystem functions we recommend a combined observing strategy.

Description: Terrestrial productivity in semi-arid woodlands is strongly susceptible to changes in precipitation, and semi-arid woodlands constitute an important element of the global water and carbon cycles. Here, we use the Carbon Cycle Data Assimilation System (CCDAS) to investigate the key parameters controlling ecological and hydrological activities for a semi-arid savanna woodland site in Maun, Botswana. Twenty-four eco-hydrological process parameters of a terrestrial ecosystem model are optimized against two data streams separately and simultaneously: daily averaged latent heat flux (LHF) derived from eddy covariance measurements, and decadal fraction of absorbed photosynthetically active radiation (FAPAR) derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Assimilation of both data streams LHF and FAPAR for the years 2000 and 2001 leads to improved agreement between measured and simulated quantities not only for LHF and FAPAR, but also for photosynthetic CO2 uptake. The mean uncertainty reduction (relative to the prior) over all parameters is 14.9% for the simultaneous assimilation of LHF and FAPAR, 8.5% for assimilating LHF only, and 6.1% for assimilating FAPAR only. The set of parameters with the highest uncertainty reduction is similar between assimilating only FAPAR or only LHF. The highest uncertainty reduction for all three cases is found for a parameter quantifying maximum plant-available soil moisture. This indicates that not only LHF but also satellite-derived FAPAR data can be used to constrain and indirectly observe hydrological quantities.

Description: The terrestrial biosphere is currently a strong sink for anthropogenic CO2 emissions. Through the radiative properties of CO2, the strength of this sink has a direct influence on the radiative budget of the global climate system. The accurate assessment of this sink and its evolution under a changing climate is, hence, paramount for any efficient management strategies of the terrestrial carbon sink to avoid dangerous climate change. Unfortunately, simulations of carbon and water fluxes with terrestrial biosphere models exhibit large uncertainties. A considerable fraction of this uncertainty reflects uncertainty in the parameter values of the process formulations within the models. This paper describes the systematic calibration of the process parameters of a terrestrial biosphere model against two observational data streams: remotely sensed FAPAR (fraction of absorbed photosynthetically active radiation) provided by the MERIS (ESA's Medium Resolution Imaging Spectrometer) sensor and in situ measurements of atmospheric CO2 provided by the GLOBALVIEW flask sampling network. We use the Carbon Cycle Data Assimilation System (CCDAS) to systematically calibrate some 70 parameters of the terrestrial BETHY (Biosphere Energy Transfer Hydrology) model. The simultaneous assimilation of all observations provides parameter estimates and uncertainty ranges that are consistent with the observational information. In a subsequent step these parameter uncertainties are propagated through the model to uncertainty ranges for predicted carbon fluxes. We demonstrate the consistent assimilation at global scale, where the global MERIS FAPAR product and atmospheric CO2 are used simultaneously. The assimilation improves the match to independent observations. We quantify how MERIS data improve the accuracy of the current and future (net and gross) carbon flux estimates (within and beyond the assimilation period). We further demonstrate the use of an interactive mission benefit analysis tool built around CCDAS to support the design of future space missions. We find that, for long-term averages, the benefit of FAPAR data is most pronounced for hydrological quantities, and moderate for quantities related to carbon fluxes from ecosystems. The benefit for hydrological quantities is highest for semi-arid tropical or sub-tropical regions. Length of mission or sensor resolution is of minor importance.

Description: Human impact on wildfires, a major earth system component, remains poorly understood. While local studies have found more fires close to settlements and roads, assimilated charcoal records and analyses of regional fire patterns from remote-sensing observations point to a decline in fire frequency with increasing human population. Here, we present a global analysis using three multi-year satellite-based burned-area products combined with a parameter estimation and uncertainty analysis with a non-linear model. We show that at the global scale, the impact of increasing population density is mainly to reduce fire frequency. Only for areas with up to 0.1 people per km2, we find that fire frequency increases by 10 to 20% relative to its value at no population. The results are robust against choice of burned-area data set, and indicate that at only very few places on earth, fire frequency is limited by human ignitions. Applying the results to historical population estimates results in a moderate but accelerating decline of global burned area by around 14% since 1800, with most of the decline since 1950.

Description: We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian plume diffusion model (LPDM). The tangent and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com), with additional manual pre- and post-processing aimed at improving the performance of the computing, including MPI (Message Passing Interface). As results, the adjoint of Eulerian model is discrete. Construction of the adjoint of the Lagrangian component did not require any code modification, as LPDMs are able to track a significant number of particles back in time and thereby calculate the sensitivity of observations to the neighboring emissions areas. Eulerian and Lagrangian adjoint components were coupled at the time boundary in the global domain.The results are verified using a series of test experiments. The forward simulation shown the coupled model is effective in reproducing the seasonal cycle and short-term variability of CO2 even in the case of multiple limiting factors, such as high uncertainty of fluxes and the low resolution of the Eulerian model. The adjoint model demonstrates the high accuracy compared to direct forward sensitivity calculations and fast performance. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation.

Description: Terrestrial ecosystem models are employed to calculate the sources and sinks of carbon dioxide between land and atmosphere. These models may be heavily parameterised. Where reliable estimates are unavailable for a parameter, it remains highly uncertain; uncertainty of parameters can substantially contribute to overall model output uncertainty. This paper builds on the work of the terrestrial Carbon Cycle Data Assimilation System (CCDAS), which, here, assimilates atmospheric CO2 concentrations to optimise 19 parameters of the underlying terrestrial ecosystem model (Biosphere Energy Transfer and Hydrology scheme, BETHY). Previous experiments have shown that the identified minimum may contain non-physical parameter values. One way to combat this problem is to use constrained optimisation and so avoid the optimiser searching non-physical regions. Another technique is to use penalty terms in the cost function, which are added when the optimisation searches outside of a specified region. The use of parameter transformations is a further method of avoiding this problem, where the optimisation is carried out in a transformed parameter space, thus ensuring that the optimal parameters at the minimum are in the physical domain. We compare these different methods of achieving meaningful parameter values, finding that the parameter transformation method shows consistent results and the other two provide no useful results.