Description: This paper provides a comprehensive description of the newest version of the Dynamic Global Vegetation Model with managed Land, LPJmL4. This model simulates – internally consistently – the growth and productivity of both natural and agricultural vegetation as coherently linked through their water, carbon, and energy fluxes. These features render LPJmL4 suitable for assessing a broad range of feedbacks within and impacts upon the terrestrial biosphere as increasingly shaped by human activities such as climate change and land use change. Here we describe the core model structure, including recently developed modules now unified in LPJmL4. Thereby, we also review LPJmL model developments and evaluations in the field of permafrost, human and ecological water demand, and improved representation of crop types. We summarize and discuss LPJmL model applications dealing with the impacts of historical and future environmental change on the terrestrial biosphere at regional and global scale and provide a comprehensive overview of LPJmL publications since the first model description in 2007. To demonstrate the main features of the LPJmL4 model, we display reference simulation results for key processes such as the current global distribution of natural and managed ecosystems, their productivities, and associated water fluxes. A thorough evaluation of the model is provided in a companion paper. By making the model source code freely available at https://gitlab.pik-potsdam.de/lpjml/LPJmL, we hope to stimulate the application and further development of LPJmL4 across scientific communities in support of major activities such as the IPCC and SDG process.

Description: The dynamic global vegetation model LPJmL4 is a process-based model that simulates climate and land use change impacts on the terrestrial biosphere, agricultural production, and the water and carbon cycle. Different versions of the model have been developed and applied to evaluate the role of natural and managed ecosystems in the Earth system and the potential impacts of global environmental change. A comprehensive model description of the new model version, LPJmL4, is provided in a companion paper (Schaphoff et al., 2018c). Here, we provide a full picture of the model performance, going beyond standard benchmark procedures and give hints on the strengths and shortcomings of the model to identify the need for further model improvement. Specifically, we evaluate LPJmL4 against various datasets from in situ measurement sites, satellite observations, and agricultural yield statistics. We apply a range of metrics to evaluate the quality of the model to simulate stocks and flows of carbon and water in natural and managed ecosystems at different temporal and spatial scales. We show that an advanced phenology scheme improves the simulation of seasonal fluctuations in the atmospheric CO2 concentration, while the permafrost scheme improves estimates of carbon stocks. The full LPJmL4 code including the new developments will be supplied open source through https://gitlab.pik-potsdam.de/lpjml/LPJmL. We hope that this will lead to new model developments and applications that improve the model performance and possibly build up a new understanding of the terrestrial biosphere.

Description: Invasive ecosystem engineers (IEE) are potentially one of the most influential types of biological invaders. They are expected to have extensive ecological impacts by altering the physical–chemical structure of ecosystems, thereby changing the rules of existence for a broad range of resident biota. To test the generality of this expectation, we used a global systematic review and meta-analysis to examine IEE effects on the abundance of individual species and communities, biodiversity (using several indices) and ecosystem functions, focusing on marine and estuarine environments. We found that IEE had a significant effect (positive and negative) in most studies testing impacts on individual species, but the overall (cumulative) effect size was small and negative. Many individual studies showed strong IEE effects on community abundance and diversity, but the direction of effects was variable, leading to statistically non-significant overall effects in most categories. In contrast, there was a strong overall effect on most ecosystem functions we examined. IEE negatively affected metabolic functions and primary production, but positively affected nutrient flux, sedimentation and decomposition. We use the results to develop a conceptual model by highlighting pathways whereby IEE impact communities and ecosystem functions, and identify several sources of research bias in the IEE-related invasion literature. Only a few of the studies simultaneously quantified IEE effects on community/diversity and ecosystem functions. Therefore, understanding how IEE may alter biodiversity–ecosystem function relationships should be a primary focus of future studies of invasion biology. Moreover, the clear effects of IEE on ecosystem functions detected in our study suggest that scientists and environmental managers ought to examine how the effects of IEE might be manifested in the services that marine ecosystems provide to humans.

