Transient Behaviour of Tropospheric Ozone Precursors in a Global 3-D CTM and Their Indirect Greenhouse Effects

The global three-dimensional Lagrangian chemistry-transport model STOCHEM has been used to follow the changes in the tropospheric distributions of the two major radiatively-active trace gases, methane and tropospheric ozone, following the emission of pulses of the short-lived tropospheric ozone precursor species, methane, carbon monoxide, NO_x and hydrogen. The radiative impacts of NO_x emissionswere dependent on the location chosen for the emission pulse, whether at the surface or in the upper troposphere or whether in the northern or southern hemispheres. Global warming potentials were derived for each of the short-lived tropospheric ozone precursor species by integrating the methane and tropospheric ozone responses over a 100 year time horizon. Indirect radiative forcing due to methane and tropospheric ozone changes appear to be significant for all of the tropospheric ozone precursor species studied. Whereas the radiative forcing from methane changes is likely to be dominated by methane emissions, that from tropospheric ozone changes is controlled by all the tropospheric ozone precursor gases, particularly NO_xemissions. The indirect radiative forcing impacts of tropospheric ozone changes may be large enough such that ozone precursors should be considered in the basket of trace gases through which policy-makers aim to combat global climate change.     

An evaluation of methodologies for the generation of stochastic hydraulic conductivity fields in highly heterogeneous aquifers

 Stochastic techniques, such as Monte Carlo experiments, are more and more frequently used for the study of flow and transport in heterogeneous aquifers. When the aquifer is composed of distinct hydrofacies, a common way to model heterogeneity is to first generate equally-possible hydrofacies fields, and then convert these hydrofacies fields into hydraulic conductivity (K) fields by assigning a single K value to each facies. This technique assumes relative homogeneity of K within each facies but may not be appropriate for the most conductive facies that often exhibits substantial variability. In this paper, we assessed the impacts of assigning multiple random K, rather than a uniform K value, to the highly conductive facies on the results of a flow and transport model. A set of fifty stochastic hydrofacies maps depicting an environment similar to the Snake River Plain aquifer (SRPA) in south-east Idaho were generated. Simulations demonstrated that a uniform K value, if carefully chosen, can reasonably reproduce the specific discharges and early particle arrival times produced by multiple K values. Yet, the results obtained with a uniform K value are dramatically less variable than those obtained with multiple K values. It is therefore concluded that stochastic simulations with uniform K assigned to the most conductive and variable facies do not necessarily portray the entire uncertainty in the analyses.     

Dolomitization of the Waulsortian Limestone (Lower Carboniferous) in the Irish Midlands

The Waulsortian Limestone (Lower Carboniferous) of the southern Irish Midlands is dolomitized pervasively over a much larger region than previous studies have documented. This study indicates a complex, multistage, multiple fluid history for regional dolomitization. Partially and completely dolomitized sections of Waulsortian Limestones are characterized by finely crystalline (0·01–0·3 mm) planar dolomite. Planar replacive dolomite is commonly followed by coarse (≥0·5 mm) nonplanar replacive dolomite, and pervasive void‐filling saddle dolomite cement is frequently associated with Zn–Pb mineralization. Planar dolomite has average δ¹⁸O and δ¹³C values (‰ PDB) of –4·8 and 3·9 respectively. These are lower oxygen and slightly higher carbon isotope values than averages for marine limestones in the Waulsortian (δ¹⁸O=–2·2, δ¹³C=3·7). Mean C and O isotope values of planar replacive dolomite are also distinct from those of nonplanar and saddle dolomite cement (–7·0 and 3·3; –7·4 and 2·4 respectively). Fluid inclusions indicate a complex history involving at least three chemically and thermally distinct fluids during dolomite cementation. The petrography and geochemistry of planar dolomites are consistent with an early diagenetic origin, possibly in equilibrium with modified Carboniferous sea water. Where the Waulsortian was exposed to hydrothermal fluids (70–280 °C), planar dolomite underwent a neomorphic recrystallization to a coarser crystalline, planar and nonplanar dolomite characterized by lower δ¹⁸O values. Void‐filling dolomite cement is isotopically similar to nonplanar, replacive dolomite and reflects a similar origin from hydrothermal fluids. This history of multiple stages of dolomitization is significantly more complex than earlier models proposed for the Irish Midlands and provides a framework upon which to test competing models of regional vs. localized fluid flow.     

