Global Net Primary Production Research Articles

Ocean Biogeochemical Fingerprints of Fast‐Sinking Tunicate and Fish Detritus

AbstractPelagic tunicates (salps, pyrosomes) and fishes generate jelly falls and/or fecal pellets that sink roughly 10 times faster than bulk oceanic detritus, but their impacts on biogeochemical cycles in the ocean interior are poorly understood. Using a coupled physical‐biogeochemical model, we find that fast‐sinking detritus decreased global net primary production and surface export, but increased deep sequestration and transfer efficiency in much of the extratropics and upwelling zones. Fast‐sinking detritus generally decreased total suboxic and hypoxic volumes, reducing a “large oxygen minimum zone (OMZ)” bias common in global biogeochemical models. Newly aerobic regions at OMZ edges exhibited reduced transfer efficiencies in contrast with global tendencies. Reductions in water column denitrification resulting from improved OMZs improved simulated nitrate deficits relative to phosphate. The carbon flux to the benthos increased by 11% with fast‐sinking detritus from fishes and pelagic tunicates, yet simulated benthic fluxes remained on the lower end of observation‐based estimates.

Stressors of primary productivity in the north Indian ocean revealed by satellite, reanalysis and CMIP6 data

Phytoplankton contribute to nearly half of the global net primary productivity (NPP) and are critical for marine ecosystems and global food security. Significant changes in ocean temperature, sea surface salinity (SSS), Mixed Layer Depth (MLD), ocean heat content (OHC), pH and dissolved oxygen (O2) at upper 200 m and that averaged between 200 and 600 m to represent the Oxygen Minimum Zones (OMZ) are identified as the key stressors of marine ecosystem. Here, we examine the changes in stressors and their impact on NPP in the north Indian Ocean (NIO) and its major upwelling zones. There is an increase in SST and freshening in the high saline regions as analysed from the Coupled Model Intercomparison Project 6 (CMIP6) Shared Socioeconomic Pathways (SSP) future (2015–2100) projections, which would reduce NPP in NIO in future. We identify MLD, O2 in 0–200 m and OMZ as the main factors controlling NPP in the present (2003–2014) period, but SSS and O2 in 0–200 m and OMZ in the future SSP scenarios. Contribution analysis indicates the dominant role of O2 in modulating NPP in the future scenarios in the upwelling regions (e.g. Northern AS, Somali Coast and South of Sri Lanka). Future projections also suggest a decreasing trend in pH across different SSP scenarios and its highest negative trend (-0.003 ± 11.12 × 10-5) is observed in SSP 3–7.0. These changes are more evident in Arabian Sea and its adjoining upwelling regions, but a dominant role of SSS in changing NPP is present in BoB and its upwelling regions. The interaction of stressors is manifested well in West Coast and East coast of India, where circulation and currents have dominant roles in the exchanges of temperature and salt fluxes between the regions. Changes in NPP in these dynamically active regions are controlled by all stressors unlike in other regions. Henceforth, our study provides new insights on the impact of spatio-temporal changes in stressors on the primary productivity of NIO.

Dynamic of Grassland Degradation and Its Driving Forces from Climate Variation and Human Activities in Central Asia

Central Asia is one of the most sensitive regions to climate changes in the world and the grassland degradation of this region has attracted considerable concern. Quantifying the driving force of grassland degradation is important for understanding the effects of climate variation and human activities on grassland. In this study, net primary productivity (NPP) was selected as an indicator to quantitatively evaluate the relative role of climate variation and human activities in Central Asia from 2000 to 2020. This study used the global NPP product MOD17A3 as actual NPP and estimated the potential NPP using the Thornthwaite memorial model. The potential NPP and the difference between the potential NPP and actual NPP were used to represent the influence of climate variation and human activities. The grassland degradation or restoration can be demonstrated by the slope of actual NPP (SA). A positive slope value (SA) suggested that restoration occurs, whereas a negative slope value suggested that degradation occurred. The results showed that 23.08% of the total grassland area experienced grassland degradation, whereas 2.51% of the whole grassland underwent grassland restoration. Furthermore, 53.8% of the degraded grassland areas were influenced by climate variation, and 14.5% were caused by human activities. By contrast, the relative roles of climate variation and human activities in grassland restoration were 25% and 47.9%, respectively. The NPP variation also could be calculated by assessing the effects of these factors and the results showed that 55.7% of the NPP decrease was caused by climate variation, whereas 9.6% was a result of human activities. On the contrary, climate variation and human activities resulted in 19.8% and 37.3% of grassland restoration, respectively. Therefore, climate variation was the dominant factor of grassland degradation, and human activities were the main driver of grassland restoration in Central Asia.

