Indentation of a free‐falling lance penetrometer into a poroelastic seabed

A solution is developed for the build‐up, steady and post‐arrest dissipative pore fluid pressure fields that develop around a blunt penetrometer that self‐embeds from freefall into the seabed. Arrest from freefall considers deceleration under undrained conditions in a purely cohesive soil, with constant shear strength with depth. The resulting decelerating velocity field is controlled by soil strength, geometric bearing capacity factors, and inertial components. At low impact velocities the embedment process is controlled by soil strength, and at high velocities by inertia. With the deceleration defined, a solution is evaluated for a point normal dislocation penetrating in a poroelastic medium with a prescribed decelerating velocity. Dynamic steady pressures, P_D, develop relative to the penetrating tip geometry with their distribution conditioned by the non‐dimensional penetration rate, U_D, incorporating impacting penetration rate, consolidation coefficient and penetrometer radius, and the non‐dimensional strength, N_D, additionally incorporating undrained shear strength of the sediment. Pore pressures develop to a steady peak magnitude at the penetrometer tip, and drop as P_D=1/x_D with distance x_D behind the tip and along the shaft. Peak induced pressure magnitudes may be correlated with sediment permeabilities, post‐arrest dissipation rates may be correlated with consolidation coefficients, and depths of penetration may be correlated with shear strengths. Together, these records enable strength and transport parameters to be recovered from lance penetrometer data. Penetrometer data recorded off La Palma in the Canary Islands (J. Volcanol. Geotherm. Res. 2000; 101 :253) are used to recover permeabilities and consolidation coefficients from peak pressure and dissipation response, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite element analysis of consolidation of layered clay soils using an elastic visco-plastic model

This paper presents a general one-dimensional (1-D) finite element (FE) procedure for a highly non-linear 1-D elastic visco-plastic (1-D EVP) model proposed by Yin and Graham for consolidation analysis of layered clay soils. In formulating the 1-D FE procedure, a trapezoidal formula is used to avoid the unsymmetry of the stiffness matrix for a Newton (modified Newton) iteration scheme. Unlike many other 1-D FE approaches in which the initial in situ stresses (or stress/strain states) are considered indirectly or even not considered, the initial in situ stress/strain states are taken into account directly in this paper. The proposed FE procedure is used for analysis of 1-D consolidation of a clay with published test results in the literature. The FE modelling results are in good agreement with the measured results. The FE model and procedure is then used to analyse the consolidation of a multi-layered clay soils with a parametric study on the effects of the variations of creep parameters in Yin and Graham's 1-D EVP model. It is found that the creep parameters ψ/V and t₀ have significant influence on the compression and porewater pressure dissipation. For some boundary conditions, changes of parameters in one layer will have some effects on the consolidation behaviour of another layer due to the different consolidation rates. Finally, the importance of initial stress/strain states is illustrated and discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Effect of soil layering on the ground response parameters: a parametric study

Soil layering has a pivotal role on the behavior and propagation of seismic waves, hence the ground response during seismic loading. Parametric study to estimate the effect of soil layering on ground response parameters is of prime importance considering the engineering significance of structure founded in seismic zone; nonetheless, it is yet to be well understood. The objective of this study was to investigate the effect of soil layering and soil properties on the ground response parameters. One-dimensional linear ground response analysis was conducted with variation in soil layer parameters including impedance ratio (?? _z ), layer thickness (d _r ), and damping ratio (D). The acceleration time history of the Bhuj Earthquake (M _w ?=?7.7), India, was used in the analysis. The results obtained from the analysis were presented as variation of ground response parameters such as spectral acceleration (SA), amplification ratio (M) with the soil layer parameters. Results showed higher values of SA at lower D and then decreased with increase in D, that in fact depict the resistance offered to the particle oscillation at comparatively higher values of D. Similarly, variation in SA and M was very less or negligible when the ?? _z was varied from 1 to 3 and the d _r equal to 0.2, while for d _r greater than 0.2 the variation increased with ?? _z and d _r . The outcome from the parametric study presented in this paper clearly demonstrates the significance of ?? _z , D, and d _r of the soil layers on the ground response parameters.     

Decline curve analyses revisited

