Seafloor change detection using multibeam echosounder backscatter: case study on the Belgian part of the North Sea

To characterize seafloor substrate type, seabed mapping and particularly multibeam echosounding are increasingly used. Yet, the utilisation of repetitive MBES-borne backscatter surveys to monitor the environmental status of the seafloor remains limited. Often methodological frameworks are missing, and should comprise of a suite of change detection procedures, similarly to those developed in the terrestrial sciences. In this study, pre-, ensemble and post-classification approaches were tested on an eight km² study site within a Habitat Directive Area in the Belgian part of the North Sea. In this area, gravel beds with epifaunal assemblages were observed. Flourishing of the fauna is constrained by overtopping with sand or increased turbidity levels, which could result from anthropogenic activities. Monitoring of the gravel to sand ratio was hence put forward as an indicator of good environmental status. Seven acoustic surveys were undertaken from 2004 to 2015. The methods allowed quantifying temporal trends and patterns of change of the main substrate classes identified in the study area; namely fine to medium homogenous sand, medium sand with bioclastic detritus and medium to coarse sand with gravel. Results indicated that by considering the entire study area and the entire time series, the gravel to sand ratio fluctuated, but was overall stable. Nonetheless, when only the biodiversity hotspots were considered, net losses and a gradual trend, indicative of potential smothering, was captured by ensemble and post-classification approaches respectively. Additionally, a two-dimensional morphological analysis, based on the bathymetric data, suggested a loss of profile complexity from 2004 to 2015. Causal relationships with natural and anthropogenic stressors are yet to be established. The methodologies presented and discussed are repeatable and can be applied to broad-scale geographical extents given that broad-scale time series datasets become available.     

The Creston,California, meteorite fall and the origin of L chondrites

It has been proposed that all L chondrites resulted from an ongoing collisional cascade of fragments that originated from the formation of the ~500 Ma old asteroid family Gefion, located near the 5:2 mean‐motion resonance with Jupiter in the middle Main Belt. If so, L chondrite pre‐atmospheric orbits should be distributed as expected for that source region. Here, we present contradictory results from the orbit and collisional history of the October 24, 2015, L6 ordinary chondrite fall at Creston, CA (here reclassified to L5/6). Creston's short 1.30 ± 0.02 AU semimajor axis orbit would imply a long dynamical evolution if it originated from the middle Main Belt. Indeed, Creston has a high cosmic ray exposure age of 40–50 Ma. However, Creston's small meteoroid size and low 4.23 ± 0.07° inclination indicate a short dynamical lifetime against collisions. This suggests, instead, that Creston originated most likely in the inner asteroid belt and was delivered via the ν₆ resonance. The U‐Pb systematics of Creston apatite reveals a Pb‐Pb age of 4,497.1 ± 3.7 Ma, and an upper intercept U‐Pb age of 4,496.7 ± 5.8 Ma (2σ), circa 70 Ma after formation of CAI, as found for other L chondrites. The K‐Ar (age ~4.3 Ga) and U,Th‐He (age ~1 Ga) chronometers were not reset at ~500 Ma, while the lower intercept U‐Pb age is poorly defined as 770 ± 320 Ma. So far, the three known L chondrites that impacted on orbits with semimajor axes a <2.0 AU all have high (>3 Ga) K‐Ar ages. This argues for a source of some of our L chondrites in the inner Main Belt. Not all L chondrites originate in a continuous population of Gefion family debris stretching across the 3:1 mean‐motion resonance.     

Spatial distribution of morphometric parameters of glacial cirques in the Central Pyrenees(Aran and Boí valleys)

Glacial cirques are typical landscape features of mid-latitude mountain environments like the Central Pyrenees. Their morphology as well as their spatial distribution provides insights about past glaciers and climates. In this study, we examine the distribution, morphometrical and topographical characteristics of glacial cirques in two U-shaped glacial valleys located in the Central Pyrenees – the Aran and the Boí valleys. They are located in different aspects of this mountain range(north vs south) under different climatic influences that promoted distinct glaciation patterns during the late Pleistocene. The spatial mapping of these landforms was carried out using high-resolution imagery and field observations. We analysed the data of the morphometrical and topographical variables of the glacial cirques by using different statistical and geospatial methods in order tounveil the factors controlling their formation and development. A total of 186 glacial cirques were mapped in the study area, including 119 in the Aran and 67 in the Boí valleys. The local topography and microclimate conditions lead to substantial differences in both areas in terms of the morphology and dimensions of the cirques. Glacial cirques in Boí are distributed at slightly higher elevations than in Aran and they are also larger, though their dimensions decrease with elevation in both valleys. Aran cirques are mostly oriented NE, while Boí landforms do not show any prevailing aspect. Even though lithology does not control the distribution of the glacial cirques, some specific lithological settings may favour the development of larger cirques. In general, glacial cirques in the Aran and the Boí valleys show morphometrical properties similar to those reported in other mid-latitude mountain ranges.     

The location of infinite electrodes in pole–pole electrical surveys: consequences for 2D imaging

