Interannual Variability In Net Ecosystem Production Research Articles

Drought in the middle growing season inhibited carbon uptake more critical in an anthropogenic shrub ecosystem of Northwest China

Large-scale ecological restoration efforts of the degraded desert steppe ecosystem in northwestern China have effectively combated desertification but have also resulted in the encroachment of shrubs. However, the function of the anthropogenic shrub ecosystem is still unknown under climate change, especially the response of carbon fluxes to drought, which hinders our further understanding of its sustainability in the future. Hence, based on the eddy covariance and environment sensors, continuous carbon and water fluxes and biophysical factors were monitored during 2019–2022 to investigate the characteristics of carbon fluxes and their response to seasonal drought in an anthropogenic shrub ecosystem dominated by Caragana liouana in the southern Mu Us sandy land, northwestern China. The results showed the shrub ecosystem was a stable carbon sink over the 4 years despite large interannual variability, with annual net ecosystem production (NEP) ranging from 54.37 to 150.35 g C·m−2·y−1. Interannual variability in NEP was mainly due to variation in gross primary productivity (GPP), which is driven by the fact that interannual variability in GPP was higher than interannual variability in ecosystem respiration (Reco). While drought inhibited GPP more than Reco, thus weakening the carbon sequestration ability of the ecosystem during the middle growing season, and thereby carbon fluxes in the shrub ecosystem were significantly suppressed by drought in the middle growing season but slightly in the early growing season. Furthermore, the magnitude and variation of carbon fluxes were determined by the maximum apparent photosynthetic capacity of the canopy (Pmax), leaf area index (LAI), and canopy conductance (gc), all of which are limited by water stress during drought. The work suggests that to keep the carbon sink of anthropogenic shrub ecosystems in arid and semi-arid areas, managers should focus on drought-relief measures in the middle growing season that are better adapted to future climate change.

Interannual and Seasonal Drivers of Carbon Cycle Variability Represented by the Community Earth System Model (CESM2).

Earth system models are intended to make long‐term projections, but they can be evaluated at interannual and seasonal time scales. Although the Community Earth System Model (CESM2) showed improvements in a number of terrestrial carbon cycle benchmarks, relative to its predecessor, our analysis suggests that the interannual variability (IAV) in net terrestrial carbon fluxes did not show similar improvements. The model simulated low IAV of net ecosystem production (NEP), resulting in a weaker than observed sensitivity of the carbon cycle to climate variability. Low IAV in net fluxes likely resulted from low variability in gross primary productivity (GPP)—especially in the tropics—and a high covariation between GPP and ecosystem respiration. Although lower than observed, the IAV of NEP had significant climate sensitivities, with positive NEP anomalies associated with warmer and drier conditions in high latitudes, and with wetter and cooler conditions in mid and low latitudes. We identified two dominant modes of seasonal variability in carbon cycle flux anomalies in our fully coupled CESM2 simulations that are characterized by seasonal amplification and redistribution of ecosystem fluxes. Seasonal amplification of net and gross carbon fluxes showed climate sensitivities mirroring those of annual fluxes. Seasonal redistribution of carbon fluxes is initiated by springtime temperature anomalies, but subsequently negative feedbacks in soil moisture during the summer and fall result in net annual carbon losses from land. These modes of variability are also seen in satellite proxies of GPP, suggesting that CESM2 appropriately represents regional sensitivities of photosynthesis to climate variability on seasonal time scales.

Multi-year eddy covariance measurements of net ecosystem exchange in tropical dry deciduous forest of India