A new global oceanic multi-model net primary productivity data product

Abstract. Net primary production of the oceans contributes approximately half of the total global net primary production, and long-term observational records are required to assess any climate-driven changes. The Ocean Colour Climate Change Initiative (OC-CCI) has proven to be robust whilst also being one of the longest records of ocean colour. However, to date, only one primary production algorithm has been applied to this data product, with other algorithms typically applied to single-sensor missions. The data product presented here addresses this issue by applying five algorithms to the OC-CCI data product, which allows the user to interrogate the range of distributions across multiple models and to identify consensus or outliers for their specific region of interest. Outputs are compared to single-sensor data missions, highlighting good overall global agreement, with some small regional discrepancies. Inter-model assessments address the source of these discrepancies, highlighting the choice of the mixed-layer data product as a vital component for accurate primary production estimates. The datasets are published in the Zenodo repository at https://doi.org/10.5281/zenodo.7849935, https://doi.org/10.5281/zenodo.7858590, https://doi.org/10.5281/zenodo.7860491 and https://doi.org/10.5281/zenodo.7861158 (Ryan-Keogh et al., 2023a, b, c, d).

Conserved cobalamin acquisition protein 1 is essential for vitamin B12 uptake in both Chlamydomonas and Phaeodactylum.

Microalgae play an essential role in global net primary productivity and global biogeochemical cycling. Despite their phototrophic lifestyle, over half of algal species depend for growth on acquiring an external supply of the corrinoid vitamin B12 (cobalamin), a micronutrient produced only by a subset of prokaryotic organisms. Previous studies have identified protein components involved in vitamin B12 uptake in bacterial species and humans. However, little is known about its uptake in algae. Here, we demonstrate the essential role of a protein, cobalamin acquisition protein 1 (CBA1), in B12 uptake in Phaeodactylum tricornutum using CRISPR-Cas9 to generate targeted knockouts and in Chlamydomonas reinhardtii by insertional mutagenesis. In both cases, CBA1 knockout lines could not take up exogenous vitamin B12. Complementation of the C. reinhardtii mutants with the wild-type CBA1 gene restored B12 uptake, and regulation of CBA1 expression via a riboswitch element enabled control of the phenotype. When visualized by confocal microscopy, a YFP-fusion with C. reinhardtii CBA1 showed association with membranes. Bioinformatics analysis found that CBA1-like sequences are present in all major eukaryotic phyla. In algal taxa, the majority that encoded CBA1 also had genes for B12-dependent enzymes, suggesting CBA1 plays a conserved role. Our results thus provide insight into the molecular basis of algal B12 acquisition, a process that likely underpins many interactions in aquatic microbial communities.

Biodiversity and Stoichiometric Plasticity Increase Pico‐Phytoplankton Contributions to Marine Net Primary Productivity and the Biological Pump

AbstractEarth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low‐latitude oligotrophic systems harboring pico‐sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon‐climate system. We have added three pico‐phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico‐eukaryote populations that contribute significantly to productivity and export within the subtropical to mid‐latitude transition zone, with the western subtropical regions of each basin supporting the most P‐poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C‐rich/P‐poor stoichiometry, preferentially inhabited by pico‐eukaryotes and diatoms. Collectively, pico‐phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico‐phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.