In 2D-multielectrode electrical surveys using the pole–pole array, the distance to ‘infinite electrodes' is actually finite. As a matter of fact, the available cable length generally imposes a poor approximation of theoretical location of these electrodes at infinity. This study shows that in most of the cases, the resulting apparent resistivity pseudosection is strongly distorted. Numerical simulation validated by field test also shows that a particular finite array provides results that are as close as possible to the ones of the ideal pole–pole array. This is achieved when two conditions that are weaker than an infinite location are fulfilled: (i) the ‘infinite electrodes' are placed symmetrically on both sides of the in-line electrodes with a spread angle of 30° and (ii) the length of ‘infinite lines' is at least 20 times the greatest distance between in-line electrodes. The electrical 2D image obtained with this enhanced array is the least distorted one with respect to the pole–pole image. The apparent resistivities are generally underestimated, but this deviation is almost homogeneous. Though the shift cannot be determined a priori, the interpretation of such an image with direct or inverse software designed for pole–pole data provides an accurate interpretation of the ground geometry.     

Hydrogeomorphological mapping in assessing ground water by using remote sensing data — A case study in lehra gaga block, sangrur district, Punjab

The paper presents the results of hydrogeomorphological mapping using IRS-IB LISS II data and evaluation of ground water prospects of each hydrogeomorphological unit in the Lehra Gaga block of Sangrur district, Punjab. The major geomorphic units identified in the area are, alluvial plain, sand dunes, palaeo channels and the Ghagar flood plain. The study area being part of alluvial plain has good to excellent ground water prospects. Field observations showed that ground water occurs under both confined and unconfined conditions with water table at shallow depth. The area on either side of the Ghagar river and along the major canals (about 46% of the total geographical area in the block) have good quality of ground water and is suitable for irrigation, whereas the water quality is marginal (sodic) in 52 per cent area of the block..     

The composition and thickness of the crust of Mars estimated from rare earth elements and neodymium-isotopic compositions of Martian meteorites

Abstract— Isotopic and trace element compositions of Martian meteorites show that early differentiation of Mars produced complementary crustal and mantle reservoirs that were sampled by later magmatic events. This paper describes a mass balance model that estimates the rare earth element (REE) content and thickness of the crust of Mars from the compositions of shergottites. The diverse REE and Nd isotopic compositions of shergottites are most easily explained by variable addition of light rare earth element (LREE)–enriched crust to basaltic magmas derived from LREE-depleted mantle source regions. Antarctic shergottites EET 79001, ALH 77005, LEW 88516, and QUE 94201 all have strongly LREE-depleted patterns and positive initial ¹⁴³Nd isotopic compositions, which is consistent with the generation of these magmas from depleted mantle sources and little or no interaction with enriched crust. In contrast, Shergotty and Zagami have negative initial ¹⁴³Nd isotopic compositions and less pronounced depletions of the LREE, which have been explained by incorporation of enriched crustal components into mantle-derived magmas (Jones, 1989; Longhi, 1991; Borg et al., 1997). The mass balance model presented here derives the REE composition of the crustal component in Shergotty by assuming it represents a mixture between a mantle-derived magma similar in composition to EET 79001A and a LREE-enriched crustal component. The amount of crust in Shergotty is constrained by mixing relations based on Nd-isotopic compositions, which allows the REE pattern of the crustal component to be calculated by mass balance. The effectiveness of this model is demonstrated by the successful recovery of important characteristics of the Earth's continental crust from terrestrial Columbia River basalts. Self-consistent results for Nd-isotopic compositions and REE abundances are obtained if Shergotty contains ～10–30% of LREE-enriched crust with >10 ppm Nd. This crustal component would have moderately enriched LREE (Sm/Nd = 0.25–0.27; ¹⁴⁷Sm/¹⁴⁴Nd = 0.15–0.17; La/Yb = 2.7–3.8), relatively unfractionated heavy rare earth elements (HREE), and no Eu anomaly. Crust with these characteristics can be produced from a primitive lherzolitic Martian mantle by modest amounts (2–8%) of partial melting, and it would have a globally averaged thickness of <45 km, which is consistent with geophysical estimates. Mars may serve as a laboratory to investigate planetary differentiation by extraction of a primary basaltic crust.     

Meta-igneous (non-gneissic) tonalites and quartz-diorites from an extensive ca. 3800 Ma terrain south of the Isua supracrustal belt, southern West Greenland: constraints on early crust formation

