Modelling close encounters with Öpik''s theory

The results obtained by numerical integration of the equations of motion of fictitious comets, in the restricted circular three-dimensional three-body problem, are compared with those obtained with Öpik's theory of close encounters, for an experimental set-up similar to that used by Froeschlé and Rickman (1980, 1981) to model both the infeed of comets from the trans-jovian region into the Jupiter family and their subsequent orbital evolution within the family. The distributions of perturbations in orbital energy E, eccentricity e and inclination i are well reproduced by Öpik's theory, as long as the comparison is made on the outcomes of encounters only up to a certain unperturbed distance b_max; several values of the latter are experimented with and it is found that, surprisingly, Öpik's theory seems to be still working reasonably well for values of b_max in excess of several times the Hill's radius of the planet.     

3D segmentation of single trees exploiting full waveform LIDAR data

This paper highlights a novel segmentation approach for single trees from LIDAR data and compares the results acquired both from first/last pulse and full waveform data. In a first step, a conventional watershed-based segmentation procedure is set up, which robustly interpolates the canopy height model from the LIDAR data and identifies possible stem positions of the tallest trees in the segments calculated from the local maxima of the canopy height model. Secondly, this segmentation approach is combined with a special stem detection method. Stem positions in the segments of the watershed segmentation are detected by hierarchically clustering points below the crown base height and reconstructing the stems with a robust RANSAC-based estimation of the stem points. Finally, a new three-dimensional (3D) segmentation of single trees is implemented using normalized cut segmentation. This tackles the problem of segmenting small trees below the canopy height model. The key idea is to subdivide the tree area in a voxel space and to set up a bipartite graph which is formed by the voxels and similarity measures between the voxels. Normalized cut segmentation divides the graph hierarchically into segments which have a minimum similarity with each other and whose members (= voxels) have a maximum similarity. The solution is found by solving a corresponding generalized eigenvalue problem and an appropriate binarization of the solution vector. Experiments were conducted in the Bavarian Forest National Park with conventional first/last pulse data and full waveform LIDAR data. The first/last pulse data were collected in a flight with the Falcon II system from TopoSys in a leaf-on situation at a point density of 10 points/m². Full waveform data were captured with the Riegl LMS-Q560 scanner at a point density of 25 points/m² (leaf-off and leaf-on) and at a point density of 10 points/m² (leaf-on). The study results prove that the new 3D segmentation approach is capable of detecting small trees in the lower forest layer. So far, this has been practically impossible if tree segmentation techniques based on the canopy height model were applied to LIDAR data. Compared to a standard watershed segmentation procedure, the combination of the stem detection method and normalized cut segmentation leads to the best segmentation results and is superior in the best case by 12%. Moreover, the experiments show clearly that using full waveform data is superior to using first/last pulse data.     

Assimilation of carbonate country rock by the parent magma of the Panzhihua Fe-Ti-V deposit（SW China）：Evidence from stable isotopes

The Panzhihua intrusion in southwest China is part of the Emeishan Large Igneous Province and host of a large Fe-Ti-V ore deposit.During emplacement of the main intrusion,multiple generations of mafic dykes invaded carbonate wall rocks,producing a large contact aureole.We measured the oxygen-isotope composition of the intrusions,their constituent minerals,and samples of the country rock.Magnetite and plagioclase from Panzhihua intrusion haveδ¹⁸O values that are consistent with magmatic equilibrium, and formed from magmas withδ¹⁸O values that were 1-2‰higher than expected in a mantle-derived magma.The unmetamorphosed country rock has highδ¹⁸O values,ranging from 13.2‰（sandstone） to 24.6-28.6‰（dolomite）.The skarns and marbles from the aureole have lowerδ¹⁸O andδ¹³C values than their protolith suggesting interaction with fluids that were in exchange equilibrium with the adjacent mafic magmas and especially the numerous mafic dykes that intruded the aureole.This would explain the alteration ofδ¹⁸O of the dykes which have significantly higher values than expected for a mantle-derived magma.Depending on the exactδ¹⁸O values assumed for the magma and contaminant, the amount of assimilation required to produce the elevatedδ¹⁸O value of the Panzhihua intrusion was between 8 and 13.7 wt.%,assuming simple mixing.The exact mechanism of contamination is unclear but may involve a combination of assimilation of bulk country rock,mixing with a melt of the country rock and exchange with CO₂-rich fluid derived from decarbonation of the marls and dolomites.These mechanisms,particularly the latter,were probably involved in the formation of the Fe-Ti-V ores.     

