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.     

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.     

Geological controls on submarine groundwater discharge in Long Bay,South Carolina (USA)

A combination of geophysical methods including continuous electrical resistivity and high-resolution Chirp sub-bottom profiling were utilized to characterize geologic controls on pore fluid salinity in the nearshore of Long Bay, SC. Resistivity values ranged from less than 1 Ω m to greater than 40 Ω m throughout the bay. Areas of elevated electrical resistivity suggest the influence of relatively fresher water on pore water composition. Geophysical evidence alone does not eliminate all ambiguity associated with lithological and porosity variations that may also contribute to electrical structure of shallow marine sediments. The anomalous field is of sufficient magnitude that lithological variation alone does not control the spatial distribution of elevated electrical resistivity zones. Geographical distribution of electrical anomalies and structures interpreted from nearby sub-bottom profiles indicates abrupt changes in shallow geologic units control preferential pathways for discharge of fresh water into the marine environment. Shore parallel resistivity profiles show dramatic decreases in magnitude with increasing distance from shore, suggesting a significant portion of the terrestrially driven fresh SGD in Long Bay is occurring via the surficial aquifer within a few hundred meters of shore.     

Resistivity imaging of Pleistocene alluvial aquifers in a contractional tectonic setting: A case history from the Po plain (Northern Italy)

In this work we present the hydrogeophysical imaging of a key sector of the Quaternary Po foreland basin (northern Italy), focussing on the reconstruction of clastic aquifers and aquitards in a complex tectono-sedimentary subsurface architecture. The study area includes the relic reliefs of Casalpusterlengo and Zorlesco, two smooth morphological features involving uplifted and gently folded Pleistocene marine to alluvial sediments, plausibly linked to the buried Northern Apennines thrust and fold belt. The geophysical data include 35 Direct Current Vertical Electrical Soundings collected over a 37 km² wide area, acquired with Schlumberger array and maximum half-spacing of 500 m. 1-D resistivity-depth profiles were computed for each VES. An integrated hydrostratigraphic approach was applied, to constrain the interpretation of the geophysical data along several cross-sections, including the comparison of resistivity soundings to stratigraphic logs, borehole electric logs and the pore-water properties.The resistivity interfaces, traceable with the same laterally continuous vertical polarity, were used to develop an electrostratigraphic model in order to portray the stacking of electrostratigraphic units down to 200 m below ground surface. Their vertical associations show a general upward increase of electrical resistivity. This assemblage mimics the regional coarsening upwards depositional trend, from the conductive units of the Plio-Pleistocene marine-to-transitional depositional systems to the resistive units of the Middle–Late Pleistocene fluvial and alluvial plain depositional systems. Middle Pleistocene depositional systems host an alternation of North-dipping, high-to-intermediate permeability aquifer systems (70–180 Ωm, thickness of 5–70 m) separated by low permeability aquitards (20–50 Ωm, thickness up to 40 m). These units pinch out against the Casalpusterlengo and Zorlesco relic reliefs, where they cover the uplifted and folded regional aquitard (20–50 Ωm) formed by Pliocene-Lower Pleistocene clays to sandy silts with gravel lenses in agreement with borehole data. In the deepest part of the local stratigraphy, a broad low-resistivity anomaly (< 10 Ωm) was clearly mapped through the study area. By comparison with electrical borehole logs in deep oil-wells, it could be interpreted as the fresh–saltwater interface due to the presence of connate waters and brines hosted by the marine-to-transitional shales.     

Resistivity monitoring of the tephra barrier at Crater Lake,Mount Ruapehu,New Zealand