Tropical deciduous forests are unique in terms of their geographical distribution, strong seasonality and their contribution to global carbon dynamics, but yet their carbon sequestration potential is poorly sampled. In the current study, we report the carbon balance of 65-year-old tropical dry deciduous forest in central India using the multi-year eddy covariance measurements from November 2011 to May 2019. Over the study period, the forest site was observed to be a net sink of atmospheric CO2 with a mean annual net ecosystem productivity (NEP) of 524 (± 40; ±1 SD across different years) g C m−2 yr−1 with a 233 (±15) day growing season. The NEP was partitioned into gross primary productivity (GPP) of 3358 (± 167) g C m−2 yr−1 and annual carbon loss due to respiration and decomposition (Reco) of 2834 (± 157) g C m−2 yr−1. The ecosystem showed a strong seasonality as a source of carbon during the leaf-off season (March to June), and as a carbon sink during the rest of the year with significant sequestration during the winter season (October to December). The intra-annual analysis suggested that CO2 flux is primarily controlled by canopy greenness with mean Reco/GPP ratio varying from 1.90 (± 0.12) to 0.79 (± 0.04) during leaf-off to leaf-on seasons. The monthly C fluxes in the growing season are found to be strongly correlated to the environmental variables rather than in the leaf-off season, while variability in monthly ecosystem respiration was better explained by air temperature during the leaf-off season than the growing season. Further, the inter-annual variability of NEP was mainly dependent on growing season length and mean annual temperature. The variations in annual GPP and Reco were directly dependent on the increase in mean annual temperature. The eddy covariance-based NEP was complemented by close independent biometric estimates, imparting confidence in our measurements. In conclusion, this tropical dry deciduous forest site is a substantial carbon sink and would add crucial information regarding carbon budgeting of tropical forests in global carbon balance models.

Understanding the role of phenology and summer physiology in controlling net ecosystem production: a multiscale comparison of satellite, PhenoCam and eddy covariance data

Understanding the temporal and spatial variability (SV) of net ecosystem productivity (NEP) is critical for coupling ecosystem carbon (C) cycle and climate system. Previous studies have shown responses of NEP to changes of plant phenology, but impacts of summer physiological status on annual NEP and how this may vary across different ecosystems and spatial scales were largely unknown. Combining large regional satellite derived indices (MODIS), 676 site-year local data (FLUXDATA) covering seven vegetation types, and 57 site-year regional data (PhenoCam), we found that phenological metrics and summer physiological indicators were significantly correlated with their respective gross primary production-based estimates. The interannual variability of NEP was mostly explained by summer physiology than phenology for most ecosystems, while phenology showed a better performance than summer physiology only for grassland sites. In comparison, inconsistent results were derived from three scales when explaining the SV of NEP. Summer physiology showed more potential in interpreting the SV of NEP at large regional scale, but both physiology and the length of growing season exhibited similar performances at local scale. Observations from regional scale were not able to explain NEP, given that the green chromatic coordinate signal cannot track photosynthesis in summer with a high canopy closure. The finding highlights the important role of summer physiology in controlling the C accumulation of terrestrial ecosystems and understanding the responses of summer physiology to environmental drivers is of great significance for improving the simulation accuracy of C sequestration under the global climate change.

Seasonal and Interannual Variability in Net Ecosystem Production of a Subtropical Coastal Lagoon Inferred from Monthly Oxygen Surveys

A central organizing concept in estuarine biogeochemistry is net ecosystem production (NEP). However, estimates of seasonal and interannual variability of whole-system estuarine NEP, which provide insight into how estuaries respond to climatic and anthropogenic forcing, are rare. The main objectives of this study are to (1) determine the seasonal and interannual variability in whole-system NEP of Biscayne Bay, a subtropical, shallow estuary located in southeastern Florida, USA, and (2) determine the potential driving mechanisms of NEP in this estuary. We compute monthly NEP over more than 12 years by constructing the dissolved oxygen budget for the estuary from monthly snapshot oxygen survey data (i.e., collected once per month). High-frequency observations of oxygen in similar subtropical estuaries were used to quantify the error associated with the snapshot monthly sampling. Oxygen air–water exchange and NEP closely balanced each other, with the long-term mean NEP (± 2 standard errors) equal to − 5.3 ± 0.3 mol O2 m−2 year−1, indicating net heterotrophy. Significant seasonality was found, with lowest NEP in September. At monthly time scales, NEP was significantly positively correlated with chlorophyll and total phosphorus and significantly negatively correlated with canal flow. Interannual variability in NEP was substantial, and the bay temporarily shifted from net heterotrophy to net autotrophy after the 2005 Atlantic hurricane season, probably in response to increases in nutrients from runoff and resuspension. These findings show that monthly oxygen surveys can be used to quantify whole-system estuarine NEP and that Biscayne Bay NEP is sensitive to climatic and anthropogenic forcing on seasonal and interannual timescales.