Tectonometamorphic evolution of the Atbashi high‐P units (Kyrgyz CAOB,Tien Shan): Implications for the closure of the Turkestan Ocean and continental subduction–exhumation of the South Kazakh continental margin

The South Tien Shan (STS) belt results from the last collision event in the western Central Asian Orogenic Belt (CAOB). Understanding its formation is of prime importance in the general framework of the CAOB. The Atbashi Range preserves high‐P (HP) rocks along the STS suture, but still, its global metamorphic evolution remains poorly constrained. Several HP units have been identified: (a) a HP tectonic mélange including boudins of mafic eclogites in a sedimentary matrix, (b) a large (>100 km long) high‐P metasedimentary unit (HPMU) and (c) a lower blueschist facies accretionary prism. Raman Spectroscopy on carbonaceous material combined with phengite and chlorite multiequilibria and isochemical phase diagram modelling indicates that the HPMU recorded homogeneous P–T conditions of 23–25 kbar and 560–570°C along the whole unit. ⁴⁰Ar/³⁹Ar dating on phengite from the HPMU ranges between 328 and 319 Ma at regional scale. These ages are interpreted as (re‐) crystallization ages of phengite during T_max conditions at a pressure range of 20–25 kbar. Thermobarometry on samples from the HP tectonic mélange provides similar metamorphic peak conditions. Thermobarometry on the blueschist to lower greenschist facies accretionary prism indicates that it underwent P–T conditions of 5–6 kbar and 290–340°C, highlighting a 17–20 kbar pressure gap between the HPMU‐tectonic mélange units and the accretionary prism. Comparison with available geochronological data suggests a very short time span between the prograde path (340 Ma), HP metamorphic peak (330 Ma), the T_max (328–319 Ma) and the final exhumation of the HPMU (303–295 Ma). Extrusion of the HPMU, accommodated by a basal thrust and an upper detachment, was driven by buoyant forces from 70–75 km up to 60 km depth, which directly followed continental subduction and detachment of the HPMU. At crustal depths, extrusion was controlled by collisional tectonics up to shallow levels. Lithological homogeneity of the HPMU and its continental‐derived character from the North Tien Shan suggest this unit corresponds to the hyper‐extended continental margin of the Kazakh continent, subducted southward below the north continental active margin of the Tarim craton. Integration of the available geological data allows us to propose a general geodynamic scenario for Tien Shan during the Carboniferous with a combination of (a) N‐dipping subduction below the Kazakh margin of Middle Tien Shan until 390–340 Ma and (b) S‐dipping subduction of remaining Turkestan marginal basins between 340 and 320 Ma.     

Influence of shallow infiltration on time-lapse ERT: Experience of advanced interpretation

Previous time-lapse Electrical Resistivity Tomography (ERT) studies have experienced difficulties in reconstructing reliable calculated resistivity changes in the subsurface. Increases or decreases of resistivity appear in the calculated ERT image where no changes were noted in the subsurface, leading to erroneous hydrological interpretations of the geophysical results. In this article, we investigate how a variation of actual resistivity with time and at shallow depth can influence time-lapse ERT results and produce resistivity artefacts at depth. We use 1 and 2-D numerical modelling to simulate infiltration scenarios. Using a standard time-lapse inversion, we demonstrate the resistivity artefact production according to the electrode spacing parameter. We used an advanced inversion methodology with a decoupling line at shallow depth to attenuate or remove resistivity artefacts. We also applied this methodology to a field data set obtained in a semi-arid environment in Burkina Faso, West Africa. Here, time-lapse ERT shows several resistivity artefacts of calculated resistivity if a standard inversion is used. We demonstrate the importance of a dense sampling of shallow resistivity variations at shallow depth. Advanced interpretation allows us to significantly attenuate or remove the resistivity artefact production at intermediate depth and produce reliable interpretation of hydrological processes.     

Geology and geochronology of the Betara region in south-southeastern Brazil: Evidence for possible Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extension events