Reducing spatial resolution increased net primary productivity prediction of terrestrial ecosystems: A Random Forest approach

Net primary production (NPP) is a pivotal component of the terrestrial carbon dynamic, as it directly contributes to the sequestration of atmospheric carbon by vegetation. However, significant variations and uncertainties persist in both the total amount and spatiotemporal patterns of terrestrial NPP, primarily stemming from discrepancies among datasets, modeling approaches, and spatial resolutions. In order to assess the influence of different spatial resolutions on global NPP, we employed a random forest (RF) model using a global observational dataset to predict NPP at 0.05°, 0.25°, and 0.5° resolutions. Our results showed that (1) the RF model performed satisfactorily with modeling efficiencies of 0.53–0.55 for the three respective resolutions; (2) NPP exhibited similar spatial patterns and interannual variation trends at different resolutions; (3) intriguingly, total global NPP varied greatly across different spatial resolutions, amounting 57.3 ± 3.07 for 0.05°, 61.46 ± 3.27 for 0.25°, and 66.5 ± 3.42 Pg C yr−1 for 0.5°. Such differences may be associated with the resolution transformation of the input variables when resampling from finer to coarser resolution, which significantly increased the spatial and temporal variation characteristics, particularly in regions within the southern hemisphere such as Africa, South America, and Australia. Therefore, our study introduces a new concept emphasizing the importance of selecting an appropriate spatial resolution when modeling carbon fluxes, with potential applications in establishing benchmarks for global biogeochemical models.

A paleothermometer for the northern Andes based on C3–C4 grass phytoliths

AbstractGrass-dominated ecosystems cover ~40% of Earth's terrestrial surface, with tropical grasses accounting for ~20% of global net primary productivity. C3 (cool/temperate) and C4 (tropical and subtropical) grass distribution is driven primarily by temperature. In this work, we used phytolith assemblages collected from vegetation plots along an elevation and temperature gradient in the northern Andes (Colombia and Ecuador) to develop a paleothermometer for the region. To accomplish this, we created a transfer function based on the inverse relationship between mean annual temperature (MAT) and the phytolith-based climatic index (Ic), which is the proportion of C3 over C4 grass phytoliths (GSSCP). To evaluate how accurately the index reflects C4–C3 grass abundance in vegetation plots, we compared it with semiquantitative floristic estimates of C4–C3 grass abundance. To further evaluate the 1 − Ic index as a proxy for C4–C3 grass abundance, we compared it with corresponding δ13C values (an independent proxy for C4–C3 vegetation). Results indicate that (1) GSSCP assemblages correctly estimate C4–C3 grass abundance in vegetation plots; (2) the Ic index outperforms the δ13C record in estimating C4–C3 grass abundance, even in open-vegetation types; and (3) our Ic index–based model accurately predicts MAT. This new calibrated proxy will help improve paleotemperature reconstructions in the northern Andes since at least the emergence and spread of C4 grasses in the region during the late Miocene.

GLOBAL ARIDIFICATION AND THE DECLINE OF NPP: A COMMENTARY on Projected increases in global terrestrial net primary productivity loss caused by drought under climate change by Dan Cao, Jiahua Zhang, Jiaqi Han, Tian Zhang, Shanshan Yang, Jingwen Wang, Foyez Ahmed Prodhan, and Fengmei Yao

AbstractNet primary production (NPP) is an important metric of the global carbon cycle, critical for monitoring climate stabilization and for societal economic products. Global NPP should be monitored continuously to adapt societal policy to emerging declines.

Diazotrophy as a key driver of the response of marine net primary productivity to climate change

