scholarly journals Simulating the atmospheric CO<sub>2</sub> concentration across the heterogeneous landscape of Denmark using a coupled atmosphere–biosphere mesoscale model system

2019 ◽  
Vol 16 (7) ◽  
pp. 1505-1524
Author(s):  
Anne Sofie Lansø ◽  
Thomas Luke Smallman ◽  
Jesper Heile Christensen ◽  
Mathew Williams ◽  
Kim Pilegaard ◽  
...  

Abstract. Although coastal regions only amount to 7 % of the global oceans, their contribution to the global oceanic air–sea CO2 exchange is proportionally larger, with fluxes in some estuaries being similar in magnitude to terrestrial surface fluxes of CO2. Across a heterogeneous surface consisting of a coastal marginal sea with estuarine properties and varied land mosaics, the surface fluxes of CO2 from both marine areas and terrestrial surfaces were investigated in this study together with their impact in atmospheric CO2 concentrations by the usage of a high-resolution modelling framework. The simulated terrestrial fluxes across the study region of Denmark experienced an east–west gradient corresponding to the distribution of the land cover classification, their biological activity and the urbanised areas. Annually, the Danish terrestrial surface had an uptake of approximately −7000 GgC yr−1. While the marine fluxes from the North Sea and the Danish inner waters were smaller annually, with about −1800 and 1300 GgC yr−1, their sizes are comparable to annual terrestrial fluxes from individual land cover classifications in the study region and hence are not negligible. The contribution of terrestrial surfaces fluxes was easily detectable in both simulated and measured concentrations of atmospheric CO2 at the only tall tower site in the study region. Although, the tower is positioned next to Roskilde Fjord, the local marine impact was not distinguishable in the simulated concentrations. But the regional impact from the Danish inner waters and the Baltic Sea increased the atmospheric concentration by up to 0.5 ppm during the winter months.

2014 ◽  
Vol 11 (3) ◽  
pp. 735-747 ◽  
Author(s):  
T. L. Smallman ◽  
M. Williams ◽  
J. B. Moncrieff

Abstract. The coupled numerical weather model WRF-SPA (Weather Research and Forecasting model and Soil-Plant-Atmosphere model) has been used to investigate a 3 yr time series of observed atmospheric CO2 concentrations from a tall tower in Scotland, UK. Ecosystem-specific tracers of net CO2 uptake and net CO2 release were used to investigate the contributions to the tower signal of key land covers within its footprint, and how contributions varied at seasonal and interannual timescales. In addition, WRF-SPA simulated atmospheric CO2 concentrations were compared with two coarse global inversion models, CarbonTrackerEurope and the National Oceanic and Atmospheric Administration's CarbonTracker (CTE-CT). WRF-SPA realistically modelled both seasonal (except post harvest) and daily cycles seen in observed atmospheric CO2 at the tall tower (R2 = 0.67, rmse = 3.5 ppm, bias = 0.58 ppm). Atmospheric CO2 concentrations from the tall tower were well simulated by CTE-CT, but the inverse model showed a poorer representation of diurnal variation and simulated a larger bias from observations (up to 1.9 ppm) at seasonal timescales, compared to the forward modelling of WRF-SPA. However, we have highlighted a consistent post-harvest increase in the seasonal bias between WRF-SPA and observations. Ecosystem-specific tracers of CO2 exchange indicate that the increased bias is potentially due to the representation of agricultural processes within SPA and/or biases in land cover maps. The ecosystem-specific tracers also indicate that the majority of seasonal variation in CO2 uptake for Scotland's dominant ecosystems (forests, cropland and managed grassland) is detectable in observations within the footprint of the tall tower; however, the amount of variation explained varies between years. The between years variation in detectability of Scotland's ecosystems is potentially due to seasonal and interannual variation in the simulated prevailing wind direction. This result highlights the importance of accurately representing atmospheric transport used within atmospheric inversion models used to estimate terrestrial source/sink distribution and magnitude.