Empiricaldecline curve forecasts of future production rates and cumulative production to be expected from certain classes of petroleum reservoirs has been a methodology practiced and relied on by petroleum geologists for many years. New ways to consider the subject are described in this paper. Thereby it is seen that these forecasts sometimes can be made even when some of the record of the past performance history is sparse (say, because it has been lost, or never recorded).SI Units D, D _i initial reference production rate,s ^–1 - n dimensionless constant [cf. Eq. (1)] - N cumulative production, m³ - Q, Q _i production rate att and initially, m³/s - R dimensionless time variable [cf. Eq. (9)] - t time, s     

Fast diffusion along mobile grain boundaries in calcite

Experimental measurements of grain boundary diffusion are usually conducted on static boundaries, despite the fact that grain boundaries deep in the Earth are frequently mobile. In order to explore the possible effect of boundary mobility on grain boundary diffusion rates we have measured the uptake of ⁴⁴Ca from a layer of ⁴⁴Ca-enriched calcite powder during the static recrystallization of a single crystal of calcite at 900°C. A region about 500 μm wide adjacent to the powder layer is heterogeneously enriched in ⁴⁴Ca, and complex zoning patterns, including sharp steps in composition and continuous increases and decreases in ⁴⁴Ca content, are developed. In metamorphic rocks, these would normally be interpreted in terms of changes in pressure or temperature, Rayleigh fractionation, or episodic fluid infiltration. These explanations cannot apply to our experiments, and instead the zoning patterns are interpreted as being due to variations in grain boundary migration rate. We have applied an analytical model which allows the product of grain boundary diffusion coefficient and grain boundary width (D _GB δ) to be calculated from the grain boundary migration rate and the compositional gradient away from the powder layer. The value of D _GB δ in the mobile grain boundaries is at least five orders of magnitude greater than the published value for static boundaries under the same conditions. In order to allow the scale of chemical equilibrium (and hence textural evolution) to be predicted under both experimental and geological conditions, we present quantitative diffusion-regime maps for static and mobile boundaries in calcite, using both published values and our new values for grain boundary diffusion in mobile boundaries. Enhanced diffusion in mobile boundaries has wide implications for the high temperature rheology of Earth materials, for geochronology, and for interpretations of the length- and time-scales of chemical mass-transport. Moreover, zones of anomalously high electrical conductivity in the crust and mantle could be regions undergoing recrystallization such as active shear zones, rather than regions of anomalous mineralogy, water- or melt-content as is generally suggested.     

Simplified DMT-based methods for evaluating liquefaction resistance of soils

This paper presents simplified dilatometer test (DMT)-based methods for evaluation of liquefaction resistance of soils, which is expressed in terms of cyclic resistance ratio (CRR). Two DMT parameters, horizontal stress index (K_D) and dilatometer modulus (E_D), are used as an index for assessing liquefaction resistance of soils. Specifically, CRR–K_D and CRR–E_D boundary curves are established based on the existing boundary curves that have already been developed based on standard penetration test (SPT) and cone penetration test (CPT). One key element in the development of CRR–K_D and CRR–E_D boundary curves is the correlations between K_D (or E_D) and the blow count (N) in the SPT or cone tip resistance (q_c) from the CPT. In this study, these correlations are established through regression analysis of the test results of SPT, CPT, and DMT conducted side-by-side at each of five sites selected. The validity of the developed CRR–K_D and CRR–E_D curves for evaluating liquefaction resistance is examined with published liquefaction case histories. The results of the study show that the developed DMT-based models are quite promising as a tool for evaluating liquefaction resistance of soils.     

An integrated law of the wall for hydrodynamically transitional flow over plane beds

The velocity profile for unidirectional currents is usually determined by two forms of the law of the wall, which are valid for smooth and rough boundaries, respectively. The law of the wall for smooth boundaries outside the viscous sublayer and buffer layer applies where the boundary Reynolds number (Re_*) is less than 5. The velocity in this case depends only on the dimensionless distance (y_d) from the bottom. For rough boundaries where Re_* is more than 65, the law of the wall is independent of the boundary Reynolds number, but the bottom roughness (k) has to be taken into account. The mean or median particle grain size (D) is commonly substituted for bottom roughness. At present, there is no generally accepted law of the wall for the transitional flow regime between hydrodynamically smooth and rough conditions. In this paper, an integrated law of the wall is proposed for the transitional regime, yielding velocity profiles that correspond well to observed velocity profiles under these conditions. The result is applied to a set of flume data, using an improved procedure to determine the bedload transport rate.     

