A 450-ka long record of glaciation in Northern Mongolia based on studies at Lake Khubsugul: high-resolution reflection seismic data and grain-size variations in cored sediments

We use high-resolution reflection seismic data and detailed grain-size analysis of a drill core (KDP-01) from Lake Khubsugul (northern Mongolia) to provide an improved reconstruction of the glacial history of the area for the last 450 ka. Grain-size analysis of suspended sediment load in modern rivers draining into the lake and of moraine material from the northern part of the catchment shows that the silt fraction is transported to the central part of the lake mainly by river suspension, whereas the clay fraction is mainly transported by glacial meltwater during deglaciation. The changes in of the clay/silt ratio in Lake Khubsugul sediments correlates well with the standard global paleoclimate records: low clay/silt ratios indicate warm climates, while a high clay/silt ratio reflects glacial erosion and cold climates. Pulses of clay input into the lake occur at the final stages of glacial periods (i.e., glacial maxima and subsequent onsets of deglaciation). The periodicity in glacial clay input in Lake Khubsugul is in tune with global periods of deglaciation, but there are differences in the intensity of the deglacial events for MIS-12 and MIS-2. These differences are attributed to specific conditions in regional distribution of moisture during glaciation, glacial ice volumes, and solar insolation intensity at the onset of deglaciation. Deglaciation of the Khubsugul glaciers occurred in response to an increase in summer solar insolation above a threshold value of 490 W/m². Two types of deglaciation can be distinguished: (1) slow melting during several tens of 1,000 years during weak increases in summer insolation, and (2) short and fast melting during several thousands of years in response to strong increases in summer insolation. The maximum ice volume in the area of Lake Khubsugul during the past 450 ka occurred during the period of 373–350 ka BP (MIS 11a-10) and was caused by high levels of moisture in the region, whereas the MIS-2 and MIS-12 glacial periods were characterized by minima in ice volume, due to the strong aridity in the region.     

Phase relations on the diopside-jadeite-hedenbergite join up to 24 GPa and stability of Na-bearing majoritic garnet

Phase relations on the diopside (Di)-hedenbergite (Hd)-jadeite (Jd) system modeling mineral associations of natural eclogites were studied for the compositions (mol %) Di₇₀Jd₃₀, Di₅₀Jd₅₀, Di₃₀Jd₇₀, Di₂₀Hd₈₀, and Di₄₀Hd₁₀Jd₅₀ using a toroidal anvil-with-hole (7 GPa) and a Kawai-type 6-8 multianvil apparatus (12-24 GPa). We established that Di, Hd, and Jd form complete series of solid solutions at 7 GPa, and melting temperatures of pure Di (1980 °C) and Jd (1870 °C) for that pressure were estimated experimentally. The melting temperature for the Di₅₀Jd₅₀ composition at 15.5 GPa is 2270 °C. The appearance of garnet is clearly dependent on initial clinopyroxene composition: at 1600 °C the first garnet crystals are observed at 13.5 GPa in the jadeite-rich part of the system (Di₃₀Jd₇₀), whereas diopside-rich starting material (Di₇₀Jd₃₀) produces garnet only above 17 GPa. The proportion of garnet increases rapidly above 18 GPa as pyroxene dissolves in the garnet structure and pyroxene-free garnetites are produced from diopside-rich starting materials. In all experiments, garnet coexists with stishovite (St). At a pressure above 18 GPa, pyroxene is completely replaced by an assemblage of majorite (Maj) + St + CaSiO₃-perovskite (Ca-Pv) in Ca-rich systems, whereas Maj is associated with almost pure Jd up to a pressure of 21.5 GPa. Above ∼22 GPa, Maj, and St are associated with NaAlSiO₄ with calcium ferrite structure (Cf). We established that an Hd component also spreads the range of pyroxene stability up to 20 GPa. In the Di₇₀Jd₃₀ system at 24 GPa an assemblage of Maj + Ca-Pv + MgSiO₃ with ilmenite structure (Mg-Il) was obtained. The experimentally established correlation between Na, Si, and Al contents in Maj and pressure in Grt(Maj)-pyroxene assemblages, may be the basis for a “majorite” geobarometer. The results of our experiments are applicable to the upper mantle and the transition zone of the Earth (400-670 km), and demonstrate a wide range of transformations from eclogite to perovskite-bearing garnetite. In addition, the mineral associations obtained from the experiments allowed us to simulate parageneses of inclusions in diamonds formed under the conditions of the transition zone and the lower mantle.     

Toward Improving Coastal-Fog Prediction (C-FOG)

Boundary-Layer Meteorology -     

Prediction of the vapor-liquid distribution constants for volatile nonelectrolytes in water up to its critical temperature

The distribution of solutes between coexisting liquid and vapor phases of water can be expressed by the distribution constant, K_D, defined as

     

Grid point surface air temperature calculations with a fast turnaround: Combining the results of IMAGE and a GCM

