A practical review of geostatistical processing applied to geophysical data: methods and applications

Nowadays, geostatistics is commonly applied for numerous gridding or modelling tasks. However, it is still under used and unknown for classical geophysical applications. This paper highlights the main geostatistical methods relevant for geophysical issues, for instance to improve the quality of seismic data such as velocity cubes or interpreted horizons. These methods are then illustrated through four examples. The first example, based on a gravity survey presents how a geostatistical interpolation can be used to filter out a global trend, in order to better define real anomalies. In the second case study, dedicated to refraction surveying, geostatistical filtering is used to filter out acquisition artefacts and identify the main geological structures. The third one is an example of porosity being integrated geostatistically with a seismic acoustic impedance map. The last example deals with geostatistical time to depth conversion; the interest of performing geostatistical simulations is finally discussed.     

Integration of Local Observations into the One Dimensional Fog Model PAFOG

The numerical prediction of fog requires a very high vertical resolution of the atmosphere. Owing to a prohibitive computational effort of high resolution three dimensional models, operational fog forecast is usually done by means of one dimensional fog models. An important condition for a successful fog forecast with one dimensional models consists of the proper integration of observational data into the numerical simulations. The goal of the present study is to introduce new methods for the consideration of these data in the one dimensional radiation fog model PAFOG. First, it will be shown how PAFOG may be initialized with observed visibilities. Second, a nudging scheme will be presented for the inclusion of measured temperature and humidity profiles in the PAFOG simulations. The new features of PAFOG have been tested by comparing the model results with observations of the German Meteorological Service. A case study will be presented that reveals the importance of including local observations in the model calculations. Numerical results obtained with the modified PAFOG model show a distinct improvement of fog forecasts regarding the times of fog formation, dissipation as well as the vertical extent of the investigated fog events. However, model results also reveal that a further improvement of PAFOG might be possible if several empirical model parameters are optimized. This tuning can only be realized by comprehensive comparisons of model simulations with corresponding fog observations.     

Magnetic anisotropy of HED and Martian meteorites and implications for the crust of Vesta and Mars

We investigate the petrofabric of crustal rocks from Mars and Vesta through the measurement of the anisotropy of the magnetic susceptibility (AMS) of achondrites. Previous data are integrated with new measurements to obtain a dataset that provide macroscopic information about the magnetic fabric of 41 meteorites of the howardite–eucrite–diogenite clan (HED, falls only) and 16 Martian meteorites. The interpretation takes into account the large contribution of paramagnetism to the magnetic susceptibility of these meteorites. We use a model that allows the computation of the anisotropy degree of the population of ferromagnetic grains and provides a quantitative proxy for the degree of shape preferential orientation of these grains in HED and Martian meteorites. The results also provide quantitative information about the shape of the magnetic fabric (prolate, oblate).In HED achondrites, the ductile FeNi grains are sensitive strain recorders and our magnetic fabric data provide the first quantitative insights to the strain history of the crustal rocks of Vesta. Most HED achondrites are breccias but display a strong and spatially coherent magnetic anisotropy, indicating that intense deformation of FeNi grains took place after brecciation. The average fabric of eucrites, howardites is oblate (i.e. the texture is foliated) whereas the fabric of diogenites is more neutral. The howardite results suggest the existence of an isotropic fraction of ferromagnetic minerals that can be ascribed to the presence of carbonaceous chondrite clasts that have preserved their original magnetic fabric. In this hypothesis, howardites have an intensity of petrofabric very similar to eucrites and diogenites. Thermal metamorphism (itself possibly impact-related) plus lithostatic compaction occurring after brecciation appears as the best candidate to explain the observed petrofabric in eucrites and diogenites, whereas compaction by hypervelocity impacts may be reponsible for the fabric of howardites.Martian meteorites may still possess their primary magmatic fabric. Among Martian meteorites, basaltic shergottites and nakhlites display an oblate fabric (foliated texture) with only limited variations among each group. Olivine–phyric shergottites have a neutral fabric that points to a different petrogenesis. Nakhlites have weaker fabric intensity than shergottites. The fabric intensity is comparable to what is classically observed in terrestrial volcanic and plutonic rocks.     

