Search

In this paper we present the potential of a new compact superconducting gravimeter (GWR iGrav) designed for groundwater monitoring. At first, 3 yr of continuous gravity data are evaluated and the performance of the instrument is investigated. With repeated absolute gravity measurements using a Micro-g Lacoste FG5, the calibration factor (-894.8 nm s-2 V-1) and the long-term drift of this instrument (45 nm s-2 yr-1) are estimated for the first time with a high precision and found to be respectively constant and linear for this particular iGrav. The low noise level performance is found similar to those of previous superconducting gravimeters and leads to gravity residuals coherent with local hydrology. The iGrav is located in a fully instrumented hydrogeophysical observatory on the Durzon karstic basin (Larzac plateau, south of France). Rain gauges and a flux tower (evapo-transpiration measurements) are used to evaluate the groundwater mass balance at the local scale. Water mass balance demonstrates that the karst is only capacitive: all the rainwater is temporarily stored in the matrix and fast transfers to the spring through fractures are insignificant in this area. Moreover, the upper part of the karst around the observatory appears to be representative of slow transfer of the whole catchment. Indeed, slow transfer estimated on the site fully supports the low-flow discharge at the only spring which represents all groundwater outflows from the catchment. In the last part of the paper, reservoir models are used to characterize the water transfer and storage processes. Particular highlights are done on the advantages of continuous gravity data (compared to repeated campaigns) and on the importance of local accurate meteorological data to limit misinterpretation of the gravity observations. The results are complementary with previous studies at the basin scale and show a clear potential for continuous gravity time-series assimilation in hydrological simulations, even on heterogeneous karstic systems. [ABSTRACT FROM AUTHOR] more...

Sensitive instruments like strainmeters and tiltmeters are necessary for measuring slowly varying low amplitude Earth deformations. Nonetheless, laser and fibre interferometers are particularly suitable for interrogating such instruments due to their extreme precision and accuracy. In this paper, a practical design of a simple pendulum borehole tiltmeter based on laser fibre interferometric displacement sensors is presented. A prototype instrument has been constructed using welded borosilicate with a pendulum length of 0.85 m resulting in a main resonance frequency of 0.6 Hz. By implementing three coplanar extrinsic fibre Fabry-Perot interferometric probes and appropriate signal filtering, our instrument provides tilt measurements that are insensitive to parasitic deformations caused by temperature and pressure variations. This prototype has been installed in an underground facility (Rustrel, France) where results show accurate measurements of Earth strains derived from Earth and ocean tides, local hydrologic effects, as well as local and remote earthquakes. The large dynamic range and the high sensitivity of this tiltmeter render it an invaluable tool for numerous geophysical applications such as transient fault motion, volcanic strain and reservoir monitoring. [ABSTRACT FROM AUTHOR] more...

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995-1997, respectively, and remeasured in 2008 and 2010.We obtain a deformation field with velocities less than 1mmyr-1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so,we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N--S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults. [ABSTRACT FROM AUTHOR] more...

Iran is an active continental domain accommodating the convergence between the Arabia and Eurasia plates. In northwestern Iran, deformation between the Central Iranian block and the Caucasus domain is mainly accommodated by right lateral strike-slip on the Tabriz fault. Cities and villages, including the city of Tabriz, have been destroyed by several strong historical earthquakes (M ∼ 7). In this study, we compare the slip-rates estimated from geodetic measurements (radar interferometry and GPS) with those determined by dating a geomorphological offset of an alluvial fan along the Tabriz fault.The GPS measurements along two profiles normal to the Tabriz fault suggest a slip-rate of 7.3 ± 1.3 mm yr−1. The persistent scatterer radar interferometry analysis of Envisat satellite archives from 2003 to 2010 shows a velocity gradient (6 ± 3 mm yr−1) across the Tabriz fault in agreement with GPS results. Moreover, it reveals that most of the area located south of the Tabriz fault is affected by subsidence, and that some sections of the fault probably act as barriers to fluid migration which may have an impact on its mechanical behaviour. West of Tabriz morphotectonic investigations on an alluvial fan surface show a right-lateral cumulative offset of 320 ± 40 m. Luminescence analyses of the coarse matrix alluvial fan deposits provide an age of 46 ± 3 ka. This yields a slip-rate comprised between 6.5 and 7.3 mm yr−1 along this segment. These results suggest that the Late Quaternary slip-rate is in agreement with the present-day slip-rate estimated by geodetic measurements, showing no slip-rate changes during the past 45 000 yr. Short-term variations within the 45 000 yr related to temporal earthquake clustering over few seismic cycles cannot be ruled out, but if they exist, they do not affect the geodetic and the geomorphological estimates. This study is in agreement with previous ones suggesting that long-term slip-rates (i.e. averaged over several tens of seismic cycles) are consistent with geodetic estimated slip-rates (i.e. extrapolated from few years of interseismic observations), and suggests that perturbations of fault slip-rates are related to variations over few seismic cycles. [ABSTRACT FROM AUTHOR] more...

