Present-day Simulations Research Articles

Evaluation of present-day rainfall simulations over West Africa in CORDEX regional climate models

The objective of this study is to evaluate the ability of seven CORDEX regional climate models (RCMs), driven by ERA-Interim reanalysis dataset to simulate the observed rainfall characteristics over West Africa during the period of 1990–2008. The seasonal climatology, annual rainfall cycles, interannual variability, 850 hPa specific humidity, and wind fields of the RCMs outputs were assessed over a number of spatial scales covering three climatically homogenous subregions (Guinea Coast, Savannah, and Sahel) and the entire West Africa domain. The ability of the RCMs to simulate the response to El Nino and La Nina events were further assessed. Results indicate that the RCMs captured the spatial pattern of rainfall and the three distinctive phases of the West African monsoon reasonably. It is worth noting that RCA and CRCM5 failed to distinctively reproduce the monsoon jump while CCLM, HIRHAM, and REMO largely overestimated the amount of the pre- and the post-monsoon rainfall. The analysis also showed significant biases in individual models depending on the subregion and season under consideration. These biases appear to be linked to the model’s failure to resolve convective processes and topography accurately. The majority of the RCMs used were consistent with the ground observation in capturing the dry (wet) conditions associated with the El Nino (La Nina) events. Statistical analysis conclusively revealed that the RCMs performance varies over the subregions and seasons, implying that no single model is best at all time. In general, REGCM3 was found to be the most outstanding of all the RCMs and is therefore recommended for use in rainfall assessment over West Africa.

The Present-Day Simulation and Twenty-First-Century Projection of the Climatology of Extratropical Transition in the North Atlantic

This study explores the simulations and twenty-first-century projections of extratropical transition (ET) of tropical cyclones (TCs) in the North Atlantic, with a newly developed global climate model: the Forecast-Oriented Low Ocean Resolution (FLOR) version of the Geophysical Fluid Dynamics Laboratory (GFDL) Coupled Model version 2.5 (CM2.5). FLOR exhibits good skill in simulating present-day ET properties (e.g., cyclone phase space parameters). A version of FLOR in which sea surface temperature (SST) biases are artificially corrected through flux-adjustment (FLOR-FA) shows much improved simulation of ET activity (e.g., annual ET number). This result is largely attributable to better simulation of basinwide TC activity, which is strongly dependent on larger-scale climate simulation. FLOR-FA is also used to explore changes of ET activity in the twenty-first century under the representative concentration pathway (RCP) 4.5 scenario. A contrasting pattern is found in which regional TC density increases in the eastern North Atlantic and decreases in the western North Atlantic, probably due to changes in the TC genesis location. The increasing TC frequency in the eastern Atlantic is dominated by increased ET cases. The increased density of TCs undergoing ET in the eastern subtropics of the Atlantic shows two propagation paths: one moves northwest toward the northeast coast of the United States and the other moves northeast toward western Europe, implying increased TC-related risks in these regions. A more TC-favorable future climate, evident in the projected changes of SST and vertical wind shear, is hypothesized to favor the increased ET occurrence in the eastern North Atlantic.

Transient Climate Sensitivity Depends on Base Climate Ocean Circulation

Abstract There is large uncertainty in the simulation of transient climate sensitivity. This study aims to understand how such uncertainty is related to the simulation of the base climate by comparing two simulations with the same model but in which CO2 is increased from either a preindustrial (1860) or a present-day (1990) control simulation. This allows different base climate ocean circulations that are representative of those in current climate models to be imposed upon a single model. As a result, the model projects different transient climate sensitivities that are comparable to the multimodel spread. The greater warming in the 1990-start run occurs primarily at high latitudes and particularly over regions of oceanic convection. In the 1990-start run, ocean overturning circulations are initially weaker and weaken less from CO2 forcing. As a consequence, there are smaller reductions in the poleward ocean heat transport, leading to less tropical ocean heat storage and less moderated high-latitude surface warming. This process is evident in both hemispheres, with changes in the Atlantic meridional overturning circulation and the Antarctic Bottom Water formation dominating the warming differences in each hemisphere. The high-latitude warming in the 1990-start run is enhanced through albedo and cloud feedbacks, resulting in a smaller ocean heat uptake efficacy. The results highlight the importance of improving the base climate ocean circulation in order to provide a reasonable starting point for assessments of past climate change and the projection of future climate change.