The eruptions of Mt Ruapehu in the North Island of New Zealand in 1995 and 1996 caused a tephra barrier to be formed across the outlet of Crater Lake. By 2005 seepage from the refilled lake into the barrier raised the possibility of an eventual collapse of the barrier, releasing a catastrophic lahar down the mountain.As part of an extensive monitoring programme of the tephra barrier, direct current (dc) resistivity surveys were carried out on a number of lines along and across it in order to test whether the extent of the seepage could be measured (and monitored) by geophysical means. Two dimensional inversion of measured apparent resistivity data showed that between the initial measurements, made in January 2005, and February 2006, there was a gradual decrease in resistivity above the old outlet from ~ 50–60 Ωm to ~ 30 Ωm. This gave the first indication that lake water was seeping into the barrier. Between October and December 2006 there was a rapid rise in lake level to only 2 m below the top of the barrier, and a further resistivity survey in January 2007 showed that there had been a further decrease in resistivity throughout the entire barrier with values dropping to < 10 Ωm. The extent of this low resistivity indicated that the barrier was now saturated. At this stage lake water was penetrating the barrier and starting to cause erosion on its downstream side. Catastrophic collapse occurred on 18 March 2007, accompanied by a lahar in the Whangaehu river valley.Subsequent forward 3D numerical modelling of the resistivity structure of the barrier has confirmed that the observed changes in measured resistivity were directly related to the progress of seepage of lake water into the barrier.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

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).     

Three-dimensional electrical resistivity model of a nuclear waste disposal site

A three-dimensional (3D) modeling study was completed on a very large electrical resistivity survey conducted at a nuclear waste site in eastern Washington. The acquisition included 47 pole–pole two-dimensional (2D) resistivity profiles collected along parallel and orthogonal lines over an area of 850 m × 570 m. The data were geo-referenced and inverted using EarthImager3D (EI3D). EI3D runs on a Microsoft 32-bit operating system (e.g. WIN-2K, XP) with a maximum usable memory of 2 GB. The memory limits the size of the domain for the inversion model to 200 m × 200 m, based on the survey electrode density. Therefore, a series of increasing overlapping models were run to evaluate the effectiveness of dividing the survey area into smaller subdomains. The results of the smaller subdomains were compared to the inversion results of a single domain over a larger area using an upgraded form of EI3D that incorporates multi-processing capabilities and 32 GB of RAM memory. The contours from the smaller subdomains showed discontinuity at the boundaries between the adjacent models, which do not match the hydrogeologic expectations given the nature of disposal at the site. At several boundaries, the contours of the low resistivity areas close, leaving the appearance of disconnected plumes or open contours at boundaries are not met with a continuance of the low resistivity plume into the adjacent subdomain. The model results of the single large domain show a continuous monolithic plume within the central and western portion of the site, directly beneath the elongated trenches. It is recommended that where possible, the domain not be subdivided, but instead include as much of the domain as possible given the memory of available computing resources.     

Combined three-dimensional electric and seismic tomography study on the Åknes rockslide in western Norway

We present a combined 3-D geoelectric and seismic tomography study conducted on the large Åknes rockslide in western Norway. Movements on the slope are strongly influenced by water infiltration, such that the hydrogeological regime is considered as a critical factor affecting the slope stability. The aim of our combined geophysical study was to identify and visualize the main shallow tension fractures and to determine their effect on hydraulic processes by comparing the geophysical results with information from borehole logging and tracer tests. To resolve the complex subsurface conditions of the highly fractured rock mass, a three-dimensional set-up was chosen for our seismic survey. To map the water distribution within the rock mass, a pattern of nine intersecting 2-D geoelectric profiles covered the complete unstable slope. Six of them that crossed the seismic survey area were considered as a single data set in a 3-D inversion. For both methods, smoothing-constraint inversion algorithms were used, and the forward calculations and parameterizations were based on unstructured triangular meshes. A pair of parallel shallow low-velocity anomalies (< 1400 m/s) observed in the final seismic tomogram was immediately underlain by two anomalies with resistivities <13 kΩm in the resistivity tomogram. In combination with borehole logging results, the low-velocity and resistivity anomalies could be associated with the drained and water-filled part of the tension fractures, respectively. There were indications from impeller flowmeter measurements and tracer tests that such tension fractures intersected several other water-filled fractures and were responsible for distinct changes of the main groundwater flow paths.     

Cenozoic building and deformational processes in the North Patagonian Andes

