Current states of stress in the northern Andes as indicated by focal mechanisms of earthquakes

Systematic inversion of double couple focal mechanisms of shallow earthquakes in the northern Andes reveals relatively homogeneous patterns of crustal stress in three main regions. The first region, presently under the influence of the Caribbean plate, includes the northern segment of the Eastern Cordillera of Colombia and the western flank of the Central Cordillera (north of 4°N). It is characterized by WNW–ESE compression of dominantly reverse type that deflects to NW–SE in the Merida Andes of Venezuela, where it becomes mainly strike–slip in type. A major bend of the Eastern thrust front of the Eastern Cordillera, near its junction with the Merida Andes, coincides with a local deflection of the stress regime (SW–NE compression), suggesting local accommodation of the thrust belt to a rigid indenter in this area. The second region includes the SW Pacific coast of Colombia and Ecuador, currently under the influence of the Nazca plate. In this area, approximately E–W compression is mainly reverse in type. It deflects to WSW–ENE in the northern Andes south of 4°N, where it is accommodated by right-lateral displacement of the Romeral fault complex and the Eastern front of the northern Andes. The third, and most complex, region is the area of the triple junction between the South American, Nazca and Caribbean plates. It reveals two major stress regimes, both mainly strike–slip in type. The first regime involves SW–NE compression related to the interaction between the Nazca and Caribbean plates and the Panama micro-plate, typically accommodated in an E–W left-lateral shear zone. The second regime involves NW–SE compression, mainly related to the interaction between the Caribbean plate and the North Andes block which induces left-lateral displacement on the Uramita and Romeral faults north of 4°N.Deep seismicity (about 150–170 km) concentrates in the Bucaramanga nest and Cauca Valley areas. The inversion reveals a rather homogeneous attitude of the minimum stress axis, which dips towards the E. This extension is consistent with the present plunge of the Nazca and Caribbean slabs, suggesting that a broken slab may be torn under gravitational stresses in the Bucaramanga nest. This model is compatible with current blocking of the subduction in the western northern Andes, inhibiting the eastward displacement of slabs, which are forced to break and sink in to the asthenosphere under their own weight.     

Analysis of the interannual variability of annual daily extreme water levels in the St Lawrence River and Lake Ontario from 1918 to 2010

We compared the interannual variability of annual daily maximum and minimum extreme water levels in Lake Ontario and the St Lawrence River (Sorel station) from 1918 to 2010, using several statistical tests. The interannual variability of annual daily maximum extreme water levels in Lake Ontario is characterized by a positive long‐term trend showing two shifts in mean (1929–1930 and 1942–1943) and a single shift in variance (in 1958–1959). In contrast, for the St Lawrence River, this interannual variability is characterized by a negative long‐term trend with a single shift in mean, which occurred in 1955–1956. As for annual daily minimum extreme water levels, their interannual variability shows no significant long‐term change in trend. However, for Lake Ontario, the interannual variability of these water levels shows two shifts in mean, which are synchronous with those for maximum water levels, and a single shift in variance, which occurred in 1965–1966. These changes in trend and stationarity (mean and variance) are thought to be due to factors both climatic (the Great Drought of the 1930s) and human (digging of the Seaway and construction of several dams and locks during the 1950s). Despite this change in means and variance, the four series are clearly described by the generalized extreme value distribution. Finally, annual daily maximum and minimum extreme water levels in the St Lawrence and Lake Ontario are negatively correlated with Atlantic multidecadal oscillation over the period from 1918 to 2010. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Unraveling the Peruvian Phase of the Central Andes: stratigraphy,sedimentology and geochronology of the Salar de Atacama Basin (22°30–23°S), northern Chile

The Salar de Atacama Basin holds important information regarding the tectonic activity, sedimentary environments and their variations in northern Chile during Cretaceous times. About 4000 m of high‐resolution stratigraphic columns of the Tonel, Purilactis and Barros Arana Formations reveal braided fluvial and alluvial facies, typical of arid to semi‐arid environments, interrupted by scarce intervals with evaporitic, aeolian and lacustrine sedimentation, displaying an overall coarsening‐upward trend. Clast‐count and point‐count data evidence the progressive erosion from Mesozoic volcanic rocks to Palaeozoic basement granitoids and deposits located around the Cordillera de Domeyko area, which is indicative of an unroofing process. The palaeocurrent data show that the source area was located to the west. The U/Pb detrital zircon geochronological data give maximum depositional ages of 149 Ma for the base of the Tonel Formation (Agua Salada Member), and 107 Ma for its middle member (La Escalera Member); 79 Ma for the lower Purilactis Formation (Limón Verde Member), and 73 Ma for the Barros Arana Formation. The sources of these zircons were located mainly to the west, and comprised from the Coastal Cordillera to the Precordillera. The ages and pulses record the tectonic activity during the Peruvian Phase, which can be split into two large events; an early phase, around 107 Ma, showing uplift of the Coastal Cordillera area, and a late phase around 79 Ma indicating an eastward jump of the deformation front to the Cordillera de Domeyko area. The lack of internal deformation and the thicknesses measured suggest that deposition of the units occurred in the foredeep zone of an eastward‐verging basin. This sedimentation would have ended with the K‐T phase, recognized in most of northern Chile.     

