Understanding fracture distribution and its relation to knickpoint evolution in the Rio Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) using landscape-scale hydrogeophysics

The Rio Icacos watershed in the Luquillo Mountains (Puerto Rico) is unique due to its extremely rapid weathering rates. The watershed is incised into a quartz diorite that has developed a large knickzone defining the river profile. Regolith thickness within the watershed generally decreases from 20 to 30 m at the ridges to several meters in the quartz diorite-dominated valley to tens of centimeters near the major river knickpoint, as determined from previous studies. Above the knickzone, we observe spheroidal corestones, but below this weathering is much less apparent. Measured erosion rates from previous studies are also high in the knickzone compared with upper elevations within the river profile. A suite of near-surface geophysical methods (i.e. ground penetrating radar and terrain conductivity) capable of fast data acquisition in rugged landscapes, was deployed at kilometer scales to characterize critical zone structure. Concentrations of chaotic ground penetrating radar (GPR) reflections and diffraction hyperbolas with low electrical conductivity were observed in vertical zones that outcrop at the land surface as areas of intense fracturing and spheroidally weathered corestones. The width of these fractured and weathered zones showed an increase with proximity to the knickpoint, and was attributed to dilation of these sub-vertical fractures near the knickpoint, as postulated theoretically by a stress model calculated for the topographic variability across the knickzone in the Rio Icacos, and that shows a release of compressive stress near the knickpoint. We hypothesize that erosion rates increase in the knickzone because of this inferred dilation of fractures. Specifically, opened fractures could enhance access of water and in turn promote spalling, erosion, and spheroidal weathering. This study shows that ground-based hydrogeophysical methods used at the landscape-scale (traditionally applied at smaller scales) can be used to explore critical zone architecture at the scales needed to explain the extreme variability in erosion rates across river profiles. © 2018 John Wiley & Sons, Ltd.     

Plurigaussian modeling of geological domains based on the truncation of non-stationary Gaussian random fields

The plurigaussian model is used in mining engineering, oil reservoir characterization, hydrology and environmental sciences to simulate the layout of geological domains in the subsurface, while reproducing their spatial continuity and dependence relationships. However, this model is well-established only in the stationary case, when the spatial distribution of the domains is homogeneous in space, and suffers from theoretical and practical impediments in the non-stationary case. To overcome these limitations, this paper proposes extending the model to the truncation of intrinsic random fields of order k with Gaussian generalized increments, which allows reproducing spatial trends in the distribution of the geological domains. Methodological tools and algorithms are presented to infer the model parameters and to construct realizations of the geological domains conditioned to existing data. The proposal is illustrated with the simulation of rock type domains in an ore deposit in order to demonstrate its applicability. Despite the limited number of conditioning data, the results show a remarkable agreement between the simulated domains and the lithological model interpreted by geologists, while the conventional stationary plurigaussian model turns out to be unsuccessful.     

Benthic indicators and index approaches in the three main estuaries along the French Atlantic coast (Seine, Loire and Gironde)

The three main estuaries in the French Atlantic coast – the Seine, Loire and Gironde, all with high tidal regimes – are interfaces between the continental and the coastal ecosystems. The Seine and Gironde are highly contaminated, whereas the Loire remains in a more natural state. Both the Seine and Gironde have suffered from harbour construction, and as a result, their biological units are extremely compartmentalized. Benthic species and communities have adapted to tolerate temporal physical and chemical changes ( e.g. salinity, substrata, depth, and levels of fine particles and oxygen) and human activities ( e.g. dredging, shipping traffic, and habitat reduction). Although numerous bio-indicators and indices are used to define the ecological quality status (EcoQS) of coastal waters, very few of them were developed specifically for environments with a mosaic of conditions and salinity levels, including freshwater. The main problem appears to be that all the indices for determining anthropogenic stress examine the abundances of stress-tolerant species, which may also be able to tolerate natural stressors such as those occurring in estuaries. This paper takes a look at the development status of the benthic indicators and index approaches used in the three main estuaries along the French Atlantic coast. In addition, it examines the adaptation of the different benthic indicators to the taxonomic sufficiency principle, and the adaptation of the Benthic Opportunistic Polychaetes Amphipods index (BOPA) and the Benthic Opportunistic Annelida Amphipod index (BO2A) for use in the freshwater zones of transitional waters ( i.e. up to the upper limit of the tidal range). Several perspectives are discussed in light of the diverse disturbances and the heterogeneity of such stressed zones, such as the use of multimetric and scoring approaches.     