Evolution of a carbonate delta generated by gateway‐funnelling of episodic currents

Cool‐water carbonate sedimentation has dominated Mediterranean shelves since the Early Pliocene. Skeletal sand and gravel herein consist of remains of heterozoan organisms, which are susceptible to reworking due to weak early cementation in non‐tropical waters. This study documents the Lower Pleistocene carbonate wedge of Favignana Island (Italy), which prograded from a 5   km wide passage between two palaeo‐islands into a perpendicular, 10 to 15   km wide strait between the palaeo‐islands at one side and Sicily at the other during the Emilian highstand (1·6   Ma to 1·1   Ma). The clinoformed carbonate wedge, which is 50   m thick and 6   km long, formed by east/south‐east progradation of a platform on the submarine sill by currents that were funnelled between the two palaeo‐islands. Platform‐slope clinoforms evolved from initial aggradation (thin and low‐angle) into a progradation phase (thick and high‐angle). Both clinoform types are characterized by a bimodal facies stacking pattern defined by sedimentary structures created by: (i) subaqueous dunes associated with dilute subcritical currents; and (ii) upper‐flow‐regime bedforms associated with sediment‐laden supercritical turbidity currents. Focusing of episodic currents on the platform by funnelling between the islands controlled the downstream formation of a sediment body, here named carbonate delta. The carbonate delta interfingers with subaqueous dune deposits formed in the perpendicular strait. This study uses a reconstruction of bedform dynamics to unravel the evolution of this gateway‐related carbonate accumulation.     

Distributed hydrologic and hydraulic modelling with radar rainfall input: Reconstruction of the 8–9 September 2002 catastrophic flood event in the Gard region,France

On 8–9 September 2002, an extreme rainfall event caused by a stationary mesoscale convective system (MCS) occurred in the Gard region, France. Distributed hydrologic and hydraulic modelling has been carried out to assess and compare the various sources of data collected operationally and during the post-event field surveys. Distributed hydrological modelling was performed with n-TOPMODELs and assessed for ungauged basins with the discharge estimates of the post-event surveys. A careful examination of the occurrence in time and space of the flash floods over the head watersheds indicates that flooding was controlled by the trajectory of the convective part of the MCS. Stationarity of the MCS over the Gardon watershed (1858 km² at Remoulins) for 28 h was responsible for the exceptional magnitude of the flood at this scale. The flood dynamics were characterized by an extensive inundation of the Gardonnenque plain upstream of the Gardon Gorges resulting in a significant peak flow reduction downstream. One-dimensional unsteady-flow hydraulic modelling was found to be required to reproduce these dynamics. Hydraulic modelling also proved to be potentially useful for the critical analysis and extrapolation of operational discharge rating curves.     

Experimental investigation of flow pattern and sediment deposition in rectangular shallow reservoirs

This paper reports the experimental investigation of flow pattern, preferential regions of deposition and trap efficiency as a function of the length of rectangular shallow reservoirs. Four flow patterns were identified (from longer to shorter reservoirs): an asymmetric flow with two reattachment points, an asymmetric flow with one reattachment point, an unstable flow, and a symmetric flow without any reattachment point. Using dye visualizations, the median value and the temporal variability of the reattachment lengths were precisely measured for the asymmetric flows. For each stable flow, sediment tests with plastic particles were carried out. The regions of deposition on the bed of the reservoir were clearly a function of the flow pattern. The transition from an asymmetric flow pattern to a symmetric flow pattern was responsible for an abrupt decrease of the trap efficiency; a number of regression laws were discussed to take it into account.     

Tachylyte in Cenozoic basaltic lavas from the Czech Republic and Iceland: contrasting compositional trends

Mineralogy and Petrology - Tachylytes from rift-related volcanic rocks were recognized as: (i) irregular veinlets in host alkaline lava flows of the Kozákov volcano, Czech Republic, (ii)...     

Multiarchive paleoseismic record of late Pleistocene and Holocene strong earthquakes in Switzerland

