An open and extensible framework for spatially explicit land use change modelling in R: the lulccR package (0.1.0)

Land use change has important consequences for biodiversity and the sustainability of ecosystem services, as well as for global environmental change. Spatially explicit land use change models improve our understanding of the processes driving change and make predictions about the quantity and location of future and past change. Here we present the lulccR package, an object-oriented framework for land use change modelling written in the R programming language. The contribution of the work is to resolve the following limitations associated with the current land use change modelling paradigm: (1) the source code for model implementations is frequently unavailable, severely compromising the reproducibility of scientific results and making it impossible for members of the community to improve or adapt models for their own purposes; (2) ensemble experiments to capture model structural uncertainty are difficult because of fundamental differences between implementations of different models; (3) different aspects of the modelling procedure must be performed in different environments because existing applications usually only perform the spatial allocation of change. The package includes a stochastic ordered allocation procedure as well as an implementation of the widely used CLUE-S algorithm. We demonstrate its functionality by simulating land use change at the Plum Island Ecosystems site, using a dataset included with the package. It is envisaged that lulccR will enable future model development and comparison within an open environment

Effects of Uncertainty in Climate Inputs on Simulated Evapotranspiration and Runoff in the Western Arctic

Hydrological models require accurate precipitation and air temperature inputs in order to adequately depict water fluxes and storages across Arctic regions. Biases such as gauge undercatch, as well as uncertainties in numerical weather prediction reanalysis data that propagate through water budget models, limit the ability to accurately model the terrestrial arctic water cycle. A hydrological model forced with three climate datasets and three methods of estimating potential evapotranspiration (PET) was used to better understand the impact of these processes on simulated water fluxes across the Western Arctic Linkage Experiment (WALE) domain. Climate data were drawn from the NCEP–NCAR reanalysis (NNR) (NCEP1), a modified version of the NNR (NCEP2), and the Willmott–Matsuura (WM) dataset. PET methods applied in the model were Hamon, Penman–Monteith, and Penman–Monteith using adjusted vapor pressure data. High vapor pressures in the NNR lead to low simulated evapotranspiration (ET) in model runs using the Penman–Monteith PET method, resulting in increased runoff. Annual ET derived from simulations using Penman–Monteith PET was half the magnitude of ET simulated when the Hamon method was used. Adjustments made to the reanalysis vapor pressure data increased the simulated ET flux, reducing simulated runoff. Using the NCEP2 or WM climate data, along with the Penman–Monteith PET function, results in agreement to within 7% between the simulated and observed runoff across the Yukon River basin. The results reveal the high degree of uncertainty present in climate data and the range of water fluxes generated from common model drivers. This suggests the need for thorough evaluations of model requirements and potential biases in forcing data, as well as corroborations with observed data, in all efforts to simulate arctic water balances

Modelling understorey dynamics in temperate forests under global change : challenges and perspectives

The understorey harbours a substantial part of vascular plant diversity in temperate forests and plays an important functional role, affecting ecosystem processes such as nutrient cycling and overstorey regeneration. Global change, however, is putting these understorey communities on trajectories of change, potentially altering and reducing their functioning in the future. Developing mitigation strategies to safeguard the diversity and functioning of temperate forests in the future is challenging and requires improved predictive capacity. Process-based models that predict understorey community composition over time, based on first principles of ecology, have the potential to guide mitigation endeavours but such approaches are rare. Here, we review fourteen understorey modelling approaches that have been proposed during the last three decades. We evaluate their inclusion of mechanisms that are required to predict the impact of global change on understorey communities. We conclude that none of the currently existing models fully accounts for all processes that we deem important based on empirical and experimental evidence. Based on this review, we contend new models are needed to project the complex impacts of global change on forest understoreys. Plant functional traits should be central to such future model developments, as they drive community assembly processes and provide valuable information on the functioning of the understorey. Given the important role of the overstorey, a coupling of understorey models to overstorey models will be essential to predict the impact of global change on understorey composition and structure, and how it will affect the functioning of temperate forests in the future

The tension between fire risk and carbon storage: evaluating U.S. carbon and fire management strategies through ecosystem models

Fire risk and carbon storage are related environmental issues because fire reduction results in carbon storage through the buildup of woody vegetation, and stored carbon is a fuel for fires. The sustainability of the U.S. carbon sink and the extent of fire activity in the next 100 yr depend in part on the type and effectiveness of fire reduction employed. Previous studies have bracketed the range of dynamics from continued fire reduction to the complete failure of fire reduction activities. To improve these estimates, it is necessary to explicitly account for fire reduction in terrestrial models. A new fire reduction submodel that estimates the spatiotemporal pattern of reduction across the United States was developed using gridded data on biomass, climate, land-use, population, and economic factors. To the authors’ knowledge, it is the first large-scale, gridded fire model that explicitly accounts for fire reduction. The model was calibrated to 1° × 1° burned area statistics [Global Burnt Area 2000 Project (GBA-2000)] and compared favorably to three important diagnostics. The model was then implemented in a spatially explicit ecosystem model and used to analyze 1620 scenarios of future fire risk and fire reduction strategies. Under scenarios of climate change and urbanization, burned area and carbon emissions both increased in scenarios where fire reduction efforts were not adjusted to match new patterns of fire risk. Fuel reducing management strategies reduced burned area and fire risk, but also limited carbon storage. These results suggest that to promote carbon storage and minimize fire risk in the future, fire reduction efforts will need to be increased and spatially adjusted and will need to employ a mixture of fuel-reducing and non-fuel-reducing strategies

