The Hamburg Ocean Carbon Cycle Model HAMOCC5.1 - Technical Description Release 1.1.

The report is intended to introduce the reader to the model structure of the oceanic biogeochemical model HAMOCC5, and to assist in setting-up and running HAMOCC5 driven by the ocean general circulation model MPI-OM. First, an overview over the simulated biogeochemical processes is given. Second, the modules that compute biogeochemistry are described in detail. Third, the interface to MPI-OM and the coupling between the ocean physics and the biogeochemical tracer model are described. Finally, the input and output files are described and an appendix provides some tablesand cross references for the usage of the model code. HAMOCC5 simulates biogeochemical tracers in the oceanic water column and in the sediment. The model can be used with any Ocean General Circulation Model (OGCM). Here, it is set up as a subroutine of the Ocean Model of the Max-Planck-Institute for Meteorology, MPI-OM (Marsland et al., 2003). HAMOCC5 is driven by the same radiation as the OGCM to compute photosynthesis. Temperature and salinity from MPI-OM are used to calculate chemical transformation rates and constants e.g., for solubility of carbon dioxide. The flux of carbon dioxide between atmosphere and ocean is computed depending on the local concentrations and the rates for air-sea gas exchange. With only few modifications any field of atmospheric trace gases or wind stress can be used to drive the fluxes. The biogeochemistry of HAMOCC5 is based on that of HAMOCC3.1 (Six and Maier-Reimer, 1996). Modifications have been made to account for tracers in addition to phosphorous, namely nitrogen, nitrous oxide, DMS, dissolved iron and dust. Additional simulated processes are denitrification and N-fixation, formation of calcium carbonate and opaline shells by phytoplankton, aggregation and size dependent sinking of detritus and marine snow, DMS production, uptake and release of dissolved iron by biogenic particles, and dust deposition and sinking. Biological production is now co-limited by phosphorous, nitrogen, and iron. For IPCC-type experiments additional tracers, anthropogenic DIC and alkalinity have been introduced. The model now also features a sediment module based on Heinze and Maier-Reimer (1999) and Heinze et al. (1999). The sediment model basically calculates the same tracers as the water column model

Observation of a stronger-than-adiabatic change of light trapped in an ultrafast switched GaAs-AlAs microcavity

We study the time-resolved reflectivity spectrum of a switched planar GaAs-AlAs microcavity. Between 5 and 40 ps after the switching (pump) pulse we observe a strong excess probe reflectivity and a change of the frequency of light trapped in the cavity up to 5 linewidths away from the cavity resonance. This frequency change does not adiabatically follow the fast-changing cavity resonance. The frequency change is attributed to an accumulated phase change due to the time-dependent refractive index. An analytical model predicts dynamics in qualitative agreement with the experiments, and points to crucial parameters that control future applications.Comment: Discussed effect of probe bandwidth. Included functional forms of n(z) and R(z

Climate-induced interannual variability of marine primary and export production in three global coupled climate carbon cycle models

© 2008 Author(s). This article is distributed under the terms of the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 5 (2008): 597-614, doi:10.5194/bg-5-597-2008Fully coupled climate carbon cycle models are sophisticated tools that are used to predict future climate change and its impact on the land and ocean carbon cycles. These models should be able to adequately represent natural variability, requiring model validation by observations. The present study focuses on the ocean carbon cycle component, in particular the spatial and temporal variability in net primary productivity (PP) and export production (EP) of particulate organic carbon (POC). Results from three coupled climate carbon cycle models (IPSL, MPIM, NCAR) are compared with observation-based estimates derived from satellite measurements of ocean colour and results from inverse modelling (data assimilation). Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) with stronger stratification (higher sea surface temperature; SST) being associated with negative PP anomalies. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for anomalies in globally integrated PP, but only one model (IPSL) also reproduces the inverse relationship between stratification (SST) and PP. An adequate representation of iron and macronutrient co-limitation of phytoplankton growth in the tropical ocean has shown to be the crucial mechanism determining the capability of the models to reproduce observed interactions between climate and PP.This work was supported by the EU grants 511106-2 (FP6 RTD project EUR-OCEANS) and GOCE-511176 (FP6 RTP project CARBOOCEAN) by the European Commission. TLF and FJ also acknowledge support from the Swiss National Science Foundations. SCD and MJB received support from NASA NNG06G127G

Projected 21st century decrease in marine productivity : a multi-model analysis

© Authors, 2010. This work is distributed under the Creative Commons Attribution 3.0 License. The definitive version was published in Biogeosciences 7 (2010): 979-1005, doi: 10.5194/bg-7-979-2010Changes in marine net primary productivity (PP) and export of particulate organic carbon (EP) are projected over the 21st century with four global coupled carbon cycle-climate models. These include representations of marine ecosystems and the carbon cycle of different structure and complexity. All four models show a decrease in global mean PP and EP between 2 and 20% by 2100 relative to preindustrial conditions, for the SRES A2 emission scenario. Two different regimes for productivity changes are consistently identified in all models. The first chain of mechanisms is dominant in the low- and mid-latitude ocean and in the North Atlantic: reduced input of macro-nutrients into the euphotic zone related to enhanced stratification, reduced mixed layer depth, and slowed circulation causes a decrease in macro-nutrient concentrations and in PP and EP. The second regime is projected for parts of the Southern Ocean: an alleviation of light and/or temperature limitation leads to an increase in PP and EP as productivity is fueled by a sustained nutrient input. A region of disagreement among the models is the Arctic, where three models project an increase in PP while one model projects a decrease. Projected changes in seasonal and interannual variability are modest in most regions. Regional model skill metrics are proposed to generate multi-model mean fields that show an improved skill in representing observation-based estimates compared to a simple multi-model average. Model results are compared to recent productivity projections with three different algorithms, usually applied to infer net primary production from satellite observations.This work was funded by the European Union projects CARBOOCEAN (511176-2) and EUROCEANS (511106-2) and is a contribution to the “European Project on Ocean Acidification” (EPOCA) which received funding from the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreement no. 211384. Additional support was received from the Swiss National Science Foundation. SCD acknowledges support from the NASA Ocean Biology and Biogeochemistry Program (NNX07AL80G). LB aknowledges support from the EU Project MEECE (Marine Ecosystem Evolution in a Changing Environnement, grant agreement 212085)