In the Itsaq Gneiss Complex south of the Isua supracrustal belt (West Greenland) some areas of early Archaean tonalite and quartz-diorite are non-gneissic, free of pegmatite veins, and in rarer cases are undeformed with relict igneous textures and hence were little modified by heterogeneous ductile deformation under amphibolite facies conditions in several Archaean events. Such well-preserved early Archaean rocks are extremely rare. Tonalites are high Al, and have bulk compositions close to experimental liquids. Trace element abundances and modelling suggest that they probably originated as melts derived from basaltic compositions at sufficiently high pressures to require residual garnet + amphibolites ± clinopyroxene in the source. The major element characteristics of the quartz-diorites suggest these were derived from more mafic magmas than the tonalites, and underwent either igneous differentiation or mixing with crustal material. As in modern arc magmas, high relative abundances of Sr, Ba, Pb, and alkali elements cannot be generated simply from a basaltic source formed by large degrees of melting of a depleted mantle. This may indicate an important role for fluids interacting with mafic rocks in generating the earliest preserved continental crust. The high Ba/Th, Ba/Nb, La/Nb and low Nb/Th, Ce/Pb, and Rb/Cs ratios of these tonalites are also observed in modern arc magmas. SHRIMP U-Pb zircon geochronology was undertaken on seven tonalites, one quartz-diorite, a thin pegmatitic vein and a thin diorite dyke. Cathodoluminescence images show the zircon populations of the quartz-diorite and tonalites are dominated by single-component oscillatory-zoned prismatic grains, which gave ages of 3806 ± 5 to 3818 ± 8 Ma (2σ) (quartz-diorite and 5 tonalites) and 3795 ± 3 Ma (1 tonalite). Dating of recrystallised domains cutting oscillatory-zoned zircon indicates disturbance as early as 3800–3780 Ma. There are rare ca. 3600 Ma and 3800–3780 Ma (very high U and low Th/U) ≤ 20 μm wide partial overgrowths on the prismatic grains. Given likely Zr-undersaturation of precursor melts and evidence of zircon recrystallisation and metamorphic regrowth as early as 3800–3780 Ma, the age determinations on the prismatic oscillatory-zoned zircon populations give the igneous crystallisation age of the tonalite and quartz-diorite protoliths. When the coherency of the geochemistry is considered, these samples represent the best preserved suites of ca. 3800 Ma felsic igneous rocks yet documented. Received: 1 December 1998 / Accepted: 23 July 1999     

Statistics for Modeling Heavy Tailed Distributions in Geology: Part II. Applications

In this paper, we present three diverse types of applications of extreme value statistics in geology, namely: earthquakes magnitudes, diamond values, and impact crater size distribution on terrestrial planets. Each of these applications has a different perspective toward tail modeling, yet many of these phenomena exhibit heavy or long tails which can be modeled by power laws. It is shown that the estimation of important tail characteristics, such as the extreme value index, is directly linked to the interpretation of the underlying geological process. Only the most extreme data are useful for studying such phenomena, so thresholds must be selected above which the data become power laws. In the case of earthquake magnitudes, we investigate the use of extreme value statistics in predicting large events on the global scale and for shallow intracontinental earthquakes in Asia. Large differences are found between estimates obtained from extreme value statistics and the usually applied standard statistical techniques. In the case of diamond deposits, we investigate the impact of the most precious stones in the global valuation of primary deposits. It is shown that in the case of Pareto-type behavior, the expected value of few extreme stones in the entire deposit has considerable influence on the global valuation. In the case of impact crater distributions, we study the difference between craters distributions on Earth and Mars and distributions occurring on other planets or satellites within the solar system. A striking result is that all planets display the same distributional tail except for Earth and Mars. In a concluding account, we demonstrate the apparent loghyperbolic variation in all of the above-mentioned examples.     

The Correlation Structure of Matheron's Classical Variogram Estimator Under Elliptically Contoured Distributions

The classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. When data have an elliptically contoured distribution with constant mean, the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix, the covariance matrix, and the kurtosis. Several specific cases are studied closely. A subclass of elliptically contoured distributions with a particular family of covariance matrices is shown to possess exactly the same correlation structure for the classical variogram estimator as the multivariate independent Gaussian distribution. The consequences on variogram fitting by generalized least squares are discussed.     

Park Forest (L5) and the asteroidal source of shocked L chondrites

The Park Forest (L5) meteorite fell in a suburb of Chicago, Illinois (USA) on March 26, 2003. It is one of the currently 25 meteorites for which photographic documentation of the fireball enabled the reconstruction of the meteoroid orbit. The combination of orbits with pre‐atmospheric sizes, cosmic‐ray exposure (CRE), and radiogenic gas retention ages (“cosmic histories”) is significant because they can be used to constrain the meteoroid's “birth region,” and test models of meteoroid delivery. Using He, Ne, Ar, ¹⁰Be, and ²⁶Al, as well as a dynamical model, we show that the Park Forest meteoroid had a pre‐atmospheric size close to 180 g cm^?2, 0–40% porosity, and a pre‐atmospheric mass range of ~2–6 tons. It has a CRE age of 14 ± 2 Ma, and (U, Th)‐He and K‐Ar ages of 430 ± 90 and 490 ± 70 Ma, respectively. Of the meteorites with photographic orbits, Park Forest is the second (after Novato) that was shocked during the L chondrite parent body (LCPB) break‐up event approximately 470 Ma ago. The suggested association of this event with the formation of the Gefion family of asteroids has recently been challenged and we suggest the Ino family as a potential alternative source for the shocked L chondrites. The location of the LCPB break‐up event close to the 5:2 resonance also allows us to put some constraints on the possible orbital migration paths of the Park Forest meteoroid.