Turbulent dispersion in the ocean

The mathematical framework for turbulent transport in the ocean is reasonably well established. It may be applied to large-scale fields of scalars in the ocean and to the instantaneous or continuous discharge from a point. The theory and its physical basis can also provide an interpretation of passive scalar spectra. Spatial variations in the rate of turbulent transfer can be related to the movement of the center of mass of a scalar and to a formulation in terms of entrainment. The relative dispersion of a scalar with respect to its center of mass and the streakiness of the concentration field within the relative dispersion domain need to be considered. In many of these problems it is valuable to think in terms of simple models for individual streaks, as well as overall statistical properties.     

The Creation of a National Agricultural Land Use Dataset: Combining Pycnophylactic Interpolation with Dasymetric Mapping Techniques

Agricultural Census data is summarised over spatially coarse reporting units for reasons of farm confidentiality. This is problematic for research at a local level. This article describes an approach combining dasymetric and volume preserving techniques to create a national land use dataset at 1 km² resolution. The results for an English county are compared with contemporaneous aggregated habitat data. The results show that the accurate estimates of local agricultural land use (Arable and Grass) patterns can be estimated when individual 1 km squares are combined into blocks of > 9 squares, thereby providing local estimates of agricultural land use. This in turn allows more detailed modelling of land uses related to specific livestock and cropping activities. The dataset created by this work has been subject to extensive external validation through its incorporation into a number of other national models: nitrate leaching (e.g. MAGPIE, NEAP‐N), waste, and pathogen modelling related to agricultural activity.     

Kinetics of peptide hydrolysis and amino acid decomposition at high temperature

Dipeptide hydrolysis and amino acid decomposition appear to follow a first-order rate law. The hydrolysis rate increases exponentially with increasing temperature in aqueous solution at both 265 atm and water steam pressures over the temperature range of 100 to 220 degrees C. Dipeptide hydrolysis has a lower apparent activation energy at 265 atm (44.1 KJ/mol) than at water steam pressure (98.9 KJ/mol). At lower temperatures (<200-220 degrees C), the rate of peptide bond hydrolysis is faster at 265 atm than at water steam pressure. At higher temperatures (>200-220 degrees C), however, peptide bond hydrolysis is slower at 265 atm than at water steam pressure. In aqueous solution, amino acid decomposition rates also increase exponentially with increasing temperature. Amino acid decomposition rates are much higher at 265 atm than at water steam pressure over the entire temperature range investigated.     

Calibration of Al/Si order variations in anorthite

New single crystal diffraction data for natural and heat-treated anorthite crystals (Angel et al. 1990) allow the determination of their states of Al/Si order in terms of a macroscopic order parameter,Q _OD , for the transition. Numerical values ofQ _OD obtained from estimates of site occupancies are shown to vary with the scalar spontaneous strain, _s , as _s Q _OD ² , and with the ratio of the sums of typeb (superlattice) reflections and typea (sublattice) reflections asI _b/I _a Q _OD ² . An empirical calibration for pure anorthite is obtained giving varies between 0.92 and 0.87 in samples equilibrated at T1300° C, but then falls off relatively rapidly with increasing temperature, reaching 0.7 near the melting point ( 1557° C). The observed temperature dependence does not conform to the predictions of the simplest single order parameter models; coupling ofQ _OD withQ of the transition is suspeeted.     

Variability in Early Amazonian Tharsis stress state based on wrinkle ridges and strike-slip faulting

We present evidence for a decrease in the magnitude of Tharsis-circumferential compressive stress during the Late Hesperian to the Middle Amazonian based on chronologic changes in the predominant style of faulting in southern Amazonis Planitia. Using high-resolution MOLA topography, we identify a population of strike-slip faults that exhibit Middle Amazonian-aged displacements of regional chrono-stratigraphic units. These strike-slip faults are adjacent to an older population of previously documented Late Hesperian-aged thrust faults (wrinkle ridges). Along-strike orientations of these thrust and strike-slip faults reveal the Tharsis-radial stress to be the area's most compressive remote principal stress and that this stress orientation and magnitude persisted throughout the Late Hesperian to the Middle Amazonian. We show that the change in the predominant style of faulting from thrust faulting to strike-slip faulting during this time requires a decrease of the Tharsis-circumferential compressive stress to a magnitude less than lithostatic load, with negligible change in stress orientation.     

Displacement-length scaling relations for faults on the terrestrial planets