We use 3 yr of surface absolute gravity measurements at three sites on the Larzac plateau (France) to quantify the changes induced by topography and the building on gravity time-series, with respect to an idealized infinite slab approximation. Indeed, local topography and buildings housing ground-based gravity measurement have an effect on the distribution of water storage changes, therefore affecting the associated gravity signal. We first calculate the effects of surrounding topography and building dimensions on the gravity attraction for a uniform layer of water. We show that a gravimetric interpretation of water storage change using an infinite slab, the so-called Bouguer approximation, is generally not suitable. We propose to split the time varying gravity signal in two parts (1) a surface component including topographic and building effects (2) a deep component associated to underground water transfer. A reservoir modelling scheme is herein presented to remove the local site effects and to invert for the effective hydrological properties of the unsaturated zone. We show that effective time constants associated to water transfer vary greatly from site to site. We propose that our modelling scheme can be used to correct for the local site effects on gravity at any site presenting a departure from a flat topography. Depending on sites, the corrected signal can exceed measured values by 5–15 μGal, corresponding to 120–380 mm of water using the Bouguer slab formula. Our approach only requires the knowledge of daily precipitation corrected for evapotranspiration. Therefore, it can be a useful tool to correct any kind of gravimetric time-series data. [ABSTRACT FROM PUBLISHER] more...

SUMMARY This paper presents an inversion method using the interseismic velocity field to determine effective rigidity of the lithosphere. The method is based on the minimization of a cost function defined as the quadratic measure of the difference between measured and modelled velocity fields on a discrete set of points. The continuous mapping of the rigidity is fulfilled with a limited set of parameters and the forward solution is achieved using a plane stress finite element code. The computation of the cost function gradient in the parameters' space allows one to iteratively find the best parameters set through a suitable optimization algorithm. We first design a benchmark including an abrupt rigidity variation that cannot be described by a continuous function. For such a case, we show that increasing the number of parameters is a way to accurately describe sharp variations of the rigidity map. Then, we use a dense GPS velocity field over the southwestern United States to estimate the corresponding rigidity variations for different spatial resolutions of the parameters' grid. We analyse the conceptual and practical difficulties associated with our methodology. Finally, rigidity maps obtained by our inversion method in southwestern United States and particularly across the San Andreas Fault System are reviewed and compared to current plate rigidity estimates and geophysical data over this area. [ABSTRACT FROM AUTHOR] more...

The Strait of Hormuz area is a transition zone between the continental collision of the Zagros (west) and the subduction of an oceanic part of the Arabian Plate beneath the Makran wedge (east). Geology and recent GPS measurements indicate that about 15 mm yr−1 of relative motion in N10°E direction is accommodated by two major fault systems: (1) the NNW-trending Minab-Zendan-Palami (MZP) fault system that connects the Main Zagros Thrust (MZT) to the inner Makran thrust system and the Frontal subduction thrust and (2) the N-trending Sabzevaran-Kahnuj-Jiroft (SKJ) fault system that bounds the Jazmurian depression to the west. We use dense GPS measurements along four transects across these fault systems in order to determine the strains spatial distribution. The northern GPS transect confirms the total fault slip rates for both fault systems estimated by the tectonic analyses (about 10 and 7.3 mm yr−1 in N10° direction across the MZP and SKJ fault systems, respectively). For both fault systems, the elastic deformation spreads over shear zones that are several tens of kilometres wide. However, transects located close to latitude 27°N reveal a much narrower shear zone (∼10 km) for the MZP fault system. Moreover, we confirm that most of the present-day strain is transferred towards the frontal subduction thrust rather than towards the inner Makran thrusts. In order to complement this new GPS velocity field with spatially dense measurements, we processed a set of ERS radar images by the radar interferometry (InSAR) technique. We used both a ‘stacking’ and a ‘persistant-scatterers’ approach to differentiate the ground deformation signal which spatial gradient is expected to be very low, from the atmospheric signal. Results from these interferograms appear to be relatively in agreement with the GPS-determined strain distribution. Nevertheless, they confirm the absence of any superficial creep behaviour since no sharp discontinuity on interferometric phase can be noted on any interferogram. Finally, we use a purely kinematic ‘block model’ inversion process to calculate slip rates and locking depths for each fault system from our GPS measurements. These models suggest that the relative quiescence over the last 200 yr has certainly produced a slip deficit as high as 2 m. So, we may wonder if the MZP fault system is not late in the interseismic phase of its earthquake cycle. [ABSTRACT FROM AUTHOR] more...

