Climate change will cause a substantial future greenhouse gas release from warming and thawing permafrost-affected soils to the atmosphere enabling a positive feedback mechanism. Increasing the population density of big herbivores in northern high-latitude ecosystems will increase snow density and hence decrease the insulation strength of snow during winter. As a consequence, theoretically 80% of current permafrost-affected soils (<10 m) is projected to remain until 2100 even when assuming a strong warming using the Representative Concentration Pathway 8.5. Importantly, permafrost temperature is estimated to remain below −4 °C on average after increasing herbivore population density. Such ecosystem management practices would be therefore theoretically an important additional climate change mitigation strategy. Our results also highlight the importance of new field experiments and observations, and the integration of fauna dynamics into complex Earth System models, in order to reliably project future ecosystem functions and climate.

Effects of the short-term temporal variability of meteorological variables on soil temperature in northern high-latitude regions have been investigated. For this, a process-oriented land surface model has been driven using an artificially manipulated climate dataset. Short-term climate variability mainly impacts snow depth, and the thermal diffusivity of lichens and bryophytes. These impacts of climate variability on insulating surface layers together substantially alter the heat exchange between atmosphere and soil. As a result, soil temperature is 0.1 to 0.8 degrees C higher when climate variability is reduced. Earth system models project warming of the Arctic region but also increasing variability of meteorological variables and more often extreme meteorological events. Therefore, our results show that projected future increases in permafrost temperature and active-layer thickness in response to climate change will be lower (i) when taking into account future changes in short-term variability of meteorological variables and (ii) when representing dynamic snow and lichen and bryophyte functions in land surface models.

Non-vascular vegetation has been shown to capture considerable quantities of rainfall, which may affect the hydrological cycle and climate at continental scales. However, direct measurements of rainfall interception by non-vascular vegetation are confined to the local scale, which makes extrapolation to the global effects difficult. Here we use a process-based numerical simulation model to show that non-vascular vegetation contributes substantially to global rainfall interception. Inferred average global water storage capacity including non-vascular vegetation was 2.7 mm, which is consistent with field observations and markedly exceeds the values used in land surface models, which average around 0.4 mm. Consequently, we find that the total evaporation of free water from the forest canopy and soil surface increases by 61% when non-vascular vegetation is included, resulting in a global rainfall interception flux that is 22% of the terrestrial evaporative flux (compared with only 12% for simulations where interception excludes non-vascular vegetation). We thus conclude that non-vascular vegetation is likely to significantly influence global rainfall interception and evaporation with consequences for regional-to continental-scale hydrologic cycling and climate.

It is important that climate models can accurately simulate the terrestrial carbon cycle in the Arctic due to the large and potentially labile carbon stocks found in permafrost-affected environments, which can lead to a positive climate feedback, along with the possibility of future carbon sinks from northward expansion of vegetation under climate warming. Here we evaluate the simulation of tundra carbon stocks and fluxes in three land surface schemes that each form part of major Earth system models (JSBACH, Germany; JULES, UK; ORCHIDEE, France). We use a site-level approach in which comprehensive, high-frequency datasets allow us to disentangle the importance of different processes. The models have improved physical permafrost processes and there is a reasonable correspondence between the simulated and measured physical variables, including soil temperature, soil moisture and snow. We show that if the models simulate the correct leaf area index (LAI), the standard C3 photosynthesis schemes produce the correct order of magnitude of carbon fluxes. Therefore, simulating the correct LAI is one of the first priorities. LAI depends quite strongly on climatic variables alone, as we see by the fact that the dynamic vegetation model can simulate most of the differences in LAI between sites, based almost entirely on climate inputs. However, we also identify an influence from nutrient limitation as the LAI becomes too large at some of the more nutrient-limited sites. We conclude that including moss as well as vascular plants is of primary importance to the carbon budget, as moss contributes a large fraction to the seasonal CO2 flux in nutrient-limited conditions. Moss photosynthetic activity can be strongly influenced by the moisture content of moss, and the carbon uptake can be significantly different from vascular plants with a similar LAI. The soil carbon stocks depend strongly on the rate of input of carbon from the vegetation to the soil, and our analysis suggests that an improved simulation of photosynthesis would also lead to an improved simulation of soil carbon stocks. However, the stocks are also influenced by soil carbon burial (e.g. through cryoturbation) and the rate of heterotrophic respiration, which depends on the soil physical state. More detailed below-ground measurements are needed to fully evaluate biological and physical soil processes. Furthermore, even if these processes are well modelled, the soil carbon profiles cannot resemble peat layers as peat accumulation processes are not represented in the models. Thus, we identify three priority areas for model development: (1) dynamic vegetation including (a) climate and (b) nutrient limitation effects; (2) adding moss as a plant functional type; and an (3) improved vertical profile of soil carbon including peat processes.