The Oligocene to present evolution of the North Patagonian Andes is analyzed linking geological and geophysical data in order to decipher the deformational processes that acted through time and relate them to basin formation processes. Seismic reflection profiles reveal the shallow structure of the retroarc area where contractional structures, associated with Oligocene to early Miocene inverted extensional depocenters, are partially onlapped by early to late Miocene synorogenic deposits. From the construction of five structural cross sections along the retroarc area between 40° and 43°30′ S, constrained by surface, gravity and seismic data, a shortening gradient is observed along Andean strike. The highest shortening of 18.7 km (15.34%) is determined near 41°30′ S coincidentally with maximum mean topographic values on the eastern Andean slope, where basement blocks were uplifted in the orogenic front area, and the deepest and broadest synorogenic depocenters were formed towards the foreland. Additionally, eastward shifting of Miocene calc-alkaline rocks occurred at these latitudes, which is interpreted as indicative of a change in the subduction parameters at this time. Deep crustal retroarc structure is evaluated through inversion of gravity models that made possible to infer Moho attenuated zones. These coincide with the occurrence of younger than 5 Ma within-plate volcanics as well as with crustal thermal anomalies suggested by shallowing of the Curie isotherm calculated from magnetic data. Younger volcanism and thermal anomalies are explained by slab steepening since early Pliocene, after a mild-shallow subduction setting in the middle to late Miocene, age of the main compressive event.     

Characterization of aquifers in the Vientiane Basin,Laos, using Magnetic Resonance Sounding and Vertical Electrical Sounding

The aim of this study is to define and characterize water bearing geological formation and to test the possibility of using geophysical techniques to determine the hydrogeological parameters in three areas in the Vientiane basin, Laos. The investigated areas are part of the Khorat Plateau where halite is naturally occurring at depths as shallow as 50 m in the Thangon Formation. Magnetic Resonance Sounding (MRS) has been used in combination with Vertical Electrical Sounding (VES) in different geological environments. In total, 46 sites have been investigated and the MRS and VES recognized the stratigraphic unit N₂Q_1–3, consisting of alluvial unconsolidated sediments, as the main water bearing unit. The aquifer thickness varies usually between 10 and 40 m and the depth to the main aquifer range from 5 to 15 m. The free water content is here up to 30%, and the decay times vary between 100 and 400 ms, suggesting a mean pore size equivalent to fine sand to gravel. The resistivity is highly variable, but usually around 10–1500 Ω-m, except for some sites in areas 1 and 2, where the aquifer is of low resistivity, probably related to salt water. Hydraulic and storage-related parameters such as transmissivity, hydraulic column, have been estimated from the MRS. The MRS together with VES has been shown to be a useful and important tool for identifying and distinguishing freshwater from possible salt-affected water as well as the salt-related clay layer of the Thangon Formation. This clay layer is characterized by very low free water content and a resistivity lower than 5 Ω-m and can be found in all 3 areas at depths from 15 to 50 m.     

Role of Lithology and Subsurface structures detected by potential field data in controlling the radioactive mineral accumulation at Natash area,Eastern Desert,Egypt

Wadi Natash area is located in the southern part of the Eastern desert of Egypt. It has a great importance for containing accumulations from the radioactive minerals of Uranium, Thorium and Potassium. An integrated potential study was carried out on the study area with the aim of locating depths to causative bodies with sufficient magnetic susceptibility that may represent magmatic intrusions with relation to the radioactivity location and delineate the subsurface structures affecting the area. Both magnetic and Bouguer data as well as radiometric data were interpreted rapidly for source positions and depths using Euler deconvolution, Werner deconvolution and 3D modeling techniques. The results deduced from the trend analyses show that the major fault trend affecting the area have NNW-SSE (Red Sea–Gulf of Suez trend) direction intersected by the less predominant NNE-SSW(The Gulf of Aqaba–Dead Sea trend) and WNW-ESE (Najd Fault System) fault trends. The causative bodies were imaged at depths ranging from 0.3 km to about 1.5 km. The depths along the interpreted profiles display discontinuities in potential field markers due to presence of the NNW-SSE fault trends act as pass channels for the hydrothermal solutions.It can be stated that the radioactive mineral accumulations were caused by the hydrothermal solutions rich with radioactive minerals as a result of intruding Natash volcanic to the granitic rocks. The Qouseir clastics and the Nudian sand stone were affected by these solutions and show a positive response for the radioactive minerals.     

