Natural and man-made induced hydrological signals,detected by high resolution tilt observations at the Geodynamic Observatory Moxa/Germany

It is well known, that high resolution borehole tiltmeters are able to observe deformations, caused by hydrological variations. The quantitative coherence is often unexplained, especially if the sources of deformation can be based on both natural as well as man-made hydrological variations. Since 1999 tilt observations have been taken at the Geodynamic Observatory Moxa in Thuringia/Germany. In two 50 m and one 100 m deep boreholes the ASKANIA tiltmeters are installed. The high quality of the recorded tilt data can be proved by the analysis of well known geodynamic signals like the tides of the solid Earth and the free modes of the Earth. Here we focus on investigations of induced tilt signals caused by pore pressure changes due to precipitation and/or ground water level changes and, in addition, on man-made induced pore pressure variations. The correlation of natural ground water level changes with the observed tilt data can be shown by different events of precipitation and snow melting. However, also the load effect of a big truck yields a small elastic deformation which is clearly detectable in the ground water level recording. The correlated tilt effect is discussed regarding changes of the tilt amplitude and the orientation of the induced pendulum tip movement during the load phase.     

Combining Electrical Resistivity Tomography and Ground Penetrating Radar to study geological structuring of karst Unsaturated Zone

This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).     

GPR to constrain ERT data inversion in cavity searching: Theoretical and practical applications in archeology

I used theoretical forward models to show that a cavity embedded in a stratified sedimentary sequence can induce an equivalence problem in the ERT data inversion. Conductive top soil increases the misfit between the ground feature and the ERT model. The misfit depends on array and stratigraphy sequences. The latter induce an equivalence problem that manifests itself as wrong cavity depth positioning. The misfit is greater in the data acquired with Schlumberger array than with dipole–dipole.The ambiguity of ERT data inversion problems was tested in the detection of cavities linked to an 8th–6th century B.C. Sabine tomb, 3 m wide × 3 m long × 2 m high, excavated from a shaly gray volcanic ash (cinerite) layer covered by semi-lithoid tuff and top soil layers. In the real study I reduced the ambiguity in the inverse problem of ERT data using a priori information on geometry and resistivity of the cavity. The constrains were carried out from georadar data acquired with 80 and 200 MHz antenna. I demonstrate that this procedure has a practical application in cavity detection, and is a key to the reduction of the uncertainty inherent in the inversion process of ERT data.     

Combined electrical and electromagnetic imaging of hot fluids within fractured rock in rugged Himalayan terrain

Integrated electrical resistance tomography (ERT) and short-offset transient electromagnetic (TEM) measurements were carried out to investigate a geothermal area in the Main Central Thrust (MCT) zone of Garhwal Himalayan region, India. The study area is located around Helang on either side of Alaknanda River and it is dotted with hot water springs with water temperature of 45°–55 °C emerging at the surface.To assess the geothermal potential and its lateral and vertical extension in and around the hot water springs in the study area, 7 ERT profiles and 21 TEM stations on 7 profiles were established around the hot water spring and at far distant locations. The 2D inversion of ERT data indicates a low resistivity (< 50 Ωm) zone in the vicinity of hot springs, which appears to be associated with an underground water channel through the fractured rock. The bedrock resistivity is very high (> 1000 Ωm) whereas the resistivity of the weathered near surface soil at a far distant location from the hot spring is low (< 100 Ωm) again. A common feature of all TEM data is the sign reversal observed at roughly 10 μs. The consistent sign reversal in all TEM data indicates the existence of the multi-dimensionality of the geoelectrical structure. Therefore, the TEM data were treated by using the SLDM (Spectral Lanczos Decomposition Method) 2D/3D forward modeling code based on the finite difference algorithm. The resistivity structure obtained from ERT data was used as an input for the modeling of TEM data. Based on the joint analysis of the ERT and TEM data it can be inferred that geothermal anomalies associated with the hot spring in the MCT zone are a local feature appearing as a low resistivity zone (< 50 Ωm) at shallow depth (< 100 m) in the vicinity of the hot spring region.     

