Fracture mechanics analysis of multiple cracks in anisotropic media

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.     

A multiphase approach for evaluating the horizontal and rocking impedances of pile group foundations

This paper advocates the use of a multiphase model, already developed for static or quasi‐static geotechnical engineering problems, for simulating the behaviour of piled raft foundations subject to horizontal as well as rocking dynamic solicitations. It is shown that such a model, implemented in a FEM code, yields appropriate predictions for the foundation impedance characteristics, provided that shear and bending effects in the piles are taken into account, thus corroborating the findings of the asymptotic homogenization theory. Besides, it is notably pointed out that such a multiphase‐based computational tool makes it possible to assess the dynamic behaviour of pile groups in a much quicker way than when using direct numerical simulations, which may face oversized problems when large pile groups are concerned. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Age of Gimli beach of Lake Agassiz based on new OSL dating

Gimli beach in Manitoba is one of the lowest elevation beaches in the southern Lake Agassiz basin, and is a distinct ridge composed of bedded sand and gravel that rises above the lake plain and extends for more than 40 km. Ten new optically stimulated luminescence (OSL) ages from Gimli beach yield ages mostly ranging from 9.7 ± 0.7 to 10.5 ± 0.8 ka (average 10.3 ± 0.5 ka), which is older by 0.6 to >1.0 ka than age estimates of previous researchers. Two of our new OSL ages are notably older than the others, dating to ~11.3 ± 0.8 and 13.9 ± 1.0 ka, which we attribute to poorly bleached sands. We ascribe an age of about 10 ka to Gimli beach, which is several centuries before overflow from Lake Agassiz and its vast drainage basin shifted from the western Great Lakes to glacial Lake Ojibway and the St. Lawrence Valley.     

Deciphering late Devonian–early Carboniferous P–T–t path of mylonitized garnet-mica schists from Prins Karls Forland,Svalbard

Quartz-in-garnet inclusion barometry integrated with trace element thermometry and calculated phase relations is applied to mylonitized schists of the Pinkie unit cropping out on the island of Prins Karls Forland, western part of the Svalbard Archipelago. This approach combines conventional and novel techniques and allows deciphering of the pressure–temperature (P–T) evolution of mylonitic rocks, for which the P–T conditions could not have been easily deciphered using traditional methods. The results obtained suggest that rocks of the Pinkie unit were metamorphosed under amphibolite facies conditions at 8–10 kbar and 560–630°C and mylonitized at ~500 to 550°C and 9–11 kbar. The P–T results are coupled with in-situ Th–U-total Pb monazite dating, which records amphibolite facies metamorphism at c. 359–355 Ma. This is the very first evidence of late Devonian–early Carboniferous metamorphism in Svalbard and it implies that the Ellesmerian Orogeny on Svalbard was associated with metamorphism up to amphibolite facies conditions. Thus, it can be concluded that the Ellesmerian collision between the Franklinian margin of Laurentia and Pearya and Svalbard caused not only commonly accepted brittle deformation and weak greenschist facies metamorphism, but also a burial and deformation of rock complexes at much greater depths at elevated temperatures.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Deducing the spatial variability of exchange within a longitudinal channel water balance

Developing an appropriate data collection scheme to infer stream–subsurface interactions is not trivial due to the spatial and temporal variability of exchange flowpaths. Within the context of a case study, this paper presents the results from a number of common data collection techniques ranging from point to reach scales used in combination to better understand the spatial complexity of subsurface exchanges, infer the hydrologic conditions where individual influences of hyporheic and groundwater exchange components on stream water can be characterized, and determine where gaps in information arise. We start with a tracer‐based, longitudinal channel water balance to quantify hydrologic gains and losses at a sub‐reach scale nested within two consecutive reaches. Next, we look at groundwater and stream water surface levels, shallow streambed vertical head gradients, streambed and aquifer hydraulic conductivities, water chemistry, and vertical flux rates estimated from streambed temperatures to provide more spatially explicit information. As a result, a clearer spatial understanding of gains and losses was provided, but some limitations in interpreting results were identified even when combining information collected over various scales. Due to spatial variability of exchanges and areas of mixing, each technique frequently captured a combination of groundwater and hyporheic exchange components. Ultimately, this study provides information regarding technique selection, emphasizes that care must be taken when interpreting results, and identifies the need to apply or develop more advanced methods for understanding subsurface exchanges. Copyright © 2013 John Wiley & Sons, Ltd.