Bioavailable atmospheric phosphorous supply to the global ocean: a 3-D global modeling study

Abstract. The atmospheric cycle of phosphorus (P) is parameterized here in a state-of-the-art global 3-D chemistry transport model, taking into account primary emissions of total P (TP) and soluble P (DP) associated with mineral dust, combustion particles from natural and anthropogenic sources, bioaerosols, sea spray and volcanic aerosols. For the present day, global TP emissions are calculated to be roughly 1.33 Tg-P yr−1, with the mineral sources contributing more than 80 % to these emissions. The P solubilization from mineral dust under acidic atmospheric conditions is also parameterized in the model and is calculated to contribute about one-third (0.14 Tg-P yr−1) of the global DP atmospheric source. To our knowledge, a unique aspect of our global study is the explicit modeling of the evolution of phosphorus speciation in the atmosphere. The simulated present-day global annual DP deposition flux is 0.45 Tg-P yr−1 (about 40 % over oceans), showing a strong spatial and temporal variability. Present-day simulations of atmospheric P aerosol concentrations and deposition fluxes are satisfactory compared with available observations, indicating however an underestimate of about 70 % on current knowledge of the sources that drive the P atmospheric cycle. Sensitivity simulations using preindustrial (year 1850) anthropogenic and biomass burning emission scenarios showed a present-day increase of 75 % in the P solubilization flux from mineral dust, i.e., the rate at which P is converted into soluble forms, compared to preindustrial times, due to increasing atmospheric acidity over the last 150 years. Future reductions in air pollutants due to the implementation of air-quality regulations are expected to decrease the P solubilization flux from mineral dust by about 30 % in the year 2100 compared to the present day. Considering, however, that all the P contained in bioaerosols is readily available for uptake by marine organisms, and also accounting for all other DP sources, a total bioavailable P flux of about 0.17 Tg-P yr−1 to the oceans is derived. Our calculations further show that in some regions more than half of the bioavailable P deposition flux to the ocean can originate from biological particles, while this contribution is found to maximize in summer when atmospheric deposition impact on the marine ecosystem is the highest due to ocean stratification. Thus, according to this global study, a largely unknown but potentially important role of terrestrial bioaerosols as suppliers of bioavailable P to the global ocean is also revealed. Overall, this work provides new insights to the atmospheric P cycle by demonstrating that biological materials are important carriers of bioavailable P, with very important implications for past and future responses of marine ecosystems to global change.

Present-day constraint for tropical Pacific precipitation changes due to global warming in CMIP5 models

The sensitivity of the precipitation responses to greenhouse warming can depend on the present-day climate. In this study, a robust linkage between the present-day precipitation climatology and precipitation change owing to global warming is examined in inter-model space. A model with drier climatology in the present-day simulation tends to simulate an increase in climatological precipitation owing to global warming. Moreover, the horizontal gradient of the present-day precipitation climatology plays an important role in determining the precipitation changes. On the basis of these robust relationships, future precipitation changes are calibrated by removing the impact of the present-day precipitation bias in the climate models. To validate this result, the perfect model approach is adapted, which treats a particular model’s precipitation change as an observed change. The results suggest that the precipitation change pattern can be generally improved by applying the present statistical approach.

Arctic climate changes in the 21st century: Ensemble model estimates accounting for realism in present-day climate simulation

In the course of forecasting future climate changes in the Arctic Region based on calculations and an ensemble of the state-of-the-art global climate models, the results depend on the method of construction the statistics from the models.

Separating climate change signals into thermodynamic, lapse-rate and circulation effects: theory and application to the European summer climate