The interpretation of coseismic surface deformation measurements through inversion techniques is of major importance to understand the mechanical behaviour of a seismic fault. Dense geodetic data sets in the vicinity of the ruptured fault provide unique constraints on detailed fault geometry and slip distribution at depth, making them complementary to seismological data. Bam earthquake (Mw 6.6, 2003 December 26) induced surface deformation has been precisely mapped by Envisat ASAR interferometry and by subpixel correlation techniques applied to Spot-5 and ASAR amplitude images. These oblique and horizontal estimations of deformation have been completed with one levelling profile along the main road crossing the rupture from west to east. We process these data (separately and jointly) in a two-step inversion technique, within the elastic half-space theory framework. Our objective is to determine the dislocation model at depth that satisfies simultaneously all the geodetic constraints. Also, we estimate the relative contribution of each geodetic data set to this inversion process. We first use a stochastic direct approach called neighbourhood algorithm in order to estimate the average characteristics of the rupture, and their relative uncertainty. Constraining in this way the geometry of the ruptured fault, we then linearize the inverse problem and compute the slip distribution on the fault using a standard weighted least-square technique, assuming the solution is smooth to some degree. At each step, we discuss the optimal models, their stability as well as the relative influence of each data set on the derived models parameters. Our preferred model reveals a shallow dislocation on a quasi-vertical fault, slightly dipping towards east. The slip vector has a strike-slip component as high as 2 m, while the dip-slip component seems negligible. However, the estimation of the resolution matrices emphasizes the fact that the details of deep fault slip distribution remain out of the scope of this ill-conditioned inverse problem. Yet, our preferred model suggests a main dislocation limited at depth between 1 and 6 km. By contrast, the aftershocks observed in the months following the earthquake are located just beneath the estimated main shock. [ABSTRACT FROM AUTHOR] more...

Measurements on either side of the Kazerun fault system in the Zagros Mountain Belt, Iran, show that the accommodation of the convergence of the Arabian and Eurasian Plates differs across the region. In northwest Zagros, the deformation is partitioned as 3–6 mm yr−1 of shortening perpendicular to the axis of the mountain belt, and 4–6 mm yr−1 of dextral strike-slip motion on northwest–southeast trending faults. No individual strike-slip fault seems to slip at a rate higher than ∼2 mm yr−1. In southeast Zagros, the deformation is pure shortening of occurring perpendicular to the simple folded belt and restricted to the Persian Gulf shore. The fact that most of the deformation is located in front of the simple folded belt, close to the Persian Gulf, while seismicity is more widely spread across the mountain belt, confirms the decoupling of the surface sedimentary layers from the seismogenic basement. A comparison with the folding and topography corroborates a southwestward propagation of the surface deformation. The difference in deformation between the two regions suggests that right-lateral shear cumulates on the north–south trending Kazerun strike-slip fault system to . [ABSTRACT FROM AUTHOR] more...

A combined analysis of the geodetic strain-rate field and the strain-rate field deduced from the seismicity allows us to define the style of deformation and to distinguish seismic from aseismic deformation. We perform this analysis in Iran where the present-day tectonics results from the north–south convergence between the plates of Arabia to the south and Eurasia to the north. The data consist of velocities measured with a GPS network of 28 benchmarks and of instrumental and historical earthquake catalogues. The axes of the seismic strain-rate tensor have similar orientations to those deduced from the GPS velocity field. This indicates that the seismicity can be used to improve GPS information on the style and the orientation of the deformation. Comparison of seismic and geodetic strain rates indicates that highly strained zones experience mainly aseismic deformation in southern Iran and seismic deformation in northern Iran. A large contrast is observed between the Zagros (less than 5 per cent seismic deformation) and the Alborz–Kopet-Dag regions (more than 30–100 per cent seismic deformation). The distribution of the seismic/geodetic ratio correlates with the distribution of large earthquakes: intensive, low-magnitude seismicity is observed in the Zagros whereas the largest earthquakes occur in northern Iran. The contrast of seismic deformation between the Zagros and peri-Caspian mountains is confirmed considering 300 or 1000 yr of seismicity rather than 100 or 200 yr. [ABSTRACT FROM AUTHOR] more...