Abstract. The impact of anthropogenic climate change on marine net primary production (NPP) is a reason for concern because changing NPP will have widespread consequences for marine ecosystems and their associated services. Projections by the current generation of Earth system models have suggested decreases in global NPP in response to future climate change, albeit with very large uncertainties. Here, we make use of two versions of the Institut Pierre-Simon Laplace Climate Model (IPSL-CM) that simulate divergent NPP responses to similar high-emission scenarios in the 21st century and identify nitrogen fixation as the main driver of these divergent NPP responses. Differences in the way N fixation is parameterised in the marine biogeochemical component PISCES (Pelagic Interactions Scheme for Carbon and Ecosystem Studies) of the IPSL-CM versions lead to N-fixation rates that are either stable or double over the course of the 21st century, resulting in decreasing or increasing global NPP, respectively. An evaluation of these two model versions does not help constrain future NPP projection uncertainties. However, the use of a more comprehensive version of PISCES, with variable nitrogen-to-phosphorus ratios as well as a revised parameterisation of the temperature sensitivity of N fixation, suggests only moderate changes in globally averaged N fixation in the 21st century. This leads to decreasing global NPP, in line with the model-mean changes of a recent multi-model intercomparison. Lastly, despite contrasting trends in NPP, all our model versions simulate similar and significant reductions in planktonic biomass. This suggests that projected plankton biomass may be a more robust indicator than NPP of the potential impact of anthropogenic climate change on marine ecosystems across models.

Modeling Global Carbon Costs of Plant Nitrogen and Phosphorus Acquisition.

Most Earth system models (ESMs) do not explicitly represent the carbon (C) costs of plant nutrient acquisition, which leads to uncertainty in predictions of the current and future constraints to the land C sink. We integrate a plant productivity‐optimizing nitrogen (N) and phosphorus (P) acquisition model (fixation & uptake of nutrients, FUN) into the energy exascale Earth system (E3SM) land model (ELM). Global plant N and P uptake are dynamically simulated by ELM‐FUN based on the C costs of nutrient acquisition from mycorrhizae, direct root uptake, retranslocation from senescing leaves, and biological N fixation. We benchmarked ELM‐FUN with three classes of products: ILAMB, a remotely sensed nutrient limitation product, and CMIP6 models; we found significant improvements in C cycle variables, although the lack of more observed nutrient data prevents a comprehensive level of benchmarking. Overall, we found N and P co‐limitation for 80% of land area, with the remaining 20% being either predominantly N or P limited. Globally, the new model predicts that plants invested 4.1 Pg C yr−1 to acquire 841.8 Tg N yr−1 and 48.1 Tg P yr−1 (1994–2005), leading to significant downregulation of global net primary production (NPP). Global NPP is reduced by 20% with C costs of N and 50% with C costs of NP. Modeled and observed nutrient limitation agreement increases when N and P are considered together (r 2 from 0.73 to 0.83).

Projected Increases in Global Terrestrial Net Primary Productivity Loss Caused by Drought Under Climate Change

AbstractUnderstanding present and future impacts of drought on the terrestrial carbon budget is of great significance to the evaluation of terrestrial ecosystem disturbance and terrestrial carbon sink. Here, we evaluate the effect of vegetation net primary productivity (NPP) associated with drought through the difference between the mean NPP in the drought and normal years during a specific time period (30 years). Then, the NPP effects in different vegetation types and climatic zones under baseline stage (1981–2010) and future climate scenarios (RCP2.6, RCP4.5, and RCP8.5) is assessed. The results indicate that the negative NPP extremes are captured in most regions, except for the high‐latitude in the Northern Hemisphere. The NPP loss caused by extreme droughts in 2071–2100 is largest under RCPs, followed by the effects of severe and moderate droughts. Regionally, central United States, southern Africa, central Asia, India, Amazon tropical rainforest, and Australia are projected to experience a significant increase in negative NPP extremes and most of these regions are in arid and semi‐arid and tropical rain forest areas. In contrast, tropical Asia suffers little drought effects. For different vegetation, Evergreen Broadleaf Forest, Closed Shrubland, Open Shrubland, Croplands, and Grassland are the most affected by drought. The largest NPP loss occurs in most part of regions under RCP4.5 scenario, not RCP8.5. Climate change is projected to play the largest role in aggravating the risk of drought‐induced NPP reduction. And meanwhile, the adverse effects of drought on vegetation may be resisted through rational fertilizer utilization and land management in future.