Climate controls over the net carbon uptake period and amplitude of net ecosystem production in temperate and boreal ecosystems

The seasonal and interannual variability of the terrestrial carbon cycle is regulated by the interactions of climate and ecosystem function. However, the key factors and processes determining the interannual variability of net ecosystem productivity (NEP) in different biomes are far from clear. Here, we quantified yearly anomalies of seasonal and annual NEP, net carbon uptake period (CUP), and the maximum daily NEP (NEPmax) in response to climatic variables in 24 deciduous broadleaf forest (DBF), evergreen forest (EF), and grassland (GRA) ecosystems that include at least eight years of eddy covariance observations. Over the 228 site-years studied, interannual variations in NEP were mostly explained by anomalies of CUP and NEPmax. CUP was determined by spring and autumn net carbon uptake phenology, which were sensitive to annual meteorological variability. Warmer spring temperatures led to an earlier start of net carbon uptake activity and higher spring and annual NEP values in DBF and EF, while warmer autumn temperatures in DBF, higher autumn radiation in EF, and more summer and autumn precipitation in GRA resulted in a later ending date of net carbon uptake and associated higher autumn and annual NEP. Anomalies in NEPmax s were determined by summer precipitation in DBF and GRA, and explained more than 50% of variation in summer NEP anomalies for all the three biomes. Results demonstrate the role of meteorological variability in controlling CUP and NEPmax, which in turn help describe the seasonal and interannual variability of NEP.

CO2 exchange and evapotranspiration across dryland ecosystems of southwestern North America.

Global-scale studies suggest that dryland ecosystems dominate an increasing trend in the magnitude and interannual variability of the land CO2 sink. However, such analyses are poorly constrained by measured CO2 exchange in drylands. Here we address this observation gap with eddy covariance data from 25 sites in the water-limited Southwest region of North America with observed ranges in annual precipitation of 100-1000mm, annual temperatures of 2-25°C, and records of 3-10years (150 site-years in total). Annual fluxes were integrated using site-specific ecohydrologic years to group precipitation with resulting ecosystem exchanges. We found a wide range of carbon sink/source function, with mean annual net ecosystem production (NEP) varying from -350 to +330 gCm-2 across sites with diverse vegetation types, contrasting with the more constant sink typically measured in mesic ecosystems. In this region, only forest-dominated sites were consistent carbon sinks. Interannual variability of NEP, gross ecosystem production (GEP), and ecosystem respiration (Reco ) was larger than for mesic regions, and half the sites switched between functioning as C sinks/C sources in wet/dry years. The sites demonstrated coherent responses of GEP and NEP to anomalies in annual evapotranspiration (ET), used here as a proxy for annually available water after hydrologic losses. Notably, GEP and Reco were negatively related to temperature, both interannually within site and spatially across sites, in contrast to positive temperature effects commonly reported for mesic ecosystems. Models based on MODIS satellite observations matched the cross-site spatial pattern in mean annual GEP but consistently underestimated mean annual ET by ~50%. Importantly, the MODIS-based models captured only 20-30% of interannual variation magnitude. These results suggest the contribution of this dryland region to variability of regional to global CO2 exchange may be up to 3-5 times larger than current estimates.

The Australian terrestrial carbon budget

