Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Aprubt climate change events in the Northeastern Atlantic Ocean

The western Iberian margin has been one of the key locations to study abrupt glacial climate change and associated interhemispheric linkages. The regional variability in the response to those events is being studied by combining a multitude of published and new records. Looking at the trend from Marine Isotope Stage (MIS) 10 to 2, the planktic foraminifer data, conform with the alkenone record of Martrat et al. [2007], shows that abrupt climate change events, especially the Heinrich events, became more frequent and their impacts in general stronger during the last glacial cycle. However, there were two older periods with strong impacts on the Atlantic meridional overturning circulation (AMOC): the Heinrich-type event associated with Termination (T) IV and the one occurring during MIS 8 (269 to 265 ka). During the Heinrich stadials of the last glacial cycle, the polar front reached the northern Iberian margin (ca. 41°N), while the arctic front was located in the vicinity of 39°N. During all the glacial periods studied, there existed a boundary at the latter latitude, either the arctic front during extreme cold events or the subarctic front during less strong coolings or warmer glacials. Along with these fronts sea surface temperatures (SST) increased southward by about 1°C per one degree of latitude leading to steep temperature gradients in the eastern North Atlantic and pointing to a close vicinity between subpolar and subtropical waters. The southern Iberian margin was always bathed by subtropical water masses - surface and/ or subsurface ones -, but there were periods when these waters also penetrated northward to 40.6°N. Glacial hydrographic conditions were similar during MIS 2 and 4, but much different during MIS 6. MIS 6 was a warmer glacial with the polar front being located further to the north allowing the subtropical surface and subsurface waters to reach at minimum as far north as 40.6°N and resulting in relative stable conditions on the southern margin. In the vertical structure, the Greenland-type climate oscillations during the last glacial cycle were recorded down to 2465 m during the Heinrich stadials, i.e. slightly deeper than in the western basin. This deeper boundary is related to the admixing of Mediterranean Outflow Water, which also explains the better ventilation of the intermediate-depth water column on the Iberian margin. This compilation revealed that latitudinal, longitudinal and vertical gradients existed in the waters along the Iberian margin, i.e. in a relative restricted area, but sufficient paleo-data exists now to validate regional climate models for abrupt climate change events in the northeastern North Atlantic Ocean.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Measured daily precipitation sums over Central Sulawesi for the period 1981-1999

Data compiled within the IMPENSO project. The Impact of ENSO on Sustainable Water Management and the Decision-Making Community at a Rainforest Margin in Indonesia (IMPENSO), http://www.gwdg.de/~impenso, was a German-Indonesian research project (2003-2007) that has studied the impact of ENSO (El Nino-Southern Oscillation) on the water resources and the agricultural production in the PALU RIVER watershed in Central Sulawesi. ENSO is a climate variability that causes serious droughts in Indonesia and other countries of South-East Asia. The last ENSO event occurred in 1997. As in other regions, many farmers in Central Sulawesi suffered from reduced crop yields and lost their livestock. A better prediction of ENSO and the development of coping strategies would help local communities mitigate the impact of ENSO on rural livelihoods and food security. The IMPENSO project deals with the impact of the climate variability ENSO (El Niño Southern Oscillation) on water resource management and the local communities in the Palu River watershed of Central Sulawesi, Indonesia. The project consists of three interrelated sub-projects, which study the local and regional manifestation of ENSO using the Regional Climate Models REMO and GESIMA (Sub-project A), quantify the impact of ENSO on the availability of water for agriculture and other uses, using the distributed hydrological model WaSiM-ETH (Sub-project B), and analyze the socio-economic impact and the policy implications of ENSO on the basis of a production function analysis, a household vulnerability analysis, and a linear programming model (Sub-project C). The models used in the three sub-projects will be integrated to simulate joint scenarios that are defined in collaboration with local stakeholders and are relevant for the design of coping strategies.

Modeled monthly local topographic lapse rates, Pacific Northwest, USA

An emerging approach to downscaling the projections from General Circulation Models (GCMs) to scales relevant for basin hydrology is to use output of GCMs to force higher-resolution Regional Climate Models (RCMs). With spatial resolution often in the tens of kilometers, however, even RCM output will likely fail to resolve local topography that may be climatically significant in high-relief basins. Here we develop and apply an approach for downscaling RCM output using local topographic lapse rates (empirically-estimated spatially and seasonally variable changes in climate variables with elevation). We calculate monthly local topographic lapse rates from the 800-m Parameter-elevation Regressions on Independent Slopes Model (PRISM) dataset, which is based on regressions of observed climate against topographic variables. We then use these lapse rates to elevationally correct two sources of regional climate-model output: (1) the North American Regional Reanalysis (NARR), a retrospective dataset produced from a regional forecasting model constrained by observations, and (2) a range of baseline climate scenarios from the North American Regional Climate Change Assessment Program (NARCCAP), which is produced by a series of RCMs driven by GCMs. By running a calibrated and validated hydrologic model, the Soil and Water Assessment Tool (SWAT), using observed station data and elevationally-adjusted NARR and NARCCAP output, we are able to estimate the sensitivity of hydrologic modeling to the source of the input climate data. Topographic correction of regional climate-model data is a promising method for modeling the hydrology of mountainous basins for which no weather station datasets are available or for simulating hydrology under past or future climates.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.

Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