Climate models robustly project a strong overall summer warming across Europe showing a characteristic north-south gradient with enhanced warming and drying in southern Europe. However, the processes that are responsible for this pattern are not fully understood. We here employ an extended surrogate or pseudo-warming approach to disentangle the contribution of different mechanisms to this response pattern. The basic idea of the surrogate technique is to use a regional climate model and apply a large-scale warming to the lateral boundary conditions of a present-day reference simulation, while maintaining the relative humidity (and thus implicitly increasing the specific moisture content). In comparison to previous studies, our approach includes two important extensions: first, different vertical warming profiles are applied in order to separate the effects of a mean warming from lapse-rate effects. Second, a twin-design is used, in which the climate change signals are not only added to present-day conditions, but also subtracted from a scenario experiment. We demonstrate that these extensions provide an elegant way to separate the full climate change signal into contributions from large-scale thermodynamic (TD), lapse-rate (LR), and circulation and other remaining effects (CO). The latter in particular include changes in land-ocean contrast and spatial variations of the SST warming patterns. We find that the TD effect yields a large-scale warming across Europe with no distinct latitudinal gradient. The LR effect, which is quantified for the first time in our study, leads to a stronger warming and some drying in southern Europe. It explains about 50 % of the warming amplification over the Iberian Peninsula, thus demonstrating the important role of lapse-rate changes. The effect is linked to an extending Hadley circulation. The CO effect as inherited from the driving GCM is shown to further amplify the north-south temperature change gradient. In terms of mean summer precipitation the TD effect leads to a significant overall increase in precipitation all across Europe, which is compensated and regionally reversed by the LR and CO effects in particular in southern Europe.

Human-Driven Atmospheric Deposition of N and P Controls on the East Mediterranean Marine Ecosystem

Abstract The historical and future impacts of atmospheric deposition of inorganic nitrogen (N) and phosphorus (P) on the marine ecosystem in the east Mediterranean Sea are investigated by using a 1D coupled physical– biogeochemical model, set up for the Cretan Sea as a representative area of the basin. For the present-day simulation (2010), the model is forced by observations of atmospheric deposition fluxes at Crete, while for the hindcast (1860) and forecast (2030) simulations, the changes in atmospheric deposition calculated by global chemistry–transport models are applied to the present-day observed fluxes. The impact of the atmospheric deposition on the fluxes of carbon in the food chain is calculated together with the contribution of human activities to these impacts. The results show that total phytoplanktonic biomass increased by 16% over the past 1.5 centuries. Small fractional changes in carbon fluxes and planktonic biomasses are predicted for the near future. Simulations show that atmospheric deposition of N and P may be the main mechanism responsible for the anomalous N:P ratio observed in the Mediterranean Sea.

Future U.S. ozone projections dependence on regional emissions, climate change, long-range transport and differences in modeling design

Abstract A consistent modeling framework with nested global and regional chemical transport models (CTMs) is used to separate and quantitatively assess the relative contributions to projections of future U.S. ozone pollution from the effects of emissions changes, climate change, long-range transport (LRT) of pollutants, and differences in modeling design. After incorporating dynamic lateral boundary conditions (LBCs) from a global CTM, a regional CTM's representation of present-day U.S. ozone pollution is notably improved, especially relative to results from the regional CTM with fixed LBCs or from the global CTM alone. This nested system of global and regional CTMs projects substantial surface ozone trends for the 2050's: 6–10 ppb decreases under the ‘clean’ A1B scenario and ∼15 ppb increases under the ‘dirty’ A1Fi scenario. Among the total trends of future ozone, regional emissions changes dominate, contributing negative 25–60% in A1B and positive 30–45% in A1Fi. Comparatively, climate change contributes positive 10–30%, while LRT effects through changing chemical LBCs account for positive 15–20% in both scenarios, suggesting introducing dynamic LBCs could influence projections of the U.S. future ozone pollution with a magnitude comparable to effects of climate change alone. The contribution to future ozone projections due to differences in modeling design, including model formulations, emissions treatments, and other factors between the global and the nested regional CTMs, is regionally dependent, ranging from negative 20% to positive 25%. It is shown that the model discrepancies for present-day simulations between global and regional CTMs can propagate into future U.S. ozone projections systematically but nonlinearly, especially in California and the Southeast. Therefore in addition to representations of emissions change and climate change, accurate treatment of LBCs for the regional CTM is essential for projecting the future U.S. ozone pollution.

