Distribution of icy particles across Enceladus' surface as derived from Cassini-VIMS measurements

The surface of Enceladus consists almost completely of water ice. As the band depths of water ice absorptions are sensitive to the size of particles, absorptions can be used to map variations of icy particles across the surface. The Visual and Infrared Mapping Spectrometer (VIMS) observed Enceladus with a high spatial resolution during three Cassini flybys in 2005 (orbits EN 003, EN 004 and EN 011). Based on these data we measured the band depths of water ice absorptions at 1.04, 1.25, 1.5, and 2 μm. These band depths were compared to water ice models that represent theoretically calculated reflectance spectra for a range of particle diameters between 2 μm and 1 mm. The agreement between the experimental (VIMS) and model values supports the assumption that pure water ice characterizes the surface of Enceladus and therefore that variations in band depth correspond to variations in water ice particle diameters. Our measurements show that the particle diameter of water ice increases toward younger tectonically altered surface units with the largest particles exposed in relatively “fresh” surface material. The smallest particles were generally found in old densely cratered terrains. The largest particles (∼0.2 mm) are concentrated in the so called “tiger stripes” at the south pole. In general, the particle diameters are strongly correlated with geologic features and surface ages, indicating a stratigraphic evolution of the surface that is caused by cryovolcanic resurfacing and impact gardening.     

Partitioning of trace elements between rutile and silicate melts: Implications for subduction zones

Mineral/melt trace element partition coefficients were determined for rutile (TiO₂) for a large number of trace elements (Zr, Hf, Nb, Ta, V, Co, Cu, Zn, Sr, REE, Cr, Sb, W, U, Th). Whilst the high field strength elements (Zr, Hf, Nb, Ta) are compatible in rutile, other studied trace elements are incompatible (Sr, Th, REE). In all experiments we found D_Ta > D_Nb, D_Hf > D_Zr and D_U > D_Th. Partition coefficients for some polyvalent elements (Sb, W, and Co) were sensitive to oxygen fugacity. Melt composition exerts a strong influence on HFSE partition coefficients. With increasing polymerization of the melt, rutile/melt partition coefficients for the high field strength elements Zr, Hf, Nb and Ta increase about an order of magnitude. However, D_Nb/D_Ta and D_Hf/D_Zr are not significantly affected by melt composition. Because D_U ? D_Th, partial melting of rutile-bearing eclogite in subducted lithosphere may cause excesses of ²³⁰Th over ²³⁸U in some island arc lavas, whereas dehydration of subducted lithosphere may cause excesses of ²³⁸U over ²³⁰Th. From our partitioning results we infer partition coefficients for protactinium (Pa) which we predict to be much lower than previously anticipated. Contrary to previous studies, our data imply that rutile should not significantly influence observed ²³¹Pa-²³⁵U disequilibria in certain volcanic rocks.     

The reaction MgCr2O4 + SiO2 = Cr2O3 + MgSiO3 and the free energy of formation of magnesiochromite (MgCr2O4)