2021 ◽  
Author(s):  
Dimitri Bandou ◽  
Patrick Schläfli ◽  
Michael Schwenk ◽  
Guilhem A. Douillet ◽  
Edi Kissling ◽  
...  

&lt;p&gt;Interpretations of the processes leading to the formation of overdeepened valleys, where the bedrock lies well below sea level today, are contested as the overdeepenings have been filled by sediments or host lakes making observations difficult. Here, we combine gravimetric, GNSS (Global Navigation Satellite System) and borehole data within a 3D forward modelling framework (Gravi3D) to assess the 3-D subsurface geometry of such overdeepenings in the Swiss plateau, to the North of the Alps. Gravi3D has two components (PRISMA and BGPoly), which allow to obtain analytically the gravity effect of prisms and polygons (Nagy (1966) and Talwani &amp; Ewing (1960)). PRISMA allows first to estimate the spatial extent of an overdeepening and the density contrast between the overdeepening fill and the bedrock. In contrast, BGPoly is designed to disclose the details of a complex 3-D geometry of an overdeepening fill through an approximation of its shape with polygons. Gravi3D will be open access and is designed for a larger scientific community.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;br&gt;Here, we focus on overdeepenings beneath two valleys, the Aare valley and the G&amp;#252;rbe valley to the South of Bern. In this region, the occurrence of overdeepenings has already been disclosed through drilling, but the details about the geometry have not been elaborated yet. The study region is characterized by three mountain ranges oriented North-South and comprises Burdigalian Upper Marine Molasse bedrock. The G&amp;#252;rbe and Aare valleys in-between are c. 300 m-deep and c. 1 km-wide, where overdeepenings with a &gt;100 m-thick Quaternary fill have already been identified by drilling. The gravity data collected along an 8 km-long profile with stations spaced between 100 and 300 m yield a Bouguer anomaly that ranges from c. -99 to -106 mGal. We relate this anomaly to the regional trend (c. 2 mGal over 8 km) and to the effect of the overdeepenings&amp;#8217; sedimentary fillings (2 &amp;#8211; 4 mGal/km), disclosing a sharp anomaly pattern over the inferred tunnel valleys. The removal of the signal related to the regional trend results in a residual anomaly of c. 1 mGal for the bedrock ridge in-between the valleys (Belpberg mountain), and of -2.65 and -3.56 mGal for the G&amp;#252;rbe and Aare valley overdeepenings, respectively. We observe a steeper gravity gradient for the Eastern flank of both overdeepenings. The use of Nettleton method to model the residual gravity anomaly across Belpberg yields a density of 2.5 g/cm&lt;sup&gt;3&lt;/sup&gt; for the Molasse bedrock. In addition, the estimation of the largest gravity response through the overdeepening fill, calculated with Prisma yields a density value of c. 2.0 &amp;#8211; 2.2 g/cm3 for the Quaternary sediments. As a further information, Prisma predicts a maximum thickness of 140 m for the Quaternary suite beneath the G&amp;#252;rbe valley and at least 200 m beneath the Aare valley. This yields a minimum slope of approximately 18&amp;#176; for the G&amp;#252;rbe overdeepening.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;REFERENCES&lt;/p&gt;&lt;p&gt;Nagy, D.: The gravitational attraction of a right rectangular prism. Geophysics 31, 362&amp;#8211;371, 1966.&lt;/p&gt;&lt;p&gt;Talwani, M., Ewing, M.: Rapid computation of gravitational attraction of three&amp;#8208;dimensional bodies of arbitrary shape. Geophysics 25, 203&amp;#8211;225, 1960.&lt;/p&gt;


2011 ◽  
Vol 4 (3) ◽  
pp. 2047-2080 ◽  
Author(s):  
A. Ganshin ◽  
T. Oda ◽  
M. Saito ◽  
S. Maksyutov ◽  
V. Valsala ◽  
...  