Electrical imaging of Narmada–Son Lineament Zone,Central India from magnetotellurics

The Narmada–Son Lineament (NSL) Zone is the second most important tectonic feature after Himalayas, in the Indian geology. Magnetotelluric (MT) studies were carried out in the NSL zone along a 130 km long NNE-SSW trending profile. The area of investigation extends from Edlabad (20°46′16″; 75°59′05″) in the South to Khandwa (21°53′51″; 76°18′05″) in the North. The data shows in general the validity of a two-dimensional (2D) approach. Besides providing details on the shallow crustal section, the 2D modeling results resolved four high conductive zones extending from the middle to deep crust, spatially coinciding with the major structural features in the area namely the Gavligarh, Tapti, Barwani-Sukta and Narmada South faults. The model for the shallow section has brought out a moderately resistive layer (30–150 Ω m) representing the exposed Deccan trap layer, overlying a conductive layer (10–30 Ω m) inferred to be the subtrappean Gondwana sediments, the latter resting on a high resistive basement/upper crust. The Deccan trap thickness varies from around a few hundred meters to as much as 1.5 km along the traverse. A subtrappean sedimentary basin like feature is delineated in the northern half of the traverse where a sudden thickening of subtrappean sediments amounting to as much as 2 km is noticed. The high resistive upper crust is relatively thick towards the southern end and tends to become thinner towards the middle and northern part of the traverse. The lower crustal segment is conductive over a major part of the profile. Considering the generally enhanced heat flow values in the NSL region, coupled with characteristic gravity highs and enhanced seismic velocities coinciding with the mid to lower crustal conductors delineated from MT, presence of zones of high density mafic bodies/intrusives with fluids, presumably associated with magmatic underplating of the crust in the zone of major tectonic faults in NSL region are inferred.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.     

Systematic deflection and offset of the Yangtze River drainage system along the strike-slip Ganzi-Yushu-Xianshuihe Fault Zone,Tibetan Plateau

The Ganzi-Yushu-Xianshuihe Fault Zone (GYXFZ) is a typical active strike-slip fault that has triggered many large historic earthquakes, including the 2010 M_w 6.9 Yushu earthquake in the central Tibetan Plateau. This fault zone extends for ca. 800 km from the central Tibetan Plateau to its southeastern margin and varies in trend from WNW-ESE in the northwestern segment of the fault zone to NNW-SSE in the southeastern segment, having the geometry of an arc projecting northeastwards. In this study, we present evidence for the systematical sinistral deflection and/or offset of the Yangtze River and its branch stream channels and valleys along the GYXFZ. Topographic analysis of three-dimensional (3D) perspective images constructed using digital elevation model (DEM) data, 0.5 m-resolution WorldView and GeoEye images, and 15 m-resolution Landsat-Enhanced Thematic Mapper (ETM+) images, together with analysis of geological structures, reveals the following: (i) the main river channels and valleys of the Yangtze River drainage system show systematic sinistral deflections and/or offsets along the GYXFZ; (ii) various amounts of sinistral offset have accumulated on the tributary stream channels, valleys, and gullies of the Yangtze River along the fault, with a linear relation, D = aL, between the upstream length L from the deflected point and the offset amount D with a certain coefficient a; (iii) the maximum amount of sinistral offset is up to ca. 60 km, which was accumulated in the past 13–5 Ma; and (iv) the long-term average strike-slip rate is ca. 4.6–12 mm/year. Geological and geomorphic evidence, combined with geophysical data, demonstrates that the GYXFZ is currently active as one of the major seismogenic faults in the Tibetan Plateau, dominated by left-lateral strike-slip motion. Our findings supply important evidence for the tectonic evolution of strike-slip faults in the Tibetan Plateau since the Eurasia-India continental collision.     