Geochemical Characteristics of Early Permian Pyroclastic Rocks in the Jimunai Basin,West Junggar,Xinjiang(NW China): Implications for Provenance and Tectonic Setting

During the Late Paleozoic, West Junggar(Xinjiang, NW China) experienced a shift in tectonic setting from compression to extension. Ha'erjiao is an important area for investigating collisional structures, post–collisional structures, and magmatic activities. Based on the petrological and geochemical characteristics of pyroclastic and other volcanic rocks in the Permian Kalagang Formation from the borehole ZKH1205 in the Jimunai Basin, the main types of source rock for the pyroclastic rocks deposited in the basin are identified and their implications for the Early Permian tectonic setting examined. The abundance of basalt and andesite lithic fragments in the pyroclastic rocks, together with the REE characteristics and the contents of transition and high field strength elements show that the source rocks were chiefly intermediate–basic volcanic rocks. High ICV values, low CIA values, low Rb/Sr ratios, low Th/U ratios and the mineralogical features suggest weak chemical weathering of the source rocks; the geochemical patterns of the pyroclastic rocks might not only have been impacted by crustal contamination but also might be related to the nature of the magma from the source area. The geochemical properties of the pyroclastic rocks distinguish them from arc-related ones, and such samples plot in the within-plate basalt(WPB) field in some diagrams. This is consistent with the formation background of the Early Permian volcanic rocks in this region.     

Spatial Analysis of Airborne Geophysical Data Applied to Geological Mapping and Mineral Prospecting in the Serra Leste Region, Carajás Mineral Province, Brazil

The Serra Leste region, located in the northeast portion of Carajás Mineral Province (PA), hosts important chrome-PGE, Au–Pd–Pt, Cu–Mo–Au, iron, and manganese deposits and occurrences related to different geological settings. We present a regional prospective model with focus on PGE mineralizations based on airborne geophysics (magnetometry, gamma-ray spectrometry and time domain electromagnetometry) and geological data. Geological and geophysics data analysis indicate a plausible correspondence between areas with meta-mafic–ultramafic rocks and low values in gamma spectrometry, mainly for the thorium channel. The meta-ultramafic units of the complexes, the main host rock of the PGE mineralization, are well defined by analytic signal from the anomalous magnetic field and by the apparent conductance anomalies derived from the GEOTEM system. The combination of the derivative products with updated geological data of the Serra Leste region and PGE-Luanga deposit enable to establish a regional prospective model for PGE and the identification of new targets within the studied area.     

Conditional Simulation of Nongaussian Random Functions

This paper presents a conditional simulation procedure that overcomes the limits of gaussian models and enables one to simulate regionalized variables with highly asymmetrical histograms or with partial or total connectivity of extreme values. The philosophy of the method is similar to that of sequential indicator technique, but it is more accurate because it is based on a complete bivariate model by means of an isofactorial law. The resulting simulations, which can be continuous or categorical, not only honor measured values at data points, but also reproduce the mono and bivariate laws of the random function associated to the regionalized variable, that is, every one or two-point statistic: histogram, variogram, indicator variograms. The sequential isofactorial method can also be adapted to conditional simulation of block values, without resorting to point–support simulations.     

Interaction of multiple ice streams on the Malin Shelf during deglaciation of the last British–Irish Ice Sheet