Abstract. We designed a method to simulate atmospheric CO2 concentrations at several continuous observation sites around the globe using surface fluxes at a very high spatial resolution. The simulations presented in this study were performed using a Lagrangian particle dispersion model coupled to a global atmospheric tracer transport model with prescribed global surface CO2 flux maps at a 1 × 1 km resolution. The surface fluxes used in the simulations were prepared by assembling the individual components of terrestrial, oceanic and fossil fuel CO2 fluxes. This experimental setup (i.e., a transport model running at a medium resolution, coupled to a high-resolution Lagrangian particle dispersion model together with global surface fluxes at a very high resolution), which was designed to represent high-frequency variations in atmospheric CO2 concentration, has not been reported at a global scale previously. Two sensitivity experiments were performed: (a) using the global transport model without coupling to the Lagrangian dispersion model, and (b) using the coupled model with a reduced resolution of surface fluxes, in order to evaluate the performance of Eulerian-Lagrangian coupling and the role of high-resolution fluxes in simulating high-frequency variations in atmospheric CO2 concentrations. A correlation analysis between observed and simulated atmospheric CO2 concentrations at selected locations revealed that the inclusion of both Eulerian-Lagrangian coupling and high-resolution fluxes improves the high-frequency simulations of the model. The results highlight the potential of a coupled Eulerian-Lagrangian model in simulating high-frequency atmospheric CO2 concentrations at many locations worldwide. The model performs well in representing observations of atmospheric CO2 concentrations at high spatial and temporal resolutions, especially for coastal sites and sites located close to sources of large anthropogenic emissions. While this study focused on simulations of CO2 concentrations, the model could be used for other atmospheric compounds with known estimated emissions.


Ocean Science ◽  
2010 ◽  
Vol 6 (1) ◽  
pp. 77-89 ◽  
Author(s):  
A. M. Omar ◽  
A. Olsen ◽  
T. Johannessen ◽  
M. Hoppema ◽  
H. Thomas ◽  
...  

Abstract. Data from two Voluntary Observing Ship (VOS) (2005–2007) augmented with data subsets from ten cruises (1987–2005) were used to investigate the spatiotemporal variations of the CO2 fugacity in seawater (fCO2sw) in the North Sea at seasonal and inter-annual time scales. The observed seasonal fCO2sw variations were related to variations in sea surface temperature (SST), biology plus mixing, and air-sea CO2 exchange. Over the study period, the seasonal amplitude in fCO2sw induced by SST changes was 0.4–0.75 times those resulting from variations in biology plus mixing. Along a meridional transect, fCO2sw normally decreased northwards (−12 μatm per degree latitude), but the gradient disappeared/reversed during spring as a consequence of an enhanced seasonal amplitude of fCO2sw in southern parts of the North Sea. Along a zonal transect, a weak gradient (−0.8 μatm per degree longitude) was observed in the annual mean fCO2sw. Annually and averaged over the study area, surface waters of the North Sea were CO2 undersaturated and, thus, a sink of atmospheric CO2. However, during summer, surface waters in the region 55.5–54.5° N were CO2 supersaturated and, hence, a source for atmospheric CO2. Comparison of fCO2sw data acquired within two 1°×1° regions in the northern and southern North Sea during different years (1987, 2001, 2002, and 2005–2007) revealed large interannual variations, especially during spring and summer when year-to-year fCO2sw differences (≈160–200 μatm) approached seasonal changes (≈200–250 μatm). The springtime variations resulted from changes in magnitude and timing of the phytoplankton bloom, whereas changes in SST, wind speed and total alkalinity may have contributed to the summertime interannual fCO2sw differences. The lowest interannual variation (10–50 μatm) was observed during fall and early winter. Comparison with data reported in October 1967 suggests that the fCO2sw growth rate in the central North Sea was similar to that in the atmosphere.


