Simulation of the development of gypsum maze caves

The development of gypsum maze caves under artesian conditions has been simulated. The numerical model simulations show that the evolution of maze caves in this type of setting requires structural preferences such as laterally extended fissure networks in a horizon of the gypsum layer. Without any structural preferences vertical shafts rather than maze caves are predicted to develop. The most important stage for the development of horizontal caves under artesian conditions is found to be the initial karstification period. During this period the structure of the mature conduit system is established. The solutional enlargement of conduits is spatially extended, total dissolution rates are higher than the later ones.     

Electrical conductivity measurements on olivines and pyroxenes under defined thermodynamic activities as a function of temperature and pressure

Electrical conductivities of Ni₂SiO₄, Fe₂SiO₄, and MgSiO₃ were measured on synthetic powders in the temperature range 340° to 1,100° C and at pressures up to 20 kbars. For ternary compounds such as olivines and pyroxenes the control of two further variables, like the chemical activities of two components are needed, besides temperature and pressure. The activities of the corresponding binary oxides were controlled by equilibrating the samples with their neighbour-phases. Control of the oxygen partial pressure was achieved by buffer techniques. From the slopes of the lg σ vs. 1/T lines the activation energies were calculated for 10 kbar: 0.56 eV and 2.7 eV for Ni₂SiO₄ in equilibrium with SiO₂ and Ni/NiO-buffer for the temperature range 500°–800°C and 800°–1,000°C resp. 0.52 eV for Fe₂SiO₄ in equilibrium with SiO₂ and metallic iron, and 0.38 eV in equilibrium with SiO₂ and magnetite; 1.11 eV for MgSiO₃ in equilibrium with SiO₂, and 1.25 eV in equilibrium with Mg₂SiO₄.     

Improved GRACE kinematic orbit determination using GPS receiver clock modeling

Kinematic positions of Low Earth Orbiters based on GPS tracking are frequently used as pseudo-observations for single satellite gravity field determination. Unfortunately, the accuracy of the satellite trajectory is partly limited because the receiver synchronization error has to be estimated along with the kinematic coordinates at every observation epoch. We review the requirements for GPS receiver clock modeling in Precise Point Positioning (PPP) and analyze its impact on kinematic orbit determination for the two satellites of the Gravity Recovery and Climate Experiment (GRACE) mission using both simulated and real data. We demonstrate that a piecewise linear parameterization can be used to model the ultra-stable oscillators that drive the GPS receivers on board of the GRACE satellites. Using such a continuous clock model allows position estimation even if the number of usable GPS satellites drops to three and improves the robustness of the solution with respect to outliers. Furthermore, simulations indicate a potential accuracy improvement of the satellite trajectory of at least 40 % in the radial direction and up to 7 % in the along-track and cross-track directions when a 60-s piecewise linear clock model is estimated instead of epoch-wise independent receiver clock offsets. For PPP with real GRACE data, the accuracy evaluation is hampered by the lack of a reference orbit of significantly higher accuracy. However, comparisons with a smooth reduced-dynamic orbit indicate a significant reduction of the high-frequency noise in the radial component of the kinematic orbit.     

Calculated greenschist facies mineral equilibria in the system CaO−FeO−MgO−Al2O3−SiO2−CO2−H2O

The kinetic problems associated with the experimental determination of reactions among complex solidsolution phases at low temperatures have hindered our understanding of the phase relations in greenschist facies rocks. In the absence of reliable experimental data, we have used the new, expanded internally-consistent thermodynamic dataset of Holland and Powell (1990), to present calculated phase equilibria for the system CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ (CaFMASCH) with quartz in eccess, in the range 400°–500°C at low to intermediate pressures, involving the minerals amphibole, chlorite, anorthite, clinozoisite, dolomite, chloritoid, garnet, margarite, andalusite, and calcite. By solving independent sets of non-linear equations formed from equilibrium relationships, we calculate not only the loci of reactions in pressuretemperature-x(CO₂) space, but also the compositions of coexisting minerals in terms of the substitutions, FeMg_-1 and (Fe,Mg)SiAl_-1Al_-1. Invariant, univariant and divariant equilibria are calculated and discussed in relation to naturally-occurring greenschist facies metabasic and siliceous dolomitic mineral assemblages. We thus avoid the use of activity-corrected curves so commonly presented in the literature as a substitute for genuine univariant phase diagram boundaries.     

The pace of Holocene vegetation change – testing for synchronous developments

Mid to high latitude forest ecosystems have undergone several major compositional changes during the Holocene. The temporal and spatial patterns of these vegetation changes hold potential information to their causes and triggers. Here we test the hypothesis that the timing of vegetation change was synchronous on a sub-continental scale, which implies a common trigger or a step-like change in climate parameters. Pollen diagrams from selected European regions were statistically divided into assemblage zones and the temporal pattern of the zone boundaries analysed. The results show that the temporal pattern of vegetation change was significantly different from random. Times of change cluster around 8.2, 4.8, 3.7, and 1.2 ka, while times of higher than average stability were found around 2.1 and 5.1 ka. Compositional changes linked to the expansion of Corylus avellana and Alnus glutinosa centre around 10.6 and 9.5 ka, respectively. A climatic trigger initiating these changes may have occurred 0.5 to 1 ka earlier, respectively. The synchronous expansion of C. avellana and A. glutinosa exemplify that dispersal is not necessarily followed by population expansion. The partly synchronous, partly random expansion of A. glutinosa in adjacent European regions exemplifies that sudden synchronous population expansions are not species specific traits but vary regionally.     

Nonlinear Flow Process: A New Package to Compute Nonlinear Flow in MODFLOW

A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW‐2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells.     

Deep solid-state equilibration and deep melting of plagioclase-free spinel peridotite from the slow-spreading Mid-Atlantic Ridge, ODP Leg 153

A petrological investigation of abyssal, plagioclase-free spinel peridotite drilled during ODP cruise 153 in the North Atlantic revealed that the peridotite represent refractory, partial residual mantle material that experienced depletion of incompatible trace elements during upper mantle melting. The degree of partial melting as estimated from spinel compositions was c. 12%. Fractionated middle and heavy rare earth elements imply polybaric melting, with c. 1–4% initial melting in the garnet peridotite stability field and subsequent partial melting of ~7–10% in the spinel peridotite stability field. Geothermobarometric investigations revealed that the solid-state equilibration of the spinel peridotite occurred at some 1,100–1,150°C and c. 20–23 kbar, corresponding to an equilibration depth of c. 70?±?5 km and an unusually low thermal gradient of some 11–17°C/km. A thermal re-equilibration of the peridotite occurred at ~850–1,000°C at similar depths. Naturally, the initial mantle melting in the garnet-peridotite stability field must have commenced at depths greater than 70?±?5 km. It is likely that the residual peridotite rose rapidly through the lithospheric cap towards the ridge axis. The exhumation of the abyssal peridotite occurred, at least in parts, via extensional detachment faulting. Given the shallow to moderate dip angles of the fault surfaces, the exhumation of the peridotite from its equilibration depth would imply an overall ridge-normal horizontal displacement of c. 50–160 km if tectonic stretching and detachment faulting were the sole exhumation mechanism.