Global-scale validation of model-based load deformation of the Earth's crust from continental watermass and atmospheric pressure variations using GPS

Temporal mass variations in the continental hydrosphere and in the atmosphere lead to changes in the gravitational potential field that are associated with load-induced deformation of the Earth’s crust. Therefore, models that compute continental water storage and atmospheric pressure can be validated by measured load deformation time series. In this study, water mass variations as computed by the WaterGAP Global Hydrology Model (WGHM) and surface pressure as provided by the reanalysis product of the National Centers for Environmental Prediction describe the hydrological and atmospheric pressure loading, respectively. GPS observations from 14 years at 208 stations world-wide were reprocessed to estimate admittance factors for the associated load deformation time series in order to determine how well the model-based deformation fits to real data. We found that such site-specific scaling factors can be identified separately for water mass and air pressure loading. Regarding water storage variation as computed by WGHM, weighted global mean admittances are 0.74 ± 0.09, 0.66 ± 0.10, 0.90 ± 0.06 for the north, east and vertical component, respectively. For the dominant vertical component, there is a rather good fit to the observed displacements, and, averaged over all sites, WGHM is found to slightly overestimate temporal variations of water storage. For Europe and North America, with a dense GPS network, site-specific admittances show a good spatial coherence. Regarding regional over- or underestimation of WGHM water storage variations, they agree well with GRACE gravity field data. Globally averaged admittance estimates of pre-computed atmospheric loading displacements provided by the Goddard Geodetic VLBI Group were determined to be 0.88 ± 0.04, 0.97 ± 0.08, 1.13 ± 0.01 for the north, east and vertical, respectively. Here, a relatively large discrepancy for the dominant vertical component indicates an underestimation of corresponding loading predictions.     

Mapping the spatial distribution of subsurface saline material in the Darling River valley

In the Australian landscape larg stores of soluble salt are present naturally. In many cases it is attributable to salts entrapped as marine sediment in earlier geological time. At the district level, the need for information on the presence of saline subsurface material is increasing, particularly for its application to salinity hazard assessment and environmental management. This is the case in irrigated areas, where changes in hydrology can result in secondary salinisation. To reduce the expense, environmental studies use a regression relationship to make use of more readily observed measurements (e.g. electromagnetic (EM) data) which are strongly correlated with the variable of interest. In this investigation a methodology is outlined for mapping the spatial distribution of average subsurface (6–12 m) salinity (EC_e — mS m^? 1) using an environmental correlation with EM34 survey data collected across the Bourke Irrigation District (BID) in the Darling River valley. The EM34 is used in the horizontal dipole mode at coil configurations of 10 (EM34-10), 20 (EM34-20), and 40 (EM34-40). A multiple-linear regression (MLR) relationship is established between average subsurface EC_e and the three EM34 signal data using a forward modeling stepwise linear modeling approach. The spatial distribution of average subsurface salinity generally reflects the known surface expression of point-source salinisation and provides information for future environmental monitoring and natural resource management. The generation of EM34 data on various contrived grids (i.e. 1, 1.5, 2. 2.5 and 3 km) indicates that in terms of accuracy, the data available on the 0.5 (RMSE = 188) and 1 km (RMSE = 283) grid are best, with the least biased predictions achieved using 1 (ME = ? 1) and 2 km (ME = 12) grids. Viewing the spatial distribution of subsurface saline material showed that the 0.5 km spacing is optimal, particularly in order to account for short-range spatial variation between various physiographic units. The Relative Improvement (RI) shows that increasing EM grids from 1, 1.5, 2, 2.5 to 3 km gave RI of ? 53, ? 100%, ? 107%, ? 128% and ? 140%, respectively. We conclude that at a minimum a 1 km grid is needed for reconnaissance EM34 surveying.     

Expected temporal absolute gravity change across the Taiwanese Orogen,a modeling approach