Nitrous oxide is a strong greenhouse gas and atmospheric ozone-depleting agent which is largely emitted by soils. Recently, lichens and bryophytes have also been shown to release significant amounts of nitrous oxide. This finding relies on ecosystem-scale estimates of net primary productivity of lichens and bryophytes, which are converted to nitrous oxide emissions by empirical relationships between productivity and respiration, as well as between respiration and nitrous oxide release. Here we obtain an alternative estimate of nitrous oxide emissions which is based on a global process-based non-vascular vegetation model of lichens and bryophytes. The model quantifies photosynthesis and respiration of lichens and bryophytes directly as a function of environmental conditions, such as light and temperature. Nitrous oxide emissions are then derived from simulated respiration assuming a fixed relationship between the two fluxes. This approach yields a global estimate of 0.27 (0.19-0.35) (TgN(2)O) year(-1) released by lichens and bryophytes. This is lower than previous estimates but corresponds to about 50% of the atmospheric deposition of nitrous oxide into the oceans or 25% of the atmospheric deposition on land. Uncertainty in our simulated estimate results from large variation in emission rates due to both physiological differences between species and spatial heterogeneity of climatic conditions. To constrain our predictions, combined online gas exchange measurements of respiration and nitrous oxide emissions may be helpful.

The progressive oxygenation of the Earth's atmosphere was pivotal to the evolution of life, but the puzzle of when and how atmospheric oxygen (O-2) first approached modern levels (similar to 21%) remains unresolved. Redox proxy data indicate the deep oceans were oxygenated during 435-392 Ma, and the appearance of fossil charcoal indicates O-2 > 15-17% by 420-400 Ma. However, existing models have failed to predict oxygenation at this time. Here we show that the earliest plants, which colonized the land surface from similar to 470 Ma onward, were responsible for this mid-Paleozoic oxygenation event, through greatly increasing global organic carbon burialthe net long-term source of O-2. We use a trait-based ecophysiological model to predict that cryptogamic vegetation cover could have achieved similar to 30% of today's global terrestrial net primary productivity by similar to 445 Ma. Data from modern bryophytes suggests this plentiful early plant material had a much higher molar C:P ratio (similar to 2,000) than marine biomass (similar to 100), such that a given weathering flux of phosphorus could support more organic carbon burial. Furthermore, recent experiments suggest that early plants selectively increased the flux of phosphorus (relative to alkalinity) weathered from rocks. Combining these effects in a model of long-term biogeochemical cycling, we reproduce a sustained +2% increase in the carbonate carbon isotope (delta C-13) record by similar to 445 Ma, and predict a corresponding rise in O-2 to present levels by 420-400 Ma, consistent with geochemical data. This oxygen rise represents a permanent shift in regulatory regime to one where fire-mediated negative feedbacks stabilize high O-2 levels.

Bryophyte and lichen cover on the forest floor at high latitudes exerts an insulating effect on the ground. In this way, the cover decreases mean annual soil temperature and can protect permafrost soil. Climate change, however, may change bryophyte and lichen cover, with effects on the permafrost state and related carbon balance. It is, therefore, crucial to predict how the bryophyte and lichen cover will react to environmental change at the global scale. To date, current global land surface models contain only empirical representations of the bryophyte and lichen cover, which makes it impractical to predict the future state and function of bryophytes and lichens. For this reason, we integrate a process-based model of bryophyte and lichen growth into the global land surface model JSBACH (Jena Scheme for Biosphere-Atmosphere Coupling in Hamburg). The model simulates bryophyte and lichen cover on upland sites. Wetlands are not included. We take into account the dynamic nature of the thermal properties of the bryophyte and lichen cover and their relation to environmental factors. Subsequently, we compare simulations with and without bryophyte and lichen cover to quantify the insulating effect of the organisms on the soil. We find an average cooling effect of the bryophyte and lichen cover of 2.7K on temperature in the topsoil for the region north of 50 degrees N under the current climate. Locally, a cooling of up to 5.7K may be reached. Moreover, we show that using a simple, empirical representation of the bryophyte and lichen cover without dynamic properties only results in an average cooling of around 0.5 K. This suggests that (a) bryophytes and lichens have a significant impact on soil temperature in high-latitude ecosystems and (b) a process-based description of their thermal properties is necessary for a realistic representation of the cooling effect. The advanced land surface scheme, including a dynamic bryophyte and lichen model, will be the basis for an improved future projection of land-atmosphere heat and carbon exchange.

It has been hypothesized that predecessors of today's bryophytes significantly increased global chemical weathering in the Late Ordovician, thus reducing atmospheric CO2 concentration and contributing to climate cooling and an interval of glaciations. Studies that try to quantify the enhancement of weathering by non-vascular vegetation, however, are usually limited to small areas and low numbers of species, which hampers extrapolating to the global scale and to past climatic conditions. Here we present a spatially explicit modelling approach to simulate global weathering by non-vascular vegetation in the Late Ordovician. We estimate a potential global weathering flux of 2.8 (km(3) rock) yr(-1), defined here as volume of primary minerals affected by chemical transformation. This is around three times larger than today's global chemical weathering flux. Moreover, we find that simulated weathering is highly sensitive to atmospheric CO2 concentration. This implies a strong negative feedback between weathering by non-vascular vegetation and Ordovician climate.