ALBERT

Description: Wetlands are the largest and most uncertain natural sources of atmospheric methane (CH4). Several process-based models have been developed to quantify the magnitude and estimate spatial and temporal variations in CH4 emissions from global wetlands. Reliable models are required to estimate global wetland CH4 emissions. This study aimed to test two process-based models, CH4MODwetland and Terrestrial Ecosystem Model (TEM), against the CH4 flux measurements of marsh, swamp, peatland and coastal wetland sites across the world; specifically, model accuracy and generality were evaluated for different wetland types and in different continents, and then the global CH4 emissions from 2000 to 2010 were estimated. Both models showed similar high correlations with the observed seasonal/annual total CH4 emissions, and the regression of the observed versus computed total seasonal/annual CH4 emissions resulted in R2 values of 0.81 and 0.68 for CH4MODwetland and TEM, respectively. The CH4MODwetland produced accurate predictions for marshes, peatlands, swamps and coastal wetlands, with model efficiency (EF) values of 0.22, 0.52, 0.13 and 0.72, respectively. TEM produced good predictions for peatlands and swamps, with EF values of 0.69 and 0.74, respectively, but it could not accurately simulate marshes and coastal wetlands (EF 

Description: We develop a new large-scale hydrological and water resources model, the Community Water Model (CWatM), which can simulate hydrology both globally and regionally at different resolutions from 30 arcmin to 30 arcsec at daily time steps. CWatM is open source in the Python programming environment and has a modular structure. It uses global, freely available data in the netCDF4 file format for reading, storage, and production of data in a compact way. CWatM includes general surface and groundwater hydrological processes but also takes into account human activities, such as water use and reservoir regulation, by calculating water demands, water use, and return flows. Reservoirs and lakes are included in the model scheme. CWatM is used in the framework of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP), which compares global model outputs. The flexible model structure allows for dynamic interaction with hydro-economic and water quality models for the assessment and evaluation of water management options. Furthermore, the novelty of CWatM is its combination of state-of-the-art hydrological modeling, modular programming, an online user manual and automatic source code documentation, global and regional assessments at different spatial resolutions, and a potential community to add to, change, and expand the open-source project. CWatM also strives to build a community learning environment which is able to freely use an open-source hydrological model and flexible coupling possibilities to other sectoral models, such as energy and agriculture.

Description: Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics.

Description: Chemistry plays an indispensable role in investigations of the atmosphere; however, many climate models either ignore or greatly simplify atmospheric chemistry, limiting both their accuracy and their scope. We present the development and evaluation of the online global atmospheric chemical model BCC-GEOS-Chem v1.0, coupling the GEOS-Chem chemical transport model (CTM) as an atmospheric chemistry component in the Beijing Climate Center atmospheric general circulation model (BCC-AGCM). The GEOS-Chem atmospheric chemistry component includes detailed tropospheric HOx–NOx–volatile organic compounds–ozone–bromine–aerosol chemistry and online dry and wet deposition schemes. We then demonstrate the new capabilities of BCC-GEOS-Chem v1.0 relative to the base BCC-AGCM model through a 3-year (2012–2014) simulation with anthropogenic emissions from the Community Emissions Data System (CEDS) used in the Coupled Model Intercomparison Project Phase 6 (CMIP6). The model captures well the spatial distributions and seasonal variations in tropospheric ozone, with seasonal mean biases of 0.4–2.2 ppbv at 700–400 hPa compared to satellite observations and within 10 ppbv at the surface to 500 hPa compared to global ozonesonde observations. The model has larger high-ozone biases over the tropics which we attribute to an overestimate of ozone chemical production. It underestimates ozone in the upper troposphere which is likely due either to the use of a simplified stratospheric ozone scheme or to biases in estimated stratosphere–troposphere exchange dynamics. The model diagnoses the global tropospheric ozone burden, OH concentration, and methane chemical lifetime to be 336 Tg, 1.16×106 molecule cm−3, and 8.3 years, respectively, which is consistent with recent multimodel assessments. The spatiotemporal distributions of NO2, CO, SO2, CH2O, and aerosol optical depth are generally in agreement with satellite observations. The development of BCC-GEOS-Chem v1.0 represents an important step for the development of fully coupled earth system models (ESMs) in China.