The island of Taiwan is located on the convergent boundary between the Philippine Sea plate and the Chinese continental margin. It offers very active mountain building and collapsing processes well illustrated by the rugged topography, rapid uplift and denudation, young tectonic landforms, active faulting and numerous earthquakes. In this paper, using simple models, we have estimated vertical movements and associated absolute gravity variations which can be expected along a profile crossing the southern part of the island and probably suffering the highest rates of rising. The two different tectonic styles proposed for the island, thin-skinned and thick-skinned, were taken into account. Horizontal and vertical movements were modeled by an elastic deformation code. Gravity variations due to these deformations are then modeled at a second step. They are dominated by plate and free-air effects, i.e. elevation of the topography, with several μGal yr⁻¹. By comparison, gravity changes generated by mass transfers are weak: maximum 0.1 μGal yr⁻¹ with the thin-skinned tectonic and 0.3 μGal yr⁻¹ with the thick-skinned tectonic. Though elastic rheology has limitations, this modeling offers interesting results on what gravity signal can be expected from the AGTO project (Absolute Gravity in the Taiwanese Orogen), which proposes to study the dynamic of these mountain ranges using absolute gravimetry (AG) and also including relative gravimetry (RG) and GPS measurements.     

Investigation of groundwater resources using controlled-source radio magnetotellurics (CSRMT) in glacial deposits in Heby,Sweden

We have combined tensor radio magnetotelluric- (RMT, 15–250 kHz) and controlled source tensor magnetotelluric (CSTMT, 1–12 kHz) data for the mapping of aquifers in gravel formations lying in between crystalline bedrock and clay rich sediments in the Heby area some 40 km west of Uppsala in Sweden. The estimated transfer functions, the impedance tensor and the tipper vector generally satisfy 1D or 2D necessary conditions except for the lowest CSTMT frequencies where near field effects become more dominant.The data measured from 8 profiles were inverted with the Rebocc code of Siripunvaraporn and Egbert (2000) assuming plane wave conditions. This meant that only 12 frequencies in the range of 4–180 kHz could be used. The four lowest frequencies of CSTMT in the range of 1–2.8 kHz were excluded because of source effects. Data from all profiles were inverted with a starting model of 100 Ω-m and a relative error floor of 0.02 on apparent resistivity, corresponding to less than 1^° on phase. Tipper vectors are generally small except when source effects become dominant in the lowest frequencies of CSTMT and were therefore not used for inversion. Comparing with models derived from vertical electrical soundings, refraction and reflection seismic data as well as ground truth from exploration wells assessed the reliability of the deep part of the models. Furthermore we carried out a non-linear resolution analysis to better quantify the depth extent of the aquifers.The inverted models from the Heby area show well the thickness variations of glacial deposits overlying crystalline bedrock. Generally, the upper 20 m of the models have resistivities below 40 Ω-m, taken to represent clay rich formations. Below the clay layer resistivities increase to about 40–400 Ω-m, interpreted to represent sand/gravel formations with a maximum thickness of about 40 m and a width of several hundred metres. This is a potential aquifer that extends in approximately N–S direction for some kilometres.     

Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements

The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500 km², as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.We derive a set of local Green's functions defined globally on a 1° × 1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0 mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.     

Geological,tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs

We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.     

Artificial earthquake test of gas supply networks

In this paper, seismic behavior of gas supply networks is investigated by constructing a 24 m×24 m buried network and conducting an artificial earthquake test by detonating TNT explosives. The pipe network consists of typical welded steel pipes used in gas supply networks. Ground motion is produced by detonating 30 kg of TNT explosives buried 5 m below ground surface and 15 m away from the network. The test preparations are presented, including site selection, network layout, ground motion generation, and arrangements of different sensors. The measured ground acceleration, strain, and acceleration of welded steel pipes, and pipe–soil relative slippage are discussed. The deformation patterns of welded steel pipes are also analyzed and explained, including axial and bending deformations, systematic network characteristics, deformation relationship between ground and pipes, and dynamic response of buried pipes.     