Transient Holocene experiments under orbital forcing using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) have been suggested to exert a critical influence on global climate through wind-driven upwelling of deep water in the Southern Ocean and the potentially resulting atmospheric CO2 variations. The investigation of the temporal and spatial evolution of the SWW along with forcings and feedbacks remains a significant challenge in climate research. In this study, the evolution of the SWW under orbital forcing from the early Holocene (9 kyr BP) to pre-industrial modern times is examined with transient experiments using the comprehensive coupled global climate model CCSM3. Analyses of the model results suggest that the annual and seasonal mean SWW were subject to an overall strengthening and poleward shifting trend during the course of the early-to-late Holocene under the influence of orbital forcing, except for the austral spring season, where the SWW exhibited an opposite trend of shifting towards the equator.

Geochemical and palynological measurements from the ODP Leg 189 at Site 1172 on the East Tasman Plateau

A brief (~150 kyr) period of widespread global average surface warming marks the transition between the Paleocene and Eocene epochs, ~56 million years ago. This so-called "Paleocene-Eocene thermal maximum" (PETM) is associated with the massive injection of 13C-depleted carbon, reflected in a negative carbon isotope excursion (CIE). Biotic responses include a global abundance peak (acme) of the subtropical dinoflagellate Apectodinium. Here we identify the PETM in a marine sedimentary sequence deposited on the East Tasman Plateau at Ocean Drilling Program (ODP) Site 1172 and show, based on the organic paleothermometer TEX86, that southwest Pacific sea surface temperatures increased from ~26 °C to ~33°C during the PETM. Such temperatures before, during and after the PETM are >10 °C warmer than predicted by paleoclimate model simulations for this latitude. In part, this discrepancy may be explained by potential seasonal biases in the TEX86 proxy in polar oceans. Additionally, the data suggest that not only Arctic, but also Antarctic temperatures may be underestimated in simulations of ancient greenhouse climates by current generation fully coupled climate models. An early influx of abundant Apectodinium confirms that environmental change preceded the CIE on a global scale. Organic dinoflagellate cyst assemblages suggest a local decrease in the amount of river run off reaching the core site during the PETM, possibly in concert with eustatic rise. Moreover, the assemblages suggest changes in seasonality of the regional hydrological system and storm activity. Finally, significant variation in dinoflagellate cyst assemblages during the PETM indicates that southwest Pacific climates varied significantly over time scales of 103 - 104 years during this event, a finding comparable to similar studies of PETM successions from the New Jersey Shelf.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Climate and total tree cover reconstructions based on the 43.7-kyr fossil pollen record from Khoe (NW Sakhalin Island) and modern climate variables of 236 modern surface pollen sampling sites

A late Quaternary pollen record from northern Sakhalin Island (51.34°N, 142.14°E, 15 m a.s.l.) spanning the last 43.7 ka was used to reconstruct regional climate dynamics and vegetation distribution by using the modern analogue technique (MAT). The long-term trends of the reconstructed mean annual temperature (TANN) and precipitation (PANN), and total tree cover are generally in line with key palaeoclimate records from the North Atlantic region and the Asian monsoon domain. TANN largely follows the fluctuations in solar summer insolation at 55°N. During Marine Isotope Stage (MIS) 3, TANN and PANN were on average 0.2 °C and 700 mm, respectively, thus very similar to late Holocene/modern conditions. Full glacial climate deterioration (TANN = -3.3 °C, PANN = 550 mm) was relatively weak as suggested by the MAT-inferred average climate parameters and tree cover densities. However, error ranges of the climate reconstructions during this interval are relatively large and the last glacial environments in northern Sakhalin could be much colder and drier than suggested by the weighted average values. An anti-phase relationship between mean temperature of the coldest (MTCO) and warmest (MTWA) month is documented during the last glacial period, i.e. MIS 2 and 3, suggesting more continental climate due to sea levels that were lower than present. Warmest and wettest climate conditions have prevailed since the end of the last glaciation with an optimum (TANN = 1.5 °C, PANN = 800 mm) in the middle Holocene interval (ca 8.7-5.2 cal. ka BP). This lags behind the solar insolation peak during the early Holocene. We propose that this is due to continuous Holocene sea level transgression and regional influence of the Tsushima Warm Current, which reached maximum intensity during the middle Holocene. Several short-term climate oscillations are suggested by our reconstruction results and correspond to Northern Hemisphere Heinrich and Dansgaard-Oeschger events, the Bølling-Allerød and the Younger Dryas. The most prominent fluctuation is registered during Heinrich 4 event, which is marked by noticeably colder and drier conditions and the spread of herbaceous taxa.