Abstract. This paper reports a study of the full carbon (C-CO2) budget of the Australian continent, focussing on 1990–2011 in the context of estimates over two centuries. The work is a contribution to the RECCAP (REgional Carbon Cycle Assessment and Processes) project, as one of numerous regional studies. In constructing the budget, we estimate the following component carbon fluxes: net primary production (NPP); net ecosystem production (NEP); fire; land use change (LUC); riverine export; dust export; harvest (wood, crop and livestock) and fossil fuel emissions (both territorial and non-territorial). Major biospheric fluxes were derived using BIOS2 (Haverd et al., 2012), a fine-spatial-resolution (0.05°) offline modelling environment in which predictions of CABLE (Wang et al., 2011), a sophisticated land surface model with carbon cycle, are constrained by multiple observation types. The mean NEP reveals that climate variability and rising CO2 contributed 12 ± 24 (1σ error on mean) and 68 ± 15 TgC yr−1, respectively. However these gains were partially offset by fire and LUC (along with other minor fluxes), which caused net losses of 26 ± 4 TgC yr−1 and 18 ± 7 TgC yr−1, respectively. The resultant net biome production (NBP) is 36 ± 29 TgC yr−1, in which the largest contributions to uncertainty are NEP, fire and LUC. This NBP offset fossil fuel emissions (95 ± 6 TgC yr−1) by 38 ± 30%. The interannual variability (IAV) in the Australian carbon budget exceeds Australia's total carbon emissions by fossil fuel combustion and is dominated by IAV in NEP. Territorial fossil fuel emissions are significantly smaller than the rapidly growing fossil fuel exports: in 2009–2010, Australia exported 2.5 times more carbon in fossil fuels than it emitted by burning fossil fuels.

Interannual variability of net ecosystem productivity in forests is explained by carbon flux phenology in autumn

AbstractAimTo investigate the importance of autumn phenology in controlling interannual variability of forest net ecosystem productivity (NEP) and to derive new phenological metrics to explain the interannual variability of NEP.LocationNorth America and Europe.MethodFlux data from nine deciduous broadleaf forests (DBF) and 13 evergreen needleleaf forests (ENF) across North America and Europe (212 site‐years) were used to explore the relationships between the yearly anomalies of annual NEP and several carbon flux based phenological indicators, including the onset/end of the growing season, onset/end of the carbon uptake period, the spring lag (time interval between the onset of growing season and carbon uptake period) and the autumn lag (time interval between the end of the carbon uptake period and the growing season). Meteorological variables, including global shortwave radiation, air temperature, soil temperature, soil water content and precipitation, were also used to explain the phenological variations.ResultsWe found that interannual variability of NEP can be largely explained by autumn phenology, i.e. the autumn lag. While variation in neither annual gross primary productivity (GPP) nor in annual ecosystem respiration (Re) alone could explain this variability, the negative relationship between annual NEP and autumn lag was due to a larger Re/GPP ratio in years with a prolonged autumn lag. For DBF sites, a longer autumn lag coincided with a significant decrease in annual GPP but showed no correlation with annual Re. However, annual GPP was insensitive to a longer autumn lag in ENF sites but annual Re increased significantly.Main conclusionsThese results demonstrate that autumn phenology plays a more direct role than spring phenology in regulating interannual variability of annual NEP. In particular, the importance of respiration may be potentially underestimated in deriving phenological indicators.

Ecosystem carbon exchanges of a subtropical evergreen coniferous plantation subjected to seasonal drought, 2003–2007

Abstract. The behavior of whole ecosystems in humid climate zone in response to drought is less well documented. Continuous measurements of carbon dioxide and water vapor exchange using the eddy covariance technique during 2003–2007 were made above a subtropical evergreen coniferous plantation in Southeast China. Seasonal distribution of precipitation and temperature was inconsistent during the summer, which resulted in droughts frequently. Gross ecosystem photosynthesis (GEP), ecosystem respiration (RE) and net ecosystem production (NEP) showed clear seasonality, with low rates during droughts and in winter. With statistical significance, the monthly GEP and RE increased exponentially with air temperature, but the monthly NEP was less sensitive and only increased linearly with air temperature. However, the function response of monthly NEP and its components (GEP and RE) disappeared under high temperature and low soil moisture. Monthly and yearly RE showed a statistically linear and positive dependence on the corresponding GEP. Annual NEP varied among years from 306 to 430 g C m−2 (mean 395, SD 50). Interannual variability in NEP was controlled primarily by air temperature and the ratio of precipitation and evapotranspiration that affected the RE/GEP ratio, which varied between 0.76 and 0.81 (mean 0.78, SD 0.02).

Carbon sequestration in boreal jack pine stands following harvesting