Provenance of buried esker groundwater: the case of Vars-Winchester esker aquifer,Eastern Ontario,Canada

An innovative mode of groundwater recharge to a buried esker aquifer is considered. The current conceptual model affords a natural safeguard to underlying aquifers from the overlying muds. A hypothesis of groundwater recharge to a buried esker aquifer via preferential pathways across its overlying muds is tested here by heuristic numerical one-dimensional and two-dimensional modeling simulations. The hypothesis has been tested against two other conventionally accepted scenarios involving: (1) distal esker outcrop areas and (2) remote shallow-bedrock recharge areas. The main evidence comes from documented recharge pressure pulses in the overlying mud aquitard and in the underlying esker hydraulic-head time series for the Vars-Winchester esker aquifer in Eastern Ontario, Canada. These perturbations to the potentiometric surface are believed to be the aquifer response to recharge events. The migration rate of these pressure pulses is directly related to the hydraulic diffusivity of the formation. The measured response time and response amplitude between singular radar precipitation events and well hydrographs constituted the heuristic model calibration targets. The main evidence also includes mud-layering deformation (water escape features) which was observed in seismic surveys of the over-esker muds. These disturbed stratigraphic elements provide a realistic mechanism for migrating water to transit through the muds. The effective hydraulic conductivities of these preferential pathways in the muds were estimated to be between 2?×?10^?6 and 7?×?10^?6 m/s. The implications of these findings relate to the alleged natural safeguard of these overlying muds.     

Ba同位素揭示上地幔中广泛存在再循环的沉积物

正Ba同位素由于具有在低温环境中分馏的特性,相比放射性同位素与微量元素,可以更好地示踪地球深部动力学过程中的再循环物质。研究表明沉积物与蚀变洋壳具有明显不同的Ba同位素组成(沉积物富集来自重晶石的Ba),使得Ba同位素成为研究表层与地幔储库间物质交换的有力工具。最近美国伍兹霍尔海洋研究所报道了来自全球洋中脊的21个MORB样品的Ba、Sr、Nd同位素和微量元素组成以及两个钻井中沉积物和蚀变洋壳     

苏鲁超高压变质岩中地震波速随围压的变化规律

我们在10~800MPa的静水压力下测量了中国大陆科学钻探(CCSD)主孔岩心和采自苏鲁地区地表露头上的共68块典型超高压变质岩试样的地震波速,并建立了波速和围压的定量关系,成功地解释了描述这一定量关系方程中4个参数的物理意义。笔者相信,本文将为超高压变质岩地区以及新老俯冲带地震波资料的解释提供必不可少的理论和实验支撑。     

Access routes to high mountain huts facing climate-induced environmental changes and adaptive strategies in the Western Alps since the 1990s

ABSTRACT

In the European Alps, high mountain environments are subject to major impacts resulting from climate change, which strongly affect human activities such as mountaineering. The purpose of the study was to examine changes in access routes to 30 high mountain huts in the Western Alps since the 1990s. Data were derived from the use of two different methods, geo-historical studies and a questionnaire, and were used to identify both the climate-related processes affecting the climbing routes and the strategies implemented by public entities, Alpine clubs, guide companies, and hut keepers to maintain acceptable safety and technical conditions. The case studies revealed issues affecting three access routes and the results from the questionnaire showed that the main processes affecting access routes were loss of ice thickness and retreat from the front of the glaciated areas. Commonly, in situ equipment was installed to facilitate access for mountaineers and/or a part of a route was relocated to a safer area. The authors conclude that in most cases, the measures were effective but they were limited by financial, ethical and legal issues, especially in protected or classified areas that could jeopardise their durability and effectiveness.     

Characteristics of developing waves as a function of atmospheric conditions,water properties,fetch and duration

For any specific wind speed, waves grow in period, height and length as a function of the wind duration and fetch until maximum values are reached, at which point the waves are considered to be fully developed. Although equations and nomograms exist to predict the parameters of developing waves for shorter fetch or duration conditions at different wind speeds, these either do not incorporate important variables such as the air and water temperature, or do not consider the combined effect of fetch and duration. Here, the wind conditions required for a fully developed sea are calculated from maximum wave heights as determined from the wind speed, together with a published growth law based on the friction velocity. This allows the parameters of developing waves to be estimated for any combination of wind velocity, fetch and duration, while also taking account of atmospheric conditions and water properties.