Analytical Model for Permeability Evolution in Microcracking Rock

— Analytical expressions to predict the enhancement of permeability due to stress-induced microcracking in initially low porosity rock are presented. A fracture mechanical microcrack model is employed to derive integrated effective hydraulic variables as a function of stress, which are then used to calculate the evolution of permeability using the statistically-based Dienes model. The model enables determination of permeability enhancement as a function of two loading parameters and three material parameters. Results are in reasonable agreement with experimental measurements and indicate that appreciable increases in permeability can be anticipated during brittle failure. The analytical nature of the model makes it easily incorporatable into numerical models that require quantification of the permeability evolution as a function of stress, for which there is currently no law.     

Development of large‐scale seismites in Upper Cretaceous fluvial sandstones in a fault‐proximal setting

Large‐scale soft‐sediment deformation structures occur within fluvial sandstone bodies of the Upper Cretaceous Wahweap Formation in the Kaiparowits basin, southern Utah, USA. These structures represent an exceptional example of metre‐scale fault‐proximal, seismogenic load structures in nearly homogenous sandstones. The load structures consist of two types: large‐scale load casts and wedge‐shaped load structures. Large‐scale load casts penetrate up to 4·5 m into the underlying sandstone bed. Wedge‐shaped load structures include metre‐scale, parallel, sub‐vertical features and decimetre‐scale features along the periphery of the large‐scale load casts or other wedge‐shaped load structures. Wedge‐shaped load structures contain well‐developed, medial cataclastic shear deformation bands. All load structures contain pervasive well‐defined millimetre‐thick to centimetre‐thick internal laminae, oriented parallel to the outside form of the load structures and asymptotic to deformation bands. Both types of load structures formed because of an inverted density profile, earthquake‐triggered liquefaction and growth of irregularities (a Rayleigh–Taylor instability) on the sandstone–sandstone erosional contact. The internal laminae and deformation bands formed during deformation and clearly demonstrate polyphase deformation, recording a transition from liquefied to hydroplastic to brittle modes of deformation. Decimetre‐scale wedge‐shaped load structures on the edge of the large‐scale load casts probably formed towards the end of a seismic event after the sediment dewatered and increased the frictional contact of grains enough to impart strength to the sands. Metre‐scale wedge‐shaped load structures were created as the tips of downward foundering sediments were driven into fractures, which widened incrementally with seismic pulsation. With each widening of the fracture, gravity and a suction effect would draw additional sediment into the fracture. Superimposed laminae indicate a secondary syndeformational origin for internal laminae, probably by flow‐generated shearing and vibrofluidization mechanisms. Large‐scale and wedge‐shaped load structures, polyphase deformation and secondary laminae may characterize soft‐sediment deformation in certain fault‐proximal settings.     

Sublittoral seaweed communities on natural and artificial substrata in a high-latitude coral community in South Africa

Coral mortality may result in macroalgal proliferation or a phase shift into an alga-dominated state. Subtidal, high-latitude western Indian Ocean coral communities at Sodwana Bay on the KwaZulu-Natal coast, South Africa, have experienced some mortality because of warm-water anomalies, storms and other causes, but the response of the macroalgae is unknown. We investigated the abundance and diversity of benthic algae on different hard natural substrata (dead digitate, brain and plate corals and beach rock) on Two-Mile Reef, Sodwana Bay. We also compared algal communities colonising ceramic, marble and pretreated ceramic tiles placed on the reef for six months. We identified 95 algae (14 Chlorophyta, 11 Phaeophyceae, 69 Rhodophyta and one cyanobacterium). Assemblages on natural and artificial substrata were dominated by the brown alga Lobophora variegata (Lamouroux) Womersley ex Oliveira and non-geniculate corallines (Rhodophyta, Corallinaceae). Cluster and ordination analyses revealed that the algae showed no affinity for particular substrata, whether natural or artificial. Algal cover was occasionally higher on rougher tiles and crustose corallines were significantly more abundant on marble than ceramic tiles. Two-Mile Reef had 23.1% dead and 48.4% live scleractinian coral cover, where dead corals were colonised indiscriminately by many small algal species, but there was no evidence of algal proliferation. The results provide a baseline for monitoring this high-latitude reef system.     

