Three-dimensional linear solution for wave propagation with slopingbottom

This paper presents a three-dimensional analytic linear wave solution for surface gravity wave propagation over a sloping bottom that is valid for small, but realistic, slopes. The sloping-bottom linear model is compared to published laboratory data and to predictions of two-dimensional, constant-bottom nonlinear theories. The model is shown to describe the measured wave-height growth in the wave transformation region up to a limiting local Ursell number U_r of 0.35-1.0, depending on the wave type, although, as a linear model, it does not predict the harmonics observed in that range. For U_r<0.35, the harmonics can generally be neglected and the sloping-bottom linear theory agrees closely with both the published wave-height data and third-order Stokes nonlinear theory. As a three-dimensional linear model, superposition can be invoked to synthesize and relate wave structure in the transformation region to complex incident ocean spectra with both wind wave and swell components that arrive with a range of incidence angles. As such, the sloping-bottom linear model presented here should be a convenient useful tool for ocean modeling through a significant portion of the wave transformation region     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.     

Depositional and post‐depositional history of warm stage deposits at Knocknacran,Co. Monaghan,Ireland: implications for preservation of Irish last interglacial deposits

Organic‐rich deposits, uncovered during overburden removal from mantled gypsum karst at Knocknacran opencast gypsum mine, Co. Monaghan, are the best candidate to date for a last interglacial record in Ireland. The two till and organic‐rich deposits (preserved at different quarry elevations) were emplaced on to a Tertiary dolerite surface during high‐energy flood events and subsequently folded and faulted by movement towards sinkholes in underlying gypsum. Uranium–thorium disequilibrium dating suggests that the organic‐rich deposits in the upper section were hydrologically isolated at ca. 41 ka and those in the lower section at ca. 86 ka. Interpretation of the pollen content, although tentative because of the depositional and post‐depositional history of the material, suggests that the organic material originated in a warm stage possibly warmer than the post‐Eemian interstadials. The unusual setting of preservation may indicate that in situ, last interglacial deposits have generally been removed by erosion in Ireland. Copyright © 2004 John Wiley & Sons, Ltd.     

Stable isotope geochemistry and structural elements of the Sabina region (Central Apennines, Italy)

Various tectonic features are present in the Meso-Cenozoic basin units of the Sabina region (Central Apennines, Italy): Mio-Pliocene northeasterly verging thrusts are followed by Plio-Pleistocene, N-S oriented right-lateral strike-slip faults. Stable isotope geochemistry and examination of meso- and microstructures show that strain conditions differed through the course of tectonic history. Carbon and oxygen isotope analyses of the calate-filled extensional fractures, the sigmoidal veins present between stylolitic cleavage surfaces, and fault plane surfaces with differing motion, demonstrate those different geneses.
The "C/'" C of the older calcite-filling fractures (present both in the thrust and the strike-slip systems) suggests a deposition from shallow, fresh water circulation. Furthermore, the calcite fill of en echelon systems, that occur in the southernmost Plio-Pleistocene units, is clearly the result of a more recent, right-lateral strike-slip movement, connected with shallow water circulation within Mesozoic limestones.
The sigmoidal vein fills are derived from solid-state pressure solution processes which were the result of strike-slip movement. The deformation pattern related to the older thrust system is similar, but less intense; this also demonstrates general recrystallization processes in a closed system.
This suggests that the total shortening of the deformed sections is lower than that obtained on the basis of solution on stylolitic planes, because a sigruficant volume of dissolved carbonates remained in the system.
Stable isotope analysis also confirms that the deformational history of strongly cleaved rocks in the Sabina region took place in two phases and that extensional fractures formed before stylolithic planes, as suggested by structural and field observations.     

No major thermal event on the mid-Cretaceous Côte d'Ivoire–Ghana Transform Margin

Transient mid-Cretaceous thermal uplift induced by lateral heating from passing oceanic lithosphere is often invoked as a mechanism for the formation of the Côte d'Ivoire–Ghana basement ridge in the Equatorial Atlantic. This heating event should have affected mid-Cretaceous sedimentary rocks along the ridge. However, organic maturity and clay mineral data on the thermal evolution of these rocks suggest that burial temperatures did not exceed 80 °C and that palaeo–geothermal gradients are not anomalous. Optical petrography and the stratigraphic pattern of temperature-sensitive parameters indicate that higher palaeotemperature estimates are related to admixtures of preheated, detrital organic and inorganic matter. Erosion brought the sediments to their present shallow burial depths. Lack of evidence for significant thermal alteration implies that either thermal exchange between oceanic and continental lithosphere along the Côte d'Ivoire–Ghana Transform Margin was negligible, or that lateral heating by oceanic lithosphere was not strong enough to affect the sedimentary cover of the basement ridge.     

40Ar/39Ar and oxygen isotope studies of polymetamorphism from Tinos Island, Cycladic blueschist belt, Greece

Abstract Petrological, oxygen isotope and ⁴⁰Ar/³⁹Ar studies were used to constrain the Tertiary metamorphic evolution of the lower tectonic unit of the Cyclades on Tinos. Polyphase high-pressure metamorphism reached pressures in excess of 15 kbar, based on measurements of the Si content in potassic white mica. Temperatures of 450–500° C at the thermal peak of high-pressure metamorphism were estimated from critical metamorphic assemblages, the validity of which is confirmed by a quartz–magnetite oxygen isotope temperature of 470° C. Some ⁴⁰Ar/³⁹Ar spectra of white mica give plateau ages of 44–40 Ma that are considered to represent dynamic recrystallization under peak or slightly post-peak high-pressure metamorphic conditions. Early stages in the prograde high-pressure evolution may be documented by older apparent ages in the high-temperature steps of some spectra. Eclogite to epidote blueschist facies mineralogies were partially or totally replaced by retrograde greenschist facies assemblages during exhumation. Oxygen isotope thermometry of four quartz–magnetite pairs from greenschist samples gives temperatures of 440–470° C which cannot be distinguished from those deduced for the high-pressure event. The exhumation and overprint is documented by decreasing ages of 32–28 Ma in some greenschists and late-stage blueschist rocks, and ages of 30–20 Ma in the lower temperature steps of the Ar release patterns of blueschist micas. Almost flat parts of Ar–Ar release spectra of some greenschist micas gave ages of 23–21 Ma which are assumed to represent incomplete resetting caused by a renewed prograde phase of greenschist metamorphism. Oxygen isotope compositions of blueschist and greenschist facies minerals show no evidence for the infiltration of a δ¹⁸O-enriched fluid. Rather, the compositions indicate that fluid to rock ratios were very low, the isotopic compositions being primarily controlled by those of the protolith rocks. We assume that the fundamental control catalysing the transformation of blueschists into greenschists and the associated resetting of their isotopic systems was the selective infiltration of metamorphic fluid. A quartz–magnetite sample from a contact metamorphic skarn, taken near the Miocene monzogranite of Tinos, gave an oxygen isotope temperature of 555° C and calculated water composition of 9.1%. The value of δ¹⁸O obtained from this water is consistent with a primary magmatic fluid, but is lower than that of fluids associated with the greenschist overprint, which indicates that the latter event cannot be directly related to the monozogranite intrusion.     

Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.