Numerical modeling of three-dimensional contaminant migration from dissolution of multicomponent NAPL pools in saturated porous media

A three-dimensional model for contaminant transport resulting from the dissolution of multicomponent nonaqueous phase liquid (NAPL) pools in three-dimensional saturated subsurface formations is developed. The solution is obtained numerically by a finite-difference scheme, and it is suitable for homogeneous porous media with unidirectional interstitial velocity. Each dissolved component may undergo first-order decay and may sorb under local equilibrium conditions. It is also assumed that the dissolution process is mass transfer limited. The nonaqueous phase activity coefficients of the NAPL pool components are evaluated at each time step. The model behavior is illustrated through a synthetic example with a NAPL pool consisting of a mixture of TCA (1,1,2-trichloroethane) and TCE (trichloroethylene). The numerical solution presented in this work is in good agreement with a recently developed analytical solution for the special case of a single component NAPL pool. The results indicate the importance of accounting for the necessary changes in the organic phase activity which significantly affects the equilibrium aqueous solubility.Notation C liquid phase solute concentration (solute mass/liquid volume) (M L^–3) - C _s single component aqueous saturation concentration (solubility) (M L^–3) - C ^w equilibrium aqueous solubility (M L^–3) - D molecular diffusion coefficient (L² t ^–1) - D _e effective molecular diffusion coefficient (L² t ^–1) - D _x longitudinal hydrodynamic dispersion coefficient (L² t ^–1) - D _y lateral hydrodynamic dispersion coefficient (L² t ^–1) - D _z hydrodynamic dispersion coefficient in the vertical direction (L² t ^–1) - I() integer mode arithmetic operator - k local mass transfer coefficient (Lt ^–1) - k ^* average mass transfer coefficient (Lt ^–1) - L length - l _x ,l _y pool dimensions inx andy directions (L) - ll _x ,l _y x andy Cartesian coordinates of the pool origin (L) - M number of moles remaining in a pool (moles) - M initial number of moles (moles) - n finite-difference scheme time level - R retardation factor (dimensionless) - t time (t) - U _x average interstitial velocity (Lt ^–1) - x, y, z spatial Cartesian coordinates (L) - X dimensionless mole fraction - dimensionless activity coefficient - _w viscosity of water (=0.8904 cp at 25°C) - decay coefficient (t ^–1) - ^* tortuosity ( 1) - i,j, k finite-difference scheme grid indicators - p component number indicator - P total number of components - s pure single component - o nonaqueous phase - w aqueous phase     

Unravelling the tectonothermal evolution of reworked Archean granulite facies metapelites using in situ geochronology: an example from the Gawler Craton,Australia

The use of in situ geochronological techniques allows for direct age constraints to be placed on fabric development and the metamorphic evolution of polydeformed and reworked terranes. The Shoal Point region of the southern Gawler Craton consists of a series of reworked granulite facies metapelitic and metaigneous units which belong to the Late Archean Sleaford Complex. Structural evidence indicates three phases of fabric development with D₁ retained within boudins, D₂ consisting of a series of upright open to isoclinal folds producing an axial planar fabric and D₃ composed of a highly planar vertical high‐strain fabric which overprints the D₂ fabric. Th–U–total Pb EPMA monazite and garnet Sm–Nd geochronology constrain the D₁ event to the c. 2450 Ma Sleaford Orogeny, whereas the D₂ and D₃ events are constrained to the 1730–1690 Ma Kimban Orogeny. P–T pseudosections constrain the metamorphic conditions for the Sleafordian Orogeny to between 4.5 and 6 kbar and between 750 and 780 °C. Subsequent Kimban‐aged reworking reached peak metamorphic conditions of 8–9 kbar at 820–850 °C during the D₂ event, followed by high‐temperature decompression to metamorphic conditions <6 kbar and 790–850 °C associated with the development of the D₃ high‐strain fabric. The P–T–t evolution of the Shoal Point rocks reflects the transpressional exhumation of lower crustal rocks during the Kimban Orogeny and the development of a regional ‘flower structure’.     

