Latent heat effects of the major mantle phase transitions on low-angle subduction

Very low to zero shallow dip angles are observed at several moderately young subduction zones with an active trenchward moving overriding plate. We have investigated the effects of latent heat for this situation, where mantle material is pushed through the major mantle phase transitions during shallow low-angle subduction below the overriding plate. The significance of the buoyancy forces, arising from the latent heat effects, on the dynamics of the shallowly subducting slab is examined by numerical modeling. When a 32-Ma-old slab is overridden with 2.5 cm/yr by a continent, flat subduction occurs with a 4–5 cm/yr convergence rate. When latent heat is included in the model, forced downwellings cause a thermal anomaly and consequently thermal and phase buoyancy forces. Under these circumstances, the flat slab segment subducts horizontally about 350 km further and for about 11 Ma longer than in the case without latent heat, before it breaks through the 400-km phase transition. The style of subduction strongly depends on the mantle rheology: increasing the mantle viscosity by one order of magnitude can change the style of subduction from steep to shallow. Similarly, an overriding velocity of less than 1 cm/yr leads to steep subduction, which gradually changes to flat subduction when increasing the overriding velocity. However, these model parameters do not change the aforementioned effect of the latent heat, provided that low-angle subduction occurs. In all models latent heat resulted in a substantial increase of the flat slab length by 300–400 km. Varying the olivine–spinel transition Clapeyron slope γ from 1 to 6 MPa/K reveals a roughly linear relation between γ and the horizontal length of the slab. Based on these results, we conclude that buoyancy forces due to latent heat of phase transitions play an important role in low-angle subduction below an overriding plate.     

Global spherical harmonic computation by two-dimensional Fourier methods

A method is presented for performing global spherical harmonic computation by two-dimensional Fourier transformations. The method goes back to old literature (Schuster 1902) and tackles the problem of non-orthogonality of Legendre-functions, when discretized on an equi-angular grid. Both analysis and synthesis relations are presented, which link the spherical harmonic spectrum to a two-dimensional Fourier spectrum. As an alternative, certain functions of co-latitude are introduced, which are orthogonal to discretized Legendre functions. Several independent Fourier approaches for spherical harmonic computation fit into our general scheme.     

On the Existence of Coverage and Integration Time Regimes in Bistatic SAR Configurations

Bistatic synthetic aperture radar (SAR) is an extension of traditional monostatic SAR, which increases the flexibility in designing SAR missions. We describe a scheme for the computation of integration time and azimuth coverage of bistatic SARs based on space-time diagrams. A classification of bistatic SAR configurations is introduced in terms of size and velocity on the ground of antenna footprints. Bistatic SAR regimes are also identified.     

Late Pleistocene–Holocene rise and collapse of Lake Suguta,northern Kenya Rift

The Late Pleistocene to Middle Holocene African Humid Period (AHP) was characterized by dramatic hydrologic fluctuations in the tropics. A better knowledge of the timing, spatial extent, and magnitude of these hydrological fluctuations is essential to decipher the climate-forcing mechanisms that controlled them. The Suguta Valley (2°N, northern Kenya Rift) has recorded extreme environmental changes during the AHP. Extensive outcrops of lacustrine sediments, ubiquitous wave-cut notches, shorelines, and broad terrace treads along the valley margins are the vestiges of Lake Suguta, which once filled an 80 km long and 20 km wide volcano–tectonic depression. Lake Suguta was deep between 16.5 and 8.5 cal ka BP. During its maximum highstand, it attained a water depth of ca 300 m, a surface area of ca 2150 km², and a volume of ca 390 km³. The spatial distribution of lake sediments, the elevation of palaeo-shorelines, and other geomorphic evidences suggest that palaeo-Lake Suguta had an overflow towards the Turkana basin to the north. After 8.5 cal ka BP, Lake Suguta abruptly disappeared. A comparison of the Lake Suguta water-level curve with other reconstructed water levels from the northern part of the East African Rift System shows that local insolation, which is dominated by precessional cycles, may have controlled the timing of lake highstands in this region. Our data show that changes of lake levels close to the Equator seem to be driven by fluctuations of spring insolation, while fluctuations north of the Equator are apparently related to variations in summer insolation. However, since these inferred timings of lake-level changes are mostly based on the radiocarbon dating of carbonate shells, which may have been affected by a local age reservoir, alternative dating methods are needed to support this regional synthesis. Between 12.7 and 11.8 cal ka BP, approximately during the Northern Hemisphere high-latitude Younger Dryas, the water level of Lake Suguta fell by ca 50 m, suggesting that remote influences also affected local hydrology.     

Structure, thermodynamics, and diffusion in CaAl2Si2O8 liquid from first-principles molecular dynamics

We have performed first-principles molecular dynamics simulations of CaAl₂Si₂O₈ (anorthite) liquid at pressures up to 120 GPa and temperatures of 3000, 4000 and 6000 K. At the lowest degrees of compression the liquid is seen to accommodate changes in density through decreasing the abundance of 3- and 4-membered rings, while increases in coordination of network forming cations take effect at somewhat higher degrees of compression. Results are fit to a fundamental thermodynamic relation with 4th order finite strain and 1st order thermal variable expansions. Upon compression by a factor of two, the Grüneisen parameter (γ) is found to increase continuously from 0.35 to 1.10. Weak temperature dependence in γ is thermodynamically consistent with a slight decrease in isochoric heat capacity (C_V), for which values of between 4.4 and 5.2 Nk_B are obtained, depending on the temperature. Pressure and temperature dependence of self-diffusivities is found to be well represented by an Arrhenius relation, except at 3000 K and pressures lower than 5 GPa, where self-diffusivities of Si, Al, and O increase with pressure. Analysis of the lifetimes of individual coordination species reveals that this phenomenon arises due to the disproportionately high stability of 4-fold coordinated Si, and to a lesser extent 4-fold coordinated Al. Our results represent a marked improvement in accuracy and reliability in describing the physics of CaAl₂Si₂O₈ liquid at deep mantle pressures, pointing the way to a general thermodynamic model of melts at extreme pressures and temperatures relevant to planetary-scale magma oceans and deep mantle partial melting.     