AbstractA large area of boreal jack pine (Pinus banksiana Lamb.) forest in Canada is recovering from clear‐cut harvesting, and the carbon (C) balance of these regenerating forests remains uncertain. Net ecosystem CO2 exchange was measured using the eddy‐covariance technique at four jack pine sites representing different stages of stand development: three postharvest sites (HJP02, HJP94, and HJP75) and one preharvest site (OJP). The four sites, located in the southern Canadian boreal forest, Saskatchewan, Canada, are typical of low productivity jack pine stands and were 2, 10, 29, and 90 years old in 2004, respectively. Mean annual net ecosystem production (NEP) for 2004 and 2005 was −137±11, 19±16, 73±28, and 22±30 g C m−2 yr−1 at HJP02, HJP94, HJP75 and OJP, respectively, showing the postharvest jack pine stands to be moderate C sources immediately after harvesting, weak sinks at 10 years, moderate C sinks at 30 years, then weak C sinks at 90 years. Mean annual gross ecosystem photosynthesis (GEP) for the 2 years was 96±10, 347±20, 576±34, and 583±35 g C m−2 yr−1 at HJP02, HJP94, HJP75, and OJP, respectively. The ratio of annual ecosystem respiration (R) to annual GEP was 2.51±0.15, 0.95±0.04, 0.87±0.03, and 0.96±0.03. Seasonally, NEP peaked in May or June at all four sites but GEP and R were highest in July. R at a reference soil temperature of 10 °C, ecosystem quantum yield and photosynthetic capacity were lowest for the 2‐year‐old stand. R was most sensitive to soil temperature for the 90‐year‐old stand. The primary source of variability in NEP over the course of succession of the jack pine ecosystem following harvesting was stand age due to the changes in leaf area index. Intersite variability in GEP and R was an order of magnitude greater than interannual variability at OJP. For both young and old stands, GEP had greater interannual variability than R and played a more important role than R in interannual variation in NEP. Based on year‐round flux measurements from 2000 to 2005, the 10‐year stand had larger interannual variability in GEP and R than the 90‐year stand. Interannual variability in NEP was driven primarily by early‐growing‐season temperature and growing‐season length. Photosynthesis played a dominant role in the rapid rise in NEP early in stand development. Late in stand development, however, the subtle decrease in NEP resulted primarily from increasing respiration.

Simulating terrestrial carbon fluxes using the new biosphere model “biosphere model integrating eco‐physiological and mechanistic approaches using satellite data” (BEAMS)

In this study we present a new biosphere model called the Biosphere model integrating Eco‐physiological And Mechanistic approaches using Satellite data (BEAMS). BEAMS provides a new method of calculating the environmental stress affecting plant growth (Stress). Stress is calculated eco‐physiologically using a photosynthesis model and stomatal conductance formulation, providing a more realistic result than previous models. Stress values are used to estimate Gross Primary Production (GPP) estimates via the light use efficiency concept. We used BEAMS, including our new Stress approach, to investigate global spatial and temporal patterns of net primary production (NPP) and net ecosystem production (NEP). BEAMS was run for the years 1982–2000 using global scale satellite and climate data. Comparison of model results with observational measurements at flux sites reveals that GPP values predicted by BEAMS agree with measured GPP. Obtained Stress values were compared with those of MOD17 and CASA; the three methods produce contrasting spatial patterns. Upon comparing predicted and observed NPP, the pattern of NPP for each plant functional type can be adequately estimated. In terms of trend analysis, NPP increased for the years 1982–2000 in most regions. Different NPP trends were observed in Europe, Russia, and northeast Canada than those proposed by Nemani et al. (2003); we attribute these differences to climate‐related processes. Simulated interannual variations in global NEP are similar to results from inverse modeling. A sensitivity study of obtained NEP shows that the interannual variability in NEP is strongly controlled by air temperature, precipitation, CO2, and the fraction of absorbed photosynthetically active radiation.

Continental-scale comparisons of terrestrial carbon sinks estimated from satellite data and ecosystem modeling 1982–1998