Microstructural evolution of the Seridó Belt, NE-Brazil: the effect of two tectonic events on development of c-axis preferred orientation in quartz

The polyphase evolution of the Seridó Belt (NE-Brazil) includes D₁ crust formation at 2.3–2.1 Ga, D₂ thrust tectonics at 1.9 Ga and crustal reworking by D₃ strike-slip shear zones at 600 Ma. Microstructural investigations within mylonites associated with D₂ and D₃ events were used to constrain the tectono-thermal evolution of the belt. D₂ shear zones commenced at deeper crustal levels and high amphibolite facies conditions (600–650 °C) through grain boundary migration, subgrain rotation and operation of quartz c-prism slip. Continued shearing and exhumation of the terrain forced the re-equilibration of high-T fabrics and the switching of slip systems from c-prism to positive and negative a-rhombs. During D₃, enhancement of ductility by dissipation of heat that came from syn-D₃ granites developed wide belts of amphibolite facies mylonites. Continued shearing, uplift and cooling of the region induced D₃ shear zones to act in ductile-brittle regimes, marked by fracturing and development of thinner belts of greenschist facies mylonites. During this event, switching from a-prism to a-basal slip indicates a thermal path from 600 to 350 °C. Therefore, microstructures and quartz c-axis fabrics in polydeformed rocks from the Seridó Belt preserve the record of two major events, which includes contrasting deformation mechanisms and thermal paths.     

Metamorphic evolution of a subduction complex, South Shetland Islands, Antarctica

A subduction complex composed of ocean floor material mixed with arc-derived metasediments crops out in the Elephant Island group and at Smith Island, South Shetland Islands, Antarctica, with metamorphic ages of 120–80 Ma and 58–47 Ma, respectively. Seven metamorphic zones (I–VII) mapped on Elephant Island delineate a gradual increase in metamorphic grade from the pumpellyite–actinolite facies, through the crossite–epidote blueschist facies, to the lower amphibolite facies. Geothermometry in garnet–amphibole and garnet–biotite pairs yields temperatures of about 350 °C in zone III to about 525 °C in zone VII. Pressures were estimated on the basis of Si content in white mica, Al₂O₃ content in alkali amphibole, Na^M4/Al^IV in sodic-calcic and calcic amphibole, Al^VI/Si in calcic amphibole, and jadeite content in clinopyroxene. Mean values vary from about 6–7.5 kbar in zone II to about 5 kbar in zone VII. Results from the other islands of the Elephant Island group are comparable to those from the main island; Smith Island yielded slightly higher pressures, up to 8 kbar, with temperatures estimated between 300 and 350 °C. Zoned minerals and other textural indications locally enable inference of P–T –t trajectories, all with a clockwise evolution. A reconstruction in space and time of these P–T –t paths allows an estimate of the thermal structure in the upper crust during the two ductile deformation phases (D₁ & D₂) that affected the area. This thermal structure is in good agreement with the one expected for a subduction zone. The arrival and collision of thickened oceanic crust may have caused the accretion and preservation of the subduction complex. In this model, D₁ represents the subduction movements expressed by the first vector of the clockwise P–T–t path, D₂ reflects the collision corresponding to the second vector with increasing temperature and decreasing pressure, and D₃ corresponds to isostatic uplift accompanied by erosion, under circumstances of decreasing temperature and pressure.     

Diffusion along interphase boundaries and its effect on retrograde zoning patterns of metamorphic minerals