Displacement-length (D/L)scaling relations for normal and thrust faults from Mars, and thrust faults from Mercury, for which sufficiently accurate measurements are available, are consistently smaller than terrestrial D/L ratios by a factor of about 5, regardless of fault type (i.e. normal or thrust). We demonstrate that D/L ratios for faults scale, to first order, with planetary gravity. In particular, confining pressure modulates: (1) the magnitude of shear driving stress on the fault; (2) the shear yield strength of near-tip rock; and (3) the Young's (or shear) modulus of crustal rock. In general, all three factors decrease with gravity for the same rock type and pore-pressure state (e.g. wet conditions). Faults on planets with lower surface gravities, such as Mars and Mercury, demonstrate systematically smaller D/L ratios than faults on larger planets, such as Earth. Smaller D/L ratios of faults on Venus and the Moon are predicted by this approach, and we infer still smaller values of D/L ratio for faults on icy satellites in the outer solar system. Collection of additional displacement-length and down-dip height data from terrestrial normal, strike-slip, and thrust faults, located within fold-and-thrust belts, plate margins, and continental interiors, is required to evaluate the influence of fault shape and progressive deformation on the scaling relations for faults from Earth and elsewhere.     

Geomechanical model for fracture deformation under hydraulic, mechanical and thermal loads

Hydraulic flow and transport (heat and solute) within crystalline rocks is dominated by the fracture systems found within them. In situ stress conditions have a significant impact on the hydraulic, mechanical and thermal coupled processes, and quantification of these processes provides a key to understanding the often transient time-dependent behaviour of crystalline rocks. In this paper, a geomechanical model is presented which describes fracture closure as a function of effective stress and the changes in parameters such as storage, permeability, porosity and aperture. Allowing the fracture closure to be defined by the change in normal effective stress provides a link to the numerical consideration of parametrical changes due to rock stress alterations caused for example by changes in fracture fluid pressure, stress release, tectonic stress, thermal stress, orientation of the natural fracture in the pervasive stress system and local changes in a rock mass due to stress alteration. The model uses geometrical considerations based on a fractal distribution of apertures on the fracture surface, and applies well-established analytical elastic deformation solutions to calculate the deformation response to changes in effective stress. Analysis of the fractal generation method allows a standard normal distribution of fracture apertures to be predicted for all common fractal dimensions relating to a 2D surface. Changes in the fracture aperture are related to hydraulic functions such as permeability, storage and porosity of the fracture. The geomechanical model is experimentally validated against laboratory scale experimental data gained from the closure of a fractured sample recovered at a depth of 3,800 m from the KTB pilot borehole. Parameters for matching the experimental data were established externally, the only fitting parameters applied were the minimum and maximum contact area between the surfaces and the number of allowable contacts. The model provides an insight into the key processes determining the closure of a fracture, and can act as a material input function for numerical models linking the effects of changes in the stress field, hydraulic or thermal conditions, to the flow and transport parameters of a fractured system.

---

Résumé L’écoulement et le transport (chaleur et soluté) dans les roches cristallines sont dominés par les systèmes de fracture. Les conditions de stress in-situ ont un impact significatif sur l’hydraulique, les processus couplés de mécanique et thermique et la quantification de ces processus apportent une clé pour comprendre le comportement transitoire des roches cristallines. Dans cet article un modèle géomécanique est présenté, modèle qui décrit la fermeture des fractures comme une fonction de la contrainte effective et des changements de paramètres tels le coefficient d’emmagasinement, la perméabilité, la porosité et l’ouverture. En s’accordant que la fermeture des fractures est définit par les changements de la contrainte effective normale, on apporte le lien avec la considération numérique des changements paramétriques dus aux altérations de la contrainte des roches, causés par exemple par des variations de la pression des fluides dans les fractures, du dégagement de la contrainte, des contraintes tectoniques et thermiques, des orientations des fractures naturelles dans le système de contraintes pénétrantes, et des changements locaux dans un massif de roches dus à l’altération des contraintes. Le modèle utilise des considérations géométriques basées sur une distribution fractale des ouvertures à la surface des fractures, et permet d’établir des solutions analytiques de la déformation élastique pour calculer la réponse de la déformation à la contrainte effective. L’analyse de la méthode par génération fractale permet de prédire une distribution normale standard de l’ouverture des fractures, pour toutes les dimensions fractales en relation avec les surfaces 2D. Les changements dans l’ouverture des fractures sont mis en relation avec les fonctions hydrauliques tels la perméabilité, l’emmagasinement et la porosité de la fracture. Le modèle géoméchanique est expérimentalement validé à l’échelle du laboratoire sur un échantillon fracturé récupéré à une profondeur de 3,800 mètres sur le puits du site pilote KTB. Les paramètres du calibrage des données expérimentales ont été établies extérieurement, les seuls paramètres utilisés étant les surfaces de contact minimum et maximum, et le nombre de contacts permis. Le modèle apporte une connaissance perspicace sur le processus clé déterminant la fermeture des fractures, et peut servir de fonction input dans les modèles numériques reliant les effets des variations de la contrainte du terrain, les conditions hydrauliques ou thermales, les paramètres de l’écoulement et du transport et les systèmes de fracture.