Description: We developed the WRF-GC model, an online coupling of the Weather Research and Forecasting (WRF) mesoscale meteorological model and the GEOS-Chem atmospheric chemistry model, for regional atmospheric chemistry and air quality modeling. WRF and GEOS-Chem are both open-source community models. WRF-GC offers regional modellers access to the latest GEOS-Chem chemical module, which is state of the science, well documented, traceable, benchmarked, actively developed by a large international user base, and centrally managed by a dedicated support team. At the same time, WRF-GC enables GEOS-Chem users to perform high-resolution forecasts and hindcasts for any region and time of interest. WRF-GC uses unmodified copies of WRF and GEOS-Chem from their respective sources; the coupling structure allows future versions of either one of the two parent models to be integrated into WRF-GC with relative ease. Within WRF-GC, the physical and chemical state variables are managed in distributed memory and translated between WRF and GEOS-Chem by the WRF-GC coupler at runtime. We used the WRF-GC model to simulate surface PM2.5 concentrations over China during 22 to 27 January 2015 and compared the results to surface observations and the outcomes from a GEOS-Chem Classic nested-China simulation. Both models were able to reproduce the observed spatiotemporal variations of regional PM2.5, but the WRF-GC model (r=0.68, bias =29 %) reproduced the observed daily PM2.5 concentrations over eastern China better than the GEOS-Chem Classic model did (r=0.72, bias =55 %). This was because the WRF-GC simulation, nudged with surface and upper-level meteorological observations, was able to better represent the pollution meteorology during the study period. The WRF-GC model is parallelized across computational cores and scales well on massively parallel architectures. In our tests where the two models were similarly configured, the WRF-GC simulation was 3 times more efficient than the GEOS-Chem Classic nested-grid simulation due to the efficient transport algorithm and the Message Passing Interface (MPI)-based parallelization provided by the WRF software framework. WRF-GC v1.0 supports one-way coupling only, using WRF-simulated meteorological fields to drive GEOS-Chem with no chemical feedbacks. The development of two-way coupling capabilities, i.e., the ability to simulate radiative and microphysical feedbacks of chemistry to meteorology, is under way. The WRF-GC model is open source and freely available from http://wrf.geos-chem.org (last access: 10 July 2020).

Description: Biomass burning emissions are a major source of trace gases and aerosols. Wildfires being highly variable in time and space, calculating emissions requires a numerical tool able to estimate fluxes at the kilometer scale and with an hourly time step. Here, the APIFLAME model version 2.0 is presented. It is structured to be modular in terms of input databases and processing methods. The main evolution compared to version 1.0 is the possibility of merging burned area and fire radiative power (FRP) satellite observations to modulate the temporal variations of fire emissions and to integrate small fires that may not be detected in the burned area product. Accounting for possible missed detection due to small fire results in an increase in burned area ranging from ∼5 % in Africa and Australia to ∼30 % in North America on average over the 2013–2017 time period based on the Moderate-Resolution Imaging Spectroradiometer (MODIS) Collection 6 fire products. An illustration for the case of southwestern Europe during the summer of 2016, marked by large wildfires in Portugal, is presented. Emissions calculated using different possible configurations of APIFLAME show a dispersion of 80 % on average over the domain during the largest wildfires (8–14 August 2016), which can be considered as an estimate of uncertainty of emissions. The main sources of uncertainty studied, by order of importance, are the emission factors, the calculation of the burned area, and the vegetation attribution. The aerosol (PM10) and carbon monoxide (CO) concentrations simulated with the CHIMERE regional chemistry transport model (CTM) are consistent with observations (good timing for the beginning and end of the events, ±1 d for the timing of the peak values) but tend to be overestimated compared to observations at surface stations. On the contrary, vertically integrated concentrations tend to be underestimated compared to satellite observations of total column CO by the Infrared Atmospheric Sounding Interferometer (IASI) instrument and aerosol optical depth (AOD) by MODIS. This underestimate is lower close to the fire region (5 %–40 % for AOD depending on the configuration and 8 %–18 % for total CO) but rapidly increases downwind. For all comparisons, better agreement is achieved when emissions are injected higher into the free troposphere using a vertical profile as estimated from observations of aerosol plume height by the Multi-angle Imaging SpectroRadiometer (MISR) satellite instrument (injection up to 4 km). Comparisons of aerosol layer heights to observations by the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) show that some parts of the plume may still be transported at too low an altitude. The comparisons of the different CTM simulations to observations point to uncertainties not only on emissions (total mass and daily variability) but also on the simulation of their transport with the CTM and mixing with other sources. Considering the uncertainty of the emission injection profile and of the modeling of the transport of these dense plumes, it is difficult to fully validate emissions through comparisons between model simulations and atmospheric observations.