2009 ◽  
Vol 6 (2) ◽  
pp. 1655-1686
Author(s):  
A. M. Omar ◽  
A. Olsen ◽  
T. Johannessen ◽  
M. Hoppema ◽  
H. Thomas ◽  
...  

Abstract. Data from two Voluntary Observing Ship (VOS) (MS Trans Carrier and MV Nuka Arctica), acquired along one zonal and one meridional transect (2005–2007) augmented with data subsets from ten cruises (1987–2005) were used to investigate the spatio-temporal variations of the CO2 fugacity in seawater (fCO2sw) in the North Sea at seasonal and inter-annual time scales. The observed seasonal fCO2sw variations were related to variations in sea surface temperature (SST), biology plus mixing, and air-sea CO2 exchange. Over the study period, the seasonal amplitude in fCO2sw induced by SST changes was 0.4–0.75 times those resulting from variations in biology plus mixing. Along the meridional transect, fCO2sw normally decreased northwards (−12 μatm per degree latitude), but the gradient disappeared/reversed during spring as a consequence of an enhanced seasonal amplitude of fCO2sw in southern parts of the North Sea. Along the zonal transect, a weak gradient (−0.8 μatm per degree longitude) was observed in the mean annual fCO2sw. Annually and averaged over the study area, surface waters of the North Sea were CO2 undersaturated and thus a sink of atmospheric CO2 throughout the year. However, during summer, surface waters in the region 55.5–54.5° N were CO2 supersaturated and, hence, a source for atmospheric CO2. Comparison of fCO2sw data acquired within two 1°×1° regions in the northern and southern North Sea during different years (1987, 2001, 2002, and 2005–2007) revealed large interannual variations, especially during spring and summer when year-to-year fCO2sw differences (≈160–200 μatm) approached seasonal changes (≈200–250 μatm). The springtime variations resulted from changes in magnitude and timing of the phytoplankton bloom, whereas changes in SST, wind speed, and total alkalinity may have contributed to the summertime interannual fCO2sw differences. The lowest interannual variation (10–50 μatm) was observed during fall and early winter. The comparison with data reported in October 1967 suggests that the fCO2sw growth rate in the central North Sea is similar to that in the atmosphere.


2014 ◽  
Vol 14 (23) ◽  
pp. 13281-13293 ◽  
Author(s):  
X. Tian ◽  
Z. Xie ◽  
Y. Liu ◽  
Z. Cai ◽  
Y. Fu ◽  
...  

Abstract. We have developed a novel framework ("Tan-Tracker") for assimilating observations of atmospheric CO2 concentrations, based on the POD-based (proper orthogonal decomposition) ensemble four-dimensional variational data assimilation method (PODEn4DVar). The high flexibility and the high computational efficiency of the PODEn4DVar approach allow us to include both the atmospheric CO2 concentrations and the surface CO2 fluxes as part of the large state vector to be simultaneously estimated from assimilation of atmospheric CO2 observations. Compared to most modern top-down flux inversion approaches, where only surface fluxes are considered as control variables, one major advantage of our joint data assimilation system is that, in principle, no assumption on perfect transport models is needed. In addition, the possibility for Tan-Tracker to use a complete dynamic model to consistently describe the time evolution of CO2 surface fluxes (CFs) and the atmospheric CO2 concentrations represents a better use of observation information for recycling the analyses at each assimilation step in order to improve the forecasts for the following assimilations. An experimental Tan-Tracker system has been built based on a complete augmented dynamical model, where (1) the surface atmosphere CO2 exchanges are prescribed by using a persistent forecasting model for the scaling factors of the first-guess net CO2 surface fluxes and (2) the atmospheric CO2 transport is simulated by using the GEOS-Chem three-dimensional global chemistry transport model. Observing system simulation experiments (OSSEs) for assimilating synthetic in situ observations of surface CO2 concentrations are carefully designed to evaluate the effectiveness of the Tan-Tracker system. In particular, detailed comparisons are made with its simplified version (referred to as TT-S) with only CFs taken as the prognostic variables. It is found that our Tan-Tracker system is capable of outperforming TT-S with higher assimilation precision for both CO2 concentrations and CO2 fluxes, mainly due to the simultaneous estimation of CO2 concentrations and CFs in our Tan-Tracker data assimilation system. A experiment for assimilating the real dry-air column CO2 retrievals (XCO2) from the Japanese Greenhouse Gases Observation Satellite (GOSAT) further demonstrates its potential wide applications.