Review: Thermal water resources in carbonate rock aquifers

The current knowledge on thermal water resources in carbonate rock aquifers is presented in this review, which also discusses geochemical processes that create reservoir porosity and different types of utilisations of these resources such as thermal baths, geothermal energy and carbon dioxide (CO₂) sequestration. Carbonate aquifers probably constitute the most important thermal water resources outside of volcanic areas. Several processes contribute to the creation of porosity, summarised under the term hypogenic (or hypogene) speleogenesis, including retrograde calcite solubility, mixing corrosion induced by cross-formational flow, and dissolution by geogenic acids from deep sources. Thermal and mineral waters from karst aquifers supply spas all over the world such as the famous bath in Budapest, Hungary. Geothermal installations use these resources for electricity production, district heating or other purposes, with low CO₂ emissions and land consumption, e.g. Germany’s largest geothermal power plant at Unterhaching near Munich. Regional fault and fracture zones are often the most productive zones, but are sometimes difficult to locate, resulting in a relatively high exploration uncertainty. Geothermal installations in deep carbonate rocks could also be used for CO₂ sequestration (carbonate dissolution would partly neutralise this gas and increase reservoir porosity). The use of geothermal installations to this end should be further investigated.     

Comparing seven candidate mission configurations for temporal gravity field retrieval through full-scale numerical simulation

The goal of this contribution is to focus on improving the quality of gravity field models in the form of spherical harmonic representation via alternative configuration scenarios applied in future gravimetric satellite missions. We performed full-scale simulations of various mission scenarios within the frame work of the German joint research project “Concepts for future gravity field satellite missions” as part of the Geotechnologies Program, funded by the German Federal Ministry of Education and Research and the German Research Foundation. In contrast to most previous simulation studies including our own previous work, we extended the simulated time span from one to three consecutive months to improve the robustness of the assessed performance. New is that we performed simulations for seven dedicated satellite configurations in addition to the GRACE scenario, serving as a reference baseline. These scenarios include a “GRACE Follow-on” mission (with some modifications to the currently implemented GRACE-FO mission), and an in-line “Bender” mission, in addition to five mission scenarios that include additional cross-track and radial information. Our results clearly confirm the benefit of radial and cross-track measurement information compared to the GRACE along-track observable: the gravity fields recovered from the related alternative mission scenarios are superior in terms of error level and error isotropy. In fact, one of our main findings is that although the noise levels achievable with the particular configurations do vary between the simulated months, their order of performance remains the same. Our findings show also that the advanced pendulums provide the best performance of the investigated single formations, however an accuracy reduced by about 2–4 times in the important long-wavelength part of the spectrum (for spherical harmonic degrees ${<}50$ ), compared to the Bender mission, can be observed. Concerning state-of-the-art mission constraints, in particular the severe restriction of heterodyne lasers on maximum range-rates, only the moderate Pendulum and the Bender-mission are beneficial options, of course in addition to GRACE and GRACE-FO. Furthermore, a Bender-type constellation would result in the most accurate gravity field solution by a factor of about 12 at long wavelengths (up to degree/order 40) and by a factor of about 200 at short wavelengths (up to degree/order 120) compared to the present GRACE solution. Finally, we suggest the Pendulum and the Bender missions as candidate mission configurations depending on the available budget and technological progress.     

Trends,rhythms and events in Plio-Pleistocene African climate

We analyzed published records of terrigenous dust flux from marine sediments off subtropical West Africa, the eastern Mediterranean Sea, and the Arabian Sea, and lake records from East Africa using statistical methods to detect trends, rhythms and events in Plio-Pleistocene African climate. The critical reassessment of the environmental significance of dust flux and lake records removes the apparent inconsistencies between marine vs. terrestrial records of African climate variability. Based on these results, major steps in mammalian and hominin evolution occurred during episodes of a wetter, but highly variable climate largely controlled by orbitally induced insolation changes in the low latitudes.     

Multi-scale characterization of a complex karst and alluvial aquifer system in southern Germany using a combination of different tracer methods

Water suppliers face major challenges such as climate change and population growth. To prepare for the future, detailed knowledge of water resources is needed. In southern Germany, the state water supplier Zweckverband Landeswasserversorgung provides 3 million people with drinking water obtained from a complex karst and alluvial aquifer system and the river Danube. In this study, a combination of different tracing techniques was used with the goal of a multi-scale characterization of the aquifer system and to gain additional knowledge about groundwater flow toward the extraction wells in the Danube Valley. For the small-scale characterization, selected groundwater monitoring wells were examined using single-borehole dilution tests. With these tests, a wide range of flow behavior could be documented, including fast outflow within just a few hours in wells with good connection to the aquifer, but also durations of many weeks in low-permeability formations. Vertical flow, caused by multiple flow horizons or uprising groundwater, was detected in 40% of the tested wells. A regional multi-tracer test with three injections was used to investigate the aquifer on a large scale. For the highly karstified connection between a swallow hole and a spring group, high flow velocities of around 80 m/h could be documented. Exceptionally delayed arrivals, 250 and 307 days after the injection, respectively showing maximum velocities of 0.44 and 0.39 m/h, were observed in an area where low-permeability sediments overlay the karst conduits. With the chosen methods, a distinct heterogeneity caused by the geological setting could be documented on both scales.