ALBERT

Description: We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

Description: The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP’s biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948– 2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, Cabio T and 14Cabio T , to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C03019, doi:10.1029/2007JC004153.

Description: Estimates of temporal trends in oceanic anthropogenic carbon dioxide (CO2) rely on the ability of empirical methods to remove the large natural variability of the ocean carbon system. A coupled carbon-climate model is used to evaluate these empirical methods. Both the ΔC* and multiple linear regression (MLR) techniques reproduce the predicted increase in dissolved inorganic carbon for the majority of the ocean and have similar average percent errors for decadal differences (24.1% and 25.5%, respectively). However, this study identifies several regions where these methods may introduce errors. Of particular note are mode and deep water formation regions, where changes in air-sea disequilibrium and structure in the MLR residuals introduce errors. These results have significant implications for decadal repeat hydrography programs, indicating the need for subannual sampling in certain regions of the oceans in order to better constrain the natural variability in the system and to robustly estimate the intrusion of anthropogenic CO2.

Description: We would like to acknowledge funding from NSF (OCE02-23869), NCAR, the WHOI Ocean Climate Institute, a Linden Earth Systems Graduate Fellowship (MIT), and a National Defense Science and Engineering Graduate Fellowship. NCAR is sponsored by the National Science Foundation. R.W. is supported by the Office of Oceanic and Atmospheric Research at NOAA.

Description: Author Posting. © American Geophysical Union, 2015. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 29 (2015): 744–759, doi:10.1002/2014GB005079.

Description: Improved constraints on carbon cycle responses to climate change are needed to inform mitigation policy, yet our understanding of how these responses may evolve after 2100 remains highly uncertain. Using the Community Earth System Model (v1.0), we quantified climate-carbon feedbacks from 1850 to 2300 for the Representative Concentration Pathway 8.5 and its extension. In three simulations, land and ocean biogeochemical processes experienced the same trajectory of increasing atmospheric CO2. Each simulation had a different degree of radiative coupling for CO2 and other greenhouse gases and aerosols, enabling diagnosis of feedbacks. In a fully coupled simulation, global mean surface air temperature increased by 9.3 K from 1850 to 2300, with 4.4 K of this warming occurring after 2100. Excluding CO2, warming from other greenhouse gases and aerosols was 1.6 K by 2300, near a 2 K target needed to avoid dangerous anthropogenic interference with the climate system. Ocean contributions to the climate-carbon feedback increased considerably over time and exceeded contributions from land after 2100. The sensitivity of ocean carbon to climate change was found to be proportional to changes in ocean heat content, as a consequence of this heat modifying transport pathways for anthropogenic CO2 inflow and solubility of dissolved inorganic carbon. By 2300, climate change reduced cumulative ocean uptake by 330 Pg C, from 1410 Pg C to 1080 Pg C. Land fluxes similarly diverged over time, with climate change reducing stocks by 232 Pg C. Regional influence of climate change on carbon stocks was largest in the North Atlantic Ocean and tropical forests of South America. Our analysis suggests that after 2100, oceans may become as important as terrestrial ecosystems in regulating the magnitude of the climate-carbon feedback.

Description: We are grateful for support from the U.S. Department of Energy Office of Science and the National Science Foundation (NSF). J.T.R. and F.H. received support from the Regional and Global Climate Modeling Program in the Climate and Environmental Sciences Division of the Biological and Environmental Research (BER) Program in the U.S. Department of Energy Office of Science. J.T.R., K.L., E.M., W.F., J.K.M., S.C.D., and N.N.M. received funding from the NSF project “Collaborative Research: Improved Regional and Decadal Predictions of the Carbon Cycle“ (AGS-1048827, AGS-1021776, and AGS-1048890). The Community Earth System Modeling project receives support from both NSF and BER.

Description: 2015-12-02

Description: © 2009 The Authors. This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 6 (2009): 2099-2120, doi:10.5194/bg-6-2099-2009

Description: Inclusion of fundamental ecological interactions between carbon and nitrogen cycles in the land component of an atmosphere-ocean general circulation model (AOGCM) leads to decreased carbon uptake associated with CO2 fertilization, and increased carbon uptake associated with warming of the climate system. The balance of these two opposing effects is to reduce the fraction of anthropogenic CO2 predicted to be sequestered in land ecosystems. The primary mechanism responsible for increased land carbon storage under radiatively forced climate change is shown to be fertilization of plant growth by increased mineralization of nitrogen directly associated with increased decomposition of soil organic matter under a warming climate, which in this particular model results in a negative gain for the climate-carbon feedback. Estimates for the land and ocean sink fractions of recent anthropogenic emissions are individually within the range of observational estimates, but the combined land plus ocean sink fractions produce an airborne fraction which is too high compared to observations. This bias is likely due in part to an underestimation of the ocean sink fraction. Our results show a significant growth in the airborne fraction of anthropogenic CO2 emissions over the coming century, attributable in part to a steady decline in the ocean sink fraction. Comparison to experimental studies on the fate of radio-labeled nitrogen tracers in temperate forests indicates that the model representation of competition between plants and microbes for new mineral nitrogen resources is reasonable. Our results suggest a weaker dependence of net land carbon flux on soil moisture changes in tropical regions, and a stronger positive growth response to warming in those regions, than predicted by a similar AOGCM implemented without land carbon-nitrogen interactions. We expect that the between-model uncertainty in predictions of future atmospheric CO2 concentration and associated anthropogenic climate change will be reduced as additional climate models introduce carbon-nitrogen cycle interactions in their land components.