Reynolds stress and TKE production in an estuary with a tidal bore

We report new measurements of the turbulent properties of the flow in a tidally energetic estuarine channel of almost uniform cross-section. A high-frequency (1.2 MHz), bottom-mounted Acoustic Doppler Current Profiler (ADCP) has been used to observe the velocity field at a sampling rate of 10 Hz in parallel with measurements of the surface elevation by tide gauges. Our data have been analysed using the Variance Method to determine turbulent kinetic energy (TKE), shear stress and TKE production over the tidal cycle with a time resolution of 60 s. During the highly energetic but brief flood period, when the surface axial velocity reaches 2 m s⁻¹, we observed large values of stress (>2 Pa) and shear production (5 W m⁻³). TKE is also input through the release of energy in the bore itself which results in a brief but intense injection of energy at the bore front with large transient TKE levels (100 J m⁻³). Subsequent input by shear production maintains TKE levels which are generally lower (20 J m⁻³) than the strong peak associated with the bore for the rest of the flood. On the ebb, the flow is relatively tranquil with maximum speeds 0.5 m s⁻¹ and peak TKE production rates of 0.1 W m⁻³.The flow and elevation data have also been used to estimate the energy fluxes into and out of the estuary. Short (1 h), intense energy inputs (8 MW at springs) on the flood flow are largely balanced by longer, less intense seaward energy flow on the ebb. The net energy input is found to be 0.1 MW at springs which is consistent with estimates of upstream dissipation. Peak dissipation in the bore itself may exceed the mean energy input but it is active only for a small fraction of the tidal cycle and its average contribution does not exceed 12% of total dissipation.     

Dynamic recrystallization and chemical evolution of clinoamphibole from Senja,Norway

Clinoamphibole from a mylonitic amphibolite exhibits microstructures characteristic of dynamic recrystallization, including porphyroclasts in a finer grained matrix of needle-shaped amphibole. The matrix amphibole defines an LS fabric and porphyroclasts have core and mantle structures with a core containing undulose to patchy extinction and (100) deformation twinning surrounded by a mantle of recrystallized grains. In addition intragranular grains also occur within the cores. TEM analyses of the porphyroclasts reveal that they contain a wide variety of lattice defects including high densities (5 × 10⁸cm^–2) of free dislocations and dislocation arrays, dissociated dislocations, stacking faults, and (100) micro-twins. TEM also shows that matrix grains and intragranular grains have relatively low defect densities, and that the intragranular new grains occur at localities in the porphyroclasts characterized by high densities of dislocations. These observations along with the chemical and orientation relationships between the recrystallized grains and porphyroclasts indicate that the new grains may have formed by heterogeneous nucleation and that further growth probably occurred by both strain assisted and chemically induced grain boundary migration or liquid film migration. This recrystallization event is interpreted to be synkinematic based on the fact that no recrystallization textures are present in the matrix grains and that the matrix grains define an LS fabric. However, the low defect densities in the matrix grains and the lack of intracrystalline strain in other phases indicate that post-kinematic recovery processes were active.     

Porphyroclast systems as kinematic indicators