Advances in upscaling of eddy covariance measurements of carbon and water fluxes

Eddy covariance flux towers provide continuous measurements of ecosystem-level net exchange of carbon, water, energy, and other trace gases between land surface and the atmosphere. The upscaling of flux observations from towers to broad regions provides a new and independent approach for quantifying these fluxes over regions, continents, or the globe. The seven contributions of this special section reflect the most recent advances in the upscaling of fluxes from towers to these broad regions. The section mainly stems from presentations at the recent North American Carbon Program (NACP), FLUXNET, and AGU meetings. These studies focus on different aspects of upscaling: (1) assessing the representativeness of flux networks; (2) upscaling fluxes from towers to broad spatial scales; (3) examining the magnitude, distribution, and interannual variability of fluxes over regions, continents, or the globe; and (4) evaluating the impacts of spatial heterogeneity and parameter variability on flux estimates. Collectively, this special issue provides a timely update on upscaling science and also generates gridded flux data that can be used for model evaluations. Future upscaling studies are expected to advance toward incorporating the impacts of disturbance on ecosystem carbon dynamics, quantifying uncertainties associated with gridded flux estimates, and comparing various upscaling methods and the resulting gridded flux fields

Heat risk assessment for the Brussels capital region under different urban planning and greenhouse gas emission scenarios

Urban residents are exposed to higher levels of heat stress in comparison to the rural population. As this phenomenon could be enhanced by both global greenhouse gas emissions (GHG) and urban expansion, urban planners and policymakers should integrate both in their assessment. One way to consider these two concepts is by using urban climate models at a high resolution. In this study, the influence of urban expansion and GHG emission scenarios is evaluated at 100 m spatial resolution for the city of Brussels (Belgium) in the near (2031-2050) and far (2081-2100) future. Two possible urban planning scenarios (translated into local climate zones, LCZs) in combination with two representative concentration pathways (RCPs 4.5 and 8.5) have been implemented in the urban climate model UrbClim. The projections show that the influence of GHG emissions trumps urban planning measures in each period. In the near future, no large differences are seen between the RCP scenarios; in the far future, both heat stress and risk values are twice as large for RCP 8.5 compared to RCP 4.5. Depending on the GHG scenario and the LCZ type, heat stress is projected to increase by a factor of 10 by 2090 compared to the present-day climate and urban planning conditions. The imprint of vulnerability and exposure is clearly visible in the heat risk assessment, leading to very high levels of heat risk, most notably for the North Western part of the Brussels Capital Region. The results demonstrate the need for mitigation and adaptation plans at different policy levels that strive for lower GHG emissions and the development of sustainable urban areas safeguarding livability in cities

Shocks in coupled socio-ecological systems: what are they and how can we model them?

Coupled socio-ecological systems (SES) are complex systems characterized by self-organization, non-linearities, interactions among heterogeneous elements within each subsystem, and feedbacks across scales and among subsystems. When such a system experiences a shock or a crisis, the consequences are difficult to predict. In this paper we first define what a shock or a crisis means for SES. Depending on where the system boundary is drawn, shocks can be seen as exogenous or endogenous. For example, human intervention in environmental systems could be seen as exogenous, but endogenous in a socio-environmental system. This difference in the origin and nature of shocks has certain consequences for coupled SES and for policies to ameliorate negative consequences of shocks. Having defined shocks, the paper then focuses on modelling challenges when studying shocks in coupled SES. If we are to explore, study and predict the responses of coupled SES to shocks, the models used need to be able to accommodate (exogenous) or produce (endogenous) a shock event. Various modelling choices need to be made. Specifically, the ‘sudden’ aspect of a shock suggests the time period over which an event claimed to be a shock occurred might be ‘quick’. What does that mean for a discrete event model? Turning to magnitude, what degree of change (in a variable or set of variables) is required for the event to be considered a shock? The ‘surprising’ nature of a shock means that none of the agents in the model should expect the shock to happen, but may need rules enabling them to generate behaviour in exceptional circumstances. This requires a certain design of the agents’ decision-making algorithms, their perception of a shock, memory of past events and formation of expectations, and the information available to them during the time the shock occurred

Agricultural Production and Externalities Simulator (APES) prototype to be used in Prototype 1 of SEAMLESS-IF

Production Economics,