Partial melting in the Inzie Head gneisses: the role of water and a petrogenetic grid in KFMASH applicable to anatectic pelitic migmatites

A sequence of prograde isograds is recognized within the Dalradian Inzie Head gneisses where pelitic compositions have undergone variable degrees of partial melting via incongruent melting reactions consuming biotite. Three leucosome types are identified. At the lowest grades, granitic leucosomes containing porphyroblasts of cordierite (CRD‐melt) are abundant. At intermediate grades, CRD‐melt mingles with garnetiferous leucosomes (GT‐melt). At the highest grades, CRD‐melt coexists with orthopyroxene‐bearing leucosomes (OPX‐melt), while garnet is conspicuously absent. The prograde metamorphic field gradient is constrained to pressures of 2–3 kbar below the CRD‐melt isograd, and no greater than 4.5 kbar at the highest grade around Inzie Head. A petrogenetic grid, calculated using thermocalc , is presented for the K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH) system for the phases orthopyroxene, garnet, cordierite, biotite, sillimanite, H₂O and melt with quartz and K‐feldspar in excess. For the implied field gradient, the reaction sequence predicted by the grid is consistent with the successive prograde development of each leucosome type. Compatibility diagrams suggest that, as anatexis proceeded, bulk compositions may have been displaced towards higher MgO content by the removal of (relatively) ferroan granitic leucosome. An isobaric (P = 4 kbar) T–a_H2O diagram shows that premigmatization fluids must have been water‐rich (a_H2O > 0.85) and suggests that, following the formation of small volumes of CRD‐melt, the system became fluid‐absent and melting reactions buffered a_H2O to lower values as temperatures rose. GT‐ and OPX‐melt formed by fluid‐absent melting reactions, but a maximum of 7–11% CRD‐melt fraction can be generated under fluid‐absent conditions, much less than the large volumes observed in the field. There is strong evidence that the CRD‐melt leucosomes could not have been derived by buoyantly aided upwards migration from levels beneath the migmatites. Their formation therefore required a significant influx of H₂O‐rich fluid, but in a quantity insufficient to have exhausted the buffering capacity of the solid assemblage plus melt. Fluid : rock ratios cannot have exceeded 1 : 30. The fluid was channelled through a regionally extensive shear zone network following melt‐induced failure. Such an influx of fluid at such depths has obvious consequences for localized crustal magma production and possibly for cordierite‐bearing granitoids in general.     

Abrupt change in magma generation processes across the Central American arc in southeastern Guatemala: flux-dominated melting near the base of the wedge to decompression melting near the top of the wedge

Lavas erupted behind the volcanic front in southeastern Guatemala have many important distinctions from lavas erupted on the volcanic front. These include: generally higher MgO, Nb, Sr, TiO₂, and rare earth element concentrations; higher La/Yb and Nb/Y ratios; and lower Ba/La, La/Nb, Ba/Zr and Zr/Nb ratios. These major and trace element distinctions are caused by reduced fractionation during ascent and storage in the crust, lower degrees of melting in the source, and greatly reduced contributions from the subducted Cocos plate in the source. In addition, because all of these important distinctions are even borne in lavas erupted within 20 km of the front, there is little apparent petrogenetic continuity between front and behind-the-front magmas. What little geochemical continuity exists is in radiogenic isotopes: ¹⁴³Nd/¹⁴⁴Nd falls across the arc, Pb isotopic ratios (except ²⁰⁶Pb/²⁰⁴Pb) rise across the arc, and ⁸⁷Sr/⁸⁶Sr rise across the arc after an initial discontinuity within 20 km of the front. These continuous across-arc changes in radiogenic isotopes are caused by increased contamination with older, more isotopically disparate rocks, away from the front. Once the effects of crustal contamination are removed, the remaining isotopic variability behind the front is non-systematic and reflects the inherent isotopic heterogeneity of the source, the mantle wedge. Geochemical disconnection in southeastern Guatemala suggests that behind-the-front magmas are produced by decompression melting near the top of the wedge, not by flux-dominated melting near the base of the wedge.     

Seismic reflection survey of an oblique aseismic basement trend on the Reykjanes Ridge