Low-temperature heat capacity measurements for MgCr₂O₄ have only been performed down to 52 K, and the commonly quoted third-law entropy at 298 K (106 J K⁻¹ mol⁻¹) was obtained by empirical extrapolation of these measurements to 0 K without considering the magnetic or electronic ordering contributions to the entropy. Subsequent magnetic measurements at low temperature reveal that the Néel temperature, at which magnetic ordering of the Cr³⁺ ions in MgCr₂O₄ occurs, is at ∼15 K. Hence a substantial contribution to the entropy of MgCr₂O₄ has been missed. We have determined the position of the near-univariant reaction MgCr₂O₄+SiO₂=MgSiO₃+Cr₂O₃. The reaction, which has a small positive slope in P-T space, has been bracketed at 100 K intervals between 1273 and 1773 K by reversal experiments. The reaction is extremely sluggish, and lengthy run times with a flux (H₂O, BaO-B₂O₃ or K₂O-B₂O₃) are needed to produce tight reversal brackets. The results, combined with assessed thermodynamic data for Cr₂O₃, MgSiO₃ and SiO₂, give the entropy and enthalpy of formation of MgCr₂O₄ spinel. As expected, our experimental results are not in good agreement with the presently available thermodynamic data. We obtain Δ _f H ^∘ ₂₉₈=−1759.2±1.5 kJ mol⁻¹ and S ^∘ ₂₉₈=122.1±1.0 J K⁻¹ mol⁻¹ for MgCr₂O₄. This entropy is some 16 J K⁻¹ mol⁻¹ more than the calorimetrically determined value, and implies a value for the magnetic entropy of MgCr₂O₄ consistent with an effective spin quantum number (S') for Cr³⁺ of 1/2 rather than the theoretical 3/2, indicating, as in other spinels, spin quenching. Received: 9 May 1997 / Accepted: 28 July 1997     

Early Palaeozoic deep subduction of continental crust in the Kyrgyz North Tianshan: evidence from Lu–Hf garnet geochronology and petrology of mafic dikes

High-pressure and ultrahigh-pressure (UHP) eclogite-bearing metamorphic assemblages in the North Tianshan of Kyrgyzstan are known from the Aktyuz and Makbal areas, where eclogites and garnet amphibolites are associated with continental rocks such as granitoid gneisses in Aktyuz and shallow-water clastic (passive margin?) metasediments in Makbal. We present the first Lu–Hf isotope data for an eclogite and two garnet amphibolite samples from the two metamorphic terranes which, combined with petrological analysis, tightly constrain the age of high-pressure metamorphism in the Kyrgyz North Tianshan. A five-point isochron for an Aktyuz eclogite sample provides a Lu–Hf age of 474.3 ± 2.2 Ma, and a four-point isochron on a Makbal sample corresponds to 470.1 ± 2.5 Ma. A prograde, subduction-related path is inferred for both samples with peak P–T conditions ranging from 1.4 to 1.6 GPa and 610–620 °C. A further Makbal sample provided a significantly older Lu–Hf age of 486 ± 5.4 Ma, most likely due to late alteration in the sample (late addition of unradiogenic Hf). We conclude that garnet growth in all three samples occurred around ca. 474 Ma and that these rocks likely experienced UHP metamorphism contemporaneously. Our results support previous geochronological evidence for an Early Ordovician collision belt in the North Tianshan and allow refinement of a tectonic model involving subduction of thinned continental crust to considerable depth along the margin of a small microcontinent.     

Observations of Atmospheric Methane and Carbon Dioxide Mixing Ratios: Tall-Tower or Mountain-Top Stations?

Mountain-top observations of greenhouse gas mixing ratios may be an alternative to tall-tower measurements for regional scale source and sink estimation. To investigate the equivalence or limitations of a mountain-top site as compared to a tall-tower site, we used the unique opportunity of comparing in situ measurements of methane (\(\hbox {CH}_{4}\)) and carbon dioxide (\(\hbox {CO}_{2}\)) mixing ratios at a mountain top (986 m above sea level, a.s.l.) with measurements from a nearby (distance 28.4 km) tall tower, sampled at almost the same elevation (1009 m a.s.l.). Special attention was given to, (i) how local wind statistics and greenhouse gas sources and sinks at the mountain top influence the observations, and (ii) whether mountain-top observations can be used as for those from a tall tower for constraining regional greenhouse gas emissions. Wind statistics at the mountain-top site are clearly more influenced by local flow systems than those at the tall-tower site. Differences in temporal patterns of the greenhouse gas mixing ratios observed at the two sites are mostly related to the influence of local sources and sinks at the mountain-top site. Major influences of local sources can be removed by applying a statistical filter (\(5{\mathrm{th}}\) percentile) or a filter that removes periods with unfavourable flow conditions. In the best case, the bias in mixing ratios between the mountain-top and the tall-tower sites after the application of the wind filter was \({-}0.0005\pm 0.0010\) ppm for methane (September, 0000–0400 UTC) and \(0.11\pm 0.18\) ppm for \(\hbox {CO}_{2}\) (February, 1200–1600 UTC). Temporal fluctuations of atmospheric \(\hbox {CH}_{4}\) and \(\hbox {CO}_{2}\) mixing ratios at both stations also showed good agreement (apart from \(\hbox {CO}_{2}\) during summertime) as determined by moving bi-weekly Pearson correlation coefficients (up to 0.96 for \(\hbox {CO}_{2}\) and 0.97 for \(\hbox {CH}_{4}\)). When only comparing mixing ratios minimally influenced by local sources (low bias and high correlation coefficients), our measurements indicate that mountain-top observations are comparable to tall-tower observations.     

Numerical Analysis of Groundwater Ridging Processes Considering Water‐Air Flow in a Hillslope

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.     

The pre-weichselian glaciations of North-West Europe

In North-West Europe no evidence has been found for glaciations older than the Elsterian. The Elsterian seems to consist of at least two separate ice advances with a minor interval in between. During the Elsterian in the area south of the Baltic Sea over 400 m deep channels were cut into the underlying substratum by a combination of glacial erosion and meltwater activity. The channels were active until the Late Elsterian and were subsequently filled with meltwater deposits and glaciolacustrine silt and clay. During the Saalian no comparative channels were being formed. Three different ice advances can be distinguished within the Saalian, the stratigraphical correlation of which is to be discussed. No convincing evidence was found for any interglacial within the Saalian sequence. There are good reasons to assume that the ice advances were only separated from each other by minor ice-free intervals.     

The Scaling Behaviour of a Turbulent Kinetic Energy Closure Model for Stably Stratified Conditions

We investigate the scaling behaviour of a turbulent kinetic energy (TKE) closure model for stably stratified conditions. The mixing length scale for stable stratification is proportional to the ratio of the square root of the TKE and the local Brunt–Väisälä frequency, which is a commonly applied formulation. We analyze the scaling behaviour of our model in terms of traditional Monin–Obukov Similarity Theory and local scaling. From the model equations, we derive expressions for the stable limit behaviour of the flux–gradient relations and other scaling quantities. It turns out that the scaling behaviour depends on only a few model parameters and that the results obey local scaling theory. The analytical findings are illustrated with model simulations for the second GABLS intercomparison study. We also investigate solutions for the case in which an empirical correction function is used to express the eddy diffusivity for momentum as a function of the Richardson number (i.e. an increasing turbulent Prandtl number with stability). In this case, it seems that for certain parameter combinations the model cannot generate a steady-state solution. At the same time, its scaling behaviour becomes unrealistic. This shows that the inclusion of empirical correction functions may have large and undesired consequences for the model behaviour.