Search Results

Now showing 1 - 9 of 9
  • Item
    Emulating atmosphere-ocean and carbon cycle models with a simpler model, MAGICC6 - Part 2: Applications
    (München : European Geopyhsical Union, 2011) Meinshausen, M.; Raper, S.C.B.; Wigley, T.M.L.
    Intercomparisons of coupled atmosphere-ocean general circulation models (AOGCMs) and carbon cycle models are important for galvanizing our current scientific knowledge to project future climate. Interpreting such intercomparisons faces major challenges, not least because different models have been forced with different sets of forcing agents. Here, we show how an emulation approach with MAGICC6 can address such problems. In a companion paper (Meinshausen et al., 2011a), we show how the lower complexity carbon cycle-climate model MAGICC6 can be calibrated to emulate, with considerable accuracy, globally aggregated characteristics of these more complex models. Building on that, we examine here the Coupled Model Intercomparison Project's Phase 3 results (CMIP3). If forcing agents missed by individual AOGCMs in CMIP3 are considered, this reduces ensemble average temperature change from pre-industrial times to 2100 under SRES A1B by 0.4 °C. Differences in the results from the 1980 to 1999 base period (as reported in IPCC AR4) to 2100 are negligible, however, although there are some differences in the trajectories over the 21st century. In a second part of this study, we consider the new RCP scenarios that are to be investigated under the forthcoming CMIP5 intercomparison for the IPCC Fifth Assessment Report. For the highest scenario, RCP8.5, relative to pre-industrial levels, we project a median warming of around 4.6 °C by 2100 and more than 7 °C by 2300. For the lowest RCP scenario, RCP3-PD, the corresponding warming is around 1.5 °C by 2100, decreasing to around 1.1 °C by 2300 based on our AOGCM and carbon cycle model emulations. Implied cumulative CO2 emissions over the 21st century for RCP8.5 and RCP3-PD are 1881 GtC (1697 to 2034 GtC, 80% uncertainty range) and 381 GtC (334 to 488 GtC), when prescribing CO2 concentrations and accounting for uncertainty in the carbon cycle. Lastly, we assess the reasons why a previous MAGICC version (4.2) used in IPCC AR4 gave roughly 10% larger warmings over the 21st century compared to the CMIP3 average. We find that forcing differences and the use of slightly too high climate sensitivities inferred from idealized high-forcing runs were the major reasons for this difference.
  • Item
    Emulating coupled atmosphere-ocean and carbon cycle models with a simpler model, MAGICC6 - Part 1: Model description and calibration
    (München : European Geopyhsical Union, 2011) Meinshausen, M.; Raper, S.C.B.; Wigley, T.M.L.
    Current scientific knowledge on the future response of the climate system to human-induced perturbations is comprehensively captured by various model intercomparison efforts. In the preparation of the Fourth Assessment Report (AR4) of the Intergovernmental Panel on Climate Change (IPCC), intercomparisons were organized for atmosphere-ocean general circulation models (AOGCMs) and carbon cycle models, named "CMIP3" and "C4MIP", respectively. Despite their tremendous value for the scientific community and policy makers alike, there are some difficulties in interpreting the results. For example, radiative forcings were not standardized across the various AOGCM integrations and carbon cycle runs, and, in some models, key forcings were omitted. Furthermore, the AOGCM analysis of plausible emissions pathways was restricted to only three SRES scenarios. This study attempts to address these issues. We present an updated version of MAGICC, the simple carbon cycle-climate model used in past IPCC Assessment Reports with enhanced representation of time-varying climate sensitivities, carbon cycle feedbacks, aerosol forcings and ocean heat uptake characteristics. This new version, MAGICC6, is successfully calibrated against the higher complexity AOGCMs and carbon cycle models. Parameterizations of MAGICC6 are provided. The mean of the emulations presented here using MAGICC6 deviates from the mean AOGCM responses by only 2.2% on average for the SRES scenarios. This enhanced emulation skill in comparison to previous calibrations is primarily due to: making a "like-with-like comparison" using AOGCM-specific subsets of forcings; employing a new calibration procedure; as well as the fact that the updated simple climate model can now successfully emulate some of the climate-state dependent effective climate sensitivities of AOGCMs. The diagnosed effective climate sensitivity at the time of CO2 doubling for the AOGCMs is on average 2.88 °C, about 0.33 °C cooler than the mean of the reported slab ocean climate sensitivities. In the companion paper (Part 2) of this study, we examine the combined climate system and carbon cycle emulations for the complete range of IPCC SRES emissions scenarios and the new RCP pathways.
  • Item
    Sensitivity of polar stratospheric ozone loss to uncertainties in chemical reaction kinetics
    (Göttingen : Copernicus GmbH, 2009) Kawa, S.R.; Stolarski, R.S.; Newman, P.A.; Douglass, A.R.; Rex, M.; Hofmann, D.J.; Santee, M.L.; Frieler, K.
    The impact and significance of uncertainties in model calculations of stratospheric ozone loss resulting from known uncertainty in chemical kinetics parameters is evaluated in trajectory chemistry simulations for the Antarctic and Arctic polar vortices. The uncertainty in modeled ozone loss is derived from Monte Carlo scenario simulations varying the kinetic (reaction and photolysis rate) parameters within their estimated uncertainty bounds. Simulations of a typical winter/spring Antarctic vortex scenario and Match scenarios in the Arctic produce large uncertainty in ozone loss rates and integrated seasonal loss. The simulations clearly indicate that the dominant source of model uncertainty in polar ozone loss is uncertainty in the Cl2O 2 photolysis reaction, which arises from uncertainty in laboratory-measured molecular cross sections at atmospherically important wavelengths. This estimated uncertainty in JCl 2O2 from laboratory measurements seriously hinders our ability to model polar ozone loss within useful quantitative error limits. Atmospheric observations, however, suggest that the Cl2O2 photolysis uncertainty may be less than that derived from the lab data. Comparisons to Match, South Pole ozonesonde, and Aura Microwave Limb Sounder (MLS) data all show that the nominal recommended rate simulations agree with data within uncertainties when the Cl2O2 photolysis error is reduced by a factor of two, in line with previous in situ ClOx measurements. Comparisons to simulations using recent cross sections from Pope et al. (2007) are outside the constrained error bounds in each case. Other reactions producing significant sensitivity in polar ozone loss include BrO + ClO and its branching ratios. These uncertainties challenge our confidence in modeling polar ozone depletion and projecting future changes in response to changing halogen emissions and climate. Further laboratory, theoretical, and possibly atmospheric studies are needed.
  • Item
    Potential yield simulated by global gridded crop models: using a process-based emulator to explain their differences
    (Katlenburg-Lindau : Copernicus, 2021-3-23) Ringeval, Bruno; Müller, Christoph; Pugh, Thomas A. M.; Mueller, Nathaniel D.; Ciais, Philippe; Folberth, Christian; Liu, Wenfeng; Debaeke, Philippe; Pellerin, Sylvain
    How global gridded crop models (GGCMs) differ in their simulation of potential yield and reasons for those differences have never been assessed. The GGCM Intercomparison (GGCMI) offers a good framework for this assessment. Here, we built an emulator (called SMM for simple mechanistic model) of GGCMs based on generic and simplified formalism. The SMM equations describe crop phenology by a sum of growing degree days, canopy radiation absorption by the Beer–Lambert law, and its conversion into aboveground biomass by a radiation use efficiency (RUE). We fitted the parameters of this emulator against gridded aboveground maize biomass at the end of the growing season simulated by eight different GGCMs in a given year (2000). Our assumption is that the simple set of equations of SMM, after calibration, could reproduce the response of most GGCMs so that differences between GGCMs can be attributed to the parameters related to processes captured by the emulator. Despite huge differences between GGCMs, we show that if we fit both a parameter describing the thermal requirement for leaf emergence by adjusting its value to each grid-point in space, as done by GGCM modellers following the GGCMI protocol, and a GGCM-dependent globally uniform RUE, then the simple set of equations of the SMM emulator is sufficient to reproduce the spatial distribution of the original aboveground biomass simulated by most GGCMs. The grain filling is simulated in SMM by considering a fixed-in-time fraction of net primary productivity allocated to the grains (frac) once a threshold in leaves number (nthresh) is reached. Once calibrated, these two parameters allow for the capture of the relationship between potential yield and final aboveground biomass of each GGCM. It is particularly important as the divergence among GGCMs is larger for yield than for aboveground biomass. Thus, we showed that the divergence between GGCMs can be summarized by the differences in a few parameters. Our simple but mechanistic model could also be an interesting tool to test new developments in order to improve the simulation of potential yield at the global scale.
  • Item
    NCIO 1.0: A simple Fortran NetCDF interface
    (München : European Geopyhsical Union, 2015) Robinson, A.; Perrette, M.
    The NetCDF (Network Common Data Form) library has become an indispensable tool for data and model output management in geoscience. However for simple tasks, particularly in Fortran, the complexity of native NetCDF functionality can be cumbersome. The NCIO (NetCDF Input/Output) module has been designed as an interface to the NetCDF library with simplicity and ease of use in mind. While this implies that some NetCDF functionality is masked from the user, the subroutines provided here are adequate for basic serial reading and writing tasks of up to 6-D data arrays along with corresponding data attributes. The code is available online via a GitHub repository (http://www.github.com/alex-robinson/ncio), which includes an example program to illustrate the approach.
  • Item
    An extended singular spectrum transformation (SST) for the investigation of Kenyan precipitation data
    (Göttingen : Copernicus GmbH, 2013) Itoh, N.; Marwan, N.
    In this paper a change-point detection method is proposed by extending the singular spectrum transformation (SST) developed as one of the capabilities of singular spectrum analysis (SSA). The method uncovers change points related with trends and periodicities. The potential of the proposed method is demonstrated by analysing simple model time series including linear functions and sine functions as well as real world data (precipitation data in Kenya). A statistical test of the results is proposed based on a Monte Carlo simulation with surrogate methods. As a result, the successful estimation of change points as inherent properties in the representative time series of both trend and harmonics is shown. With regards to the application, we find change points in the precipitation data of Kenyan towns (Nakuru, Naivasha, Narok, and Kisumu) which coincide with the variability of the Indian Ocean Dipole (IOD) suggesting its impact of extreme climate in East Africa.
  • Item
    Review: Visual analytics of climate networks
    (Göttingen : Copernicus GmbH, 2015) Nocke, T.; Buschmann, S.; Donges, J.F.; Marwan, N.; Schulz, H.-J.; Tominski, C.
  • Item
    Lévy noise improves the electrical activity in a neuron under electromagnetic radiation
    (San Francisco, CA : Public Library of Science (PLoS), 2017) Wu, J.; Xu, Y.; Ma, J.
    As the fluctuations of the internal bioelectricity of nervous system is various and complex, the external electromagnetic radiation induced by magnet flux on membrane can be described by the non-Gaussian type distribution of Levy noise. Thus, the electrical activities in an improved Hindmarsh-Rose model excited by the external electromagnetic radiation of Levy noise are investigated and some interesting modes of the electrical activities are exhibited. The external electromagnetic radiation of Levy noise leads to the mode transition of the electrical activities and spatial phase, such as from the rest state to the firing state, from the spiking state to the spiking state with more spikes, and from the spiking state to the bursting state. Then the time points of the firing state versus Levy noise intensity are depicted. The increasing of Levy noise intensity heightens the neuron firing. Also the stationary probability distribution functions of the membrane potential of the neuron induced by the external electromagnetic radiation of Levy noise with different intensity, stability index and skewness papremeters are analyzed. Moreover, through the positive largest Lyapunov exponent, the parameter regions of chaotic electrical mode of the neuron induced by the external electromagnetic radiation of Levy noise distribution are detected.
  • Item
    Stabilizing the West Antarctic Ice Sheet by surface mass deposition
    (Washington, DC [u.a.] : Assoc., 2019) Feldmann, Johannes; Levermann, Anders; Mengel, Matthias
    There is evidence that a self-sustaining ice discharge from the West Antarctic Ice Sheet (WAIS) has started, potentially leading to its disintegration. The associated sea level rise of more than 3m would pose a serious challenge to highly populated areas including metropolises such as Calcutta, Shanghai, New York City, and Tokyo. Here, we show that the WAIS may be stabilized through mass deposition in coastal regions around Pine Island and Thwaites glaciers. In our numerical simulations, a minimum of 7400 Gt of additional snowfall stabilizes the flow if applied over a short period of 10 years onto the region (−2 mm year−1 sea level equivalent). Mass deposition at a lower rate increases the intervention time and the required total amount of snow. We find that the precise conditions of such an operation are crucial, and potential benefits need to be weighed against environmental hazards, future risks, and enormous technical challenges. Copyright © 2019 The Authors,