A detailed (5 km track separation) seismic reflection survey of a portion of the upper flank of Reykjanes Ridge supports the existence of an oblique aseismic ridge, previously postulated from other data. The oblique basement ridge may have been formed by a magma center moving southwest under this portion of the Reykjanes Ridge at about 6 cm/yr between 7 and 5 mbyp. The oblique ridge is complex, being interrupted by saddles about every 30 km length. This spacing could reflect incipient, very weakly developed transverse fractures, or more probably the concentration of volcanic activity at particularly active vents, which shift southwestward every million years or so in response to the south-westward moving magma chambers entrained in the asthenosphere. Minor irregularities in the oblique ridge parallel crustal isochrons; such small features are probably elongate fissure eruptions restricted to a narrow spreading axis.     

Mechanics of growth of some laccolithic intrusions in the Henry mountains, Utah, I: Field observations, Gilbert's model, physical properties and flow of the magma

The shapes of sills and laccolithic intrusions and associated host rock deformation were studied at several locations on the flanks of the Henry Mountains. Diorite sills range from 0.5 to 10 m in thickness, are less than 1 km² in areal extent, and have blunt terminations. The laccolithic intrusions range from 10 to 200 m in thickness, and from 1 to 3 km² in areal extent. The host rock, principally sandstone and shale, is deformed along closely spaced cataclastic shear planes. This deformation is concentrated at contacts, especially near sill terminations and over laccolith peripheries. The diorite contains plagioclase phenocrysts which are usually sheared in a thin zone adjacent to each contact. Field observations suggest that sills are the forerunners of laccolithic intrusions which form only after magma has spread far enough laterally (greater than about 1 km² in the Henry Mountains) to gain leverage to bend the overburden upward. Further injection of magma results in laccolithic peripheries or terminations with one of three distinct cross-sectional forms: (1) blunt termination of the diorite accompanied by bending and minor faulting of the host rock; (2) termination at a peripheral diorite dike cutting upward across the host rock; or (3) abrupt termination of the diorite against a nearly vertical fault zone.In order to study some of the processes of sill and laccolith intrusion, mechanical models for the driving pressure, physical properties, and flow behavior of the diorite magma are derived and discussed. A static driving pressure (equal to the difference between total magma pressure and lithostatic pressure) of up to 700 bar is estimated. The rheological behavior of the magma in the Henry Mountains is unknown. However, flow behavior is calculated assuming three of the more common models for fluids: Newtonian viscous, pseudoplastic, and Bingham. Suspended crystals probably contributed to the finite strength of the magma (estimated to be at least 10³ dyn/cm² for the Henry Mountains magma) which enables it to support dense zenoliths and also fixes maximum limits on the lengths of sills or dikes. Pressure in magma flowing along tabular intrusions of uniform thickness drops linearly in the flow direction for all three rheological materials. Thickening of tabular intrusions tends to make the pressure drop less rapidly, but pressure drops more rapidly in the tapered region near a termination. Pressure distributions under these and other conditions are derived in order to use them in the models of host rock deformation presented in Part II.     

An experiment on reservoir representation schemes to improve hydrologic prediction: coupling the national water model with the HEC-ResSim

ABSTRACT

This study experiments with reservoir representation schemes to improve the ability to model active water management in the National Water Model (NWM). For this purpose, we developed an integrated water management model, NWM-ResSim, by coupling the NWM with HEC-ResSim, and two reservoir representation schemes are tested: simulation of reservoir operations and retrieval of scheduled operations. The experiments focus on a pilot reservoir domain in the Russian River basin – Lake Mendocino, California – and its contributing watershed. The evaluation results suggest that the NWM-ResSim improves the simulation performance of reservoir outflow from this managed reservoir over the NWM default level pool routing scheme. The degree of this improvement depends on the suitability of the operation guidance; the reservoir operations simulation scheme could have acceptable errors for the purposes of water resources management, but not for flood operations. Results of the retrieval scheme of scheduled operations demonstrated better performance for sub-daily flood operations.     

On the application of the mode-acceleration method to structural engineering problems

Mode-superposition has been extensively used in computing the dynamic response of complex structures. Two versions of mode-superposition, namely the mode-displacement method and the mode-acceleration method, have been employed. The present paper summarizes the results of a systematic study comparing the accuracy of the mode-displacement and mode-acceleration methods when applied to structures with various levels of damping or various excitation frequencies. The paper also discusses several details concerning the implementation of the mode-acceleration method.