scholarly journals Comparing Lagrangian and Eulerian models for CO<sub>2</sub> transport – a step towards Bayesian inverse modeling using WRF/STILT-VPRM

2012 ◽  
Vol 12 (19) ◽  
pp. 8979-8991 ◽  
Author(s):  
D. Pillai ◽  
C. Gerbig ◽  
R. Kretschmer ◽  
V. Beck ◽  
U. Karstens ◽  
...  

Abstract. We present simulations of atmospheric CO2 concentrations provided by two modeling systems, run at high spatial resolution: the Eulerian-based Weather Research Forecasting (WRF) model and the Lagrangian-based Stochastic Time-Inverted Lagrangian Transport (STILT) model, both of which are coupled to a diagnostic biospheric model, the Vegetation Photosynthesis and Respiration Model (VPRM). The consistency of the simulations is assessed with special attention paid to the details of horizontal as well as vertical transport and mixing of CO2 concentrations in the atmosphere. The dependence of model mismatch (Eulerian vs. Lagrangian) on models' spatial resolution is further investigated. A case study using airborne measurements during which two models showed large deviations from each other is analyzed in detail as an extreme case. Using aircraft observations and pulse release simulations, we identified differences in the representation of details in the interaction between turbulent mixing and advection through wind shear as the main cause of discrepancies between WRF and STILT transport at a spatial resolution such as 2 and 6 km. Based on observations and inter-model comparisons of atmospheric CO2 concentrations, we show that a refinement of the parameterization of turbulent velocity variance and Lagrangian time-scale in STILT is needed to achieve a better match between the Eulerian and the Lagrangian transport at such a high spatial resolution (e.g. 2 and 6 km). Nevertheless, the inter-model differences in simulated CO2 time series for a tall tower observatory at Ochsenkopf in Germany are about a factor of two smaller than the model-data mismatch and about a factor of three smaller than the mismatch between the current global model simulations and the data.

2012 ◽  
Vol 12 (1) ◽  
pp. 1267-1298
Author(s):  
D. Pillai ◽  
C. Gerbig ◽  
R. Kretschmer ◽  
V. Beck ◽  
U. Karstens ◽  
...  

Abstract. We present simulations of atmospheric CO2 concentrations provided by two modeling systems, run at high spatial resolution: the Eulerian-based Weather Research Forecasting (WRF) model and the Lagrangian-based Stochastic Time-Inverted Lagrangian Transport (STILT) model, both of which are coupled to a diagnostic biospheric model, the Vegetation Photosynthesis and Respiration Model (VPRM). The consistency of the simulations is assessed with special attention paid to the details of horizontal as well as vertical transport and mixing of CO2 concentrations in the atmosphere. The dependence of model mismatch (Eulerian vs. Lagrangian) on models' spatial resolution is further investigated. A case study using airborne measurements during which both models showed large deviations from each other is analyzed in detail as an extreme case. Using aircraft observations and pulse release simulations, we identified differences in the representation of details in the interaction between turbulent mixing and advection through wind shear as the main cause of discrepancies between WRF and STILT transport at a spatial resolution such as 2 and 6 km. Based on observations and inter-model comparisons of atmospheric CO2 concentrations, we show that a refinement of the parameterization of turbulent velocity variance and Lagrangian time-scale in STILT is needed to achieve a better match between the Eulerian and the Lagrangian transport at such a high spatial resolution (e.g. 2 and 6 km). Nevertheless, the inter-model differences in simulated CO2 time series for a tall tower observatory at Ochsenkopf in Germany are about a factor of two smaller than the model-data mismatch and about a factor of three smaller than the mismatch between the current global model simulations and the data. Thus suggests that it is reasonable to use STILT as an adjoint model of WRF atmospheric transport.


2017 ◽  
Vol 13 (11) ◽  
pp. 1539-1571 ◽  
Author(s):  
Natalie S. Lord ◽  
Michel Crucifix ◽  
Dan J. Lunt ◽  
Mike C. Thorne ◽  
Nabila Bounceur ◽  
...  

Abstract. Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high-resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled general circulation models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel long-term continuous projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations, which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted, which would not be possible with other modelling techniques on the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling surface air temperature – SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 kyr (thousand years) (by modelling SAT and precipitation). A range of CO2 scenarios are prescribed for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy sea surface temperature (SST) data from different sites and emulator results, finds that relatively similar CO2 concentrations are estimated based on sites at lower latitudes, whereas higher-latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ∼ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost. However, there are uncertainties associated with the approach used, including the inability of the emulator to capture deviations from a quasi-stationary response to the forcing, such as transient adjustments of the deep-ocean temperature and circulation, in addition to its limited range of fixed ice sheet configurations and its requirement for prescribed atmospheric CO2 concentrations.


2018 ◽  
Vol 242 ◽  
pp. 53-61 ◽  
Author(s):  
Romina Beleggia ◽  
Mariagiovanna Fragasso ◽  
Franco Miglietta ◽  
Luigi Cattivelli ◽  
Valeria Menga ◽  
...  

2021 ◽  
Vol 15 (3) ◽  
pp. 1627-1644
Author(s):  
Andrea J. Pain ◽  
Jonathan B. Martin ◽  
Ellen E. Martin ◽  
Åsa K. Rennermalm ◽  
Shaily Rahman

Abstract. Accelerated melting of the Greenland Ice Sheet has increased freshwater delivery to the Arctic Ocean and amplified the need to understand the impact of Greenland Ice Sheet meltwater on Arctic greenhouse gas budgets. We evaluate subglacial discharge from the Greenland Ice Sheet for carbon dioxide (CO2) and methane (CH4) concentrations and δ13C values and use geochemical models to evaluate subglacial CH4 and CO2 sources and sinks. We compare discharge from southwest (a sub-catchment of the Isunnguata Glacier, sub-Isunnguata, and the Russell Glacier) and southern Greenland (Kiattut Sermiat). Meltwater CH4 concentrations vary by orders of magnitude between sites and are saturated with respect to atmospheric concentrations at Kiattut Sermiat. In contrast, meltwaters from southwest sites are supersaturated, even though oxidation reduces CH4 concentrations by up to 50 % during periods of low discharge. CO2 concentrations range from supersaturated at sub-Isunnguata to undersaturated at Kiattut Sermiat. CO2 is consumed by mineral weathering throughout the melt season at all sites; however, differences in the magnitude of subglacial CO2 sources result in meltwaters that are either sources or sinks of atmospheric CO2. At the sub-Isunnguata site, the predominant source of CO2 is organic matter (OM) remineralization. However, multiple or heterogeneous subglacial CO2 sources maintain atmospheric CO2 concentrations at Russell but not at Kiattut Sermiat, where CO2 is undersaturated. These results highlight a previously unrecognized degree of heterogeneity in greenhouse gas dynamics under the Greenland Ice Sheet. Future work should constrain the extent and controls of heterogeneity to improve our understanding of the impact of Greenland Ice Sheet melt on Arctic greenhouse gas budgets, as well as the role of continental ice sheets in greenhouse gas variations over glacial–interglacial timescales.


Sign in / Sign up

Export Citation Format

Share Document