A Rhaetian 40Ar/39Ar age for the Rochechouart impact structure (France) and implications for the latest Triassic sedimentary record

Abstract– ⁴⁰Ar/³⁹Ar dating of potassium feldspar (primary spherulitic‐blocky and secondary idiomorphic K‐feldspar) separated from impact‐metamorphosed gneiss found near Videix in the western central part of the Rochechouart impact structure (NW Massif Central, France) yielded a Rhaetian combined age of 201 ± 2 Ma (2σ), indistinguishable within uncertainty from the age of the Triassic/Jurassic boundary. Ballen quartz intergrown with the primary K‐feldspar indicates post‐shock temperatures exceeding approximately 1000 °C that affected the precursor gneiss. Geochemically, both feldspar types represent essentially pure potassium end‐members. Apart from the approximately 15 km diameter impact deposit area, the youngest crystallization age known for basement rocks in this part of the Massif Central is approximately 300 Ma. No endogenic magmatic‐thermal events are known to have occurred later in this region. The K‐feldspar recrystallized from local feldspar melts and superimposed post‐shock hydrothermal crystallization, probably within some thousands of years after the impact. It is, therefore, suggested that the ⁴⁰Ar/³⁹Ar age for the Videix gneiss (as a potassic “impact metasomatite”) dates the Rochechouart impact, in consistence with evidence for K‐metasomatism in the Rochechouart impactites. The new age value is distinctly younger than the previously obtained Karnian–Norian age for Rochechouart and, thus, contradicts the Late Triassic multiple impact theory postulated some years ago. In agreement with the paleogeographic conditions in the western Tethys domain around the Triassic/Jurassic boundary, the near‐coastal to shallow marine Rochechouart impact is compatible with the formation of seismites and tsunami deposits in the latest Triassic of the British Isles and possible related deposits in other parts of Europe.     

Long-term soil moisture content estimation using satellite and climate data in agricultural area of Mongolia

The purpose of this study is to estimate long-term SMC and find its relation with soil moisture (SM) of climate station in different depths and NDVI for the growing season. The study area is located in agricultural regions in the North of Mongolia. The Pearson’s correlation methodology was used in this study. We used MODIS and SPOT satellite data and 14 years data for precipitation, temperature and SMC of 38 climate stations. The estimated SMC from this methodology were compared with SM from climate data and NDVI. The estimated SMC was compared with SM of climate stations at a 10-cm depth (r² = 0.58) and at a 50-cm depth (r² = 0.38), respectively. From the analysis, it can be seen that the previous month’s SMC affects vegetation growth of the following month, especially from May to August. The methodology can be an advantageous indicator for taking further environmental analysis in the region.     

Intrusion and its implication for subsidence: A case from the Baiyun Sag, on the northern margin of the South China Sea

We develop a series of simple numerical models to explain the anomalous subsidence and deposition phenomena on the northern continental margin of the South China Sea, in particular in the Baiyun Sag. The results suggest that a short-period high rate deposition of around 17 Ma is related to a rapid subsidence event, which may be due to episodic emplacement of a dense intrusion. Necking and gravity models indicate that in the basement of the Baiyun Sag, there is a dense zone that is 100–200 kg/m³ more dense than the surrounding country rock. Considering its high magnetic intensity and regional igneous activity, the dense zone is thought to be related to a phase of basalt intrusion that may have taken place around 17 Ma. Thermal and subsidence models indicate that a hot denser intrusion can cause significant subsidence immediately after the intrusion. The subsidence rate then slows down with cooling, thus becoming consistent with the observed subsidence curves at around 17 Ma. The results also indicate that the lithospheric strength under the Baiyun Sag is negligible, and that the high-velocity layer in the lowermost crust may be not an original part of the pre-rift crust. Instead, it is thought to be underplated intrusion emplaced at around 30 Ma when the continental margin broke up.     

GETEMME—a mission to explore the Martian satellites and the fundamentals of solar system physics