Intercomparison of global terrestrial carbon fluxes estimated by MODIS and Earth system models

Earth system models (ESMs) have been widely used to simulate global terrestrial carbon fluxes, including gross primary production (GPP) and net primary production (NPP). Assessment of such GPP and NPP products can be valuable for understanding the efficacy of certain ESMs in simulating the global carbon cycle and future climate impacts. In this work, we studied the model performance of 22 ESMs participating in the fifth and sixth phases of the Coupled Model Intercomparison Project (CMIP5 and CMIP6) by comparing historical GPP and NPP simulations with satellite data from MODIS and further evaluating potential model improvement from CMIP5 to CMIP6. In CMIP6, the average global total GPP and NPP estimated by the 22 ESMs are 16% and 13% higher than MODIS data, respectively. The multi-model ensembles (MME) of the 22 ESMs can fairly reproduce the spatial distribution, zonal distribution and seasonal variations of both GPP and NPP from MODIS. They perform much better in simulating GPP and NPP for grasslands, wetlands, croplands and other biomes than forests. However, there are noticeable differences among individual ESM simulations in terms of overall fluxes, temporal and spatial flux distributions, and fluxes by biome and region. The MME consistently outperforms all individual models in nearly every respect. Even though several ESMs have been improved in CMIP6 relative to CMIP5, there is still much work to be done to improve individual ESM and overall CMIP performance. Future work needs to focus on more comprehensive model mechanisms and parametrizations, higher resolution and more reasonable coupling of land surface schemes and atmospheric/oceanic schemes.

Elevated CO2 negates O3 impacts on terrestrial carbon and nitrogen cycles

<h2>Summary</h2> Increasing tropospheric concentrations of ozone (e[O₃]) and carbon dioxide (e[CO₂]) profoundly perturb terrestrial ecosystem functions through carbon and nitrogen cycles, affecting beneficial services such as their capacity to combat climate change and provide food. However, the interactive effects of e[O₃] and e[CO₂] on these functions and services remain unclear. Here, we synthesize the results of 810 studies (9,109 observations), spanning boreal to tropical regions around the world, and show that e[O₃] significantly decreases global net primary productivity and food production as well as the capacity of ecosystems to store carbon and nitrogen, which are stimulated by e[CO₂]. More importantly, simultaneous increases in [CO₂] and [O₃] negate or even overcompensate the negative effects of e[O₃] on ecosystem functions and carbon and nitrogen cycles. Therefore, the negative effects of e[O₃] on terrestrial ecosystems would be overestimated if e[CO₂] impacts are not considered, stressing the need for evaluating terrestrial carbon and nitrogen feedbacks to concurrent changes in global atmospheric composition.

Nonlinear interactions of land carbon cycle feedbacks in Earth System Models.

Carbon cycle feedbacks were often quantified through the carbon-concentration and carbon-climate feedbacks with the assumption of no significant interaction between the two feedbacks in most previous studies. Here we calculated the strength of the interactions between the two responses using simulations of models participated in the phase 6 of the Coupled Model Intercomparison Project (CMIP6). We found that the nonlinear interaction contributed 11% of the land-atmosphere carbon exchange on average with large intermodel variation (from -20% to +162%). This nonlinear interaction is largely driven by the pattern of net primary production (NPP), with shifts in heterotrophic respiration that dampen the overall positive interactions from NPP. Photosynthetic rate per unit leaf area alone cannot adequately explain a wide variation of interactions in global NPP simulated by CMIP6models. Plant respiration and processes that regulate leaf area are also important contributors to the interactions. Dominant factors that induce carbon-concentration and carbon-climate interactions are highly variable among models. One of those dominant factors is nutrient limitation. Using additional simulations of ACCESS-ESM1.5 that include both nitrogen and phosphorus limitation, we found that the estimated interactions by ACCESS-ESM1.5 with or without nutrient limitations covered the large intermodel variations among the CMIP6models. It remains largely unknown how nutrient limitation complicates ecosystem's responses to simultaneously CO2 fertilization and warming at the global scale. Our modeling results point to a potential important role of nutrients, especially phosphorus on the nonlinear interactions. Yet, more studies are needed on ecosystem responses to concurrent changes in nutrient availability, atmospheric CO2 concentration, and warming.

