|
O'Connor, F. M.
Johnson, C. E.
Morgenstern, O.
Abraham, N. L.
Braesicke, P.
Dalvi, M.
Folberth, G. A.
Sanderson, M. G.
Telford, P. J.
Voulgarakis, A.
Young, P. J.
Zeng, G.
Collins, W. J.
Pyle, J. A. (2014). Evaluation of the new UKCA climate-composition model – Part 2: The Troposphere. Geoscientific Model Development, 7(1), 41–91.
Abstract: In this paper, we present a description of the tropospheric chemistry component of the UK Chemistry and Aerosols (UKCA) model which has been coupled to the Met Office Hadley Centre's HadGEM family of climate models. We assess the model's transport and scavenging processes, in particular focussing on convective transport, boundary layer mixing, wet scavenging and inter-hemispheric exchange. Simulations with UKCA of the short-lived radon tracer suggest that modelled distributions are comparable to those of other models and the comparison with observations indicate that apart from a few locations, boundary layer mixing and convective transport are effective in the model as a means of vertically redistributing surface emissions of radon. Comparisons of modelled lead tracer concentrations with observations suggest that UKCA captures surface concentrations in both hemispheres very well, although there is a tendency to underestimate the observed geographical and interannual variability in the Northern Hemisphere. In particular, UKCA replicates the shape and absolute concentrations of observed lead profiles, a key test in the evaluation of a model's wet scavenging scheme. The timescale for inter-hemispheric transport, calculated in the model using a simple krypton tracer experiment, does appear to be long relative to other models and could indicate deficiencies in tropical deep convection and/or insufficient boundary layer mixing. We also describe the main components of the tropospheric chemistry and evaluate it against observations and other tropospheric chemistry models. In particular, from a climate forcing perspective, present-day observed surface methane concentrations and tropospheric ozone concentrations are reproduced very well by the model, thereby making it suitable for long centennial integrations as well as studies of biogeochemical feedbacks. Results from both historical and future simulations with UKCA tropospheric chemistry are presented. Future projections of tropospheric ozone vary with the Representative Concentration Pathway (RCP). In RCP2.6, for example, tropospheric ozone increases up to 2010 and then declines by 13% of its year-2000 global mean by the end of the century. In RCP8.5, tropospheric ozone continues to rise steadily throughout the 21st century, with methane being the main driving factor. Finally, we highlight aspects of the UKCA model which are undergoing and/or have undergone recent developments and are suitable for inclusion in a next-generation Earth System Model.
Programme: 416
|
|
|
Tsuruta, A., T. Aalto, L. Backman, J. Hakkarainen, I.T. van der Laan-Luijkx, M.C. Krol, R. Spahni, S. Houweling, M. Laine, E. Dlugokencky, A.J. Gomez-Pelaez, M. van der Schoot, R. Langenfelds, R. Ellul, A. J., F. Apadula, C. Gerbig, D.G. Feist, R. Kivi, Y. Yoshida, and W. Peters. (2016). Development of CarbonTracker Europe-CH4 – Part 2: global methane emission estimates and their evaluation for 2000–2012. Geoscientific model development discussions, .
Abstract: Gobal methane emissions were estimated for 2000–2012 using the CarbonTracker Europe-CH4 (CTE-CH4) data assimilation system. In CTE-CH4, the anthropogenic and biosphere emissions of CH4 are simultaneously constrained by global atmospheric in-situ methane mole fraction observations. We use three configurations developed in Tsuruta et al. (2016) to assess the sensitivity of the CH4 flux estimates to (a) the number of unknown flux scaling factors to be optimized which in turn depends on the choice of underlying land-ecosystem map, and (b) on the parametrization of vertical mixing in the 30 atmospheric transport model TM5. The posterior emission estimates were evaluated by comparing simulations to surface in-situ observation sites, to profile observations made by aircraft, to dry air total column-averaged mole fractions (XCH4) observations from the Total Carbon Column Observing Network (TCCON), and to XCH4 retrievals from the Greenhouse gases Observing SATellite (GOSAT). Our estimated posterior mean global total emissions during 2000–2012 are 516 $\pm$ 51 Tg CH4 yr\textminus1, and emission estimates during 2007–2012 are 18 Tg CH4 yr\textminus1 greater than those from 2001–2006, mainly driven by an 35 increase in emissions from the south America temperate region, the Asia temperate region and Asia tropics. The sensitivity of the flux estimates to the underlying ecosystem map was large for the Asia temperate region and Australia, but not significant in the northern latitude regions, i.e. the north American boreal region, the north American temperate region and Europe. Instead, the posterior estimates for the northern latitude regions show larger sensitivity to the choice of convection scheme in TM5. The Gregory et al. (2000) mixing scheme with faster interhemispheric exchange leads to higher estimated CH4 emissions at northern latitudes, and lower emissions in southern latitudes, compared to the estimates using Tiedtke (1989) convection scheme. Our evaluation with non-assimilated observations showed that posterior mole fractions were better matched with the 5 observations when Gregory et al. (2000) convection scheme was used.
Programme: 416
|
|
|
Tsuruta, A., T. Aalto, L. Backman, J. Hakkarainen, I.T. van der Laan-Luijkx, M.C. Krol, R. Spahni, S. Houweling, M. Laine, M. van der Schoot, R. Langenfelds, R. Ellul, and W. Peters, Development of CarbonTracker Europe-CH4 – Part 1: system set-up and sensitivity analyses, Geosci. Model Dev. Discuss., 2016a. (2016). Development of CarbonTracker Europe-CH4 – Part 1: system set-up and sensitivity analyses. Geoscientific Model Development Discussions, .
Abstract: 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. Academic paper: Development of CarbonTracker Europe-CH4 – Part 1: system set-up and sensitivity analyses. Available from: https://www.researchgate.net/publication/306338110DevelopmentofCarbonTrackerEurope-CH4-Part1systemset-upandsensitivityanalyses [accessed Apr 4, 2017].
Programme: 416
|
|
|
Pengyang Song, Dmitry Sidorenko, Patrick Scholz, Maik Thomas, Gerrit Lohmann. (2023). The tidal effects in the Finite-volumE Sea ice–Ocean Model (FESOM2.1): a comparison between parameterised tidal mixing and explicit tidal forcing (Vol. 16).
Abstract: Tides are proved to have a significant effect on the ocean and climate. Previous modelling research either adds a tidal mixing parameterisation or an explicit tidal forcing to the ocean models. However, no research compares the two approaches in the same framework. Here we implement both schemes in a general ocean circulation model and assess both methods by comparing the results. The aspects for comparison involve hydrography, sea ice, meridional overturning circulation (MOC), vertical diffusivity, barotropic streamfunction and energy diagnostics. We conclude that although the mesh resolution is poor in resolving internal tides in most mid-latitude and shelf-break areas, explicit tidal forcing still shows stronger tidal mixing at the Kuril–Aleutian Ridge and the Indonesian Archipelago than the tidal mixing parameterisation. Beyond that, the explicit tidal forcing method leads to a stronger upper cell of the Atlantic MOC by enhancing the Pacific MOC and the Indonesian Throughflow. Meanwhile, the tidal mixing parameterisation leads to a stronger lower cell of the Atlantic MOC due to the tidal mixing in deep oceans. Both methods maintain the Antarctic Circumpolar Current at a higher level than the control run by increasing the meridional density gradient. We also show several phenomena that are not considered in the tidal mixing parameterisation, for example, the changing of energy budgets in the ocean system, the bottom drag induced mixing on the continental shelves and the sea ice transport by tidal motions. Due to the limit of computational capacity, an internal-tide-resolving simulation is not feasible for climate studies. However, a high-resolution short-term tidal simulation is still required to improve parameters and parameterisation schemes in climate studies.
Programme: 688
|
|
|
Nairita Pal, Kristin N. Barton, Mark R. Petersen, Steven R. Brus, Darren Engwirda, Brian K. Arbic, Andrew F. Roberts, Joannes J. Westerink, Damrongsak Wirasaet. (2023). Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities (Vol. 16).
Abstract: Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.
Programme: 688
|
|
|
Ole Richter, David E. Gwyther, Benjamin K. Galton-Fenzi, Kaitlin A. Naughten. (2022). The Whole Antarctic Ocean Model (WAOM v1.0): development and evaluation (Vol. 15).
Abstract: The Regional Ocean Modeling System (ROMS), including an ice shelf component, has been applied on a circum-Antarctic domain to derive estimates of ice shelf basal melting. Significant improvements made compared to previous models of this scale are the inclusion of tides and a horizontal spatial resolution of 2 km, which is sufficient to resolve on-shelf heat transport by bathymetric troughs and eddy-scale circulation. We run the model with ocean–atmosphere–sea ice conditions from the year 2007 to represent nominal present-day climate. We force the ocean surface with buoyancy fluxes derived from sea ice concentration observations and wind stress from ERA-Interim atmospheric reanalysis. Boundary conditions are derived from the ECCO2 ocean state estimate; tides are incorporated as sea surface height and barotropic currents at the open boundary. We evaluate model results using satellite-derived estimates of ice shelf melting and established compilations of ocean hydrography. The Whole Antarctic Ocean Model (WAOM v1.0) qualitatively captures the broad scale difference between warm and cold regimes as well as many of the known characteristics of regional ice–ocean interaction. We identify a cold bias for some warm-water ice shelves and a lack of high-salinity shelf water (HSSW) formation. We conclude that further calibration and development of our approach are justified. At its current state, the model is ideal for addressing specific, process-oriented questions, e.g. related to tide-driven ice shelf melting at large scales.
Programme: 688
|
|
|
Charles Amory, Christoph Kittel, Louis Le Toumelin, Cécile Agosta, Alison Delhasse, Vincent Favier, Xavier Fettweis. (2021). Performance of MAR (v3.11) in simulating the drifting-snow climate and surface mass balance of Adélie Land, East Antarctica (Vol. 14).
Abstract: Drifting snow, or the wind-driven transport of snow particles originating from clouds and the surface below and above 2 m above ground and their concurrent sublimation, is a poorly documented process on the Antarctic ice sheet, which is inherently lacking in most climate models. Since drifting snow mostly results from erosion of surface particles, a comprehensive evaluation of this process in climate models requires a concurrent assessment of simulated drifting-snow transport and the surface mass balance (SMB). In this paper a new version of the drifting-snow scheme currently embedded in the regional climate model MAR (v3.11) is extensively described. Several important modifications relative to previous version have been implemented and include notably a parameterization for drifting-snow compaction of the uppermost snowpack layer, differentiated snow density at deposition between precipitation and drifting snow, and a rewrite of the threshold friction velocity above which snow erosion initiates. Model results at high resolution (10 km) over Adélie Land, East Antarctica, for the period 2004–2018 are presented and evaluated against available near-surface meteorological observations at half-hourly resolution and annual SMB estimates. The evaluation demonstrates that MAR resolves the local drifting-snow frequency and transport up to the scale of the drifting-snow event and captures the resulting observed climate and SMB variability, suggesting that this model version can be used for continent-wide applications.
Programme: 411
|
|
|
Alexandra Touzeau, Amaëlle Landais, Samuel Morin, Laurent Arnaud, Ghislain Picard. (2018). Numerical experiments on vapor diffusion in polar snow and firn and its impact on isotopes using the multi-layer energy balance model Crocus in SURFEX v8.0 (Vol. 11).
Abstract: Abstract. To evaluate the impact of vapor diffusion on isotopic composition variations in snow pits and then in ice cores, we introduced water isotopes in the detailed snowpack model Crocus. At each step and for each snow layer, (1) the initial isotopic composition of vapor is taken at equilibrium with the solid phase, (2) a kinetic fractionation is applied during transport, and (3) vapor is condensed or snow is sublimated to compensate for deviation to vapor pressure at saturation. We study the different effects of temperature gradient, compaction, wind compaction, and precipitation on the final vertical isotopic profiles. We also run complete simulations of vapor diffusion along isotopic gradients and of vapor diffusion driven by temperature gradients at GRIP, Greenland and at Dome C, Antarctica over periods of 1 or 10 years. The vapor diffusion tends to smooth the original seasonal signal, with an attenuation of 7 to 12% of the original signal over 10 years at GRIP. This is smaller than the observed attenuation in ice cores, indicating that the model attenuation due to diffusion is underestimated or that other processes, such as ventilation, influence attenuation. At Dome C, the attenuation is stronger (18%), probably because of the lower accumulation and stronger δ18O gradients.
Programme: 1110
|
|
|
Mathieu Barrere, Florent Domine, Bertrand Decharme, Samuel Morin, Vincent Vionnet, Matthieu Lafaysse. (2017). Evaluating the performance of coupled snow–soil models in SURFEXv8 to simulate the permafrost thermal regime at a high Arctic site (Vol. 10).
Abstract: Abstract. Climate change projections still suffer from a limited representation of the permafrost–carbon feedback. Predicting the response of permafrost temperature to climate change requires accurate simulations of Arctic snow and soil properties. This study assesses the capacity of the coupled land surface and snow models ISBA-Crocus and ISBA-ES to simulate snow and soil properties at Bylot Island, a high Arctic site. Field measurements complemented with ERA-Interim reanalyses were used to drive the models and to evaluate simulation outputs. Snow height, density, temperature, thermal conductivity and thermal insulance are examined to determine the critical variables involved in the soil and snow thermal regime. Simulated soil properties are compared to measurements of thermal conductivity, temperature and water content. The simulated snow density profiles are unrealistic, which is most likely caused by the lack of representation in snow models of the upward water vapor fluxes generated by the strong temperature gradients within the snowpack. The resulting vertical profiles of thermal conductivity are inverted compared to observations, with high simulated values at the bottom of the snowpack. Still, ISBA-Crocus manages to successfully simulate the soil temperature in winter. Results are satisfactory in summer, but the temperature of the top soil could be better reproduced by adequately representing surface organic layers, i.e., mosses and litter, and in particular their water retention capacity. Transition periods (soil freezing and thawing) are the least well reproduced because the high basal snow thermal conductivity induces an excessively rapid heat transfer between the soil and the snow in simulations. Hence, global climate models should carefully consider Arctic snow thermal properties, and especially the thermal conductivity of the basal snow layer, to perform accurate predictions of the permafrost evolution under climate change.
Programme: 1042
|
|
|
Ingrid T. van der Laan-Luijkx, Ivar R. van der Velde, Emma van der Veen, Aki Tsuruta, Karolina Stanislawska, Arne Babenhauserheide, Hui Fang Zhang, Yu Liu, Wei He, Huilin Chen, Kenneth A. Masarie, Maarten C. Krol, Wouter Peters. (2017). The CarbonTracker Data Assimilation Shell (CTDAS) v1.0: implementation and global carbon balance 2001–2015 (Vol. 10).
Abstract: Abstract. Data assimilation systems are used increasingly to constrain the budgets of reactive and long-lived gases measured in the atmosphere. Each trace gas has its own lifetime, dominant sources and sinks, and observational network (from flask sampling and in situ measurements to space-based remote sensing) and therefore comes with its own optimal configuration of the data assimilation. The CarbonTracker Europe data assimilation system for CO2 estimates global carbon sources and sinks, and updates are released annually and used in carbon cycle studies. CarbonTracker Europe simulations are performed using the new modular implementation of the data assimilation system: the CarbonTracker Data Assimilation Shell (CTDAS). Here, we present and document this redesign of the data assimilation code that forms the heart of CarbonTracker, specifically meant to enable easy extension and modification of the data assimilation system. This paper also presents the setup of the latest version of CarbonTracker Europe (CTE2016), including the use of the gridded state vector, and shows the resulting carbon flux estimates. We present the distribution of the carbon sinks over the hemispheres and between the land biosphere and the oceans. We show that with equal fossil fuel emissions, 2015 has a higher atmospheric CO2 growth rate compared to 2014, due to reduced net land carbon uptake in later year. The European carbon sink is especially present in the forests, and the average net uptake over 2001–2015 was 0. 17 ± 0. 11PgC yr−1 with reductions to zero during drought years. Finally, we also demonstrate the versatility of CTDAS by presenting an overview of the wide range of applications for which it has been used so far.
Programme: 416
|
|