An estimate of the influence of loading effects on tectonic velocities in the Pyrenees

Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.     

Tracking emission sources of sulfur and elemental carbon in Hong Kong/Pearl River Delta region

The Tagged Species Source Apportionment (TSSA) algorithm is applied to study contributions to pollutants PM sulfate, SO2 and elemental carbon (EC) in Hong Kong by emitting sectors as well as non-emitting sources within and beyond the Pearl River Delta (PRD) region. We look at three areas of Hong Kong: western new towns, central downtown, and eastern countryside. Pollutant transport from beyond the PRD influenced all three areas during January and October 2004 but the source sectors impacting the three areas were different. Specifically, power plant SO2 from Hong Kong and Shenzhen, and vehicle EC from Shenzhen contribute to 11?~?66?% of SO2 concentration, and 33?~?75?% of EC concentration in the western new towns, respectively. Ships in and around Hong Kong contribute 8?~?56?% to the sulfate concentration in the downtown area, while local moving vehicles take up 28?~?81?% of the EC concentration there. The study also shows that diurnal variation of planetary boundary layer (PBL) causes day-night difference in SO2 by as much as 50?%. In addition, 13?~?38?% of all SO2 emissions have been converted to PM sulfate.     

Petrological evidence for stepwise accretion of metamorphic soles during subduction infancy (Semail ophiolite,Oman and UAE)

Metamorphic soles are tectonic slices welded beneath most large‐scale ophiolites. These slivers of oceanic crust metamorphosed up to granulite facies conditions are interpreted as forming during the first million years of intraoceanic subduction following heat transfer from the incipient mantle wedge towards the top of the subducting plate. This study reappraises the formation of metamorphic soles through detailed field and petrological work on three key sections from the Semail ophiolite (Oman and United Arab Emirates). Based on thermobarometry and thermodynamic modelling, it is shown that metamorphic soles do not record a continuous temperature gradient, as expected from simple heating by the upper plate or by shear heating as proposed in previous studies. The upper, high‐T metamorphic sole is subdivided in at least two units, testifying to the stepwise formation, detachment and accretion of successive slices from the down‐going slab to the mylonitic base of the ophiolite. Estimated peak pressure–temperature conditions through the metamorphic sole, from top to bottom, are 850°C and 1 GPa, 725°C and 0.8 GPa and 530°C and 0.5 GPa. These estimates appear constant within each unit but differing between units by 100–200°C and ~0.2 GPa. Despite being separated by hundreds of kilometres below the Semail ophiolite and having contrasting locations with respect to the ridge axis position, metamorphic soles show no evidence for significant petrological variations along strike. These constraints allow us to refine the tectonic–petrological model for the genesis of metamorphic soles, formed via the stepwise stacking of several homogeneous slivers of oceanic crust and its sedimentary cover. Metamorphic soles result not so much from downward heat transfer (ironing effect) as from progressive metamorphism during strain localization and cooling of the plate interface. The successive thrusts originate from rheological contrasts between the sole, initially the top of the subducting slab, and the peridotite above as the plate interface progressively cools. These findings have implications for the thickness, the scale and the coupling state at the plate interface during the early history of subduction/obduction systems.     

Quantification of Storm-Induced Bathymetric Change in a Back-Barrier Estuary

Geomorphology is a fundamental control on ecological and economic function of estuaries. However, relative to open coasts, there has been little quantification of storm-induced bathymetric change in back-barrier estuaries. Vessel-based and airborne bathymetric mapping can cover large areas quickly, but change detection is difficult because measurement errors can be larger than the actual changes over the storm timescale. We quantified storm-induced bathymetric changes at several locations in Chincoteague Bay, Maryland/Virginia, over the August 2014 to July 2015 period using fixed, downward-looking altimeters and numerical modeling. At sand-dominated shoal sites, measurements showed storm-induced changes on the order of 5 cm, with variability related to stress magnitude and wind direction. Numerical modeling indicates that the predominantly northeasterly wind direction in the fall and winter promotes southwest-directed sediment transport, causing erosion of the northern face of sandy shoals; southwesterly winds in the spring and summer lead to the opposite trend. Our results suggest that storm-induced estuarine bathymetric change magnitudes are often smaller than those detectable with methods such as LiDAR. More precise fixed-sensor methods have the ability to elucidate the geomorphic processes responsible for modulating estuarine bathymetry on the event and seasonal timescale, but are limited spatially. Numerical modeling enables interpretation of broad-scale geomorphic processes and can be used to infer the long-term trajectory of estuarine bathymetric change due to episodic events, when informed by fixed-sensor methods.     