Updated respiration routines alter spatio-temporal patterns of carbon cycling in a global land surface model

We updated the routines used to estimate leaf maintenance respiration (MR) in the Energy Land Model (ELM) using a comprehensive global respiration data base. The updated algorithm includes a temperature acclimating base rate, an updated instantaneous temperature response, and new plant functional type specific parameters. The updated MR algorithm resulted in a very large increase in global MR of 16.1 Pg (38%), but the signal was not geographically uniform. The increase was concentrated in the tropics and humid warm-temperate forests. The increase in MR led to large but proportionally smaller decreases in global net primary production (19%) and in average global leaf area index (15%). The effect on global gross primary production (GPP) was a more modest 5.7 Pg (4%). A detailed site level analysis also demonstrated a wide range of effects the updated algorithm can have on the seasonal cycle of GPP. Output from the updated and old models did not differ markedly in how closely they matched a suite of benchmarks. Given the substantial impact on the land surface carbon cycle, a neutral influence on model benchmarks, and better alignment with empirical evidence, an MR algorithm similar to the one presented here should be adopted into ELM.

Recent Trends in SST, Chl‐a, Productivity and Wind Stress in Upwelling and Open Ocean Areas in the Upper Eastern North Atlantic Subtropical Gyre

AbstractThe global upper ocean has been warming during the last decades accompanied with a chlorophyll‐a (Chl‐a) and productivity decrease. Whereas subtropical gyres show similar trends, Eastern Boundary Upwelling Systems are thought to increase in productivity due to increased trade winds. This study analyzes recent trends in sea surface temperature (SST), Chl‐a, net primary production (NPP) and meridional wind stress in the Eastern North Atlantic subtropical gyre (NASE) in order to examine if the global trends can be detected in open ocean and upwelling areas and how the ocean biota responds. Satellite data of such variables of the last 15–40 years were analyzed to calculate mean trends in upwelling areas in the Canary upwelling system and open ocean areas around the Azores, Madeira and the Canary Islands. Our results show significant warming in the area with a maximum of 2.7°C per century for the Azores. Moreover, a general decreasing trend for Chl‐a and NPP seems to be more evident in the permanent upwelling areas, which will be responsible for a loss of 0.13% of the global NPP per century. Our results also highlight a significant expansion of the oceanic desert area of 10% with an increase in unproductive days of up to 84 days in the last 20 years. The competitive relationship of stratification and wind stress in the Canary upwelling system might be a more plausible explanation for the decrease in Chl‐a and NPP in upwelling areas linked to the increase in upwelling favorable wind stress and the surface warming.

An operational overview of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Northeast Pacific field deployment

The goal of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) field campaign is to develop a predictive understanding of the export, fate, and carbon cycle impacts of global ocean net primary production. To accomplish this goal, observations of export flux pathways, plankton community composition, food web processes, and optical, physical, and biogeochemical (BGC) properties are needed over a range of ecosystem states. Here we introduce the first EXPORTS field deployment to Ocean Station Papa in the Northeast Pacific Ocean during summer of 2018, providing context for other papers in this special collection. The experiment was conducted with two ships: a Process Ship, focused on ecological rates, BGC fluxes, temporal changes in food web, and BGC and optical properties, that followed an instrumented Lagrangian float; and a Survey Ship that sampled BGC and optical properties in spatial patterns around the Process Ship. An array of autonomous underwater assets provided measurements over a range of spatial and temporal scales, and partnering programs and remote sensing observations provided additional observational context. The oceanographic setting was typical of late-summer conditions at Ocean Station Papa: a shallow mixed layer, strong vertical and weak horizontal gradients in hydrographic properties, sluggish sub-inertial currents, elevated macronutrient concentrations and low phytoplankton abundances. Although nutrient concentrations were consistent with previous observations, mixed layer chlorophyll was lower than typically observed, resulting in a deeper euphotic zone. Analyses of surface layer temperature and salinity found three distinct surface water types, allowing for diagnosis of whether observed changes were spatial or temporal. The 2018 EXPORTS field deployment is among the most comprehensive biological pump studies ever conducted. A second deployment to the North Atlantic Ocean occurred in spring 2021, which will be followed by focused work on data synthesis and modeling using the entire EXPORTS data set.