2012 ◽  
Vol 5 (1) ◽  
pp. 231-243 ◽  
Author(s):  
A. Ganshin ◽  
T. Oda ◽  
M. Saito ◽  
S. Maksyutov ◽  
V. Valsala ◽  
...  

Abstract. We designed a method to simulate atmospheric CO2 concentrations at several continuous observation sites around the globe using surface fluxes at a very high spatial resolution. The simulations presented in this study were performed using the Global Eulerian-Lagrangian Coupled Atmospheric model (GELCA), comprising a Lagrangian particle dispersion model coupled to a global atmospheric tracer transport model with prescribed global surface CO2 flux maps at a 1 × 1 km resolution. The surface fluxes used in the simulations were prepared by assembling the individual components of terrestrial, oceanic and fossil fuel CO2 fluxes. This experimental setup (i.e. a transport model running at a medium resolution, coupled to a high-resolution Lagrangian particle dispersion model together with global surface fluxes at a very high resolution), which was designed to represent high-frequency variations in atmospheric CO2 concentration, has not been reported at a global scale previously. Two sensitivity experiments were performed: (a) using the global transport model without coupling to the Lagrangian dispersion model, and (b) using the coupled model with a reduced resolution of surface fluxes, in order to evaluate the performance of Eulerian-Lagrangian coupling and the role of high-resolution fluxes in simulating high-frequency variations in atmospheric CO2 concentrations. A correlation analysis between observed and simulated atmospheric CO2 concentrations at selected locations revealed that the inclusion of both Eulerian-Lagrangian coupling and high-resolution fluxes improves the high-frequency simulations of the model. The results highlight the potential of a coupled Eulerian-Lagrangian model in simulating high-frequency atmospheric CO2 concentrations at many locations worldwide. The model performs well in representing observations of atmospheric CO2 concentrations at high spatial and temporal resolutions, especially for coastal sites and sites located close to sources of large anthropogenic emissions. While this study focused on simulations of CO2 concentrations, the model could be used for other atmospheric compounds with known estimated emissions.


2013 ◽  
Vol 10 (8) ◽  
pp. 14301-14331
Author(s):  
T. L. Smallman ◽  
M. Williams ◽  
J. B. Moncrieff

Abstract. The Weather Research and Forecasting (WRF) meteorological model has been coupled to the Soil Plant Atmosphere (SPA) terrestrial ecosystem model, hereafter known as WRF-SPA. SPA generates realistic land-atmosphere exchanges through fully coupled hydrological, carbon and energy cycles. Here we have used WRF-SPA to investigate regional scale observations of atmospheric CO2 concentrations made over a multi-annual period from a tall tower in Scotland. WRF-SPA realistically models both seasonal and daily cycles, predicting CO2 at the tall tower (R2 = 0.67, RMSE = 3.5 ppm, bias = 0.58 ppm), indicating realistic transport, and appropriate source sink distribution and magnitude of CO2 exchange. We have highlighted a consistent post harvest increase in model-observation residuals in atmospheric CO2 concentrations. This increase in model-observation residuals post harvest is likely related to a lack of an appropriate representation of uncultivated components (~ 36% of agricultural holding in Scotland) of agricultural land (e.g., hedgerows and forest patches) which continue to photosynthesise after the crop has been harvested. Through the use of ecosystem specific CO2 tracers we have shown that tall tower observations here do not detect a representative fraction of Scotland's ecosystem CO2 uptake. Cropland CO2 uptake is the dominant ecosystem signal detected at the tall tower, consistent with the dominance of cropland in the area surrounding the tower. However cropland is over-represented in the atmospheric CO2 concentrations simulated to be at the tall tower, relative to the simulated surface cropland CO2 uptake. Observations made at the tall tower were able to detect seasonal variation in ecosystem CO2 uptake, however a majority of variation was only detected for croplands. We have found evidence that interannual variation in weather has a greater impact than interannual variation of the simulated land surface CO2 exchange on tall tower observations for the simulated years. This highlights the importance of accurately representing atmospheric transport used within atmospheric inversion models used to estimate terrestrial source/sink distribution and magnitude.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Lahouari Bounoua ◽  
Ping Zhang ◽  
Kurtis Thome ◽  
Jeffrey Masek ◽  
Abdelmounaime Safia ◽  
...  

