Ozone suppression of carbon uptake by vegetation : A model study of the effect of ozone on carbon uptake and storage in boreal forests in northern Europe

Det er kjent at ozon har skadelig effekt på både mennesker og vegetasjon, i tillegg til å ha en negativ effekt på klimaet gjennom et direkte strålingspådriv. I tillegg har nyere studier vist at ozon også har en inndirekte effekt på klimaet ved å hemme opptak av karbon til vegetasjon. Formålet med dette studiet er å kaste lys over hva dagens, tidligere og fremtidige ozonkosentrasjoner betyr for strålingsbalansen til atmosfæren ved å redusere opptak av CO2 til vegetasjon, med spesielt fokus på boreal skog i nordeuropa. Med dette målet for øyet brukes en regional klimamodell koblet med kjemi (WRF-chem) til å simulere ozonkonsentrasjoner i nordeuropa for året 2009. Disse konsentrasjonene sammenlignes med observasjoner fre EMEP-nettverket, og brukes deretter i en landmodell (NoahMP) uten tilbakekoblinger, som er endret til å inkludere ozoneffekter på planter. NoahMP modellen er validert med målinger fra SMEAR II stasjonen, og brukes til å simulere endringer i total lagret karbon i boreal skog i nordeuropa. I tillegg brukes resultater fra OsloCTM-simuleringer til å produsere konsentrasjoner som er representative for år 1900 og år 2100 under SRES A2 senarioet. Endringene i totalt karbon sammenlignet med simuleringene uten ozoneffekter viser en klar innvirkning av ozon ved dagens konsentrasjoner som resulterer i betydelig reduksjon i totalt lagret karbon. Økningen i landkarbon fra 1900 til 2009 på grunn av økt atmosfærisk CO2 blir betydelig redusert på grunn av økning i ozonkonsentrasjoner, mens 2100 simuleringene viser en redusert effekt av ozon, selv for områder med betydelig økning i ozonkonsentrasjoner som følge av redusert stomata konduktans på grunn av økende CO2 konsentrasjoner

Integration of a Frost Mortality Scheme Into the Demographic Vegetation Model FATES

Frost is damaging to plants when air temperature drops below their tolerance threshold. The set of mechanisms used by cold-tolerant plants to withstand freezing is called “hardening” and typically take place in autumn to protect against winter damage. The recent incorporation of a hardening scheme in the demographic vegetation model FATES opens up the possibility to investigate frost mortality to vegetation. Previously, the hardening scheme was used to improve hydraulic processes in cold-tolerant plants. In this study, we expand upon the existing hardening scheme by implementing hardiness-dependent frost mortality into CLM5.0-FATES to study the impacts of frost on vegetation in temperate and boreal sites from 1950 to 2015. Our results show that the original freezing mortality approach of FATES, where each plant type had a fixed freezing tolerance threshold—an approach common to many other dynamic vegetation models, was restricted to predicting plant type distribution. The main results emerging from the new scheme are a high autumn and spring frost mortality, especially at colder sites, and increasing mid-winter frost mortality due to global warming, especially at warmer sites. We demonstrate that the new frost scheme is a major step forward in dynamically representing vegetation in ESMs by for the first time including a level of frost tolerance that is responding to the environment and includes some level of cost (implicitly) and benefit. By linking hardening and frost mortality in a land surface model, we open new ways to explore the impact of frost events in the context of global warming.publishedVersio

Explicitly modelling microtopography in permafrost landscapes in a land surface model (JULES vn5.4_microtopography)