Carbon Exchange and Accumulation in an Orinoco High Plains Native Savanna Ecosystem as Measured by Eddy Covariance

Savanna ecosystems cover ∼20% of the total land surface and account for ∼30% of the terrestrial global net primary production. They are also highly sensitive to climate change, since their carbon (C) sink capacity may decline under rising temperatures and irregular rainfall. These responses, which will define the future climate role of the savanna ecosystems, are currently not well understood. The Colombian Orinoco River basin (“Llanos”) natural savannas are being rapidly converted to agriculture. The impact of this transformation on C fluxes and accumulation is not clear. It is thus urgent to understand the Llanos natural savanna ecosystem services, including their C cycle and underlying mechanisms. Here we report and analyze 2 years of measurements of carbon dioxide fluxes from a naturally-restored (secondary) Llanos High Plains savanna ecosystem, using eddy covariance. Meteorological conditions, particularly rainfall, were quite variable during the measurement period. During the first year of measurements, the savanna was a weak carbon source (35 gC m−2), while during the second year, the system was a comparatively strong carbon sink (−273 gC m−2), despite receiving less rainfall than during the first year. As expected, the savanna net ecosystem exchange (NEE) was highly dependent on global solar radiation, soil water content, and ecosystem respiration. We found that after ∼10 days of nominal drought, i.e., with less than ∼5 mm/day of precipitation, the NEE became highly dependent on drought duration. The ecosystem reached a critical condition of low photosynthetic activity after ∼60 days of nominal drought. Based on these findings, we developed and applied a simple standard meteorology-based model that properly reproduced the observations. Our results indicate that a shift to a climate with similar total precipitation but split into extreme dry and wet seasons might eventually suppress the savanna C uptake capacity.

Effects of Sediment Characteristics on Carbon Dioxide Fluxes Based on Interacting Factors in Unvegetated Tidal Flats

The contribution of unvegetated tidal flats to global net primary production is comparable to that of some vegetated coastal habitats. However, compared to carbon flux regulating factors in vegetated habitats, those in unvegetated tidal flats are not well understood, particularly in terms of their cause-effect relationships. Maximum gross primary production (GPPm), net primary production (NPP) and total respiration (TR) during emersion at noon when the irradiance level was at the saturation level for photosynthesis in nine unvegetated tidal flats across the Taiwan Strait in tropical and subtropical regions were determined in winter and summer from 2011 to 2016 to examine the direct and indirect relationships between sediment characteristics and carbon fluxes by using structural equation models (SEM). Most of the GPPm and NPP values were higher in winter than in summer. Conversely, the TR values were higher in summer than in winter. The NPP values at some sites shifted to negative values in summer, indicating the sites were carbon sources. The values of GPPm, TR and NPP for the tidal flats all increased significantly with increasing sediment mud content. The SEM results showed that the sediment mud content directly increased GPPm or indirectly increased GPPm via a compositional shift in benthic microalgae. The sediment mud content also directly increased TR or indirectly increased TR via increased organic matter content. The contribution of benthic microalgal and macrofaunal biomass to TR was relatively minor. This result suggests that primary production was stimulated mainly by the organic matter content rather than by increased microalgal biomass. With the integrated SEM framework, our results showed mechanistic evidence of how sediment mud content acted as a critical factor regulating carbon fluxes in unvegetated tidal flats.