Multiphase simulation of mine waters and aqueous leaching processes

Managing of large amounts of water in mining and mineral processing sites remains a concern in both actively operated and closed mining areas. When the mining site with its metal or concentrate producing units is operational, the challenge is to find either ways for economical processing with maximum yields, while minimizing the environmental impact of the water usage and waste salt treatments. For safe closure of the site, the environmental control of possible drainage will be needed. For both challenges, the present-day multiphase process simulations tools can be used to provide improved accuracy and better economy in controlling the smooth and environmentally sound operation of the plant. One of the pioneering studies in using the multiphase thermodynamic software in simulation of hydrometallurgical processes was that of Koukkari et al. [1]. The study covered the use of Solgasmix equilibrium software for a number of practical acid digesters. The models were made for sulfuric acid treatments in titania pigment production and in NPK fertilizer manufacturing. During the past two decades the extensive data assessment has taken place particularly in geochemistry and a new versions of geochemical multiphase equilibrium software has been developed. On the other hand, there has been some progress in development of the process simulation software in all the aforementioned fields. Thus, the thermodynamic simulation has become a tool of great importance in development of hydrometallurgical processes. The presentation will cover three example cases of either true pilot or industrial systems including a South African acid mine water drainage treatment, hydrometallurgical extraction of rare earths from uranium leachate in Russia and a multistage process simulation of a Finnish heap leaching mine with its subsequent water treatment system.

Summer Convective Afternoon Rainfall Simulation and Projection Using WRF Driven by Global Climate Model. Part I: Over Taiwan

This paper is the first of a two-part study that investigates summer convective afternoon rainfall (CAR) activity over the East Asian continent through simulation and projection. We focus on the CAR activity over Taiwan in Part I. In Part II, the changes in CAR activity over South China and Luzon are examined. Using the WRF (Weather Research and Forecasting) model driven by two super-high resolution global climate models: HiRAM (High Resolution Atmospheric Model) and MRI (Meteorological Research Institute Atmospheric General Circulation Model), this study evaluates the performance of models (HiRAM, MRI, and WRF) used to simulate the summer CAR activity over Taiwan. The evaluations focus on the spatial-temporal variations in CAR activity during two time periods: the present-day (1979 - 2003, historical run) and the future (2075 - 2099, RCP8.5 scenario). For the present-day simulations, analyses indicate that both HiRAM and MRI cannot accurately simulate the temporal evolution of diurnal rainfall. This timing shift problem can be fixed after dynamical downscaling using WRF. The WRF dynamical downscaling approach also helps generate a more realistic CAR amount simulation over western Taiwan, where most residents live. In future projections all models predict that a significant decrease in CAR amount will occur over southwest Taiwan. This decrease in CAR amount is suggested due to the decrease in CAR frequency, not a change in the CAR rate. All models suggest that the cause of decrease in CAR frequency over southwest Taiwan is weaker local afternoon surface wind convergence and thermal instability to suppress CAR genesis.

A Python tool to set up relative free energy calculations in GROMACS.

Free energy calculations based on molecular dynamics (MD) simulations have seen a tremendous growth in the last decade. However, it is still difficult and tedious to set them up in an automated manner, as the majority of the present-day MD simulation packages lack that functionality. Relative free energy calculations are a particular challenge for several reasons, including the problem of finding a common substructure and mapping the transformation to be applied. Here we present a tool, alchemical-setup.py, that automatically generates all the input files needed to perform relative solvation and binding free energy calculations with the MD package GROMACS. When combined with Lead Optimization Mapper (LOMAP;Liu et al. in J Comput Aided Mol Des 27(9):755-770, 2013), recently developed in our group, alchemical-setup.py allows fully automated setup of relative free energy calculations in GROMACS. Taking a graph of the planned calculations and a mapping, both computed by LOMAP, our tool generates the topology and coordinate files needed to perform relative free energy calculations for a given set of molecules, and provides a set of simulation input parameters. The tool was validated by performing relative hydration free energy calculations for a handful of molecules from the SAMPL4 challenge (Mobley et al. in J Comput Aided Mol Des 28(4):135-150, 2014). Good agreement with previously published results and the straightforward way in which free energy calculations can be conducted make alchemical-setup.py a promising tool for automated setup of relative solvation and binding free energy calculations.