Post-cyclic settlement and tilting potential of mat foundations

A simplified procedure for the estimation of cyclic shear stresses on foundation soil layers due to interactions of seismic excitation, foundation mat and overlying structural system, incorporating both kinematic and inertial aspects was described earlier in Cetin et al. [1], [2]. This simplified procedure is now used, along with laboratory-based cyclic shear and volumetric strain relationships for the assessment of cyclically induced settlement and tilting potential of mat foundations. The proposed methodology is calibrated by well-documented seismic foundation performance case histories of 3–6 story, relatively rigid residential structures with no basements. Immediately after 1999 Kocaeli and Duzce, Turkey earthquakes, foundation settlements of these case history buildings were carefully mapped relative to available elevation benchmarks, such as relatively rigid concrete pavements, drainage pipes and ditches, and entrance stairs, located in the immediate vicinity of the buildings. Tilt angles of these buildings were also mapped in orthogonal multi-directions. Relative settlement and tilt angle values vary in the range of 0–60 cm and 0–5.5 rad, respectively. As part of subsurface characterization studies, SPT and seismic CPTU were performed. For the retrieved disturbed and undisturbed samples, an intensive laboratory testing program including sieve analysis, Atterberg Limit, hydrometer, oedometer, static and cyclic triaxial and cyclic simple shear tests was executed. The foundation soil profiles generally consist of silty soils, sand–silt mixtures and silt–clay mixtures of SPT-N values varying in the range of 2–5 blows/30 cm in the upper 5 m and gradually increases up to a maximum value of 25 blows/30 cm beyond depths of 5–8 m. As part of the proposed framework, displacement potential indices, defined as the product of induced-cyclic strains and the thickness of soil sub-layers, were estimated, the sum of which produce settlement and tilting potentials. Then, these settlement and tilting potentials were calibrated against recorded settlement and tilt values through a statistically based calibration scheme. The proposed simplified procedure is shown to reliably and precisely capture both deviatoric and volumetric components of post-cyclic settlements, as well as tilting responses of mat foundations.     

Aminostratigraphy of surface and subsurface Quaternary sediments,North Carolina coastal plain,USA

The Quaternary stratigraphy and geochronology of the Albemarle Embayment of the North Carolina (NC) Coastal Plain is examined using amino acid racemization (AAR) in marine mollusks, in combination with geophysical, lithologic, and biostratigraphic analysis of 28 rotasonic cores drilled between 2002 and 2006. The Albemarle Embayment is bounded by structural highs to the north and south, and Quaternary strata thin westward toward the Suffolk paleoshoreline, frequently referred to as the Suffolk Scarp. The Quaternary section is up to ～90 m thick, consists of a variety of estuarine, shelf, back-barrier, and lagoonal deposits, and has been influenced by multiple sea-level cycles. The temporal resolution of the amino acid racemization method is tested statistically and with the stratigraphic control provided by this geologic framework, and it is then applied to the correlation and age estimation of subsurface units throughout the region. Over 500 specimens (primarily Mercenaria and Mulinia) from the subsurface section have been analyzed using either gas chromatographic (GC) or reverse-phase liquid chromatographic (RPLC) techniques. The subsurface stratigraphic data are compared with AAR results from numerous natural or excavated exposures from the surrounding region, as well as results from NC beach collections, to develop a comprehensive aminostratigraphic database for the entire Quaternary record within the NC coastal system. Age mixing, recognized in the beach collections, is also seen in subsurface sections, usually where major seismic reflections or core lithology indicate the presence of stratigraphic discontinuities. Kinetic models for racemization are tested within the regional stratigraphic framework, using either radiocarbon or U-series calibrations or comparison with regional biostratigraphy. Three major Pleistocene aminozones [AZ2, AZ3, and AZ4] are found throughout the region, all being found in superposition in several cores. Each can be subdivided, yielding a total of at least eight stratigraphically and statistically distinct aminozones. Kinetic modeling, supplemented with local calibration, indicates that these aminozones represent depositional events ranging from ～80 ka to nearly 2 Ma. Three prominent seismic reflections are interpreted to represent the base of the early, middle, and late Pleistocene, respectively, roughly 2 Ma, 800 ka, and 130 ka. The large number of samples and the available stratigraphic control provide new insights into the capabilities and limitations of aminostratigraphic methods in assessing relative and numerical ages of Atlantic Coastal Plain Quaternary deposits.     