Shallow fault segmentation of the Alpine fault zone,New Zealand revealed from 2- and 3-D GPR surveying

Where they are preserved, landforms that have been truncated and offset by past fault movements provide potentially valuable quantitative data that can be used to estimate slip rates. At such locations, it is important to investigate the fault zone in sufficient detail to understand how displacements are accommodated on individual fault strands. At a site along a northern section of the Alpine fault zone on the South Island of New Zealand, surface mapping of a series of faulted river terraces and channels has revealed a complicated and poorly understood paleoearthquake history. We have acquired high-resolution 2- and 3-D ground-penetrating radar (GPR) data over a large area (~ 500 × 500 m) of the terraces to map along-strike changes in shallow (<20 m) fault zone morphology. By identifying distinct reflection patterns within the topographically migrated 3-D GPR volumes and extrapolating them to the longer and more widely spaced GPR profiles, we determined the subsurface extent of two main structural/depositional facies that were juxtaposed by three left-stepping en-echelon fault strands. Two regions of warped strata are interpreted to result from transpressive folding between the overlapping strands, where displacement is transferred from one fault to the next. We suggest that diffuse deformation between the overlapping fault tips results in anomalously low estimates for horizontal and vertical fault displacements of some geomorphic features.     

Structure and development of the Andean system between 36° and 39°S

One of the main morphological changes along the Southern Central Andes occurs from 36° to 39°S. The northern portion is characterized by prominent basement structures and a thick-skinned orogenic front with relief of over 2000 m with a deep level of exhumation where more than 4 km of section has been eroded. Contrastingly, the southern part is formed by mildly inverted basement structures restricted mainly to the hinterland zone, which reaches only 1500–1700 m relief. We quantify the variable contributions of two main contractional stages through the construction of three regionally balanced sections across the Andes, constrained by field and geophysical data. Extensional re-activation described for this segment in late Oligocene-early Miocene and Pliocene to Quaternary times, after the two main contractional episodes, suggests only 3 km of stretching that represents 30–10% of the original longitude. We, therefore, conclude that while initial Late Cretaceous to Eocene compression was similar along strike (∼10–7 km), it is the contrasting degrees of Neogene shortening (∼16–6 km) that have played the largest role in the along strike differences in structure and morphology along this portion of the southern Andes. Variable Neogene arc expansion could be responsible for the contrasting contractional deformation: In the north, late Miocene arc-related rocks cover most of the retroarc zone (>200 km with respect to the late Miocene arc front in the south), presumably driven by a shallow subduction episode in the area, whereas to the south they remain restricted to the continental drainage divide. Other factors involving architecture of previous rift structures, are proposed as additional mechanisms that accommodated variable shortening magnitudes through inversion.     

Crustal structure variation along 30°N in the eastern Tibetan Plateau and its tectonic implications

We determined crustal structure along the latitude 30°N through the eastern Tibetan Plateau using a teleseismic receiver function analysis. The data came mostly from seismic stations deployed in eastern Tibet and western Sichuan region from 2004 to 2006. Crustal thickness and Vp/Vs ratio at each station were estimated by the H–k stacking method. On the profile, the mean crustal thickness and Vp/Vs ratio were found to be 62.3 km and 1.74 in the Lhasa block, 71.2 km and 1.79 near the Bangong–Nujiang suture, 66.3 km and 1.80 in the Qiangtang block, 59.8 km and 1.81 in the Songpan–Garze block, and 42.9 km and 1.76 in the Yangtze block, respectively. The estimated crustal thicknesses are consistent with predictions based on the topography and the Airy isostasy, except near the Bangong–Nujiang suture and in the Qiangtang block where the crust is 5–10 km thicker than predicted, indicating that the crust may be denser, possibly due to mafic underplating. We also inverted receiver functions for crustal velocity structure along the profile, which reveals a low S-wave velocity zone in the lower crust beneath the eastern Tibetan Plateau, although the extent of the low-velocity zone varies considerably. The low-velocity zone, together with previous results, suggests limited partial melting and localized crustal flow in the lower crust of the eastern Tibetan Plateau.     