Evaluating the climate and air quality impacts of short-lived pollutants

Evaluating the climate and air quality impacts of short-lived pollutants

Early-Holocene warming in Beringia and its mediation by sea-level and vegetation changes

Abstract. Arctic land-cover changes induced by recent global climate change (e.g., expansion of woody vegetation into tundra and effects of permafrost degradation) are expected to generate further feedbacks to the climate system. Past changes can be used to assess our understanding of feedback mechanisms through a combination of process modeling and paleo-observations. The subcontinental region of Beringia (northeastern Siberia, Alaska, and northwestern Canada) was largely ice-free at the peak of deglacial warming and experienced both major vegetation change and loss of permafrost when many arctic regions were still ice covered. The evolution of Beringian climate at this time was largely driven by global features, such as the amplified seasonal cycle of Northern Hemisphere insolation and changes in global ice volume and atmospheric composition, but changes in regional land-surface controls, such as the widespread development of thaw lakes, the replacement of tundra by deciduous forest or woodland, and the flooding of the Bering–Chukchi land bridge, were probably also important. We examined the sensitivity of Beringia's early Holocene climate to these regional-scale controls using a regional climate model (RegCM). Lateral and oceanic boundary conditions were provided by global climate simulations conducted using the GENESIS V2.01 atmospheric general circulation model (AGCM) with a mixed-layer ocean. We carried out two present-day simulations of regional climate – one with modern and one with 11 ka geography – plus another simulation for 6 ka. In addition, we performed five ~ 11 ka climate simulations, each driven by the same global AGCM boundary conditions: (i) 11 ka Control, which represents conditions just prior to the major transitions (exposed land bridge, no thaw lakes or wetlands, widespread tundra vegetation), (ii) sea-level rise, which employed present-day continental outlines, (iii) vegetation change, with deciduous needleleaf and deciduous broadleaf boreal vegetation types distributed as suggested by the paleoecological record, (iv) thaw lakes, which used the present-day distribution of lakes and wetlands, and (v) post-11 ka All, incorporating all boundary conditions changed in experiments (ii)–(iv). We find that regional-scale controls strongly mediate the climate responses to changes in the large-scale controls, amplifying them in some cases, damping them in others, and, overall, generating considerable spatial heterogeneity in the simulated climate changes. The change from tundra to deciduous woodland produces additional widespread warming in spring and early summer over that induced by the 11 ka insolation regime alone, and lakes and wetlands produce modest and localized cooling in summer and warming in winter. The greatest effect is the flooding of the land bridge and shelves, which produces generally cooler conditions in summer but warmer conditions in winter and is most clearly manifest on the flooded shelves and in eastern Beringia. By 6 ka continued amplification of the seasonal cycle of insolation and loss of the Laurentide ice sheet produce temperatures similar to or higher than those at 11 ka, plus a longer growing season.

Influence of In-Situ Oil Sands Development on Caribou (Rangifer tarandus) Movement.