Nutrients,irradiance, and mixing as factors regulating primary production in coastal waters impacted by the Mississippi River plume

Relationships among primary production, chlorophyll, nutrients, irradiance and mixing processes were examined along the salinity gradient in the Mississippi River outflow region. A series of six cruises were conducted during 1988–1992 at various times of year and stages of river discharge. Maximum values of biomass and primary production were typically observed at intermediate salinities and coincided with non-conservative decreases in nutrients along the salinity gradient. Highest values of productivity (>10 gC m⁻² d⁻¹) and biomass (>30 mg chlorophyll a m⁻³) were observed in April 1988, July–August 1990 and April–May 1992; values were lower in March and September 1991. Rates of primary production were apparently constrained by low irradiance and mixing in the more turbid, low salinity regions of the plume, and by nutrient limitation outside the plume. Highest values of primary production occurred at stations where surface nutrient concentrations exhibited large deviations from conservative mixing relationships, indicating that depletion of nutrients was related to phytoplankton uptake. Mixing and advection were important in determining the location and magnitude of primary production maxima and nutrient depletion. In addition to growth within plume surface waters, enhanced growth and/or retention of biomass may have occurred in longer residence time waters at the plume edge and/or beneath the surface plume. Vertical structure of some plume stations revealed the presence of subsurface biomass maxima in intermediate salinity water that was depleted in nutrients presumably by uptake processes. Exchange between subsurface water and the surface plume apparently contributed to the reduction in nutrients at intermediate salinities in the surface layer. DIN (=nitrate+nitrite+ammonium) : PO₄ (=phosphate) ratios in river water varied seasonally, with high values in winter and spring and low values in late summer and fall. Periods of high DIN : PO₄ ratios in river nutrients coincided with cruises when surface nutrient concentrations and their ratios indicated a high probability for P limitation. N limitation was more likely to occur at high salinities and during late summer and fall. Evidence for Si limitation was also found, particularly in spring.     

Bioavailability of heavy metals in water and sediments from a typical Mediterranean Bay (Málaga Bay,Region of Andalucía,Southern Spain)

Concentrations of heavy metals were measured in sediment and water from Málaga Bay (South Spain). In the later twentieth century, cities such as Málaga, have suffered the impact of mass summer tourism. The ancient industrial activities, and the actual urbanization and coastal development, recreation and tourism, wastewaters treatment facilities, have been sources of marine pollution. In sediments, Ni was the most disturbing metal because Ni concentrations exceeded the effects range low (ERL), concentration at which toxicity could start to be observed in 85% of the samples analyzed. The metal bioavailability decreased in the order: Cd > Ni > Pb > Cu > Cr. In the sea water samples, Cd and Pb were the most disturbing metals because they exceeded the continuous criteria concentration (CCC) of US EPA in a 22.5% and 10.0% of the samples, respectively. Statistical analyses (ANOVA, PCA, CA) were performed.     

Water mass variation in the Mediterranean and Black Seas

