Does continental crust transform during eclogite facies metamorphism?

Eclogites within exhumed continental collision zones indicate regional burial to depths of at least 60 km, and often more than 100 km in the coesite‐stable, ultra‐high pressure (UHP) eclogite facies. Garnet, omphacitic pyroxene, high‐Si mica, kyanite ± coesite should grow at the expense of low‐P minerals in most felsic compositions, if equilibrium obtained at these conditions. The quartzofeldspathic rocks that comprise the bulk of eclogite facies terranes, however, contain mainly amphibolite facies, plagioclase‐bearing assemblages. To what extent these lower‐P minerals persisted metastably during (U)HP metamorphism, or whether they grew afterwards, reflects closely upon crustal parameters such as density, strength and seismic character. The Nordfjord area in western Norway offers a detailed view into a large crustal section that was subducted into the eclogite facies. The degree of transformation in typical pelite, paragneiss, granitic and granodioritic gneiss was assessed by modelling the equilibrium assemblage, comparing it with existing parageneses in these rocks and using U/Th–Pb zircon geochronology from laser ablation ICPMS to establish the history of mineral growth. U–Pb dates define a period of zircon recrystallization and new growth accompanying burial and metamorphism lasting from 430 to 400 Ma. Eclogite facies mafic rock (~2 vol.% of crust) is the most transformed composition and records the ambient peak conditions. Rare garnet‐bearing pelitic rocks (<10 vol.% of crust) preserve a mostly prograde mineral evolution to near‐peak conditions; REE concentrations in zircon indicate that garnet was present after 425 Ma and feldspar broke down after 410 Ma. Felsic gneiss – by far the most abundant rock type – is dominated by quartz + biotite + feldspar, but minor zoisite/epidote, phengitic white mica, garnet and rutile point to a prograde HP overprint. Relict textures indicate that much of the microstructural framework of plagioclase, K‐feldspar, and perhaps biotite, persisted through at least 25 Ma of burial, and ultimately UHP metamorphism. The signature reaction of the eclogite facies in felsic rocks – jadeite/omphacite growth from plagioclase – cannot be deduced from the presence of pyroxene or its breakdown products. We conclude that prograde dehydration in orthogneiss leads to fluid absent conditions, impeding equilibration beyond ~high‐P amphibolite facies.     

Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater

Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.     

A hybrid approach to modeling territorial control in violent armed conflicts

Territorial control is central to the understanding of violent armed conflicts, yet reliable and valid measures of this concept do not exist. We argue that geospatial analysis provides an important perspective to measure the concept. In particular, measuring territorial control can be seen as an application of calculating service areas around points of control. The modeling challenge is acute for areas with limited road infrastructure, where no complete network is available to perform the analysis, and movements largely occur off road. We present a new geospatial approach that applies network analysis on a hybrid transportation network with both actual road data and hexagon‐fishnet‐based artificial road data representing on‐road and off‐road movements, respectively. Movement speed or restriction can be readily adjusted using various input data. Simulating off‐road movement with hexagon‐fishnet‐based artificial road data has a number of advantages including scalability to small or large study areas and flexibility to allow all‐directional travel. We apply this method to measuring territorial control of armed groups in Sub‐Saharan Africa where inferior transport infrastructure is the norm. Based on the Uppsala Conflict Data Program's (UCDP) Georeferenced Event Data (GED) as well as spatial data on terrain, population locations, and limited transportation networks, we enhance the delineation of the specific areas directly controlled by each warring party during civil wars within a given travel time.     

Controls on sediment evacuation from glacially modified and unmodified catchments in the eastern Sierra Nevada,California

The degree of glacial modification in small catchments along the eastern Sierra Nevada, California, controls the timing and pattern of sediment flux to the adjacent fans. There is a close relationship between the depth of fan‐head incision and the pattern and degree of Late Pleistocene catchment erosion by valley glaciers; catchments with significant glacial activity are associated with deeply incised fan heads, whereas fans emerging from glacially unmodified catchments are unincised. We suggest that the depth of fan‐head incision is controlled by the potential for sediment storage during relatively dry ice‐free periods, which in turn is related to the downstream length of the glacially modified valley and creation of accommodation through valley floor slope lowering and glacial valley overdeepening and widening. Significant storage in glacially modified basins during ice‐free periods leads to sediment supply‐limited conditions at the fan head and causes deep incision. In contrast, a lack of sediment trapping allows quasi‐continuous sediment supply to the fan and prevents incision of the fan head. Sediment evacuation rates should thus show large variations in glacially modified basins, with major peaks during glacial and lows during interglacial or ice‐free periods, respectively. In contrast, sediment removal from glacially unmodified catchments in this type of setting should be free of this effect, and will be dominated instead by short‐term variations, modulated for example by changes in vegetation cover or storm frequency. This distinction may help improve our understanding of long‐term sediment yields as a measure of erosional efficiency. Copyright © 2008 John Wiley & Sons, Ltd.     

Performance of Full‐Scale Enhanced Reductive Dechlorination in Clay Till