Description: There has been increasing interest in algae-based bioassessment, particularly, trait-based approaches are increasingly suggested. However, the main drivers, especially the contribution of hydrological variables, of species composition, trait composition, and beta diversity of algae communities are less studied. To link species and trait composition to multiple factors (i.e., hydrological variables, local environmental variables, and spatial factors) that potentially control species occurrence/abundance and to determine their relative roles in shaping species composition, trait composition, and beta diversities of pelagic algae communities, samples were collected from a German lowland catchment, where a well-proven ecohydrological modeling enabled to predict long-term discharges at each sampling site. Both trait and species composition showed significant correlations with hydrological, environmental, and spatial variables, and variation partitioning revealed that the hydrological and local environmental variables outperformed spatial variables. A higher variation of trait composition (57.0%) than species composition (37.5%) could be explained by abiotic factors. Mantel tests showed that both species and trait-based beta diversities were mostly related to hydrological and environmental heterogeneity with hydrological contributing more than environmental variables, while purely spatial impact was less important. Our findings revealed the relative importance of hydrological variables in shaping pelagic algae community and their spatial patterns of beta diversities, emphasizing the need to include hydrological variables in long-term biomonitoring campaigns and biodiversity conservation or restoration. A key implication for biodiversity conservation was that maintaining the instream flow regime and keeping various habitats among rivers are of vital importance. However, further investigations at multispatial and temporal scales are greatly needed.

Description: Marine life is controlled by multiple physical and chemical drivers and by diverse ecological processes. Many of these oceanic properties are being altered by climate change and other anthropogenic pressures. Hence, identifying the influences of multifaceted ocean change, from local to global scales, is a complex task. To guide policy-making and make projections of the future of the marine biosphere, it is essential to understand biological responses at physiological, evolutionary and ecological levels. Here, we contrast and compare different approaches to multiple driver experiments that aim to elucidate biological responses to a complex matrix of ocean global change. We present the benefits and the challenges of each approach with a focus on marine research, and guidelines to navigate through these different categories to help identify strategies that might best address research questions in fundamental physiology, experimental evolutionary biology and community ecology. Our review reveals that the field of multiple driver research is being pulled in complementary directions: the need for reductionist approaches to obtain process-oriented, mechanistic understanding and a requirement to quantify responses to projected future scenarios of ocean change. We conclude the review with recommendations on how best to align different experimental approaches to contribute fundamental information needed for science-based policy formulation.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use and land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was 9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1, SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of 0.5 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in EFF was about 1.6 % and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017, ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 405.0±0.1 ppm averaged over 2017. For 2018, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959–2017, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018, 2016, 2015a, b, 2014, 2013)

Description: In climate reanalyses for multi-decadal or longer scales with coupled atmosphere-ocean General Circulation models (CGCMs) it can be assumed that the growth of prediction errors arises chiefly from imprecisely known model parameters, which have a nonlinear relationship with the climate observations (paleoclimate proxies). Also, high-resolution CGCMs for climate analysis are extremely expensive to run, which constrains the applicability of assimilation schemes. In a model framework where we assume that model dynamic parameters account for (nearly) all forecast errors at observation times, we compare two computationally efficient iterative schemes for approximate nonlinear model parameter estimation and joint flux estimation (taking the specific shape of freshwater from melting in the Greenland ice sheet), and its physically consistent state. First, a trivial adaptation of the strong constraint incremental 4D-Var formulation leads to what we refer to as the parameter space iterative extended Kalman smoother (pIKS); a Gauss-Newton scheme. Second, a so-called parameter space fractional Kalman smoother (pFKS) is an alternative controlled-step line search, which can potentially be a more stable approach. While these iterative schemes have been used in data assimilation, we revisit them together within the context of parameter estimation in climate reanalysis, as compared to the more general 4D-Var formulation. Then, the two schemes are evaluated in numerical experiments with a simple 1D energy balance model (Ebm1D) and with a fully-coupled Community Earth System Model (CESM v1.2). Firstly, with Ebm1D the pFKS obtains a cost function similar to the adjoint method with highly reduced computational cost, while an ensemble transform Kalman filter with an m = 60 ensemble size (ETKF60) behaves slightly worse. The pIKS behaves worse than the ETKF60, but an ETKF10 (m = 10) is even worst. Accordingly, with CESM we evaluate the pKFS and the ETKF60 along with an ETKF with Gaussian Anamorphosis (ETKF-GA60). From all the options, the pFKS has the lowest cost function and seems the favored overall option under heavy computational restrictions, but the ETKF obtains better estimates of the flux term.