The mass-induced sea level variability and the net mass transport between Mediterranean Sea and Black Sea are derived for the interval between August 2002 and July 2008 from satellite-based observations and from model data. We construct in each basin two time series representing the basin mean mass signal in terms of equivalent water height. The first series is obtained from steric-corrected altimetry while the other is deduced from GRACE data corrected for the contamination by continental hydrology. The series show a good agreement in terms of annual and inter-annual signals, which is in line with earlier works, although different model corrections influence the consistency in terms of seasonal signal and trend.In the Mediterranean Sea, we obtain the best agreement using a steric correction from the regional oceanographic model MFSTEP and a continental hydrological leakage correction derived from the global continental hydrological model WaterGAP2. The inter-annual time series show a correlation of 0.85 and a root mean square (RMS) difference of 15 mm. The two estimates have similar accuracy and their annual amplitude and phase agree within 3 mm and 23 days respectively. The GRACE-derived mass-induced sea level variability yields an annual amplitude of 27 ± 5 mm peaking in December and a trend of 5.3 ± 1.9 mm/yr, which deviates within 3 mm/yr from the altimetry-derived estimate.In the Black Sea, the series are less consistent, with lower accuracy of the GRACE-derived estimate, but still show a promising agreement considering the smaller size of the basin. The best agreement is realized choosing the corrections from WaterGAP2 and from the regional oceanographic model NEMO. The inter-annual time series have a correlation and RMS differences of 0.68 and 55 mm, their annual amplitude and phase agree within 4 mm and 6 days respectively. The GRACE-derived seawater mass signal has an annual amplitude of 32 ± 4 mm peaking in April. On inter-annual time scales, the mass-induced sea level variability is stronger than in the Mediterranean Sea, with an increase from 2003 to 2005 followed by a decrease from 2006 to 2008.Based on mass conservation, the mass-induced sea level variations, river runoff and precipitation minus evaporation are combined to derive the strait flows between the basins and with the Atlantic Ocean. At the Gibraltar strait, the net inflow varies annually with an amplitude of 52 ± 10 × 10⁻³ Sv peaking end of September (1 Sv = 10⁶ m³ s⁻¹). The inflow through the Bosphorus strait displays an annual amplitude of 13 ± 3 ×10⁻³ Sv peaking in the middle of March. Additionally, an increase of the Gibraltar net inflow (3.4 ± 0.8 × 10⁻³ Sv/yr) is detected.     

Ca. 13 Ma strike-slip deformation in coastal Sonora from a large-scale,en-echelon,brittle-ductile,dextral shear indicator: implications for the evolution of the California rift

The Rancho Nuevo semi-circular structure is a geomorphological structure defined by drainage patterns in coastal Sonora, about 160 km NW of Hermosillo. The structure is about 15 by 30 km, and it is cored by felsic to intermediate plutons (granodiorite, monzogranite, quartz-porphyry) covered by Miocene volcanic rocks. This work is focused on the deformation of the intrusives which cover most of study area.The plutons are a co-magmatic suite dated between 71 ± 1.1 and 67.9 ± 1.0 Ma (U-Pb zircon, LA-ICPMS). The most voluminous unit is a granodiorite characterized by conspicuous sigmoidal fractures at the scale of high resolution satellite images, along which rhyolite dikes were emplaced about 13.2 Ma. Magnetic fabric (AMS) and paleomagnetic data were collected from 27 sites in the granodiorite. Magnetic fabrics are weak but well developed, and are characterized by steep foliation planes with strikes that follow the sigmoidal fracture pattern and suggest NE-SW to NW-SE flattening after emplacement. The characteristic magnetization is of dual polarity, but it is dominantly reverse consistent with emplacement during chron C31r. The prevalent magnetization is southwest and moderately steep negative (ten sites), a discordant direction rotated clockwise about 41 ° ± 11 with respect to the expected Late Cretaceous reference direction, also indicating gentle southward tilt. There is, however, paleomagnetic evidence suggesting that the structure did not rotate as a rigid body, but it deformed internally instead. These data are interpreted to indicate that the Rancho Nuevo semicircular structure is a large-scale, dextral, brittle-ductile shear indicator. The age of the dikes and the fact that they are covered discordantly by rocks assigned to the tuff of San Felipe indicate that northwest, srike-slip, motion of Baja California peninsula (and thus the Pacific plate relative to North America) was accommodated by faults in coastal Sonora about 13 Ma ago.     

Hydrological effects on gravity and correlations between gravitational variations and level of the Alzette River at the station of Walferdange,Luxembourg