In this study we use two dimensional chemical patterns and numerical modeling to estimate the relative rates of chemical transport along interphase boundaries (ib) and through grain (s) interiors during retrograde Fe–Mg exchange between garnet and biotite at a garnet–biotite–quartz triple junction. We demonstrate that systematic variations in garnet–rim compositions, which are frequently observed along garnet–quartz interfaces, and deviations from concentric retrograde zoning patterns start to develop when chemical transport along the interphase boundaries becomes slow during cooling. The capacities for chemical transport along an interphase boundary depend on the product D ^ib K ^ib/s a, where D ^ib is the diffusion coefficient of the exchangeable components within the interphase boundary medium, K ^ib/s is the equilibrium partitioning coefficient between the cation exchange partners and the interphase boundary medium and a is the interphase boundary width. The model is applied to analyze the retrograde zoning patterns in garnets from the Mozambique belt (SE-Kenya), which cooled from 820°C at a rate of ca. 2°C/my. It is found that non-equilibrated compositions in garnet along garnet/quartz interphase boundaries started to develop below 700°C due to insufficient rates of chemical transport along these boundaries. The transport capacities of garnet/quartz interphase boundaries was estimated to have been between about 1E-23 cm³/s (575°C) and 1E-20 cm³/s (700°C) from modeling the observed X _Fe pattern in garnet close to a garnet–quartz–biotite triple junction and relying on published data on the diffusivity of Fe²⁺ in garnet. Similar transport capacities are obtained; when the interphase boundary is assumed to be filled with a material that has the transport properties and chemical composition of a free melt in equilibrium with garnet, biotite and quartz at the respective conditions. In contrast, if the transport properties of the interphase boundary medium are related to the diffusivities and solubility of Fe²⁺ and FeOH⁺ within a free aqueous solution, chemical transport along the interphase boundaries would be much more efficient, and exchange equilibrium would have been maintained during the entire cooling history of the rocks. The observation of systematic deviations from local equilibrium along the garnet–quartz interphase boundaries leads us to exclude the presence of an aqueous fluid along the interphase boundary at any time during cooling.     

An integral-free expression for short-term changes in fault stability due to pore pressure induced when a point load is placed on the pervious boundary of a porous elastic half space containing a fault

A concentrated load with step-function time behaviour is placed normal to the planar, pervious boundary of a porous elastic half space (PEHS) with compressible constituents. A planar fault exists in the PEHS in such a way that the poroelastic behaviour of the medium is unhindered. We derive an approximate but integral-free expression for CFSCPP, i.e., changes in fault stability due to changes in pore pressure, at a point not too far off the line along which the load acts. But, in the interest of simplicity, the main discussion is focussed on a consideration of CFSCPP at a point P located on the fault at depth z directly beneath the load. It is convenient to introduce dimensionless time t _D directly proportional to real time t. The constant of proportionality is 4c/z ², where c is hydraulic diffusivity. The derived approximate expression gives results with an accuracy of greater than 99% for limited values of t _D after the load is imposed. We learn from the derived expression that, for a given z, fault stability undergoes an initial sudden decrease commensurate with the undrained pore pressure induced in the PEHS. This is followed by a more gradual decrease in fault stability with increasing t _D until a minimum is reached. The real time t to minimum fault stability increases with z. The magnitude of CFSCPP decreases with z as z ^?2 for a given t _D in the permissible range. The derived expression and the inferences based on it should be useful during earth science investigations of the possible hazards due to reactivation of a pre-existing shallow fault when a civil engineering project involving imposition of a heavy load on the earth’s surface is to be executed nearby. They should be useful also for investigations if a shallow earthquake occurs near such a project soon after its execution.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Deformation,metamorphism, and mobilization of Ni–Cu–PGE sulfide ores at Garson Mine,Sudbury