Twenty-first-century compatible co2 emissions and airborne fraction simulated by cmip5 earth system models under four representative concentration pathways

PublishedJournal ArticleThe carbon cycle is a crucial Earth system component affecting climate and atmospheric composition. The response of natural carbon uptake to CO2 and climate change will determine anthropogenic emissions compatible with a target CO2 pathway. For phase 5 of the Coupled Model Intercomparison Project (CMIP5), four future representative concentration pathways (RCPs) have been generated by integrated assessment models (IAMs) and used as scenarios by state-of-the-art climate models, enabling quantification of compatible carbon emissions for the four scenarios by complex, process-based models. Here, the authors present results from 15 such Earth system GCMs for future changes in land and ocean carbon storage and the implications for anthropogenic emissions. The results are consistent with the underlying scenarios but show substantial model spread. Uncertainty in land carbon uptake due to differences among models is comparable with the spread across scenarios. Model estimates of historical fossil-fuel emissions agree well with reconstructions, and future projections for representative concentration pathway 2.6 (RCP2.6) and RCP4.5 are consistent with the IAMs. For high-end scenarios (RCP6.0 and RCP8.5), GCMs simulate smaller compatible emissions than the IAMs, indicating a larger climate-carbon cycle feedback in the GCMs in these scenarios. For the RCP2.6 mitigation scenario, an average reduction of 50% in emissions by 2050 from 1990 levels is required but with very large model spread (14%-96%). The models also disagree on both the requirement for sustained negative emissions to achieve the RCP2.6 CO2 concentration and the success of this scenario to restrict global warming below 28C. All models agree that the future airborne fraction depends strongly on the emissions profile with higher airborne fraction for higher emissions scenarios. ©2013 American Meteorological Society.MOHC authors were supported by the JointDECC/Defra MetOffice Hadley Centre Climate Programme (GA01101), and work to performHadGEM2- ES and MPI-ESM CMIP5 simulations was supported by the EU-FP7 COMBINE project (Grant 226520). JS was supported by the EU-FP7 CARBOCHANGE project (Grant 284679). We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. JT and CR were supported by the Research Council of Norway through the EarthClim (207711/E10) project

Inconsistent strategies to spin up models in CMIP5: implications for ocean biogeochemical model performance assessment

International audienceDuring the fifth phase of the Coupled Model Inter-comparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hind-casts were compared with available modern biogeochemi-cal observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to-model uncertainty

Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon

Global carbon budget 2013

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere 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 a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr−1, ELUC 0.9 ± 0.5 GtC yr−1, GATM 4.3 ± 0.1 GtC yr−1, SOCEAN 2.5 ± 0.5 GtC yr−1, and SLAND 2.8 ± 0.8 GtC yr−1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr−1, 2.2% above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and assuming an ELUC of 1.0 ± 0.5 GtC yr−1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr−1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70% from EFF (390 ± 20 GtC) and 30% from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3)

Trajectories of the Earth System in the Anthropocene

This is the final version of the article. Available from National Academy of Sciences via the DOI in this record.We explore the risk that self-reinforcing feedbacks could push the Earth System toward a planetary threshold that, if crossed, could prevent stabilization of the climate at intermediate temperature rises and cause continued warming on a "Hothouse Earth" pathway even as human emissions are reduced. Crossing the threshold would lead to a much higher global average temperature than any interglacial in the past 1.2 million years and to sea levels significantly higher than at any time in the Holocene. We examine the evidence that such a threshold might exist and where it might be. If the threshold is crossed, the resulting trajectory would likely cause serious disruptions to ecosystems, society, and economies. Collective human action is required to steer the Earth System away from a potential threshold and stabilize it in a habitable interglacial-like state. Such action entails stewardship of the entire Earth System-biosphere, climate, and societies-and could include decarbonization of the global economy, enhancement of biosphere carbon sinks, behavioral changes, technological innovations, new governance arrangements, and transformed social values.W.S. and C.P.S. are members of the Anthropocene Working Group. W.S., J.R., K.R., S.E.C., J.F.D., I.F., S.J.L., R.W. and H.J.S. are members of the Planetary Boundaries Research Network PB.net and the Earth League’s EarthDoc Programme supported by the Stordalen Foundation. T.M.L. was supported by a Royal Society Wolfson Research Merit Award and the European Union Framework Programme 7 Project HELIX. C.F. was supported by the Erling– Persson Family Foundation. The participation of D.L. was supported by the Haury Program in Environment and Social Justice and National Science Foundation (USA) Decadal and Regional Climate Prediction using Earth System Models Grant 1243125. S.E.C. was supported in part by Swedish Research Council Formas Grant 2012-742. J.F.D. and R.W. were supported by Leibniz Association Project DOMINOES. S.J.L. receives funding from Formas Grant 2014-589. This paper is a contribution to European Research Council Advanced Grant 2016, Earth Resilience in the Anthropocene Project 743080