The results of geological mapping, chemical analysis and radiometric dating of metabasic rocks of Betara Formation, and mapping and dating of those present in the Betara basement nucleus together with mylonitic granodiorite and syenogranite are reported here. U–Pb analysis of bulk zircon fractions from the metabasic rocks of the basement nucleus yielded a Statherian age of 1790 ± 22 Ma, while the metabasic rocks from the upper part of the Betara Formation yielded a Calymmian age between 1500 and 1450 Ma. This age is a minimum for the deposition of the Betara Formation. The older metabasic rocks are associated with post-tectonic, possibly anorogenic syenogranite, while the younger ones are gabbro or very porphyritic ankaramite whose REE patterns are consistent with crystallization from an N-MORB parent magma. The observations and data point to the probable events associated with extensional processes of the end of Paleoproterozoic and early Mesoproterozoic. Similar registers of Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extensional events are recorded in other parts of the South American and African continents. The Neoproterozoic witnessed the formation and junction of the tectonic slices which formed the Apiaí domain during the assemblage of western Gondwana.     

Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure

We consider an immiscible incompressible two-phase flow in a porous medium composed of two different rocks so that the capillary pressure field is discontinuous at the interface between the rocks. This leads us to apply a concept of multivalued phase pressures and a notion of weak solution for the flow which have been introduced in Cancès and Pierre (SIAM J Math Anal 44(2):966–992, 2012). We discretize the problem by means of a numerical algorithm which reduces to a standard finite volume scheme in each rock and prove the convergence of the approximate solution to a weak solution of the two-phase flow problem. The numerical experiments show in particular that this scheme permits to reproduce the oil-trapping phenomenon.     

Evolution of a Holocene banner bank controlled by morphodynamics and structural setting of a macrotidal coast: Saint-Brieuc Bay(NW-Europe)

The morphology and internal structure of the Horaine Bank(Bay of Saint-Brieuc, NW France) are described based on multibeam echosounder and high-resolution seismic datasets coupled with vibro-core data. The Horaine Bank shows large-scale bedforms in the lee of a submerged rocky shoal, which allowed defining it as a Banner Bank.The internal structure of the sandbank reveals four seismic units(U1–U4) on a Cambrian basement(U0). The basal unit U1 is interpreted as reworked lowstand fluvial sediments those infilled micro incised valleys during a rise in sea level. This unit is overlain by paleo-coastal barrier sand-spit(U2) whose development was controlled by swell in the context of a rapid rise in sea level. The successive prograding unit(U3) is interpreted as flooding deposits in continuity with unit U2. The unit U4 is characterized by oblique reflectors oriented in two opposite directions. This last unit, dated post 3500 yr BP, corresponds to migrating dunes superimposed on the bank and observable in the high-resolution bathymetric data. The strong correlation between tidal currents and the apparent clockwise migration of dune crests suggests the presence of a tidal gyre controlling the present-day dynamics of most of the Horaine bank dunes. This study proposes a new model for the construction of banner banks characterized by the gradual transition of a sand spit to a banner bank during marine transgression and ensuing hydrodynamic variability.     

Stress Measurement by Overcoring at Shallow Depths in a Rock Slope: the Scattering of Input Data and Results

This paper describes a field experiment of stress measurement using the overcoring method performed in a rock slope, called Rochers de Valabres (located in France’s Southern Alps Region), a field laboratory site prone to rockfalls. Six measurements were conducted at shallow depths from the surface, moving deeper along a sub-horizontal borehole. The experiment was conducted in heterogeneous and anisotropic gneiss, with the overcored rock elastic properties, as evaluated by biaxial and uniaxial tests, being widely scattered. Since stress calculations are sensitive to all input data uncertainties, strain inversion was, thus, performed using an experimental device and Monte Carlo simulations. The experimental device allows the assessment of rather broad confidence intervals for both stress magnitude and orientation. The results indicate that the stress state in the surface area is quite heterogeneous and may be correlated with topography. The measurements show a nonlinear stress distribution with distance to the free surface, along with high values of principal stresses, despite the vicinity of the surface. Although influenced by local topography, orientations of the principal computed stresses are characterized by a high turnover due to local heterogeneities. The results are roughly in accordance to a 2D finite element model of the site.