Spatiotemporal patterns of fault slip rates across the Central Sierra Nevada frontal fault zone

Patterns in fault slip rates through time and space are examined across the transition from the Sierra Nevada to the Eastern California Shear Zone–Walker Lane belt. At each of four sites along the eastern Sierra Nevada frontal fault zone between 38 and 39° N latitude, geomorphic markers, such as glacial moraines and outwash terraces, are displaced by a suite of range-front normal faults. Using geomorphic mapping, surveying, and ¹⁰Be surface exposure dating, mean fault slip rates are defined, and by utilizing markers of different ages (generally, ~ 20 ka and ~ 150 ka), rates through time and interactions among multiple faults are examined over 10⁴–10⁵ year timescales.At each site for which data are available for the last ~ 150 ky, mean slip rates across the Sierra Nevada frontal fault zone have probably not varied by more than a factor of two over time spans equal to half of the total time interval (~ 20 ky and ~ 150 ky timescales): 0.3 ± 0.1 mm year^? 1 (mode and 95% CI) at both Buckeye Creek in the Bridgeport basin and Sonora Junction; and 0.4 + 0.3/?0.1 mm year^? 1 along the West Fork of the Carson River at Woodfords. Data permit rates that are relatively constant over the time scales examined. In contrast, slip rates are highly variable in space over the last ~ 20 ky. Slip rates decrease by a factor of 3–5 northward over a distance of ~ 20 km between the northern Mono Basin (1.3 + 0.6/?0.3 mm year^? 1 at Lundy Canyon site) to the Bridgeport Basin (0.3 ± 0.1 mm year^? 1). The 3-fold decrease in the slip rate on the Sierra Nevada frontal fault zone northward from Mono Basin is indicative of a change in the character of faulting north of the Mina Deflection as extension is transferred eastward onto normal faults between the Sierra Nevada and Walker Lane belt.A compilation of regional deformation rates reveals that the spatial pattern of extension rates changes along strike of the Eastern California Shear Zone-Walker Lane belt. South of the Mina Deflection, extension is accommodated within a diffuse zone of normal and oblique faults, with extension rates increasing northward on the Fish Lake Valley fault. Where faults of the Eastern California Shear Zone terminate northward into the Mina Deflection, extension rates increase northward along the Sierra Nevada frontal fault zone to ~ 0.7 mm year^? 1 in northern Mono Basin. This spatial pattern suggests that extension is transferred from more easterly fault systems, e.g., Fish Lake Valley fault, and localized on the Sierra Nevada frontal fault zone as the Eastern California Shear Zone–Walker Lane belt faulting is transferred through the Mina Deflection.     

ELF transients associated with sprites and elves in eastern Mediterranean winter thunderstorms

During the northern hemisphere winter of 2005–2006, transient luminous events (TLEs) known as ‘sprites’ and ‘elves’ were imaged over thunderstorm cells in the eastern Mediterranean. Simultaneously, extremely low frequency (ELF) data (ELF: 3–3000 Hz) were recorded at two observation stations in Israel and Hungary in order to qualify and quantify parameters of the parent lightning discharge associated with the transient optical emissions in the upper atmosphere. In this study, we found that for 87% (Israel) and 77% (Hungary) of optically observed TLEs an intense ELF transient event was recorded. These stations are located some 500 and 2100 km, respectively, from the region of the TLEs. All ELF transients that were associated with TLEs were caused by lightning discharges with positive polarity. Calculation of the charge moment change showed values between 600 and 2800 C km with a peak around 1000 C km. Additionally, the time delay between the +CG and ensuing sprite was 76±34 ms and it was displaced up to 50 km from its parent CG.One of our objectives in the present study was to characterize, based on the ELF radiation from lightning, the electromagnetic (EM) waveforms of the lightning discharges which generate TLEs in the time and frequency domains, and to compare them with other lightning discharges occurring in the same thunderstorm cell at approximately the same time, but which did not produce TLEs. The survey for a typical EM waveform showed no unique ELF signature for lightning discharges associated with either sprites or elves.