Microtopography can be a key driver of heterogeneity in the ground thermal and hydrological regime of permafrost landscapes. In turn, this heterogeneity can influence plant communities, methane fluxes, and the initiation of abrupt thaw processes. Here we have implemented a two-tile representation of microtopography in JULES (the Joint UK Land Environment Simulator), where tiles are representative of repeating patterns of elevation difference. Tiles are coupled by lateral flows of water, heat, and redistribution of snow, and a surface water store is added to represent ponding. Simulations are performed of two Siberian polygon sites, (Samoylov and Kytalyk) and two Scandinavian palsa sites (Stordalen and Iškoras). The model represents the observed differences between greater snow depth in hollows vs. raised areas well. The model also improves soil moisture for hollows vs. the non-tiled configuration (“standard JULES”) though the raised tile remains drier than observed. The modelled differences in snow depths and soil moisture between tiles result in the lower tile soil temperatures being warmer for palsa sites, as in reality. However, when comparing the soil temperatures for July at 20 cm depth, the difference in temperature between tiles, or “temperature splitting”, is smaller than observed (3.2 vs. 5.5 ∘C). Polygons display small (0.2 ∘C) to zero temperature splitting, in agreement with observations. Consequently, methane fluxes are near identical (+0 % to 9 %) to those for standard JULES for polygons, although they can be greater than standard JULES for palsa sites (+10 % to 49 %). Through a sensitivity analysis we quantify the relative importance of model processes with respect to soil moisture and temperatures, identifying which parameters result in the greatest uncertainty in modelled temperature. Varying the palsa elevation between 0.5 and 3 m has little effect on modelled soil temperatures, showing that using only two tiles can still be a valid representation of sites with a range of palsa elevations. Mire saturation is heavily dependent on landscape-scale drainage. Lateral conductive fluxes, while small, reduce the temperature splitting by ∼ 1 ∘C and correspond to the order of observed lateral degradation rates in peat plateau regions, indicating possible application in an area-based thaw model

Climate–ecosystem modelling made easy: The Land Sites Platform

Dynamic Global Vegetation Models (DGVMs) provide a state-of-the-art process-based approach to study the complex interplay between vegetation and its physical environment. For example, they help to predict how terrestrial plants interact with climate, soils, disturbance and competition for resources. We argue that there is untapped potential for the use of DGVMs in ecological and ecophysiological research. One fundamental barrier to realize this potential is that many researchers with relevant expertize (ecology, plant physiology, soil science, etc.) lack access to the technical resources or awareness of the research potential of DGVMs. Here we present the Land Sites Platform (LSP): new software that facilitates single-site simulations with the Functionally Assembled Terrestrial Ecosystem Simulator, an advanced DGVM coupled with the Community Land Model. The LSP includes a Graphical User Interface and an Application Programming Interface, which improve the user experience and lower the technical thresholds for installing these model architectures and setting up model experiments. The software is distributed via version-controlled containers; researchers and students can run simulations directly on their personal computers or servers, with relatively low hardware requirements, and on different operating systems. Version 1.0 of the LSP supports site-level simulations. We provide input data for 20 established geo-ecological observation sites in Norway and workflows to add generic sites from public global datasets. The LSP makes standard model experiments with default data easily achievable (e.g., for educational or introductory purposes) while retaining flexibility for more advanced scientific uses. We further provide tools to visualize the model input and output, including simple examples to relate predictions to local observations. The LSP improves access to land surface and DGVM modelling as a building block of community cyberinfrastructure that may inspire new avenues for mechanistic ecosystem research across disciplines.publishedVersio

Projecting circum-Arctic excess-ground-ice melt with a sub-grid representation in the Community Land Model

To address the long-standing underrepresentation of the influences of highly variable ground ice content on the trajectory of permafrost conditions simulated in Earth system models under a warming climate, we implement a sub-grid representation of excess ground ice within permafrost soils using the latest version of the Community Land Model (CLM5). Based on the original CLM5 tiling hierarchy, we duplicate the natural vegetated land unit by building extra tiles for up to three cryostratigraphies with different amounts of excess ice for each grid cell. For the same total amount of excess ice, introducing sub-grid variability in excess-ice contents leads to different excess-ice melting rates at the grid level. In addition, there are impacts on permafrost thermal properties and local hydrology with sub-grid representation. We evaluate this new development with single-point simulations at the Lena River delta, Siberia, where three sub-regions with distinctively different excess-ice conditions are observed. A triple-land-unit case accounting for this spatial variability conforms well to previous model studies for the Lena River delta and displays markedly different dynamics of future excess-ice thaw compared to a single-land-unit case initialized with average excess-ice contents. For global simulations, we prescribed a tiling scheme combined with our sub-grid representation to the global permafrost region using presently available circum-Arctic ground ice data. The sub-grid-scale excess ice produces significant melting of excess ice under a warming climate and enhances the representation of sub-grid variability of surface subsidence on a global scale. Our model development makes it possible to portray more details on the permafrost degradation trajectory depending on the sub-grid soil thermal regime and excess-ice melting, which also shows a strong indication that accounting for excess ice is a prerequisite of a reasonable projection of permafrost thaw. The modeled permafrost degradation with sub-grid excess ice follows the pathway that continuous permafrost transforms into discontinuous permafrost before it disappears, including surface subsidence and talik formation, which are highly permafrost-relevant landscape changes excluded from most land models. Our development of sub-grid representation of excess ice demonstrates a way forward to improve the realism of excess-ice melt in global land models, but further developments require substantially improved global observational datasets on both the horizontal and vertical distributions of excess ground ice.publishedVersio