A network of 27 GPS sites was implemented in Iran and northern Oman to measure displacements in this part of the Alpine–Himalayan mountain belt. We present and interpret the results of two surveys performed in 1999 September and 2001 October. GPS sites in Oman show northward motion of the Arabian Plate relative to Eurasia slower than the NUVEL-1A estimates (e.g. at instead of 30.5 mm yr−1 at N6°E at Bahrain longitude). We define a GPS Arabia–Eurasia Euler vector of . The Arabia–Eurasia convergence is accommodated differently in eastern and western Iran. East of 58°E, most of the shortening is accommodated by the Makran subduction zone (19.5 ± 2 mm yr−1) and less by the Kopet-Dag (6.5 ± 2 mm yr−1). West of 58°E, the deformation is distributed in separate fold and thrust belts. At the longitude of Tehran, the Zagros and the Alborz mountain ranges accommodate 6.5 ± 2 mm yr−1 and respectively. The right-lateral displacement along the Main Recent Fault in the northern Zagros is about , smaller than what was generally expected. By contrast, large right-lateral displacement takes place in northwestern Iran (up to 8 ± mm yr−1). The Central Iranian Block is characterized by coherent plate motion (internal deformation <2 mm yr−1). Sites east of 61°E show very low displacements relative to Eurasia. The kinematic contrast between eastern and western Iran is accommodated by strike-slip motions along the Lut Block. To the south, the transition zone between Zagros and Makran is under transpression with right-lateral displacements of 11 ± 2 mm yr−1. [ABSTRACT FROM AUTHOR] more...

Mechanisms that control seismic activity in low strain rate areas such as western Europe remain poorly understood. For example, in spite of low shortening rates of <0.5 mm/ yr, the Western Alps and the Pyrenees are underlain by moderate but frequent seismicity detectable by instruments. Beneath the elevated part of these mountain ranges, analysis of earthquake focal mechanisms indicates extension, which is commonly interpreted as the result of gravitational collapse. Here we show that erosional processes are the predominant control on present-day deformation and seismicity. We demonstrate, using finite element modeling, that erosion induces extension and rock uplift of the elevated region of mountain ranges accommodating relatively low overall convergence. Our results suggest that an erosion rate of ~1 mm/yr can lead to extension in mountain ranges accommodating significant shortening of <3 mm/yr. Based on this study, the seismotectonic framework and seismic hazard assessment for low strain rate areas need to be revisited, because erosion-related earthquakes could increase seismic hazard. [ABSTRACT FROM AUTHOR] more...

A non-acid-based, di-functional epoxide, neopentyl glycol diglycidyl ether (NPGDGE), was used to modify cotton fabrics. Direct characterization of the modified cotton was conducted by Nuclear Magnetic Resonance (NMR) without grinding the fabric into a fine powder. NaOH and MgBr 2 were compared in catalyzing the reaction between the epoxide groups of NPGDGE and the hydroxyl groups of cellulose. Possible reaction routes were discussed. Scanning electron microscopy (SEM) images showed that while the MgBr 2 -catalyzed reaction resulted in self-polymerization of NPGDGE, the NaOH-catalyzed reaction did not. Fourier transform infrared spectroscopy (FTIR) showed that at high NaOH concentration cellulose restructures from allomorph I to II. NMR studies verified the incorporation of NPGDGE into cotton fabrics with a clear NMR signal, and confirmed that at higher NaOH concentration the efficiency of grafting of NPGDGE was increased. This demonstrates that use of solid state NMR directly on woven fabric samples can simultaneously characterize chemical modification and crystalline polymorph of cotton. No loss of tensile strength was observed for cotton fabrics modified with NPGDGE., (Copyright © 2017 Elsevier Ltd. All rights reserved.) more...