This paper describes a methodology that combines the outputs of (1) the Integrated Model to Assess the Greenhouse Effect (IMAGE Version 1.0) of the Netherlands National Institute of Public Health and Environmental Protection (RIVM) (given a greenhouse gas emission policy, this model can estimate the effects such as global mean surface air temperature change for a wide variety of policies) and (2) ECHAM-1/LSG, the Global Circulation Model (GCM) of the Max-Planck Institute for Meteorology in Hamburg, Germany. The combination enables one to calculate grid point surface air temperature changes for different scenarios with a turnaround time that is much quicker than that for a GCM. The methodology is based upon a geographical pattern of the ratio of grid point temperature change to global mean values during a certain period of the simulation, as calculated by ECHAM-1/LSG for the 1990 Scenarios A and D of the Intergovernmental Panel on Climate Change (IPCC). A procedure, based upon signal-to noise ratios in the outputs, enabled us to estimate where we have confidence in the methodology; this is at about 23% to 83% of the total of 2,048 grid points, depending upon the scenario and the decade in the simulation. It was found that the methodology enabled IMAGE to provide useful estimates of the GCM-predicted grid point temperature changes. These estimates were within 0.5K (0.25K) throughout the 100 years of a given simulation for at least 79% (74%) of the grid points where we are confident in applying the methodology. The temperature ratio pattern from Scenario A enabled IMAGE to provide useful estimates of temperature change within 0.5K (0.25K) in Scenario D for at least 88% (68%) of the grid points where we have confidence; indicating that the methodology is transferable to other scenarios. Tests with the Geophysical Fluid Dynamics Laboratory GCM indicated, however, that a temperature ratio pattern may have to be developed for each GCM. The methodology, using a temperature ratio pattern from the 1990 IPCC Scenario A and involving IMAGE, gave gridded surface air temperature patterns for the 1992 IPCC radiative-forcing Scenarios C and E and the RIVM emission Scenario B; none of these scenarios has been simulated by ECHAM-1/LSG. The simulations reflect the uncertainty range of a future warming.The work reported by the authors was carried out during their stay at the project Forestry and Climate Change of the International Institute for Applied Systems Analysis, Laxenburg, Austria.     

Oxygen isotope heterogeneity and disequilibria of olivine crystals in large volume Holocene basalts from Iceland: Evidence for magmatic digestion and erosion of Pleistocene hyaloclastites

This work considers petrogenesis of the largest Holocene basaltic fissure eruptions of Iceland, which are also the largest in the world: Laki (1783-84 AD, 15 km³), Eldgjá (934 AD, 18 km³), Veidivötn (900, 1480 AD, multiple eruptions, >2 km³), Núpahraun (ca. 4000 BP, >1 km³) and Thjórsárhraun (ca 8000 BP, >20 km³). We present oxygen isotope laser fluorination analyses of 55 individual and bulk olivine crystals, coexisting individual and bulk plagioclase phenocrysts, and their host basaltic glasses with average precision of better than 0.1‰ (1SD). We also report O isotope analyses of cores and rims of 61 olivine crystals by SIMS with average precision on single spots of 0.24‰ (1SD) in 13 samples coupled with electron microprobe data for major and trace elements in these olivines. Within each individual sample, we have found that basaltic glass is relatively homogeneous with respect to oxygen isotopes, plagioclase phenocrysts exhibit crystal to crystal variability, while individual olivines span from the values in equilibrium with the low-δ¹⁸O matrix glass to those being three permil higher in δ¹⁸O than the equilibrium. Olivine cores with maximum value of 5.2‰ are found in many of these basalts and suggest that the initial magma was equilibrated with normal-δ¹⁸O mantle. No olivines or their intracrystalline domains are found with bulk or spot value higher than those found in MORB olivines. The δ¹⁸O variability of 0.3-3‰ exists for olivine grains from different lavas, and variable core-to-rim oxygen isotopic zoning is present in selected olivine grains. Many olivines in the same sample are not zoned, while a few grains are zoned with respect to oxygen isotopes and exhibit small core-to-core variations in Fe-Mg, Ni, Mn, Ca. Grains that are zoned in both Mg# and δ¹⁸O exhibit positive correlation of these two parameters. Electron microprobe analysis shows that most olivines equilibrated with the transporting melt, and thin Fe-richer rim is present around many grains, regardless of the degree of olivine-melt oxygen isotope disequilibrium.The preservation of isotopic and compositional zoning in selected grains, and subtle to severe Δ¹⁸O (melt-olivine) and Δ¹⁸O (plagioclase-olivine) disequilibria suggests rather short crystal residence times of years to centuries. Synglacially-altered upper crustal, tufaceous hyaloclastites of Pleistocene age serve as a viable source for low-δ¹⁸O values in Holocene basalts through assimilation, mechanical and thermal erosion, and devolatilization of stoped blocks. Cumulates formed in response to cooling during assimilation, and xenocrysts derived from hyaloclastites, contribute to the diverse δ¹⁸O crystalline cargo. The magma plumbing systems under each fissure are likely to include a network of interconnected dikes and sills with high magma flow rates that contribute to the efficacy of magmatic erosion of large quantities (10-60% mass) of hyaloclastites required by isotopic mass balance.Olivine diversity and the pervasive lack of phenocryst-melt oxygen isotopic equilibrium suggest that a common approach of analyzing bulk olivine for oxygen isotopes, as a proxy for the basaltic melt or to infer mantle δ¹⁸O value, needs to proceed with caution. The best approach is to analyze olivine crystals individually and demonstrate their equilibrium with matrix.