Refine
Year of publication
- 2011 (134) (remove)
Document Type
- Article (108)
- Doctoral Thesis (20)
- Review (2)
- Habilitation Thesis (1)
- Other (1)
- Postprint (1)
- Preprint (1)
Language
- English (134) (remove)
Keywords
- Holocene (6)
- Tibetan Plateau (6)
- Site effects (4)
- Holozän (3)
- Pollen (3)
- Wave propagation (3)
- climate change (3)
- Benzene (2)
- Climate change (2)
- Crustal structure (2)
Institute
- Institut für Geowissenschaften (134) (remove)
Humidity is an important determinant of the mycotoxin production (DON, ZEA) by Fusarium species in the grain ears. From a landscape perspective humidity is not evenly distributed across fields. The topographically-controlled redistribution of water within a single field rather leads to spatially heterogeneous soil water content and air humidity. Therefore we hypothesized that the spatial distribution of mycotoxins is related to these topographically-controlled factors. To test this hypothesis we studied the mycotoxin concentrations at contrasting topographic relief positions, i.e. hilltops and depressions characterized by soils of different soil moisture regimes, on ten winter wheat fields in 2006 and 2007. Maize was the preceding crop and minimum tillage was practiced in the fields. The different topographic positions were associated with moderate differences in DON and ZEA concentrations in 2006, but with significant differences in 2007, with six times higher median ZEA and two times higher median DON detected at depression sites compared to the hilltops. The depression sites correspond to a higher topographic wetness index as well as redoximorphic properties in soil profiles, which empirically supports our hypothesis at least for years showing wetter conditions in sensitive time windows for Fusarium infections.
The open source computational fluid dynamics (CFD) wind model (CFD-WEM) for wind erosion research in the Xilingele grassland in Inner Mongolia (autonomous region, China) is compared with two open source CFD models Gerris and OpenFOAM. The evaluation of these models was made according to software technology, implemented methods, handling, accuracy and calculation speed. All models were applied to the same wind tunnel data set. Results show that the simplest CFD-WEM has the highest calculation speed with acceptable accuracy, and the most powerful OpenFOAM produces the simulation with highest accuracy and the lowest calculation speed. Gerris is between CFD-WEM and OpenFOAM. It calculates faster than OpenFOAM, and it is capable to solve different CFD problems. CFD-WEM is the optimal model to be further developed for wind erosion research in Inner Mongolia grassland considering its efficiency and the uncertainties of other input data. However, for other applications using CFD technology, Gerris and OpenFOAM can be good choices. This paper shows the powerful capability of open source CFD software in wind erosion study, and advocates more involvement of open source technology in wind erosion and related ecological researches.
Wavelet modelling of the gravity field by domain decomposition methods: an example over Japan
(2011)
With the advent of satellite gravity, large gravity data sets of unprecedented quality at low and medium resolution become available. For local, high resolution field modelling, they need to be combined with the surface gravity data. Such models are then used for various applications, from the study of the Earth interior to the determination of oceanic currents. Here we show how to realize such a combination in a flexible way using spherical wavelets and applying a domain decomposition approach. This iterative method, based on the Schwarz algorithms, allows to split a large problem into smaller ones, and avoids the calculation of the entire normal system, which may be huge if high resolution is sought over wide areas. A subdomain is defined as the harmonic space spanned by a subset of the wavelet family. Based on the localization properties of the wavelets in space and frequency, we define hierarchical subdomains of wavelets at different scales. On each scale, blocks of subdomains are defined by using a tailored spatial splitting of the area. The data weighting and regularization are iteratively adjusted for the subdomains, which allows to handle heterogeneity in the data quality or the gravity variations. Different levels of approximations of the subdomains normals are also introduced, corresponding to building local averages of the data at different resolution levels.
We first provide the theoretical background on domain decomposition methods. Then, we validate the method with synthetic data, considering two kinds of noise: white noise and coloured noise. We then apply the method to data over Japan, where we combine a satellite-based geopotential model, EIGEN-GL04S, and a local gravity model from a combination of land and marine gravity data and an altimetry-derived marine gravity model. A hybrid spherical harmonics/wavelet model of the geoid is obtained at about 15 km resolution and a corrector grid for the surface model is derived.
Vertical flow filters and vertical flow constructed wetlands are established wastewater treatment systems and have also been proposed for the treatment of contaminated groundwater. This study investigates the removal processes of volatile organic compounds in a pilot-scale vertical flow filter. The filter is intermittently irrigated with contaminated groundwater containing benzene, MTBE and ammonium as the main contaminants. The system is characterized by unsaturated conditions and high contaminant removal efficiency. The aim of the present study is to evaluate the contribution of biodegradation and volatilization to the overall removal of benzene and MTBE. Tracer tests and flow rate measurements showed a highly transient flow and heterogeneous transport regime. Radon-222, naturally occurring in the treated groundwater, was used as a gas tracer and indicated a high volatilization potential. Radon-222 behavior was reproduced by numerical simulations and extrapolated for benzene and MTBE, and indicated these compounds also have a high volatilization potential. In contrast, passive sampler measurements on top of the filter detected only low benzene and MTBE concentrations. Biodegradation potential was evaluated by the analysis of catabolic genes involved in organic compound degradation and a quantitative estimation of biodegradation was derived from stable isotope fractionation analysis. Results suggest that despite the high volatilization potential, biodegradation is the predominant mass removal process in the filter system, which indicates that the volatilized fraction of the contaminants is still subject to subsequent biodegradation. In particular, the upper filter layer located between the injection tubes and the surface of the system might also contribute to biodegradation, and might play a crucial role in avoiding the emission of volatilized contaminants into the atmosphere.
The East African Plateau provides a spectacular example of geodynamic plateau uplift, active continental rifting, and associated climatic forcing. It is an integral part of the East African Rift System and has an average elevation of approximately 1,000 m. Its location coincides with a negative Bouguer gravity anomaly with a semi-circular shape, closely related to a mantle plume, which influences the Cenozoic crustal development since its impingement in Eocene-Oligocene time. The uplift of the East African Plateau, preceding volcanism, and rifting formed an important orographic barrier and tectonically controlled environment, which is profoundly influenced by climate driven processes. Its location within the equatorial realm supports recently proposed hypotheses, that topographic changes in this region must be considered as the dominant forcing factor influencing atmospheric circulation patterns and rainfall distribution. The uplift of this region has therefore often been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While the far-reaching influence of the plateau uplift is widely accepted, the timing and the magnitude of the uplift are ambiguous and are still subject to ongoing discussion. This dilemma stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. In order to quantify the amount of plateau uplift and to find evidence for the existence of significant relief along the East African Plateau prior to rifting, I analyzed and modeled one of the longest terrestrial lava flows; the 300-km-long Yatta phonolite flow in Kenya. This lava flow is 13.5 Ma old and originated in the region that now corresponds to the eastern rift shoulders. The phonolitic flow utilized an old riverbed that once drained the eastern flank of the plateau. Due to differential erosion this lava flow now forms a positive relief above the parallel-flowing Athi River, which is mimicking the course of the paleo-river. My approach is a lava-flow modeling, based on an improved composition and temperature dependent method to parameterize the flow of an arbitrary lava in a rectangular-shaped channel. The essential growth pattern is described by a one-dimensional model, in which Newtonian rheological flow advance is governed by the development of viscosity and/or velocity in the internal parts of the lava-flow front. Comparing assessments of different magma compositions reveal that length-dominated, channelized lava flows are characterized by high effusion rates, rapid emplacement under approximately isothermal conditions, and laminar flow. By integrating the Yatta lava flow dimensions and the covered paleo-topography (slope angle) into the model, I was able to determine the pre-rift topography of the East African Plateau. The modeling results yield a pre-rift slope of at least 0.2°, suggesting that the lava flow must have originated at a minimum elevation of 1,400 m. Hence, high topography in the region of the present-day Kenya Rift must have existed by at least 13.5 Ma. This inferred mid-Miocene uplift coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa. Accordingly, the combination of my results regarding the Yatta lava flow emplacement history, its location, and its morphologic character, validates it as a suitable “paleo-tiltmeter” and has thus to be considered as an important topographic and volcanic feature for the topographic evolution in East Africa.
Overland flow is an important hydrological pathway in many forests of the humid tropics. Its generation is subject to topographic controls at differing spatial scales. Our objective was to identify such controls on the occurrence of overland flow in a lowland tropical rainforest. To this end, we installed 95 overland flow detectors (OFDs) in four nested subcatchments of the Lutzito catchment on Barro Colorado Island, Panama, and monitored the frequency of overland flow occurrence during 18 rainfall events at each OFD location temporal frequency. For each such location, we derived three non-digital terrain attributes and 17 digital ones, of which 15 were based on Digital Elevation Models (DEMs) of three different resolutions. These attributes then served as input into a Random Forest ensemble tree model to elucidate the importance and partial and joint dependencies of topographic controls for overland flow occurrence.
Lutzito features a high median temporal frequency in overland flow occurrence of 0.421 among OFD locations. However, spatial temporal frequencies of overland flow occurrence vary strongly among these locations and the subcatchments of Lutzito catchment. This variability is best explained by (1) microtopography, (2) coarse terrain sloping and (3) various measures of distance-to-channel, with the contribution of all other terrain attributes being small. Microtopographic features such as concentrated flowlines and wash areas produce highest temporal frequencies, whereas the occurrence of overland flow drops sharply for flow distances and terrain sloping beyond certain threshold values.
Our study contributes to understanding both the spatial controls on overland flow generation and the limitations of terrain attributes for the spatially explicit prediction of overland flow frequencies.
Despite the importance of rhizosphere properties for water flow from soil to roots, there is limited quantitative information on the distribution of water in the rhizosphere of plants.
Here, we used neutron tomography to quantify and visualize the water content in the rhizosphere of the plant species chickpea (Cicer arietinum), white lupin (Lupinus albus), and maize (Zea mays) 12 d after planting.
We clearly observed increasing soil water contents (h) towards the root surface for all three plant species, as opposed to the usual assumption of decreasing water content. This was true for tap roots and lateral roots of both upper and lower parts of the root system. Furthermore, water gradients around the lower part of the roots were smaller and extended further into bulk soil compared with the upper part, where the gradients in water content were steeper.
Incorporating the hydraulic conductivity and water retention parameters of the rhizosphere into our model, we could simulate the gradual changes of h towards the root surface, in agreement with the observations. The modelling result suggests that roots in their rhizosphere may modify the hydraulic properties of soil in a way that improves uptake under dry conditions.
Three-dimensional hydrostratigraphic models from ground-penetrating radar and direct-push data
(2011)
Three-dimensional models of hydraulic conductivity and porosity are essential to understand and simulate groundwater flow in heterogeneous geological environments. However, considering the inherent limitations of traditional hydrogeological field methods in terms of resolution, alternative field approaches are needed to establish such 3-D models with sufficient accuracy. In this study, we developed a workflow combining 3-D structural information extracted from ground penetrating radar (GPR) images with 1-D in situ physical-property estimates from direct-push (DP) logging to construct a 3-D hydrostratigraphic model. To illustrate this workflow, we collected an similar to 70 m x 90 m 100 MHz 3-D GPR data set over a shallow sedimentary aquifer system resolving six different GPR facies down to similar to 15 m depth. DP logs of the relative dielectric permittivity, the relative hydraulic conductivity, the cone resistance, the sleeve friction and the pore pressure provided crucial data (1) to establish a GPR velocity model for 3-D depth migration and to check the time-to-depth conversion of the GPR data, and (2) to construct a 3-D hydrostratigraphic model. This model was built by assigning porosity values, which were computed from the DP relative dielectric permittivity logs, and DP relative hydraulic conductivity estimates to the identified GPR facies. We conclude that the integration of 3-D GPR structural images and 1-D DP logs of target physical parameters provides an efficient way for detailed 3-D subsurface characterization as needed, for example, for groundwater flow simulations.
The Valais units in Savoy (Zone des BrSches de Tarentaise) have been re-mapped in great detail and are subject of combined stratigraphic, structural and petrological investigations summarized in this contribution. The sediments and rare relics of basement, together with Cretaceous age mafic and ultramafic rocks of the Valais palaeogeographical domain, represent the heavily deformed relics of the former distal European margin (External Valais units) and an ocean-continent transition (Internal Valais unit or Versoyen unit) that formed during rifting. This rifting led to the opening of the Valais ocean, a northern branch of the Alpine Tethys. Post-rift sediments referred to as "Valais trilogy" stratigraphically overlie both External and Internal Valais successions above an angular unconformity formed in Barremian to Aptian times, providing robust evidence for the timing of the opening of the Valais ocean. The Valais units in Savoy are part of a second and more external mid-Eocene high-pressure belt in the Alps that sutured the Brian double dagger onnais microcontinent to Europe. Top-N D1-deformation led to the formation of a nappe stack that emplaced the largely eclogite-facies Internal Valais unit (Versoyen) onto blueschist-facies External Valais units. The latter originally consisted of, from internal to external, the Petit St. Bernard unit, the Roc de l'Enfer unit, the MoA >> tiers unit and the Quermoz unit. Ongoing top-N D2-thrusting and folding substantially modified this nappe stack. Post 35 Ma D3 folding led to relatively minor modifications of the nappe stack within the Valais units but was associated with substantial top-WNW thrusting of the Valais units over the Dauphinois units along the Roselend thrust during W-directed indentation of the Adria block contributing to the formation of the arc of the Western Alps.
Eclogites from the main borehole of the Chinese Continental Scientific Drilling project yield highly precise Lu-Hf garnet-clinopyroxene ages of 216.9 +/- 1.2 Ma (four samples) and 220.5 +/- 2.7 Ma (one sample). The spatial distribution of the rare earth elements in garnet is consistent with the preservation of primary growth zoning, unmodified by diffusion, which supports the interpretation that the Lu-Hf ages date the time of formation of garnet, the major rock forming mineral in the eclogites. The preservation of primary REE-zoning, despite peak metamorphic temperatures around 800-850 degrees C. indicates that the Lu-Hf chronometer is perfectly suitable to date garnet-forming reactions in high grade rocks. The range of Lu-Hf ages for eclogites in the Dabie-Sulu UHP terrane point to episodic rather than continuous growth of garnets and thus punctuated metamorphism during the collision of the North China Block and the Yangtze Block. The U-Pb ages and Hf-isotope systematics of zircon grains from one eclogite sample imply a protracted geologic history of the eclogite precursors that started around 2 Ga and culminated in the UHP metamorphism around 220 Ma.
Synorogenic extension has been recognized as an integral structural constituent of mountain belts and high-elevation plateaus during their evolution. In the Himalaya, both orogen-parallel and orogen-normal extension has been recognized. However, the underlying driving forces for extension and their timing are still a matter of debate. Here we present new fault kinematic data based on systematic measurements of hundreds of outcrop-scale brittle fault planes in the NW Indian Himalaya. This new data set, as well as field observations including crosscutting relationships, mineral fibers on fault planes, and correlations with deformation structures in lake sediments, allows us to distinguish different deformation styles. The overall strain pattern derived from our data reflects the large regional contractional deformation pattern very well but also reveals significant extensional deformation in a region, which is dominated by shortening. In total, we were able to identify six deformation styles, most of which are temporally and spatially linked, representing protracted shortening. Our observations also furnish the basis for a detailed overview of the younger deformation history in the NW Himalaya, which has been characterized by extension overprinting previously generated structures related to shortening. The four dominant deformation styles are (1) shortening parallel to the regional convergence direction; (2) arc-normal extension; (3) arc-parallel extension; and finally, (4) E-W extension. This is the first data set where a succession of both arc-normal and E-W extension has been documented in the Himalaya. Importantly, our observations help differentiate E-W extension triggered by processes within the Tibetan Plateau from arc-parallel and arc-normal extension originating from the curvature of the Himalayan orogen.
Shell architectures of the larger foraminiferal genera cyclopseudedomia, and Rhapydionina were studied by comparing topotypes of previously described species with new specimens retrieved from Late Cretaceous shallow-water carbonates of Pylos (Peloponnese, Greece), where the three genera are found in association. The megalospheric generation of each genus exhibits a distinctive shell shape in adult specimens (i.e., fan-shaped in Cyclopseudedomia, conical in Rhapydionina, and cylindrical in Cuvillierinella). Although their microspheric adults are similarly thin, flat, and discoidal, they can be identified at the genus level by means of a detailed structural analysis. Cavillierinella shows the septula to be interrupted by a large preseptal space, while Cyclopseudedomia and Rhapydionina exhibit continuous, non-interrupted septula. In addition, Cyclopseudedomia presents only one row of medullar chamberlets, whereas Rhapydionina shows numerous medullar chamberlets distributed in a thick basal layer.
Two new species, Cuvillierinella pylosensis and Rhapydionina fleuryi, are described. The former is a more complex taxon than the type species, C salentina, while the latter corresponds to a more primitive species, R. liburnica. Strontium-isotope stratigraphy indicates an uppermost Campanian-lowermost Maastrichtian age for these new species.
The Guarguaraz Complex in West Argentina formed during collision between the microplate Chilenia and South America. It is composed of neritic clastic metasediments with intercalations of metabasic and ultrabasic rocks of oceanic origin. Prograde garnet growth in metapelite and metabasite occurred between 1.2 GPa, 470 degrees C and 1.4 GPa, 530 degrees C, when the penetrative s(2)-foliation was formed. The average age of garnet crystallization of 390 +/- 2 Ma (2 sigma) was determined from three four-point Lu-Hf mineral isochrones from metapelite and metabasite samples and represents the time of collision. Peak pressure conditions are followed by a decompression path with slight heating at 0.5 GPa, 560 degrees C. Fluid release during decompression caused equilibration of mineral compositions at the rims and also aided Ar diffusion. An Ar-40/39 Ar plateau age of white mica at 353 +/- 1 Ma (1 sigma) indicates the time of cooling below 350-400 degrees C. These temperatures were attained at pressures of 0.2-0.3 GPa, indicative of an average exhumation rate of >= 1 mm/a for the period 390-353 Ma. Late hydrous influx at 0.1-0.3 GPa caused pervasive growth of sericite and chlorite and reset the Ar/Ar ages of earlier coarse-grained white mica. At 284-295 Ma, the entire basement cooled below 280 degrees C (fission track ages of zircon) after abundant post-collisional granitoid intrusion. The deeply buried epicontinental sedimentary rocks, the high peak pressure referring to a low metamorphic geotherm of 10-12 degrees C/km, and the decompression/heating path are characteristics of material buried and exhumed within a (micro) continent-continent collisional setting.
Rapid population growth and economic development have led to increased anthropogenic pressures on the Tibetan Plateau, causing significant land cover changes with potentially severe ecological consequences. To assess whether or not these pressures are also affecting the remote montane-boreal lakes on the SE Tibetan Plateau, fossil pollen and diatom data from two lakes were synthesized. The interplay of aquatic and terrestrial ecosystem response was explored in respect to climate variability and human activity over the past 200 years. Nonmetric multidimensional scaling and Procrustes rotation analysis were undertaken to determine whether pollen and diatom responses in each lake were similar and synchronous. Detrended canonical correspondence analysis was used to develop quantitative estimates of compositional species turnover. Despite instrumental evidence of significant climatic warming on the southeastern Plateau, the pollen and diatom records indicate very stable species composition throughout their profiles and show only very subtle responses to environmental changes over the past 200 years. The compositional species turnover (0.36-0.94 SD) is relatively low in comparison to the species reorganizations known from the periods during the mid-and early-Holocene (0.64-1.61 SD) on the SE Plateau, and also in comparison to turnover rates of sediment records from climate-sensitive regions in the circum arctic. Our results indicate that climatically induced ecological thresholds are not yet crossed, but that human activity has an increasing influence, particularly on the terrestrial ecosystem in our study area. Synergistic processes of post-Little Ice Age warming, 20th century climate warming and extensive reforestations since the 19th century have initiated a change from natural oak-pine forests to seminatural, likely less resilient pine-oak forests. Further warming and anthropogenic disturbances would possibly exceed the ecological threshold of these ecosystems and lead to severe ecological consequences.
Well-preserved primary contact relationships between a Late Proterozoic metasedimentary and the metagranitic core and Palaeozoic cover series of the Menderes Massif have been recognized in the eastern part of the Cine submassif on a regional-scale. Metaconglomerates occur as laterally discontinuous channel-fill bodies close the base of the metaquartzarenite directly above the basement. The pebbles in the metaconglomerates consist mainly of different types of tourmaline-rich leucocratic granitoids, tourmalinite and schist in a sandy matrix. Petrographic features, geochemical compositions and zircon radiometric ages (549.6 +/- 3.7-552.3 +/- 3.1 Ma) of the diagnostic clasts of the metaconglomerates (e.g. leucocratic granitoids and tourmalinites) show excellent agreement with their in situ equivalents (549.0 +/- 5.4 Ma) occurring in the Pan-African basement as stocks and veins.
The correlation between clasts in the metaconglomerates and granitoids of the basement suggests that the primary contact between the basement and cover series is a regional unconformity (supra-Pan-African Unconformity) representing deep erosion of the Pan-African basement followed by the deposition of the cover series. Hence the usage of 'core-cover' terminology in the Menderes Massif is valid. Consequently, these new data preclude the views that the granitic precursors of the leucocratic orthogneisses are Tertiary intrusions.
Ether functionalities form an important cross-linking structure within the macromolecular organic matrix of lignites and coals. To obtain a deeper insight into the complex internal structure of such macromolecules and the maturation related changes of the ether compounds within the network structure, boron tribromide (BBr3) ether cleavage was applied to a series of lignite and coal samples of different maturity (R-0 0.27-0.80%) obtained from coal mines and natural outcrops from the North and South Island of New Zealand. Terminal ether-bound alcohols rapidly decrease during diagenesis and occur only in low amounts during the catagenetic stage. Comparison between ester- and ether-bound terminal alcohols indicates a parallel decreasing trend during the diagenetic stage, suggesting that the stability differences between both linkages are not large enough to be observed in maturation processes over geological time scales:Polyether compounds were detected with chain length up to five carbon atoms. After a small decrease during the diagenetic phase these compounds occur in relatively high concentrations, even in the main catagenetic stage. This suggests that these linkage structures represent important cross-linking substructures within the macromolecular matrix of lignites and coals being sterically protected within the macromolecular network during the maturation process. Additional cross-linking substructures were (poly)ether aromatics, esters and ketones.
Controversy about the current state and future evolution of Himalayan glaciers has been stirred up by erroneous statements in the fourth report by the Intergovernmental Panel on Climate Change(1,2). Variable retreat rates(3-6) and a paucity of glacial mass-balance data(7,8) make it difficult to develop a coherent picture of regional climate-change impacts in the region. Here, we report remotely-sensed frontal changes and surface velocities from glaciers in the greater Himalaya between 2000 and 2008 that provide evidence for strong spatial variations in glacier behaviour which are linked to topography and climate. More than 65% of the monsoon-influenced glaciers that we observed are retreating, but heavily debris-covered glaciers with stagnant low-gradient terminus regions typically have stable fronts. Debris-covered glaciers are common in the rugged central Himalaya, but they are almost absent in subdued landscapes on the Tibetan Plateau, where retreat rates are higher. In contrast, more than 50% of observed glaciers in the westerlies-influenced Karakoram region in the northwestern Himalaya are advancing or stable. Our study shows that there is no uniform response of Himalayan glaciers to climate change and highlights the importance of debris cover for understanding glacier retreat, an effect that has so far been neglected in predictions of future water availability(9,10) or global sea level(11).
Dust has the potential to modify global climate by influencing the radiative balance of the atmosphere and by supplying iron and other essential limiting micronutrients to the ocean(1,2). Indeed, dust supply to the Southern Ocean increases during ice ages, and 'iron fertilization' of the subantarctic zone may have contributed up to 40 parts per million by volume (p. p. m. v.) of the decrease (80-100 p. p. m. v.) in atmospheric carbon dioxide observed during late Pleistocene glacial cycles(3-7). So far, however, the magnitude of Southern Ocean dust deposition in earlier times and its role in the development and evolution of Pleistocene glacial cycles have remained unclear. Here we report a high-resolution record of dust and iron supply to the Southern Ocean over the past four million years, derived from the analysis of marine sediments from ODP Site 1090, located in the Atlantic sector of the subantarctic zone. The close correspondence of our dust and iron deposition records with Antarctic ice core reconstructions of dust flux covering the past 800,000 years (refs 8, 9) indicates that both of these archives record large-scale deposition changes that should apply to most of the Southern Ocean, validating previous interpretations of the ice core data. The extension of the record beyond the interval covered by the Antarctic ice cores reveals that, in contrast to the relatively gradual intensification of glacial cycles over the past three million years, Southern Ocean dust and iron flux rose sharply at the Mid-Pleistocene climatic transition around 1.25 million years ago. This finding complements previous observations over late Pleistocene glacial cycles(5,8,9), providing new evidence of a tight connection between high dust input to the Southern Ocean and the emergence of the deep glaciations that characterize the past one million years of Earth history.