Description: Based on the fifth phase of the Coupled Model Intercomparison Project (CMIP5)-generation previous Institut Pierre Simon Laplace (IPSL) Earth system model, we designed a new version, IPSL-CM5A2, aiming at running multi-millennial simulations typical of deep-time paleoclimate studies. Three priorities were followed during the setup of the model: (1) improving the overall model computing performance, (2) overcoming a persistent cold bias depicted in the previous model generation and (3) making the model able to handle the specific continental configurations of the geological past. These developments include the integration of hybrid parallelization Message Passing Interface – Open Multi-Processing (MPI-OpenMP) in the atmospheric model of the Laboratoire de Météorologie Dynamique (LMDZ), the use of a new library to perform parallel asynchronous input/output by using computing cores as “I/O servers” and the use of a parallel coupling library between the ocean and the atmospheric components. The model, which runs with an atmospheric resolution of 3.75∘×1.875∘ and 2 to 0.5∘ in the ocean, can now simulate ∼100 years per day, opening new possibilities towards the production of multi-millennial simulations with a full Earth system model. The tuning strategy employed to overcome a persistent cold bias is detailed. The confrontation of a historical simulation to climatological observations shows overall improved ocean meridional overturning circulation, marine productivity and latitudinal position of zonal wind patterns. We also present the numerous steps required to run IPSL-CM5A2 for deep-time paleoclimates through a preliminary case study for the Cretaceous. Namely, specific work on the ocean model grid was required to run the model for specific continental configurations in which continents are relocated according to past paleogeographic reconstructions. By briefly discussing the spin-up of such a simulation, we elaborate on the requirements and challenges awaiting paleoclimate modeling in the next years, namely finding the best trade-off between the level of description of the processes and the computing cost on supercomputers.

Description: In the frame of the project “MOSAIK – Model-based city planning and application in climate change”, a German-wide research project within the call “Urban Climate Under Change” ([UC]2) funded by the German Federal Ministry of Education and Research (BMBF), a biometeorology module was implemented into the Parallelized Large-Eddy Simulation Model (PALM) system. The new biometeorology module is comprised of methods for the calculation of UV-exposure quantities, a human–biometeorologically weighted mean radiant temperature (Tmrt), as well as for the estimation of human thermal comfort or stress. The latter is achieved through the implementation of the three widely used thermal indices: perceived temperature (PT), Universal Thermal Climate Index (UTCI), as well as physiologically equivalent temperature (PET). Comparison calculations were performed for the PT, UTCI and PET indices based on the SkyHelios model and showing PALM calculates higher values in general. This is mostly due to a higher radiational gain leading to higher values of mean radiant temperature. For a more direct comparison, the PT, PET and UTCI indices were calculated by the biometeorology module, as well as the programs provided by the attachment to Verein Deutscher Ingenieure (VDI) guideline 3787, as well as by the RayMan model based on the very same input dataset. Results show deviations below the relevant precision of 0.1 K for PET and UTCI and some deviations of up to 2.683 K for PT caused by repeated unfavorable rounding in very rare cases (0.027 %).

Description: The complex and computationally expensive nature of landscape evolution models poses significant challenges to the inference and optimization of unknown model parameters. Bayesian inference provides a methodology for estimation and uncertainty quantification of unknown model parameters. In our previous work, we developed parallel tempering Bayeslands as a framework for parameter estimation and uncertainty quantification for the Badlands landscape evolution model. Parallel tempering Bayeslands features high-performance computing that can feature dozens of processing cores running in parallel to enhance computational efficiency. Nevertheless, the procedure remains computationally challenging since thousands of samples need to be drawn and evaluated. In large-scale landscape evolution problems, a single model evaluation can take from several minutes to hours and in some instances, even days or weeks. Surrogate-assisted optimization has been used for several computationally expensive engineering problems which motivate its use in optimization and inference of complex geoscientific models. The use of surrogate models can speed up parallel tempering Bayeslands by developing computationally inexpensive models to mimic expensive ones. In this paper, we apply surrogate-assisted parallel tempering where the surrogate mimics a landscape evolution model by estimating the likelihood function from the model. We employ a neural-network-based surrogate model that learns from the history of samples generated. The entire framework is developed in a parallel computing infrastructure to take advantage of parallelism. The results show that the proposed methodology is effective in lowering the computational cost significantly while retaining the quality of model predictions.

Description: This paper describes the implementation of a coupling between a three-dimensional ocean general circulation model (NEMO) and a wave model (WW3) to represent the interactions of upper-oceanic flow dynamics with surface waves. The focus is on the impact of such coupling on upper-ocean properties (temperature and currents) and mixed layer depth (MLD) at global eddying scales. A generic coupling interface has been developed, and the NEMO governing equations and boundary conditions have been adapted to include wave-induced terms following the approach of McWilliams et al. (2004) and Ardhuin et al. (2008). In particular, the contributions of Stokes–Coriolis, vortex, and surface pressure forces have been implemented on top of the necessary modifications of the tracer–continuity equation and turbulent closure scheme (a one-equation turbulent kinetic energy – TKE – closure here). To assess the new developments, we perform a set of sensitivity experiments with a global oceanic configuration at 1/4∘ resolution coupled with a wave model configured at 1/2∘ resolution. Numerical simulations show a global increase in wind stress due to the interaction with waves (via the Charnock coefficient), particularly at high latitudes, resulting in increased surface currents. The modifications brought to the TKE closure scheme and the inclusion of a parameterization for Langmuir turbulence lead to a significant increase in the mixing, thus helping to deepen the MLD. This deepening is mainly located in the Southern Hemisphere and results in reduced sea surface currents and temperatures.