XRF-scanned elemental concentrations from sediment cores off Peru

We present a high-resolution marine record of sediment input from the Guayas River, Ecuador, that reflects changes in precipitation along western equatorial South America during the last 18ka. We use log (Ti/Ca) derived from X-ray Fluorescence (XRF) to document terrigenous input from riverine runoff that integrates rainfall from the Guayas River catchment. We find that rainfall-induced riverine runoff has increased during the Holocene and decreased during the last deglaciation. Superimposed on those long-term trends, we find that rainfall was probably slightly increased during the Younger Dryas, while the Heinrich event 1 was marked by an extreme load of terrigenous input, probably reflecting one of the wettest period over the time interval studied. When we compare our results to other Deglacial to Holocene rainfall records located across the tropical South American continent, different modes of variability become apparent. The records of rainfall variability imply that changes in the hydrological cycle at orbital and sub-orbital timescales were different from western to eastern South America. Orbital forcing caused an antiphase behavior in rainfall trends between eastern and western equatorial South America. In contrast, millennial-scale rainfall changes, remotely connected to the North Atlantic climate variability, led to homogenously wetter conditions over eastern and western equatorial South America during North Atlantic cold spells. These results may provide helpful diagnostics for testing the regional rainfall sensitivity in climate models and help to refine rainfall projections in South America for the next century.

Sedimentological and geophysical investigation of sediment cores from the Amundsen Sea

Ice loss from the marine-based, potentially unstable West Antarctic Ice Sheet (WAIS) contributes to current sea-level rise and may raise sea level by up to 3.3 to 5 meters in the future. Over the past few decades, glaciers draining the WAIS into the Amundsen Sea Embayment (ASE) have shown accelerated ice flow, rapid thinning and grounding-line retreat. However, the long-term context of this ice-sheet retreat is poorly constrained, limiting our ability to accurately predict future WAIS behaviour. Here we present a new chronology for WAIS retreat from the inner continental shelf of the eastern ASE based on radiocarbon dates from three marine sediment cores. The ages document a retreat of the grounding line to within ~93 km of its modern position before 11.7±0.7 kyr BP (thousand years before present). This early deglaciation is consistent with ages for grounding-line retreat from the western ASE. Our new data demonstrate that, other than in the Ross Sea, WAIS retreat in the ASE has not continued progressively since the Last Glacial Maximum. Furthermore, our results suggest that the grounding-line position in the ASE was predominantly stable throughout the Holocene, and that any episodes of fast retreat similar to that observed today must have been short-lived. Alternatively, today's rapid retreat was unprecedented during the Holocene. Therefore, the current ice loss must originate in recent changes in regional climate, ocean circulation or ice-sheet dynamics. Incorporation of these results into models is essential to produce robust predictions of future ice-sheet change and its contribution to sea-level rise.

Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.

Particle size of Saharan dust in the Atlantic Ocean, from October 2012 to November 2013

Mineral dust has a large impact on regional and global climate, depending on its particle size. Especially in the Atlantic Ocean downwind of the Sahara, the largest dust source on earth, the effects can be substantial but are poorly understood. This study focuses on seasonal and spatial variations in particle size of Saharan dust deposition across the Atlantic Ocean, using an array of submarine sediment traps moored along a transect at 12° N. We show that the particle size decreases downwind with increased distance from the Saharan source, due to higher gravitational settling velocities of coarse particles in the atmosphere. Modal grain sizes vary between 4 and 33 µm throughout the different seasons and at five locations along the transect. This is much coarser than previously suggested and incorporated into climate models. In addition, seasonal changes are prominent, with coarser dust in summer, and finer dust in winter and spring. Such seasonal changes are caused by transport at higher altitudes and at greater wind velocities during summer than in winter. Also the latitudinal migration of the dust cloud, associated with the Intertropical Convergence Zone, causes seasonal differences in deposition as the summer dust cloud is located more to the north, and more directly above the sampled transect. Furthermore, increased precipitation and more frequent dust storms in summer coincide with coarser dust deposition. Our findings contribute to understanding Saharan dust transport and deposition relevant for the interpretation of sedimentary records for climate reconstructions, as well as for global and regional models for improved prediction of future climate.