Porphyroclasts of relatively strong minerals in mylonites commonly have an internal monoclinic shape symmetry defined by tails of dynamically recrystallized material. The geometry of a porphyroclast and its tails, called a ‘porhyroclast system’, can serve as a valuable indicator of the sense of vorticity. Porphyroclast systems have been divided into σ- and δ-types on the basis of the geometry of the tails. σ-Types have wedge-shaped recrystallized tails whose median lines lie on opposite sides of a reference plane parallel to the tails and containing the symmetry axis for the system. σ-Types are further subdivided into a σ_a-types, in which the porphyroclast is isolated in a relatively homogeneous matrix, and σ_b-types, in which the porphyroclast system is associated with a shear band foliation in the matrix. δ-Types typically have narrow recrystallized tails whose median lines cross the reference plane adjacent to the porphyroclast. Consequently, embayments of matrix material occur adjacent to the porphyroclasts and the tails display characteristic bends.A porphyroclast system in a mylonite develops when the relatively weak dynamically recrystallized grain aggregate in the porphyroclast mantle changes its shape due to non-coaxial flow in the adjacent matrix. This behaviour has been modelled in shear box experiments. Passive marker lines around rigid cylinders embedded in silicone putty were subjected to simple shear. The experiments were modified to simulate a change in recrystallization rate (R) with respect to rate of deformation (γ) by decreasing the diameter of the rigid cylinder during deformation at variable rates. The ratio R/γ appears to be one of the most important factors in determining which porphyroclast system will develop. At high R/γ values, flow of recrystallized material away from the porphyroclast is continuously appended by the production of new grains and wedge-shaped σ_a-type tails develop. At low R/γ values, relatively few new grains are added to the tails which become thinned and deflected by drag due to the spinning motion of the porphyroclast. In addition, most porphyroclast systems at low shear strains are of σ_a-type or lack monoclinic symmetry, whereas δ-types are only developed at high shear strain values. Complex porphyroclast systems, characterized by two generations of tails, are observed in many of the natural and model shear zones studied and may form due to fluctuating R/γ. Conditions that allow isolated σ_a- and δ-type porphyroclast systems to be used as sense of vorticity indicators are: the systems should have a monoclinic shape symmetry; matrix grain size should be small with respect to porphyroclast size; matrix fabric should be homogeneous; deformation history should be simple, and observations should be made on sections normal to the inferred bulk vorticity vector for the mylonite.     

The effects of a spillage of diesel fuel on a rocky shore in the sub-Antarctic region (Macquarie Island)

On 3 December 1987, the supply ship Nella Dan ran aground at Macquarie Island (54°29′S, 158°58′E) releasing about 270 000 I of oil, mostly light marine diesel, into the sea. At the time of the incident, many marine invertebrates were washed up dead along 2 km of shoreline. Twelve months later, the shore community was investigated using 1. algal and invertebrate populations of the littoral and sublittoral rocky shore, and 2. the invertebrate communities living in the holdfasts of the giant kelp Durvillaea antarctica, which were collected for later examination. Investigations were undertaken at both affected and control locations. Analyses of differences in community structure involved nested ANOVA and multidimensional scaling techniques. On the rocky substrate, the effect of the spill was restricted to some biota of the lower littoral and sublittoral zones—particularly echinoderms and the patellid limpet Nacella macquariensis. There were differences in cover for some algal species between locations. Within the kelp holdfasts, communities were dominated by peracarid crustaceans at control locations and by polychaetes (particularly the opportunistic groups—capitellids, cirratulids and spionids) at oil-affected locations.

The communities have recently been re-surveyed (in the summer of 1994–1995) to assist in the interpretation of the results and to gauge the extent of recovery of the affected biota.     

Measurements of Trace Gases in the Inflow of South China Sea Background Air and Outflow of Regional Pollution at Tai O, Southern China