Projecting circum-Arctic excess-ground-ice melt with a sub-grid representation in the Community Land Model

To address the long-standing underrepresentation of the influences of highly variable ground ice content on the trajectory of permafrost conditions simulated in Earth system models under a warming climate, we implement a sub-grid representation of excess ground ice within permafrost soils using the latest version of the Community Land Model (CLM5). Based on the original CLM5 tiling hierarchy, we duplicate the natural vegetated land unit by building extra tiles for up to three cryostratigraphies with different amounts of excess ice for each grid cell. For the same total amount of excess ice, introducing sub-grid variability in excess-ice contents leads to different excess-ice melting rates at the grid level. In addition, there are impacts on permafrost thermal properties and local hydrology with sub-grid representation. We evaluate this new development with single-point simulations at the Lena River delta, Siberia, where three sub-regions with distinctively different excess-ice conditions are observed. A triple-land-unit case accounting for this spatial variability conforms well to previous model studies for the Lena River delta and displays markedly different dynamics of future excess-ice thaw compared to a single-land-unit case initialized with average excess-ice contents. For global simulations, we prescribed a tiling scheme combined with our sub-grid representation to the global permafrost region using presently available circum-Arctic ground ice data. The sub-grid-scale excess ice produces significant melting of excess ice under a warming climate and enhances the representation of sub-grid variability of surface subsidence on a global scale. Our model development makes it possible to portray more details on the permafrost degradation trajectory depending on the sub-grid soil thermal regime and excess-ice melting, which also shows a strong indication that accounting for excess ice is a prerequisite of a reasonable projection of permafrost thaw. The modeled permafrost degradation with sub-grid excess ice follows the pathway that continuous permafrost transforms into discontinuous permafrost before it disappears, including surface subsidence and talik formation, which are highly permafrost-relevant landscape changes excluded from most land models. Our development of sub-grid representation of excess ice demonstrates a way forward to improve the realism of excess-ice melt in global land models, but further developments require substantially improved global observational datasets on both the horizontal and vertical distributions of excess ground ice

A Tiling Approach to Represent Subgrid Snow Variability in Coupled Land Surface–Atmosphere Models

A mosaic approach to represent subgrid snow variation in a coupled atmosphere–land surface model (WRF–Noah) is introduced and tested. Solid precipitation is scaled in 10 subgrid tiles based on precalculated snow distributions, giving a consistent, explicit representation of variable snow cover and snow depth on subgrid scales. The method is tested in the Weather Research and Forecasting (WRF) Model for southern Norway at 3-km grid spacing, using the subgrid tiling for areas above the tree line. At a validation site in Finse, the modeled transition time from full snow cover to snow-free ground is increased from a few days with the default snow cover fraction formulation to more than 2 months with the tiling approach, which agrees with in situ observations from both digital camera images and surface temperature loggers. This in turn reduces a cold bias at this site by more than 2°C during the first half of July, with the noontime bias reduced from −5° to −1°C. The improved representation of subgrid snow variation also reduces a cold bias found in the reference simulation on regional scales by up to 0.8°C and increases surface energy fluxes (in particular the latent heat flux), and it resulted in up to 50% increase in monthly (June) precipitation in some of the most affected areas. By simulating individual soil properties for each tile, this approach also accounts for a number of secondary effects of uneven snow distribution resulting in different energy and moisture fluxes in different tiles also after the snow has disappeared. This research was originally published in Journal of Hydrometeorology. © 2017 American Meteorological Societ

Regional-scale phytoplankton dynamics and their association with glacier meltwater runoff in Svalbard