The gravitational effects of water storage variations driven by local precipitation events are modeled for the Walferdange Underground Laboratory for Geodynamics (WULG) in the Grand Duchy of Luxembourg. A modified mass continuity model is implemented, which uses rainfall data from Walferdange as input. In the absence of soil moisture and groundwater level information, the model is empirically parameterized. Model outputs are compared to the gravity time series registered with the Observatory Superconducting Gravimeter CT040 located in the WULG. We find that the model explains 77% of the gravity residuals. In addition, a statistical analysis is carried out to determine the relationship between precipitation, gravity variations and water level changes in the nearby Alzette River. A time delay of 88 ± 34 min between the maximum variation rates of the water level and gravity signal has been calculated. The signals have an admittance of 45 ± 5 cm μGal^?1.     

Nighttime E region plasma irregularities over an equatorial station Trivandrum

In situ measurements of electron density were made over Trivandrum (8.5°N, 76.9°E) during nighttime to study E-region plasma density irregularities. Irregularities, with vertical scale sizes from a few km to 15 cm, were detected during rocket ascent and descent. Electron density profiles during ascent and descent of an earlier nighttime rocket flight from Trivandrum are also presented. Some of the important results are as follows: (i) horizontal gradients in electron density exist in 110–120 km region with horizontal scale size of at least 40 km, (ii) based on the presence/absence of electron density structures during ascent and descent of both flights, the horizontal distance over which the gradient drift instability operates is found to be at least 80 km and 90 km, for both the flights, (iii) observed irregularities in regions of negative density gradient are suggested to be produced through the gradient drift instability (GDI) driven by vertical polarization electric field as well as by electric field produced through wind shears and those in positive gradient regions by wind driven GDI, (iv) largest irregularity amplitude (≈30%) was associated with steepest gradients and so was the presence of smallest vertical scale sizes (12 m to 15 cm), which were absent at other altitudes, (v) the spectral index of irregularities was in the range of ?2.2±0.2 for large scales (few kilometers>λ>50 m), ?3.25±0.25 for medium scales (50 m>λ>10 m) and ?2.6±0.1 for smaller scales (10 m>λ>1 m) and (vi) irregularities in large and medium scales are expected to be produced directly through GDI and the small and sum-meter scales through non-linear GDI.     

40Ar / 39Ar geochronology constraints on late miocene weathering rates in Minas Gerais,Brazil

⁴⁰Ar / ³⁹Ar incremental heating ages for twenty one grains of cryptomelane, collected at 0, 42, 45, and 60 m depths in the Cachoeira Mine weathering profile, Minas Gerais, permit calculating long-term (10 Ma time scale) weathering rate (saprolitization rate) in SE Brazil. Pure well-crystallized cryptomelane grains with high K contents (3–5 wt.%) yield reliable geochronological results. The ⁴⁰Ar / ³⁹Ar plateau ages obtained decrease from the top to the bottom of the profile (12.7 ± 0.1 to 7.6 ± 0.1 Ma at surface; 7.6  ± 0.2 to 6.1 ± 0.2 Ma at 42 m; and 7.1 ± 0.2 to 5.9 ± 0.1 Ma at 45 m; 6.6 ± 0.1 to 5.2 ± 0.1 Ma at 60 m), yielding a weathering front propagation rate of 8.9 ± 1.1 m/m.y. From the geochronological results and the mineral transformations implicit by the current mineralogy in the weathering profiles, it is possible to calculate the saprolitization rate for the Cachoeira Mine lithologies and for adjacent weathering profiles developed on granodiorites and schists. The measured weathering front propagation rate yields a saprolitization rate of 24.9 ± 3.1 t/km²/yr. This average long-term (> 10 Ma) saprolitization rate is consistent with mass balance calculations results for present saprolitization rates in weathering watersheds. These results are also consistent with long-term saprolitization rates estimated by combining cosmogenic isotope denudation rates with mass balance calculations.