Precursors and fluid flows in the case of the 1996, ML = 5.2 Saint-Paul-de-Fenouillet earthquake (Pyrenees,France): A complete pre-, co- and post-seismic scenario

Earthquake precursors are now regularly described but often detected only after a major or moderate seismic event. Presence and influence of fluids in the seismogenic processes are often observed at the time of earthquake studies. Even today, the understanding of the physical processes involved in the source region is a real challenge for seismic hazard assessment. Here, the aftershock sequence of the M_L = 5.2, 1996 Saint-Paul-de-Fenouillet (Eastern Pyrenees, France) earthquake is first re-examined with P-wave cross-correlations, resulting in extracting three multiplets and in determining new locations. Multiplets and spatio-temporal distribution analysis of the aftershocks allow for quantifying the hydraulic diffusivity D at a maximum value of 5 m²/s and the permeability K at 10^? 15 m² in the upper Pyrenean crust. Second, a model is established in order to explain the hydrogeochemical transient anomalies, which occurred during the 15 day-period preceding the 1996 earthquake. These anomalies consist on a temporal and spatial sequence of gas emissions in the epicentral area and on chloride and lead concentration variations in a bottled mineral water 25 km north to the main shock epicenter. The proposed model processed in a standard elastic half-space, consists of creep on a low-angle crustal normal-fault, generating volumetric strain field changes over a distance of 25 km from the epicentral area. This model is able to constrain not only the mechanisms and the locations of the geochemical anomalies, but also their timing and probable casual links to the triggering of the impending major event. Also, the active extension proposed here is compatible with seismological observations in the Pyrenees. Thus, the possibility of such creep, which can be considered as a slow-slip event, is discussed in the Pyrenean tectonic and geological context. The model is discussed and compared to previous proposed models on precursor processes of earthquakes, especially concerning the preparation zone concept. Finally, a complete seismic scenario over the period beginning 15 days before the quake and ending 5 days after is proposed and discussed.     

Multi-tracer tests to evaluate the hydraulic setting of a complex aquifer system (Brévilles spring catchment, France)

For good management of groundwater resources, and to comply with European and national regulations, a detailed understanding of an aquifer’s hydraulic setting is required. In order to better characterize a sandy aquifer that is affected by diffuse pollution (Brévilles spring catchment, Val d’Oise, France), and to quantify the transfer time in the saturated zone, a multi-tracer test involving a new technique, the ‘finite volume point dilution method’, has been performed in natural flow conditions. In November 2005, injections of four different tracers took place in four piezometers involving different locations and depths in the aquifer. Recovery of the tracers was observed at two different places near the aquifer outlet. A particularly long and unusual monitoring exercise (27 months) demonstrated the existence of several different velocities within the sandy layer, which seems to be linked to the decrease of hydraulic conductivity with depth. The new insight and parameter quantification brought by interpretation of these tests contribute to a better characterization of the saturated zone. The particularly long-term monitoring exercise also gives new information to understand and forecast the trend and persistence of groundwater contamination by pesticides in the catchment.     

Typecasting prospective Au-bearing sedimentary lithologies using sedimentary geochemistry and Nd isotopes in poorly exposed Proterozoic basins of the Tanami region, Northern Australia

The development of a regional stratigraphy in Palaeoproterozoic basins within the Tanami region, Northern Australia has been hindered by the difficulty of discriminating sedimentary units and facies across this isolated and poorly exposed basin. A regional stratigraphy is important as it provides constraints on sedimentary basin evolution and assists in gold exploration. Gold is known to be more concentrated in certain rock formations. Based on Nd isotopes and whole rock geochemistry, five main sedimentary events have been identified in the Tanami region. Some sedimentary units were derived from homogeneous local sources, whereas others contain evidence of a well-mixed fine-grained remote provenance. Within the basins, major gold-bearing lithologies are characterised by mafic source indicators: (1) high Cr/Th ratios; (2) low Th/Sc ratios; (3) low (La/Yb)_PAAS ratios relative to Post-Archaean Average Shale (Taylor and McLennan 1985); (4) Eu anomalies equal to ∼1; and (5) distinctive ranges in initial ε _Nd values. Potential future exploration target areas have been identified in the Tanami region at the Cashel and Sunline prospects using these geochemical parameters.