At a low permeability clay till site contaminated with chlorinated ethenes (Gl. Kongevej, Denmark), enhanced reductive dechlorination (ERD) was applied by direct push injection of molasses and dechlorinating bacteria. The performance was investigated by long‐term groundwater monitoring, and after 4 years of remediation, the development of degradation in the clay till matrix was investigated by high‐resolution subsampling of intact cores. The formation of degradation products, the presence of specific degraders Dehalococcoides spp. with the vinyl chloride (VC) reductase gene vcrA, and the isotope fractionation of trichloroethene, cis‐dichloroethene (cis‐DCE), and VC showed that degradation of chlorinated ethenes occurred in the clay till matrix as well as in sand lenses, sand stringers, and fractures. Bioactive sections of up to 1.8 m had developed in the clay till matrix, but sections, where degradation was restricted to narrow zones around sand lenses and stringers, were also observed. After 4 years of remediation, an average mass reduction of 24% was estimated. Comparison of the results with model simulation scenarios indicate that a mass reduction of 85% can be obtained within approximately 50 years without further increase in the narrow reaction zones if no donor limitations occur at the site. Long‐term monitoring of the concentration of chlorinated ethenes in the underlying chalk aquifer revealed that the aquifer was affected by the more mobile degradation products cis‐DCE and VC generated during the remediation by ERD.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling the sediment yield from landslides in the Shihmen Reservoir watershed,Taiwan

Landslides generate enormous volumes of sediment in mountainous watersheds; however, quantifying the downstream transport of landslide‐derived sediment remains a challenge. Landslide erosion and sediment delivery to the Shihmen Reservoir watershed in Taiwan was estimated using empirical landslide frequency–area and volume–area relationships, empirical landslide runout models, and the Hydrological Simulation Program‐ FORTRAN (HSPF). Landslide erosion rates ranged from 0.4 mm yr^‐1 to 2.2 mm yr^‐1 during the period 1986–2003, but increased to 7.9 mm yr^‐1 following Typhoon Aere in 2004. The percentage of landslide sediment delivered to streams decreased from 78% during the period 1986–1997 to 55% in 2004. Although the delivery ratio was lower, the volume of landslide sediment delivered to streams was 2.81 × 10⁶ Mg yr^‐1 in 1986–1997 and 8.60 × 10⁶ Mg yr^‐1 in 2004. Model simulations indicate that only a small proportion of the landslide material was delivered downstream. An average of 13% of the landslide material delivered to rivers was moved downstream during the period 1986–1997. In 2004, the period including Typhoon Aere, the annual fluvial sediment yield accounted for approximately 23% of the landslide material delivered to streams. In general, the transfer of sediment in the fluvial system in the Shihmen Reservoir watershed is dominantly transport limited. The imbalance between sediment supply and transport capacity has resulted in a considerable quantity of landslide material remaining in the upper‐stream regions of the watershed. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

Neogene upper‐crustal cooling of the Olympus range (northern Aegean): major role of Hellenic back‐arc extension over propagation of the North Anatolia Fault Zone

The North Anatolian Fault Zone (NAFZ) is one of the most hazardous active faults on Earth, yet its Pliocene space‐time propagation across the north Aegean domain remains poorly constrained. We use low‐temperature multi‐thermochronology and inverse thermal modelling to quantify the cooling history of the upper crust across the Olympus range. This range is located in the footwall of a system of normal faults traditionally interpreted as resulting from superposed Middle–Late Miocene N–S stretching, related to the back‐arc extension of the Hellenic subduction zone, and a Pliocene‐Quaternary transtensional field, attributed to the south‐westward propagation of the NAFZ. We find that accelerated exhumational cooling occurred between 12 and 6 Ma at rates of 15–35 °C Ma^?1 and decreased to <3 °C Ma^?1 by 8–6 Ma. The absence of significant Plio‐Pleistocene cooling across Olympus suggests that crustal exhumation there is driven by late Miocene back‐arc extension, while the impact of the NAFZ remains limited.     

Optical dating of sand wedges and ice‐wedge casts from Qinghai Lake area on the northeastern Qinghai‐Tibetan Plateau and its palaeoenvironmental implications

Sand wedges and ice‐wedge casts in the Qinghai Lake area on the northeastern Qinghai‐Tibetan Plateau (QTP) occur within alluvial gravel or river terrace deposits. In this study, we report the results of quartz Optically Stimulated Luminescence (OSL) dating of the infill of five relict sand wedges and one ice‐wedge cast. Combining our dating results with previously published luminescence ages of permafrost wedges in the Qinghai Lake area, we show that sand/ice wedges formed at c. 62 ka, c. 45 ka and between 30 and 15 ka, and that the mean annual air temperature (MAAT) was depressed by at least ～3°C relative to present during the sand/ice‐wedge formation periods. This new work is partially corroborated by post‐LGM proxy records from lakes and aeolian deposits reported from the northeastern QTP. It also significantly extends the palaeoenvironmental record in the region in the period before the LGM, when other proxy records are rare, allowing a better understanding of the palaeoenvironmental conditions on the northeastern QTP.