In terms of the space cities occupy, urbanization appears as a minor land transformation. However, it permanently modifies land’s ecological functions, altering its carbon, energy, and water fluxes. It is therefore necessary to develop a land cover characterization at fine spatial and temporal scales to capture urbanization’s effects on surface fluxes. We develop a series of biophysical vegetation parameters such as the fraction of photosynthetically active radiation, leaf area index, vegetation greenness fraction, and roughness length over the continental US using MODIS and Landsat products for 2001. A 13-class land cover map was developed at a climate modeling grid (CMG) merging the 500 m MODIS land cover and the 30 m impervious surface area from the National Land Cover Database. The landscape subgrid heterogeneity was preserved using fractions of each class from the 500 m and 30 m into the CMG. Biophysical parameters were computed using the 8-day composite Normalized Difference Vegetation Index produced by the North American Carbon Program. In addition to urban impact assessments, this dataset is useful for the computation of surface fluxes in land, vegetation, and urban models and is expected to be widely used in different land cover and land use change applications.


2010 ◽  
Vol 6 (2) ◽  
pp. 489-535 ◽  
Author(s):  
A.-J. Henrot ◽  
L. François ◽  
E. Favre ◽  
M. Butzin ◽  
M. Ouberdous ◽  
...  

Abstract. The Middle Miocene was one of the last warm periods of the Neogene, culminating with the Middle Miocene Climatic Optimum (MMCO, approximatively 17–15 Ma). Several proxy-based reconstructions support warmer and more humid climate during the MMCO. The mechanisms responsible for the warming at MMCO and particulary the role of the atmospheric carbon dioxide CO2 are still highly debated. Here we carried out a series of sensitivity experiments with the model of intermediate complexity Planet Simulator, investigating the contributions of the absence of ice on the continents, the opening of the Central American and Eastern Tethys Seaways, the lowering of the topography on land, the effect of various atmospheric CO2 concentrations and the vegetation retroaction. Our results show that a higher than present-day CO2 concentration is necessary to generate a warmer climate at all latitudes at the Middle Miocene, in agreement with the terrestrial proxy reconstructions which suggest high atmospheric CO2 concentrations at MMCO. Nevertheless, the changes in sea-surface conditions and the lowering of the topography on land also produce significant local warming that may, locally, even be stronger than the CO2 induced temperature increases. The lowering of the topography leads to a more zonal atmospheric circulation and allows the westerly flow to continue over the lowered Plateaus at mid-latitudes. The reduced height of the Tibetan Plateau notably prevents the development of a monsoon-like circulation, whereas the reduction of elevations of the North American and European reliefs strongly increases precipitation from northwestern to eastern Europe. The changes in vegetation cover contributes to maintain and even to intensify the the warm and humid conditions produced by the other factors, suggesting that the vegetation-climate interactions could help to improve the model-data comparison.


Sign in / Sign up

Export Citation Format

Share Document