The Garson Ni–Cu–platinum group element deposit is a deformed, overturned, low Ni tenor contact-type deposit along the contact between the Sudbury Igneous Complex (SIC) and stratigraphically underlying rocks of the Huronian Supergroup in the South Range of the 1.85-Ga Sudbury structure. The ore bodies are coincident with steeply south-dipping, north-over-south D₁ shear zones, which imbricated the SIC, its ore zones, and underlying Huronian rocks during mid-amphibolite facies metamorphism. The shear zones were reactivated as south-over-north, reverse shear zones during D₂ at mid-greenschist facies metamorphism. Syn-D₂ metamorphic titanite yields an age of 1,849?±?6 Ma, suggesting that D₁ and D₂ occurred immediately after crystallization of the SIC during the Penokean Orogeny. The ore bodies plunge steeply to the south parallel to colinear L₁ and L₂ mineral lineations, indicating that the geometry of the ore bodies are strongly controlled by D₁ and D₂. Sulfide mineralization consists of breccia ores, with minor disseminated sulfides hosted in norite, and syn-D₂ quartz–calcite–sulfide veins. Mobilization by ductile plastic flow was the dominant mechanism of sulfide/metal mobilization during D₁ and D₂, with additional minor hydrothermal mobilization of Cu, Fe, and Ni by hydrothermal fluids during D₂. Metamorphic pentlandite overgrows a S₁ ferrotschermakite foliation in D₁ deformed ore zones. Pentlandite was exsolved from recrystallized polygonal pyrrhotite grains after cessation of D₁, which resulted in randomly distributed large pentlandite grains and randomly oriented pentlandite loops along the grain boundaries of polygonal pyrrhotite within the breccia ore. It also overgrows a S₂ chlorite foliation in D₂ shear zones. Pyrrhotite recrystallized and was flattened during D₂ deformation of breccia ore along narrow shear zones. Exsolution of pentlandite loops along the grain boundaries of these flattened grains produced a pyrrhotite–pentlandite layering that is not observed in D₁ deformed ore zones. The overprinting of the two foliations by pentlandite and exsolution of pentlandite along the grain boundaries of flattened pyrrhotite grains suggest that the Garson ores reverted to a metamorphic monosulfide solid solution at temperatures ranging between 550 and 600 °C during D₁ and continued to deform as a monosulfide solid solution during D₂.     

The evolution of a grospydite from the Roberts Victor mine,South Africa

A grospydite from Roberts Victor contains the most Ca-rich garnets yet found in South African kimberlite xenoliths and also sub-micron sized sodic nepheline in melted and quenched clinopyroxene. Three stages can be recognised in the textural evolution of the grospydite. The first is the development of a layering of large kyanite laths. Kyanite together with complex aluminous clinopyroxene precipitated and accumulated from an evolved residual eclogitic liquid which has penetrated across the garnet join so that garnet no longer precipitated. Solidus conditions for the Roberts Victor grospydite are estimated as T = 1350–1550 ° C, P = 27–39 kbars. Adjacent layers in the grospydite have slightly different mineral compositions suggesting that the small-scale layering (1–5 cm) in this, and associated rocks, may be related to varying activities of R₂O₃ components and possibly to f _{o 2}.The second stage is represented by a necklace texture in which all the garnet and some kyanite developed along grain boundaries of clinopyroxenes with triple-point textures. This is interpreted as an example of incoherent, grain-boundary exsolution resulting from large subsolidus volume changes. The conditions for subsolidus equilibration are estimated to be T= 1120–1320 ° C, P = 42–56 kbars.Moderate Ca-contents in garnet and excess Al^[6] in clinopyroxene may be subsolidus indicators of eclogite samples evolving towards grospydite at the solidus.The third stage is represented by the melting of jadeite-rich clinopyroxenes and quenching to glass, nepheline and plagioclase. Most of the glass has a composition similar to clinopyroxene, except for K₂O, though local areas of different glass, possibly the result of phase separation, also occur. The melting process seems to be a low-pressure feature involving limited addition of H₂O at temperatures between 900–1000 ° C. Water-absent melting could indicate temperatures up to 1500 ° C.The temperatures and pressures assigned to the three-stage evolution of this grospydite imply formation at moderate pressures and subsolidus equilibration at higher pressures. This is equated with downgoing mantle/asthenosphere tectonic processes. After entrainment in a kimberlite magma the grospydite fragment apparently ascended rapidly, thus allowing low-pressure melting and quenching.     

Rock mass characterization by fractal dimension

Fractal theory is used in the present study to develop a more reliable method for rock mass characterization. Field studies have been carried out in opencast mines of dolomite, limestone, fluorite; sandstone and shale in coalmines. Fractal dimension of blasted fragments (D_frag) and in situ rock blocks (D_{in situ}) is calculated using size distribution curves according to Schumann's model. Based on the co-relation between Uniaxial Compressive Strength (UCS) and D_frag, it is observed that change in fractal dimension is nominal beyond the UCS value of 20. From the co-relation between Bieniawaski's Rock Mass Rating (RMR) and D_{in situ}, it is found that there is a sharp increase in fractal dimension for RMR greater than 40. Co-relation between RMR and D_frag/D_{in situ} shows that as RMR increases, D_frag/D_{in situ} ratio decreases. Rock mass classification based on fractal geometry is suggested.     