A simulation model based on satellite observations of monthly vegetation cover was used to estimate monthly carbon fluxes in terrestrial ecosystems from 1982 to 1998. The NASA–CASA model was driven by vegetation properties derived from the Advanced Very High Resolution Radiometer (AVHRR) and radiative transfer algorithms that were developed for Moderate Resolution Imaging Spectroradiometer (MODIS). For the terrestrial biosphere, predicted net ecosystem production (NEP) flux for atmospheric CO 2 has varied widely between an annual source of −0.9 Pg C per year and a sink of +2.1 Pg C per year. The southern hemisphere tropical zones (SHT, between 0° and 30°S) have a major influence over the predicted global trends in interannual variability of NEP. In contrast, the terrestrial NEP sink for atmospheric CO 2 on the North American (NA) continent has been fairly consistent between +0.2 and +0.3 Pg C per year, except during relatively cool annual periods when continental NEP fluxes are predicted to total to nearly zero. The predicted NEP sink for atmospheric CO 2 over Eurasia (EA) increased notably in the late 1980s and has been fairly consistent between +0.3 and +0.55 Pg C per year since 1988. High correlations can be detected between the El Niño Southern Oscillation (ENSO) and predicted NEP fluxes on the EA continent and for the SHT latitude zones, whereas NEP fluxes for the North American continent as a whole do not correlate strongly with ENSO events over the same time series since 1982. These observations support the hypothesis that regional climate warming has had notable but relatively small-scale impacts on high latitude ecosystem (tundra and boreal) sinks for atmospheric CO 2.

Historical and projected carbon balance of mature black spruce ecosystems across North America: the role of carbon–nitrogen interactions

The role of carbon (C) and nitrogen (N) interactions on sequestration of atmospheric CO2 in black spruce ecosystems across North America was evaluated with the Terrestrial Ecosystem Model (TEM) by applying parameterizations of the model in which C–N dynamics were either coupled or uncoupled. First, the performance of the parameterizations, which were developed for the dynamics of black spruce ecosystems at the Bonanza Creek Long-Term Ecological Research site in Alaska, were evaluated by simulating C dynamics at eddy correlation tower sites in the Boreal Ecosystem Atmosphere Study (BOREAS) for black spruce ecosystems in the northern study area (northern site) and the southern study area (southern site) with local climate data. We compared simulated monthly growing season (May to September) estimates of gross primary production (GPP), total ecosystem respiration (RESP), and net ecosystem production (NEP) from 1994 to 1997 to available field-based estimates at both sites. At the northern site, monthly growing season estimates of GPP and RESP for the coupled and uncoupled simulations were highly correlated with the field-based estimates (coupled: R2= 0.77, 0.88 for GPP and RESP; uncoupled: R2 = 0.67, 0.92 for GPP and RESP). Although the simulated seasonal pattern of NEP generally matched the field-based data, the correlations between field-based and simulated monthly growing season NEP were lower (R2 = 0.40, 0.00 for coupled and uncoupled simulations, respectively) in comparison to the correlations between field-based and simulated GPP and RESP. The annual NEP simulated by the coupled parameterization fell within the uncertainty of field-based estimates in two of three years. On the other hand, annual NEP simulated by the uncoupled parameterization only fell within the field-based uncertainty in one of three years. At the southern site, simulated NEP generally matched field-based NEP estimates, and the correlation between monthly growing season field-based and simulated NEP (R2 = 0.36, 0.20 for coupled and uncoupled simulations, respectively) was similar to the correlations at the northern site. To evaluate the role of N dynamics in C balance of black spruce ecosystems across North America, we simulated historical and projected C dynamics from 1900 to 2100 with a global-based climatology at 0.5° resolution (latitude × longitude) with both the coupled and uncoupled parameterizations of TEM. From analyses at the northern site, several consistent patterns emerge. There was greater inter-annual variability in net primary production (NPP) simulated by the uncoupled parameterization as compared to the coupled parameterization, which led to substantial differences in inter-annual variability in NEP between the parameterizations. The divergence between NPP and heterotrophic respiration was greater in the uncoupled simulation, resulting in more C sequestration during the projected period. These responses were the result of fundamentally different responses of the coupled and uncoupled parameterizations to changes in CO2 and climate.