In-situ oil sands development (ISD) involves a network of facilities, wells, roads and pipelines to extract and transport subsurface bitumen. This technology is rapidly expanding and there is uncertainty whether ISDs restrict animal movement, leading to increased extinction probabilities for some wide-ranging species. Here we test for effects of simulated future (i.e., 50 years from now) and current ISDs on simulated movements of woodland caribou (Rangifer tarandus), a threatened species across North America. In simulations of future scenarios, we varied the spacing and permeability of ISDs and the presence/absence of protected areas. Permeability was measured as the number of times simulated caribou crossed ISDs with different levels of modelled permeability. We estimated the effects of these factors on caribou step length and annual home range size, key metrics of small and large spatiotemporal scales of movement, respectively. Current caribou crossings of above-ground pipeline features of ISDs were measured using camera traps and compared to expected caribou crossing rates based on present-day caribou movement simulations. Current crossing rates were evaluated within the context of predicted future crossing success rates necessary to maintain caribou step lengths and home ranges. With few exceptions, permeability across ISDs was the main factor affecting caribou movement, more so than spacing between developments or the presence of protected areas. However, minimal permeability (crossing rates of c. 15% to 60%, relative to an undisturbed site was needed to maintain existing home range size and step lengths. The effect of permeability on home range size and step length was non-linear, suggesting that small increases in permeability would provide a disproportionately greater benefit to caribou movement. Our predictions demonstrate that maintaining permeability across ISDs is more important than spacing between leases or including protected areas, and thus provides clear direction for mitigation efforts for features that will exist on the landscape for decades to come.

Estimating heat stress from climate-based indicators: present-day biases and future spreads in the CMIP5 global climate model ensemble

The increased exposure of human populations to heat stress is one of the likely consequences of global warming, and it has detrimental effects on health and labor capacity. Here, we consider the evolution of heat stress under climate change using 21 general circulation models (GCMs). Three heat stress indicators, based on both temperature and humidity conditions, are used to investigate present-day model biases and spreads in future climate projections. Present day estimates of heat stress indicators from observational data shows that humid tropical areas tend to experience more frequent heat stress than other regions do, with a total frequency of heat stress 250–300 d yr−1. The most severe heat stress is found in the Sahel and south India. Present-day GCM simulations tend to underestimate heat stress over the tropics due to dry and cold model biases. The model based estimates are in better agreement with observation in mid to high latitudes, but this is due to compensating errors in humidity and temperature. The severity of heat stress is projected to increase by the end of the century under climate change scenario RCP8.5, reaching unprecedented levels in some regions compared with observations. An analysis of the different factors contributing to the total spread of projected heat stress shows that spread is primarily driven by the choice of GCMs rather than the choice of indicators, even when the simulated indicators are bias-corrected. This supports the utility of the multi-model ensemble approach to assess the impacts of climate change on heat stress.

Regional climate simulations over Vietnam using the WRF model

We present an analysis of the present-day (1961–1990) regional climate simulations over Vietnam. The regional climate model Weather Research and Forecasting (WRF) was driven by the global reanalysis ERA40. The performance of the regional climate model in simulating the observed climate is evaluated with a main focus on precipitation and temperature. The regional climate model was able to reproduce the observed spatial patterns of the climate, although with some biases. The model also performed better in reproducing the extreme precipitation and the interannual variability. Overall, the WRF model was able to simulate the main regional signatures of climate variables, seasonal cycles, and frequency distributions. This study is an evaluation of the present-day climate simulations of a regional climate model at a resolution of 25 km. Given that dynamical downscaling has become common for studying climate change and its impacts, the study highlights that much more improvements in modeling might be necessary to yield realistic simulations of climate at high resolutions before they can be used for impact studies at a local scale. The need for a dense network of observations is also realized as observations at high resolutions are needed when it comes to evaluations and validations of models at sub-regional and local scales.

Changes in dissolved iron deposition to the oceans driven by human activity: a 3-D global modelling study

Abstract. The global atmospheric iron (Fe) cycle is parameterized in the global 3-D chemical transport model TM4-ECPL to simulate the proton- and the organic ligand-promoted mineral-Fe dissolution as well as the aqueous-phase photochemical reactions between the oxidative states of Fe (III/II). Primary emissions of total (TFe) and dissolved (DFe) Fe associated with dust and combustion processes are also taken into account, with TFe mineral emissions calculated to amount to ~ 35 Tg-Fe yr−1 and TFe emissions from combustion sources of ~ 2 Tg-Fe yr−1. The model reasonably simulates the available Fe observations, supporting the reliability of the results of this study. Proton- and organic ligand-promoted Fe dissolution in present-day TM4-ECPL simulations is calculated to be ~ 0.175 Tg-Fe yr−1, approximately half of the calculated total primary DFe emissions from mineral and combustion sources in the model (~ 0.322 Tg-Fe yr−1). The atmospheric burden of DFe is calculated to be ~ 0.024 Tg-Fe. DFe deposition presents strong spatial and temporal variability with an annual flux of ~ 0.496 Tg-Fe yr−1, from which about 40 % (~ 0.191 Tg-Fe yr−1) is deposited over the ocean. The impact of air quality on Fe deposition is studied by performing sensitivity simulations using preindustrial (year 1850), present (year 2008) and future (year 2100) emission scenarios. These simulations indicate that about a 3 times increase in Fe dissolution may have occurred in the past 150 years due to increasing anthropogenic emissions and thus atmospheric acidity. Air-quality regulations of anthropogenic emissions are projected to decrease atmospheric acidity in the near future, reducing to about half the dust-Fe dissolution relative to the present day. The organic ligand contribution to Fe dissolution shows an inverse relationship to the atmospheric acidity, thus its importance has decreased since the preindustrial period but is projected to increase in the future. The calculated changes also show that the atmospheric DFe supply to the globe has more than doubled since the preindustrial period due to 8-fold increases in the primary non-dust emissions and about a 3-fold increase in the dust-Fe dissolution flux. However, in the future the DFe deposition flux is expected to decrease (by about 25 %) due to reductions in the primary non-dust emissions (about 15 %) and in the dust-Fe dissolution flux (about 55 %). The present level of atmospheric deposition of DFe over the global ocean is calculated to be about 3 times higher than for 1850 emissions, and about a 30 % decrease is projected for 2100 emissions. These changes are expected to impact most on the high-nutrient–low-chlorophyll oceanic regions.

Spatial distribution of precipitation annual cycles over South Africa in 10 CORDEX regional climate model present-day simulations

This study presents an evaluation of the ability of 10 regional climate models (RCMs) participating in the COordinated Regional climate Downscaling Experiment-Africa to reproduce the present-day spatial distribution of annual cycles of precipitation over the South African region and its borders. As found in previous studies, annual mean precipitation is quasi-systematically overestimated by the RCMs over a large part of southern Africa south of about 20°S and more strongly over South Africa. The spatial analysis of precipitation over the studied region shows that in most models the distribution of biases appears to be linked to orography. Wet biases are quasi-systematic in regions with higher elevation with inversely neutral to dry biases particularly in the coastal fringes. This spatial pattern of biases is particularly obvious during summer and specifically at the beginning of the rainy season (November and December) when the wet biases are found to be the strongest across all models. Applying a k-means algorithm, a classification of annual cycles is performed using observed precipitation data, and is compared with those derived from modeled data. It is found that the in-homogeneity of the spatial and temporal distribution of biases tends to impact the modeled seasonality of precipitation. Generally, the pattern of rainfall seasonality in the ensemble mean of the 10 RCMs tends to be shifted to the southwest. This spatial shift is mainly linked to a strong overestimation of convective precipitation at the beginning of the rainy season over the plateau inducing an early annual peak and to an underestimation of stratiform rainfall in winter and spring over southwestern South Africa.

Causes of Strengthening and Weakening of ENSO Amplitude under Global Warming in Four CMIP5 Models

AbstractThe mechanisms for El Niño–Southern Oscillation (ENSO) amplitude change under global warming are investigated through quantitative assessment of air–sea feedback processes in present-day and future climate simulations of four models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5). Two models (MPI-ESM-MR and MRI-CGCM3) project strengthened ENSO amplitude, whereas the other two models (CCSM4 and FGOALS-g2) project weakened ENSO amplitude. A mixed layer heat budget diagnosis shows that the major cause of the projected ENSO amplitude difference between the two groups is attributed to the changes of the thermocline and zonal advective feedbacks. A weaker (stronger) equatorial thermocline response to a unit anomalous zonal wind stress forcing in the Niño-4 region is found in CCSM4 and FGOALS-g2 (MPI-ESM-MR and MRI-CGCM3). The cause of the different response arises from the change in the meridional scale of ENSO. A narrower (wider) meridional width of sea surface temperature (SST) and zonal wind stress anomalies causes a strengthening (weakening) of the equatorial thermocline response and thus stronger Bjerknes and zonal advective feedbacks, as the subsurface temperature and zonal current anomalies depend on the thermocline response; consequently, the ENSO amplitude increases (decreases). The change of ENSO meridional width is caused by the change in mean meridional overturning circulation in the equatorial Pacific Ocean, which depends on change of mean wind stress and SST warming patterns under global warming.