The Malin Shelf, off north-west Ireland, was an important zone of confluence for marine-based ice streams of the former British–Irish Ice Sheet (BIIS). Legacy geophysical datasets are used to construct models of the seismic character, relative age and distribution of shelf sediments and landforms. Buried and surface landform assemblages provide evidence that during deglaciation of the Late Devensian BIIS, the region was occupied not by a single Hebrides Ice Stream as previously proposed, but by four discrete ice streams, here referred to as the Sea of the Hebrides (SHIS), Inner Hebrides, North Channel and Tory Island ice streams. Our observations of stratigraphic relationships between the deposits of these ice streams indicate physical interactions between them during shelf deglaciation. We interpret an initial dominant cross-shelf flow along the SHIS impeding cross-shelf ice flow from other ice sheet sectors. Following withdrawal of the SHIS grounding line from the shelf edge to mid-shelf bathymetric highs during deglaciation, a reconfiguration of ice sheet flow paths allowed the expansion of smaller cross-shelf ice streams draining central Scotland and north-western Ireland. This internal dynamic behaviour provides a possible physical analogue for time-transgressive flow patterns reported for outlets draining the West Antarctic Ice Sheet.     

Two Tales of a City: detecting socio-economic disadvantage in an ‘advantaged’ Australian urban centre

Understanding the distribution of socio-economic status (SES) within and between communities of varying spatial scales is vital for equitable and effective resource delivery and policy development. In Australia, the Socio-economic Indexes for Areas (SEIFA) have been used extensively for these purposes, in both research and government service delivery. As area-based averages, however, their use in these areas is potentially confounded by the ecological fallacy and the modifiable areal unit problem (MAUP), statistical issues which have widespread implications for the use of spatial data. The Socio-economic Index for Individuals (SEIFI) is an individual-based index developed by the ABS to quantify SES at an individual level for 15–64 year-old Australians. Previous work investigated the potential disparities between these area- and individual-level indices, showing that considerable disparities exist, particularly in the Australian Capital Territory (ACT). Using spatial autocorrelation and pair-wise proportion analyses, we show that there was considerable diversity in individual-level SES at a very fine level of geography in 2006, a unique pattern in Australia. We also show that the SES ranking of areas within the ACT is, in many instances, not indicative of the proportions of individuals living within that area experiencing high levels of relative socio-economic disadvantage. We show that between 65–95 per cent of the ACT's disadvantaged population are masked within areas ranked as advantaged, depending on the index and level of spatial aggregation. Our results suggest that the use of SEIFA as a proxy for individual level disadvantage is highly problematic at both a theoretical and practical level, and necessitates further development of an individual- or household-level index as part of the regular suite of socio-economic indices produced by the ABS.     

Global sea-level rise and its relation to the terrestrial reference frame

We examined the sensitivity of estimates of global sea-level rise obtained from GPS-corrected long term tide gauge records to uncertainties in the International Terrestrial Reference Frame (ITRF) realization. A useful transfer function was established, linking potential errors in the reference frame datum (origin and scale) to resulting errors in the estimate of global sea level rise. Contrary to scale errors that are propagated by a factor of 100%, the impact of errors in the origin depends on the network geometry. The geometry of the network analyzed here resulted in an error propagation factor of 50% for the Z component of the origin, mainly due to the asymmetry in the distribution of the stations between hemispheres. This factor decreased from 50% to less than 10% as the geometry of the network improved using realistic potential stations that did not yet meet the selection criteria (e.g., record length, data availability). Conversely, we explored new constraints on the reference frame by considering forward calculations involving tide gauge records. A reference frame could be found in which the scatter of the regional sea-level rates was limited. The resulting reference frame drifted by 1.36 ± 0.22? mm/year from the ITRF2000 origin in the Z component and by ?0.44 ± 0.22?mm/year from the ITRF2005 origin. A bound on the rate of global sea level rise of 1.2 to 1.6?mm/year was derived for the past century, depending on the origin of the adopted reference frame. The upper bound is slightly lower than previous estimates of 1.8?mm/year discussed in the IPCC fourth report.