We present a 16-month record of ozone (O₃), carbon monoxide (CO), total reactive nitrogen (NO_y), sulphur dioxide (SO₂), methane (CH₄), C₂ – C₈ non-methane hydrocarbons (NMHCs), C₁ – C₂ halocarbons, and dimethyl sulfide (DMS) measured at a southern China coastal site. The study aimed to establish/update seasonal profiles of chemically active trace gases and pollution tracers in subtropical Asia and to characterize the composition of the `background' atmosphere over the South China Sea (SCS) and of pollution outflow from the industrialized Pearl River Delta (PRD) region and southern China. Most of the measured trace gases of anthropogenic origin exhibited a winter maximum and a summer minimum, while O₃ showed a maximum in autumn which is in contrast to the seasonal behavior of O₃ in rural eastern China and in many mid-latitude remote locations in the western Pacific. The data were segregated into two groups representing the SCS background air and the outflow of regional continental pollution (PRD plus southern China), based on CO mixing ratios and meteorological conditions. NMHCs and halocarbon data were further analyzed to examine the relationships between their variability and atmospheric lifetime and to elucidate the extent of atmospheric processing in the sampled air parcels. The trace gas variability (S) versus lifetime (τ) relationship, defined by the power law, S_lnx = Aτ^{− b}, (where X is the trace gas mixing ratio) gives a fit parameter A of 1.39 and exponent b of 0.42 for SCS air, and A of 2.86 and b of 0.31 for the regional continental air masses. An examination of ln[n-butane]/ln[ethane] versus ln[propane]/ln[ethane] indicates that their relative abundance was dominated by mixing as opposed to photochemistry in both SCS and regional outflow air masses. The very low ratios of ethyne/CO, propane/ethane and toluene/benzene suggest that the SCS air mass has undergone intense atmospheric processing since these gases were released into the atmosphere. Compared to the results from other polluted rural sites and from urban areas, the large values of these species in the outflow of PRD/southern China suggest source(s) emitting higher levels of ethyne, benzene, and toluene, relative to light alkanes. These chemical characteristics could be unique indicators of anthropogenic emissions from southern China.     

Why the Sacramento Delta area differs from other parts of the great valley: Numerical modeling of thermal structure and thermal subsidence of forearc basins

Data on present-day heat flow, subsidence history, and paleotemperature for the Sacramento Delta region, California, have been employed to constrain a numerical model of tectonic subsidence and thermal evolution of forearc basins. The model assumes an oceanic basement with an initial thermal profile dependent on its age subjected to refrigeration caused by a subducting slab. Subsidence in the Sacramento Delta region appears to be close to that expected for a forearc basin underlain by normal oceanic lithosphere of age 150 Ma, demonstrating that effects from both the initial thermal profile and the subduction process are necessary and sufficient. Subsidence at the eastern and northern borders of the Sacramento Valley is considerably less, approximating subsidence expected from the dynamics of the subduction zone alone. These results, together with other geophysical data, show that Sacramento Delta lithosphere, being thinner and having undergone deeper subsidence, must differ from lithosphere of the transitional type under other parts of the Sacramento Valley. Thermal modeling allows evaluation of the rheological properties of the lithosphere. Strength diagrams based on our thermal model show that, even under relatively slow deformation (10⁻¹⁷ s⁻¹), the upper part of the delta crystalline crust (down to 20–22 km) can fail in brittle fashion, which is in agreement with deeper earthquake occurrence. Hypocentral depths of earthquakes under the Sacramento Delta region extend to nearly 20 km, whereas, in the Coast Ranges to the west, depths are typically less than 12–15 km. The greater width of the seismogenic zone in this area raises the possibility that, for fault segments of comparable length, earthquakes of somewhat greater magnitude might occur than in the Coast Ranges to the west. The text was submitted by the authors in English.     

Seismotectonic deformation during the filling of Toktogul Reservoir,Kirghizia

Natural seismicity and induced seismic events are discussed in the context of the seismotectonic deformation method through the example of earthquakes in the Toktogul Reservoir region in the Central Tien Shan, Kirghizia. The parameters of seismotectonic deformation of various sites in the Toktogul Reservoir region are described, based on a statistical study of focal mechanisms. The relationship of induced seismicity to changes of water level in the reservoir is reviewed. The temporal stress-strain characteristics are investigated. Average focal mechanisms for the entire region, as well as areas in the immediate vicinity of Toktogul Dam are analyzed. The vertical component of seismotectonic deformation varies in time from compression to extension-opposite to what is expected from the influence of the reservoir load; strike-slip motions become oblique thrusts. Changes in the orientation of focal mechanisms coincide with the time of maximum rate of the filling of the reservoir.