Arctic amplification of global warming has accelerated mass loss of Arctic land ice over the past decades and led to increased freshwater discharge into glacier fjords and adjacent seas. Glacier freshwater discharge is typically associated with high sediment load which limits the euphotic depth but may also aid to provide surface waters with essential nutrients, thus having counteracting effects on marine productivity. In situ observations from a few measured fjords across the Arctic indicate that glacier fjords dominated by marine-terminating glaciers are typically more productive than those with only land-terminating glaciers. Here we combine chlorophyll a from satellite ocean color, an indicator of phytoplankton biomass, with glacier meltwater runoff from climatic mass-balance modeling to establish a statistical model of summertime phytoplankton dynamics in Svalbard (mid-June to September). Statistical analysis reveals significant and positive spatiotemporal associations of chlorophyll a with glacier runoff for 7 out of 14 primary hydrological regions but only within 10 km distance from the shore. These seven regions consist predominantly of the major fjord systems of Svalbard. The adjacent land areas are characterized by a wide range of total glacier coverage (35.5%to 81.2 %) and fraction of marine-terminating glacier area (40.2% to 87.4 %). We find that an increase in specific glacier-runoff rate of 10mm water equivalent per 8 d period raises summertime chlorophyll a concentrations by 5.2% to 20.0%, depending on the region. During the annual peak discharge we estimate that glacier runoff increases chlorophyll a by 13.1% to 50.2% compared to situations with no runoff. This suggests that glacier runoff is an important factor sustaining summertime phytoplankton production in Svalbard fjords, in line with findings from several fjords in Greenland. In contrast, for regions bordering open coasts, and beyond 10 km distance from the shore, we do not find significant associations of chlorophyll a with runoff. In these regions, physical ocean and sea-ice variables control chlorophyll a, pointing at the importance of a late sea-ice breakup in northern Svalbard, as well as the advection of Atlantic water masses along the West Spitsbergen Current for summertime phytoplankton dynamics. Our method allows for the investigation and monitoring of glacier-runoff effects on primary production throughout the summer season and is applicable on a pan-Arctic scale, thus complementing valuable but scarce in situ measurements in both space and time

Regional-scale phytoplankton dynamics and their association with glacier meltwater runoff in Svalbard

Arctic amplification of global warming has accelerated mass loss of Arctic land ice over the past decades and led to increased freshwater discharge into glacier fjords and adjacent seas. Glacier freshwater discharge is typically associated with high sediment load which limits the euphotic depth but may also aid to provide surface waters with essential nutrients, thus having counteracting effects on marine productivity. In situ observations from a few measured fjords across the Arctic indicate that glacier fjords dominated by marine-terminating glaciers are typically more productive than those with only land-terminating glaciers. Here we combine chlorophyll a from satellite ocean color, an indicator of phytoplankton biomass, with glacier meltwater runoff from climatic mass-balance modeling to establish a statistical model of summertime phytoplankton dynamics in Svalbard (mid-June to September). Statistical analysis reveals significant and positive spatiotemporal associations of chlorophyll a with glacier runoff for 7 out of 14 primary hydrological regions but only within 10 km distance from the shore. These seven regions consist predominantly of the major fjord systems of Svalbard. The adjacent land areas are characterized by a wide range of total glacier coverage (35.5%to 81.2 %) and fraction of marine-terminating glacier area (40.2% to 87.4 %). We find that an increase in specific glacier-runoff rate of 10mm water equivalent per 8 d period raises summertime chlorophyll a concentrations by 5.2% to 20.0%, depending on the region. During the annual peak discharge we estimate that glacier runoff increases chlorophyll a by 13.1% to 50.2% compared to situations with no runoff. This suggests that glacier runoff is an important factor sustaining summertime phytoplankton production in Svalbard fjords, in line with findings from several fjords in Greenland. In contrast, for regions bordering open coasts, and beyond 10 km distance from the shore, we do not find significant associations of chlorophyll a with runoff. In these regions, physical ocean and sea-ice variables control chlorophyll a, pointing at the importance of a late sea-ice breakup in northern Svalbard, as well as the advection of Atlantic water masses along the West Spitsbergen Current for summertime phytoplankton dynamics. Our method allows for the investigation and monitoring of glacier-runoff effects on primary production throughout the summer season and is applicable on a pan-Arctic scale, thus complementing valuable but scarce in situ measurements in both space and time