A multiarchive approach has been applied to the investigation of the late Pleistocene and Holocene record of strong earthquakes in Switzerland. The geological archives used for this study include active faults, lake deposits, slope instabilities, and caves. In the Basle area, eight trenches were opened across the Basle–Reinach fault, nearby rockfall deposits were systematically investigated, sediment cores were taken from two lakes, and nine caves were studied. In Central Switzerland, five lakes were investigated by means of high-resolution seismic lines and sediment cores. Furthermore, three caves were studied in Central Switzerland. Altogether, the investigations are based on more than 350 km of high-resolution reflection seismic lines, 450 m of core samples, 260 m of trenches, and 245 radiocarbon age determinations. The measured co-seismic displacements along the Basle–Reinach fault supply independent information for the magnitude of the AD 1356 Basle earthquake exclusively based on geological evidence. Deformation features related to three well-documented strong historic earthquake shocks were identified. Deformation features of the AD 1774 Altdorf and AD 1601 Unterwalden earthquakes can be used to calibrate paleoseismic evidence in Central Switzerland. Altogether, traces of 13 earthquakes could be found in the two study areas, all of them with magnitudes M_w  6 or greater. For the first time, the earthquake catalogue for Switzerland can be extended back beyond historic records, into the late Pleistocene, spanning 15,000 years.     

On early Solar System chronology: Implications of an heterogeneous spatial distribution of Al and Mn

Early Solar System chronology is usually built with the assumption that the distribution of short-lived radionuclides was homogeneous through the solar accretion disk. At present, there is no unambiguous evidence for a homogeneous distribution of short-lived radionuclides in the solar accretion disk, while some data point to a heterogeneous distribution of short-lived radionuclides. In this paper, we explore a possible chronology based on a heterogeneous distribution of ²⁶Al and ⁵³Mn in the accretion disk. Our basic assumption is that the different abundances of extinct short-lived radionuclides in calcium-aluminium-rich inclusions (CAIs) and chondrules are due to spatial rather than temporal differences. We develop a simple model where CAIs and chondrules form contemporaneously, in different spatial locations, and are characterised by distinct initial ²⁶Al and ⁵³Mn abundances. In this model, all evolved bodies are supposed to be originally chondritic, i.e., to be made of a mixture of CAIs, chondrules, and matrix. This mixture determines the initial content in ²⁶Al and ⁵³Mn of a chondritic parent-body as a function of its CAI and chondrule abundance fraction. This approach enables us to calculate coherent ²⁶Al and ⁵³Mn ages from the agglomeration of the parent-body precursors (CAIs and chondrules) until the isotopic closure of ²⁶Al and ⁵³Mn, thereafter called ²⁶Al-⁵³Mn age. We calculate such ²⁶Al-⁵³Mn ages for a diversity of evolved objects, with the constraint that they should be found for realistic chondritic parent-body precursors, i.e., objects having similar or identical petrograpy to the existing chondrite groups. The so defined age of the d’Orbigny angrite is 4.3 ± 1.1 Myr, for the Asuka-881394 eucrite 2.8 ± 1.0 Myr, for the H4 chondrite Sainte Marguerite ∼3 Myr, and for H4 Forest Vale ∼5 Myr. The calculated ²⁶Al-⁵³Mn ages give timescales for the evolution of the respective parent-bodies/meteorites that can be investigated in the light of further petrographic studies. We anchor the calculated relative chronology to an absolute chronology using absolute Pb-Pb ages and relative Hf-W ages of the objects under scrutiny. The precursors of Sainte Marguerite and Forest Vale agglomerated at the same time (∼4565.8 ± 1.2 Ma ago). The precursors of eucrites (Asuka-881394) agglomerated 4564.8 ± 1.2 Ma ago. The precursors of angrites agglomerated late (4561.5 ± 1.8 Ma ago). Our model provides a fully compatible Al-Mg/Mn-Cr/Pb-Pb chronology, and is shown to be robust to reasonable changes in the input parameters. The calculated initial ²⁶Al/²⁷Al ratios are high enough to have ²⁶Al as a possible heat source for differentiation.