Rainfall intensity–duration thresholds for postfire debris-flow emergency-response planning

Following wildfires, emergency-response and public-safety agencies can be faced with evacuation and resource-deployment decisions well in advance of coming winter storms and during storms themselves. Information critical to these decisions is provided for recently burned areas in the San Gabriel Mountains of southern California. A compilation of information on the hydrologic response to winter storms from recently burned areas in southern California steeplands is used to develop a system for classifying magnitudes of hydrologic response. The four-class system describes combinations of reported volumes of individual debris flows, consequences of debris flows and floods in an urban setting, and spatial extents of the hydrologic response. The range of rainfall conditions associated with different magnitude classes is defined by integrating local rainfall data with the response magnitude information. Magnitude I events can be expected when within-storm rainfall accumulations (A) of given durations (D) fall above the threshold A = 0.4D ^0.5 and below A = 0.5D ^0.6 for durations greater than 1 h. Magnitude II events will be generated in response to rainfall accumulations and durations between A = 0.4D ^0.5 and A = 0.9D ^0.5 for durations less than 1 h, and between A = 0.5D ^0.6 and A = 0.9D ^0.5 or durations greater than 1 h. Magnitude III events can be expected in response to rainfall conditions above the threshold A = 0.9D ^0.5. Rainfall threshold-magnitude relations are linked with potential emergency-response actions as an emergency-response decision chart, which leads a user through steps to determine potential event magnitudes and identify possible evacuation and resource-deployment levels. Use of this information in planning and response decision-making process could result in increased safety for both the public and emergency responders.     

Structures associated with superimposed inhomogenous shearing of basic dykes from Finnmark, Norway

Metadolerites within the Olderfjord Nappe of the Kalak Nappe Complex of North Norway show a wide range of structural relationships. During D₂ simple shear rotated the dykes from an initially orthogonal relationship towards parallelism with the sedimentary layering and boudinaged those dykes that rotated to within 25° of the sedimentary layering. During D₃ dykes which had been rotated to within 5° of the layering were folded; those lying between 5° and 25° were boudinaged whilst those at angles between 25° and 75° had folds developed against their margins. At higher angles of discordance little evidence of D₂ or D₃ deformation has been observed.In other thrust sheets within the Kalak Nappe Complex metabasites have generally been sheared into parallelism with the compositional banding. By comparison with the Olderfjord Nappe the average D₂ shear strain within the Kalak Nappe Complex ranges between 5 and 10. The anomalous preservation of high angle discordant relationships in the Olderfjord Nappe reflects the preservation of a large scale low-strain block (mega-boudin) within the nappe; samples taken from the mega-boudin are more suitable for geochemical analysis for tectonic setting discrimination.     

Grain boundary and volume diffusion experiments in yttrium aluminium garnet bicrystals at 1,723 K: a miniaturized study

Yb-Y inter-diffusion along a single grain boundary of a synthetic yttrium aluminium garnet (YAG) bicrystal has been studied using analytical transmission electron microscopy (ATEM). To investigate the diffusion, a thin-film containing Yb as the diffusant was deposited perpendicular to the bicrystal grain boundary by pulsed laser deposition (PLD). Structural properties and their change with time in both the diffusant source and the grain boundary are reported. The diffusion profiles are incorporated in a numerical diffusion model, which is applied to determine the grain boundary diffusion coefficient, D _gb , at 1.723 K it is equal to 3 × 10⁻¹⁵ m²/s. We find that grain boundary diffusion is 4.85 orders of magnitude faster than volume diffusion, which was determined from the same diffusion experiment. This result is discussed in the context of special versus general grain boundaries. Finally, we successfully tested the capability of synchrotron-based nano-X-ray fluorescence analysis to map micro-chemical patterns in two dimensions with sub-micrometre resolution.