Description: The Pacific region of Colombia, like many sparsely populated places in developing countries, has been imagined as empty in social terms, and yet full in terms of natural resources and biodiversity. These imaginaries have enabled the creation of frontiers of land and sea control, where the state as well as private and illegal actors have historically dispossessed Afro-descendant and indigenous peoples. This paper contributes to the understanding of territorialisation in the oceans, where political and legal framings of the sea as an open-access public good have neglected the existence of marine social processes. It shows how Afro-descendant communities and non-state actors are required to use the language of resources, rather than socio-cultural attachment, to negotiate state marine territorialisation processes. Drawing on a case study on the Pacific coast of Colombia, we demonstrate that Afro-descendant communities hold local aquatic epistemologies, in which knowledge and the production of space are entangled in fluid and volumetric spatio-temporal dynamics. However, despite the social importance of aquatic environments, they were excluded from Afro-descendants' collective territorial rights in the 1990s. Driven by their local aquatic epistemologies, coastal communities are reclaiming authority over the seascape through the creation of a marine protected area. We argue that they have transformed relations of authority at sea to ensure local access and control, using state institutional instruments to subvert and challenge the legal framing of the sea as an open access public good. As such, this marine protected area represents a place of resistance that ironically subjects coastal communities to disciplinary technologies of conservation.

Description: The seasonal and interannual variability of chlorophyll in the Gulf of Mexico open waters is studied using a three‐dimensional coupled physical‐biogeochemical model. A 5 years hindcast driven by realistic open‐boundary conditions, atmospheric forcings, and freshwater discharges from rivers is performed. The use of recent in situ observations allowed an in‐depth evaluation of the model nutrient and chlorophyll seasonal distributions, including the chlorophyll vertical structure. We find that different chlorophyll patterns of temporal variability coexist in the deep basin which thereby cannot be considered as a homogeneous region with respect to chlorophyll dynamics. A partitioning of the Gulf of Mexico open waters based on the winter chlorophyll concentration increase is then proposed. This partition is basically explained by the amount of nutrients injected into the euphotic layer which is highly constrained by the dynamic of the winter mixed layer. The seasonal and interannual variability appears to be affected by the variability of atmospheric fluxes and mesoscale dynamics (Loop Current eddies in particular). Finally, estimates of primary production in the deep basin are provided.

Description: In this study, high-resolution bathymetric multibeam and optical image data, both obtained within the Belgian manganese (Mn) nodule mining license area by the autonomous underwater vehicle (AUV) Abyss, were combined in order to create a predictive random forests (RF) machine learning model. AUV bathymetry reveals small-scale terrain variations, allowing slope estimations and calculation of bathymetric derivatives such as slope, curvature, and ruggedness. Optical AUV imagery provides quantitative information regarding the distribution (number and median size) of Mn nodules. Within the area considered in this study, Mn nodules show a heterogeneous and spatially clustered pattern, and their number per square meter is negatively correlated with their median size. A prediction of the number of Mn nodules was achieved by combining information derived from the acoustic and optical data using a RF model. This model was tuned by examining the influence of the training set size, the number of growing trees (ntree), and the number of predictor variables to be randomly selected at each node (mtry) on the RF prediction accuracy. The use of larger training data sets with higher ntree and mtry values increases the accuracy. To estimate the Mn-nodule abundance, these predictions were linked to ground-truth data acquired by box coring. Linking optical and hydroacoustic data revealed a nonlinear relationship between the Mn-nodule distribution and topographic characteristics. This highlights the importance of a detailed terrain reconstruction for a predictive modeling of Mn-nodule abundance. In addition, this study underlines the necessity of a sufficient spatial distribution of the optical data to provide reliable modeling input for the RF.