Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land‐atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO₂) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process‐based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century‐long (1901–2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi‐scale Synthesis and Terrestrial Model Intercomparison Project and two observation‐based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 10¹⁵ g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr⁻¹ with a median value of 51 Pg C yr⁻¹ during 2001–2010. The largest uncertainty in SOC stocks exists in the 40–65°N latitude whereas the largest cross‐model divergence in Rh are in the tropics. The modeled SOC change during 1901–2010 ranges from −70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble‐estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO₂ and nitrogen deposition over intact ecosystems increased SOC stocks—even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Quantifying the spatial and temporal patterns of the water lost to the atmosphere through land surface evapotranspiration (ET) is essential for understanding the global hydrological cycle, but remains much uncertain. In this study, we use the Dynamic Land Ecosystem Model to estimate the global terrestrial ET during 2000–2009 and project its changes in response to climate change and increasing atmospheric CO₂ under two IPCC SRES scenarios (A2 and B1) during 2010–2099. Modeled results show a mean annual global terrestrial ET of about 549 (545–552) mm yr⁻¹ during 2000–2009. Relative to the 2000s, global terrestrial ET for the 2090s would increase by 30.7 mm yr⁻¹ (5.6%) and 13.2 mm yr⁻¹ (2.4%) under the A2 and B1 scenarios, respectively. About 60% of global land area would experience increasing ET at rates of over 9.5 mm decade⁻¹ over the study period under the A2 scenario. The Arctic region would have the largest ET increase (16% compared with the 2000s level) due to larger increase in temperature than other regions. Decreased ET would mainly take place in regions like central and western Asia, northern Africa, Australia, eastern South America, and Greenland due to declines in soil moisture and changing rainfall patterns. Our results indicate that warming temperature and increasing precipitation would result in large increase in ET by the end of the 21st century, while increasing atmospheric CO₂ would be responsible for decrease in ET, given the reduction of stomatal conductance under elevated CO₂.

The Chesapeake Bay plays an important role in transforming riverine nutrients before they are exported to the adjacent continental shelf. Although the mean nitrogen budget of the Chesapeake Bay has been previously estimated from observations, uncertainties associated with interannually varying hydrological conditions remain. In this study, a land‐estuarine‐ocean biogeochemical modeling system is developed to quantify Chesapeake riverine nitrogen inputs, within‐estuary nitrogen transformation processes and the ultimate export of nitrogen to the coastal ocean. Model skill was evaluated using extensive in situ and satellite‐derived data, and a simulation using environmental conditions for 2001–2005 was conducted to quantify the Chesapeake Bay nitrogen budget. The 5 year simulation was characterized by large riverine inputs of nitrogen (154 × 10⁹ g N yr⁻¹) split roughly 60:40 between inorganic:organic components. Much of this was denitrified (34 × 10⁹ g N yr⁻¹) and buried (46 × 10⁹ g N yr⁻¹) within the estuarine system. A positive net annual ecosystem production for the bay further contributed to a large advective export of organic nitrogen to the shelf (91 × 10⁹ g N yr⁻¹) and negligible inorganic nitrogen export. Interannual variability was strong, particularly for the riverine nitrogen fluxes. In years with higher than average riverine nitrogen inputs, most of this excess nitrogen (50–60%) was exported from the bay as organic nitrogen, with the remaining split between burial, denitrification, and inorganic export to the coastal ocean. In comparison to previous simulations using generic shelf biogeochemical model formulations inside the estuary, the estuarine biogeochemical model described here produced more realistic and significantly greater exports of organic nitrogen and lower exports of inorganic nitrogen to the shelf.

The magnitude, spatiotemporal patterns, and controls of carbon flux from land to the ocean remain uncertain. Here we applied a process‐based land model with explicit representation of carbon processes in streams and rivers to examine how changes in climate, land conversion, management practices, atmospheric CO₂, and nitrogen deposition affected carbon fluxes from eastern North America to the Atlantic Ocean, specifically the Gulf of Maine (GOM), Middle Atlantic Bight (MAB), and South Atlantic Bight (SAB). Our simulation results indicate that the mean annual fluxes (±1 standard deviation) of dissolved organic carbon (DOC), particulate organic carbon (POC), and dissolved inorganic carbon (DIC) in the past three decades (1980–2008) were 2.37 ± 0.60, 1.06 ± 0.20, and 3.57 ± 0.72 Tg C yr⁻¹, respectively. Carbon export demonstrated substantial spatial and temporal variability. For the region as a whole, the model simulates a significant decrease in riverine DIC fluxes from 1901 to 2008, whereas there were no significant trends in DOC or POC fluxes. In the SAB, however, there were significant declines in the fluxes of all three forms of carbon, and in the MAB subregion, DIC and POC fluxes declined significantly. The only significant trend in the GOM subregion was an increase in DIC flux. Climate variability was the primary cause of interannual variability in carbon export. Land conversion from cropland to forest was the primary factor contributing to decreases in all forms of C export, while nitrogen deposition and fertilizer use, as well as atmospheric CO₂ increases, tended to increase DOC, POC, and DIC fluxes.