GETEMME (Gravity, Einstein??s Theory, and Exploration of the Martian Moons?? Environment), a mission which is being proposed in ESA??s Cosmic Vision program, shall be launched for Mars on a Soyuz Fregat in 2020. The spacecraft will initially rendezvous with Phobos and Deimos in order to carry out a comprehensive mapping and characterization of the two satellites and to deploy passive Laser retro-reflectors on their surfaces. In the second stage of the mission, the spacecraft will be transferred into a lower 1500-km Mars orbit, to carry out routine Laser range measurements to the reflectors on Phobos and Deimos. Also, asynchronous two-way Laser ranging measurements between the spacecraft and stations of the ILRS (International Laser Ranging Service) on Earth are foreseen. An onboard accelerometer will ensure a high accuracy for the spacecraft orbit determination. The inversion of all range and accelerometer data will allow us to determine or improve dramatically on a host of dynamic parameters of the Martian satellite system. From the complex motion and rotation of Phobos and Deimos we will obtain clues on internal structures and the origins of the satellites. Also, crucial data on the time-varying gravity field of Mars related to climate variation and internal structure will be obtained. Ranging measurements will also be essential to improve on several parameters in fundamental physics, such as the Post-Newtonian parameter ?? as well as time-rate changes of the gravitational constant and the Lense-Thirring effect. Measurements by GETEMME will firmly embed Mars and its satellites into the Solar System reference frame.     

The Greenland ice sheet and greenhouse warming

Increased melting on glaciers and ice sheets and rising sea level are often mentioned as important aspects of the anticipated greenhouse warming of the earth's atmosphere. This paper deals with the sensitivity of Greenland's ice mass budget and presents a tentative projection of the Greenland component of future sea level rise for the next few hundred years. To do this, the ‘Villach II temperature scenario’ is prescribed,output from a comprehensive mass balance model is used to drive a high-resolution 3-D thermomechanic model of the ice sheet.The mass balance model consists of two parts: the accumulation part is based on presently observed values and is forced by changes in mean anr tempeerature. The ablation model is based on the degree-day method and accounts for daily and annual temperature cycle, a different degree-day factor for ice and snow melting and superimposed ice formation. Under present-day climatic conditions, the following total mass balance results (in ice equivalent per years): 599.3 × 10⁹ m³ of accumulation, 281.7 × 10⁹m³ of runoff assuming a balanced budget, 317.6 × 10⁹m³ of iceberg calving. A 1K uniform warming is then calculated to increase the runoff by 119.5 × 10⁹ m³. Since accumulation also increases by 32 × 10⁹ m³, this leads to reduction of the total mass balance by 887.5 × 10⁹ m³ of ice, corresponding to a sea level rise of 0.22 mm/yr. For temperature increase larger than 2.7 K, runoff, exceeds accumulation, and if ice sheet dynamics were to remain unchanged, this would add an extra amount of 0.8 mmyr to the worl's oceans.Imposing the Villach II scenario (warming up to 4.23 K) and accumulating mass balance changes forward in time (static response) would then result in a global sea level rise of 7.1 cm by 2100 AD, but this figure may go up to as much as 40 cm per century in case the warming is doubled. In a subsequent dynamic model involving the ice flow, the ice sheet is found to produce a counteracting effect by dynamically producing steeper slopes at the margin, thereby reducing the area over which runoff can take place. This effect is particularly apparent in the northeastern part of the ice sheet, and is also more pronounced for the smaller temperature perturbations. Nevertheless, all these experiments certainly highlight the vulnerability of the Greenland ice sheet with respect to a climatic warming.     

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol + hy-normative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative sequence is exposed in the Dents de Bertol area (center of intrusion). P–T calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 °C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#_cpx = 83–89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment.     

Accuracy and efficiency of time integration methods for 1D diffusive wave equation

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.     

Curvature analysis of triangulated surfaces in structural geology

This paper addresses the problem of characterizing the shape of a geological surface on the basis of its principal curvatures. The surface is assumed to be modeled as a set of adjacent triangles defined by the location of their vertices and a method is proposed for estimating numerically the principal curvatures at the vertices of the triangles using a local C² interpolant. Also shown is how principal curvatures can be useful for studying the deformation of a geological surface (with application to 3D balanced unfolding), and analyzing the folding or faulting of the interface between two adjacent layers.