---

Resumen El flujo hidráulico y transporte (de calor y solutos) dentro de rocas cristalinas está dominado por los sistemas de fracturas que se encuentran en ellas. Las condiciones de esfuerzos in-situ tienen un impacto significativo en los procesos aparejados termales, mecánicos e hidráulicos y la cuantificación de estos procesos aporta una clave para entender el frecuente comportamiento transitorio dependiente de las rocas cristalinas. En este artículo se presenta un modelo geomecánico que describe el cierre de fracturas en función del esfuerzo efectivo y los cambios en parámetros tal como almacenamiento, permeabilidad, porosidad y apertura. El definir el cierre de fractura mediante el cambio en esfuerzo normal efectivo aporta un vínculo con la consideración numérica de cambios paramétricos ocasionados por alteraciones de esfuerzos en la roca causadas, por ejemplo, por cambios en presión de fluidos en fractura, liberación de esfuerzo, esfuerzo tectónico, esfuerzo termal, orientación de fracturas naturales en el sistema de esfuerzos penetrante, y cambios locales en una masa rocosa ocasionados por alteración de esfuerzos. El modelo utiliza consideraciones geométricas basadas en la distribución fractal de aperturas en la superficie de fractura y aplica soluciones analíticas bien establecidas de deformación elástica para calcular la respuesta de deformación a cambios en el esfuerzo efectivo. Los análisis del método de generación fractal permiten predecir una distribución normal standard para la distribución de aperturas de fracturas para todas las dimensiones fractales comunes que se relacionan con una superficie 2D. Los cambios en la apertura de fractura se relacionan con funciones hidráulicas tal como permeabilidad, almacenamiento y porosidad de la fractura. El modelo geomecánico se ha validado experimentalmente en contra de datos experimentales a escala de laboratorio obtenidos a partir del cierre de una muestra fracturada recuperada a una profundidad de 3,800 m en el pozo piloto KTB. Se establecieron externamente parámetros que se ajustan a los datos experimentales, con los parámetros de ajuste aplicados que fueron el área máxima y mínima de contacto entre las superficies y el número de contactos permisibles. El modelo arroja luz sobre los procesos clave que determinan el cierre de una fractura y puede actuar como un material de función de entrada para modelos numéricos que vinculan los efectos de cambios en el campo de esfuerzos, condiciones termales o hidráulicas, con los parámetros de flujo y transporte de un sistema fracturado.

     

Origin of spongy textures in clinopyroxene and spinel from mantle xenoliths, Hessian Depression, Germany

Summary Spongy textures are observed in anhydrous Group 1 mantle xenoliths (harzburgite, lherzolite and wehrlite) hosted in Tertiary alkali basaltic lavas from the Hessian Depression, Germany. These textures are developed only on clinopyroxene and spinel, and occur as rims or cross-cutting veinlets and patches showing optical continuity with the host grain. They are often associated with pools of amorphous glassy material. There is no preferential development of spongy domains against the xenolith-lava contact suggesting that the host magma did not play any significant role in their formation. Spongy clinopyroxene and spinel occur in all rock types, but, are more pervasive in wehrlite. Chemically, spongy domains of clinopyroxene and spinel are more refractory than unaffected areas, which is consistent with their formation through a partial melting event. The associated glassy material shows chemical characteristics which suggest that the melt pools are genetically related to the development of the spongy textures. The partial melting event was probably triggered by the infiltration of a low-density fluid. The fluid may have evolved from a silicate melt responsible for the metasomatic Fe-enrichment recorded in wehrlite. In this context, the more pervasive development of spongy clinopyroxene in wehrlite may be explained by a higher concentration of the evolved fluid phase at proximity to its silicate melt source. Received March 15, 2000; revised version accepted September 6, 2001