Description: This work was supported in part by NASA Earth Science Enterprise, Terrestrial Ecology Program, grant #W19,953 to P. E. Thornton. Support was provided by the National Center for Atmospheric Research (NCAR) through the NCAR Community Climate System Modeling program, and through the NCAR Biogeosciences program. Additional support was provided by the US Department of Energy, Office of Science, Office of Biological and Environmental Research. I. Fung, S. Doney, N. Mahowald, and J. Randerson acknowledge support from National Science Foundation, Atmospheric Sciences Division, through the Carbon and Water Initiative.

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 18 (2004): GB4028, doi:10.1029/2004GB002220.

Description: A global three-dimensional marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary and export production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed high nitrate, low chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 55 Tg N. Diazotrophs directly account for a small fraction of primary production (0.5%) but indirectly support 10% of primary production and 8% of sinking particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in the upper ocean is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, owing to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.

Description: This work was funded by NSF grant OCE-0222033 and the National Center for Atmospheric Research.

Description: Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 21 (2007): GB3007, doi:10.1029/2006GB002857.

Description: Results are presented of export production, dissolved organic matter (DOM) and dissolved oxygen simulated by 12 global ocean models participating in the second phase of the Ocean Carbon-cycle Model Intercomparison Project. A common, simple biogeochemical model is utilized in different coarse-resolution ocean circulation models. The model mean (±1σ) downward flux of organic matter across 75 m depth is 17 ± 6 Pg C yr−1. Model means of globally averaged particle export, the fraction of total export in dissolved form, surface semilabile dissolved organic carbon (DOC), and seasonal net outgassing (SNO) of oxygen are in good agreement with observation-based estimates, but particle export and surface DOC are too high in the tropics. There is a high sensitivity of the results to circulation, as evidenced by (1) the correlation of surface DOC and export with circulation metrics, including chlorofluorocarbon inventory and deep-ocean radiocarbon, (2) very large intermodel differences in Southern Ocean export, and (3) greater export production, fraction of export as DOM, and SNO in models with explicit mixed layer physics. However, deep-ocean oxygen, which varies widely among the models, is poorly correlated with other model indices. Cross-model means of several biogeochemical metrics show better agreement with observation-based estimates when restricted to those models that best simulate deep-ocean radiocarbon. Overall, the results emphasize the importance of physical processes in marine biogeochemical modeling and suggest that the development of circulation models can be accelerated by evaluating them with marine biogeochemical metrics.

Description: R. G. N. and J. L. S. acknowledge the support of NASA grants NAG5-6451 and NAG5-6591, respectively, as part of the JGOFS Synthesis and Modeling Program. G. K. P. and F. J. acknowledge support by the Swiss National Science Foundation. European contributions were supported by the EU GOSAC Project (ENV4-CT97- 0495).

Description: Author Posting. © American Geophysical Union, 2009. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 23 (2009): GB3016, doi:10.1029/2008GB003440.

Description: We present results from transient sensitivity studies with the Biogeochemical Elemental Cycling (BEC) ocean model to increasing anthropogenic atmospheric inorganic nitrogen (N) and soluble iron (Fe) deposition over the industrial era. Elevated N deposition results from fossil fuel combustion and agriculture, and elevated soluble Fe deposition results from increased atmospheric processing in the presence of anthropogenic pollutants and soluble Fe from combustion sources. Simulations with increasing Fe and increasing Fe and N inputs raised simulated marine nitrogen fixation, with the majority of the increase in the subtropical North and South Pacific, and raised primary production and export in the high-nutrient low-chlorophyll (HNLC) regions. Increasing N inputs alone elevated small phytoplankton and diatom production, resulting in increased phosphorus (P) and Fe limitation for diazotrophs, hence reducing nitrogen fixation (∼6%). Globally, the simulated primary production, sinking particulate organic carbon (POC) export. and atmospheric CO2 uptake were highest under combined increase in Fe and N inputs compared to preindustrial control. Our results suggest that increasing combustion iron sources and aerosol Fe solubility along with atmospheric anthropogenic nitrogen deposition are perturbing marine biogeochemical cycling and could partially explain the observed trend toward increased P limitation at station ALOHA in the subtropical North Pacific. Excess inorganic nitrogen ([NO3 −] + [NH4 +] − 16[PO4 3−]) distributions may offer useful insights for understanding changing ocean circulation and biogeochemistry.

Description: This work was supported by funding from NSF grant OCE-0452972 to J. K. Moore and C. S. Zender. Computations were supported by the Earth System Modeling Facility at UCI (NSFATMO321380) and by the Climate Simulation Laboratory at National Center for Atmospheric Research. The National Center for Atmospheric Research is sponsored by the U.S. National Science Foundation. N.M. would like to acknowledge the assistance of NSF– Carbon and Water (ATM-0628472), and N.M., S.D., and C.L. would like to acknowledge the assistance of NASA-IDS (NNX07AL80G).