Fracture sealing of quartzite under a temperature gradient: experimental results

Veined lithologies are formed by fracturing and sealing processes, with the veins representing former fluid conduits through the rock. Although detailed fieldwork and numerical simulations have provided a better understanding of vein growth, few studies have attempted to seal fractures and generate veins experimentally. In this pilot study, we subjected fractured quartzite to temperature gradients of 45–125 °C under hydrothermal conditions in a static fluid, with the aim of precipitating secondary quartz in the cooler portions of the fracture. Results show that secondary quartz precipitates due to the imposed temperature gradient, causing the initial fracture to seal locally. Although no systematic sealing pattern was observed along the fracture, samples subjected to higher temperatures exhibit a smaller fracture width and appear to have reacted more extensively. Electron microprobe mapping visualizes the spatial distribution of secondary quartz, which contains elevated concentrations of aluminium.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Tsunami Probability in the Caribbean Region

We calculated tsunami runup probability (in excess of 0.5 m) at coastal sites throughout the Caribbean region. We applied a Poissonian probability model because of the variety of uncorrelated tsunami sources in the region. Coastlines were discretized into 20 km by 20 km cells, and the mean tsunami runup rate was determined for each cell. The remarkable ~500-year empirical record compiled by O’Loughlin and Lander (2003) was used to calculate an empirical tsunami probability map, the first of three constructed for this study. However, it is unclear whether the 500-year record is complete, so we conducted a seismic moment-balance exercise using a finite-element model of the Caribbean-North American plate boundaries and the earthquake catalog, and found that moment could be balanced if the seismic coupling coefficient is c = 0.32. Modeled moment release was therefore used to generate synthetic earthquake sequences to calculate 50 tsunami runup scenarios for 500-year periods. We made a second probability map from numerically-calculated runup rates in each cell. Differences between the first two probability maps based on empirical and numerical-modeled rates suggest that each captured different aspects of tsunami generation; the empirical model may be deficient in primary plate-boundary events, whereas numerical model rates lack backarc fault and landslide sources. We thus prepared a third probability map using Bayesian likelihood functions derived from the empirical and numerical rate models and their attendant uncertainty to weight a range of rates at each 20 km by 20 km coastal cell. Our best-estimate map gives a range of 30-year runup probability from 0–30% regionally.     

Future climate in the Pacific Northwest

Climate models used in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) on the whole reproduce the observed seasonal cycle and twentieth century warming trend of 0.8°C (1.5°F) in the Pacific Northwest, and point to much greater warming for the next century. These models project increases in annual temperature of, on average, 1.1°C (2.0°F) by the 2020s, 1.8°C (3.2°F) by the 2040s, and 3.0°C (5.3°F) by the 2080s, compared with the average from 1970 to 1999, averaged across all climate models. Rates of warming range from 0.1°C to 0.6°C (0.2°F to 1.0°F) per decade. Projected changes in annual precipitation, averaged over all models, are small (+1% to +2%), but some models project an enhanced seasonal cycle with changes toward wetter autumns and winters and drier summers. Changes in nearshore sea surface temperatures, though smaller than on land, are likely to substantially exceed interannual variability, but coastal upwelling changes little. Rates of twenty-first century sea level rise will depend on poorly known factors like ice sheet instability in Greenland and Antarctica, and could be as low as twentieth century values (20 cm, 8^″) or as large as 1.3 m (50^″).     

Genesis and solvus relations of submicroscopically intergrown paragonite and phengite in a blueschist from northern California

Electron microbeam techniques have been used to examine submicroscopically intergrown paragonite, phengite and chlorite from the South Fork Mountain Schist of the Franciscan Terrane of northern California, which was subjected to blueschist facies metamorphism. The sample also contains quartz, albite, lawsonite, and rutile. The subassemblage albite-lawsonite-rutile requires metamorphic conditions on the low-temperature side of the equilibrium albite+lawsonite+rutile=paragonite+sphene+quartz+H₂O (T<200° C and P<7.4 kbars based on thermodynamic data of Holland and Powell 1990). The white micas appear to be optically homogeneous, but back-scattered electron images can distinguish two different micas by their slight difference in contrast. Electron microprobe analyses (EMPA) of micas show Na/(Na+K) ranging from 0.2 to 0.8. The two micas are resolved by transmission electron microscopy (TEM) as packets of phengite and paragonite that range from 20 to several hundred nm in thickness. The compositions, determined by analytical electron microscopy (AEM), constrain the limbs of the phengite-paragonite solvus to values of Na/(Na+K)=<0.02 and 0.97, representing less mutual solid solution than ever reported by EMPA. The textural relations imply that the sheet silicates were derived from reactions between fluids and detrital clays and that they are in an intermediate stage of textural development. We caution that microprobe analyses of apparently homogeneous sheet silicates may yield erroneous data and lead to faulty conclusions using phengite barometry and paragonite-muscovite thermometry, especially in fine-grained rocks that formed at relatively low temperatures. Contribution no. 473 from the Mineralogical Laboratory, Department of Geological Sciences, The University of Michigan, Ann Arbor, Michigan, USA     

Experimental investigation and application of the equilibrium rutile + orthopyroxene = quartz + ilmenite

Equilibria in the Sirf (Silica-Ilmenite-Rutile-Ferrosilite) system: $${\text{SiO}}_{\text{2}} + ({\text{Mg,Fe}}){\text{TiO}}_{\text{3}} {\text{ + (Mg,Fe)SiO}}_{\text{3}} $$ have been calibrated in the range 800–1100° C and 12–26 kbar using a piston-cylinder apparatus to assess the potential of the equilibria for geobarometry in granulite facies assemblages that lack garnet. Thermodynamic calculations indicate that the two end-member equilibria involving quartz + geikielite = rutile + enstatite, and quartz + ilmenite = rutile + ferrosilite, are metastable. We therefore reversed equilibria over the compositional range Fs_40–70, using Ag₈₀Pd₂₀ capsules with \(f_{{\text{O}}_{\text{2}} } \) buffered at or near iron-wüstite. Ilmenite compositions coexisting with orthopyroxene are \(X_{{\text{MgTiO}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.06 to 0.15 and \(X_{{\text{Fe}}_{\text{2}} {\text{O}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.00 to 0.01, corresponding toK _D values of 13.3, 10.2, 9.0 and 8.0 (±0.5) at 800, 900, 1000 and 1100° C, respectively, whereK _D =(XMg/XFe)^Opx/(XMg/XFe)^Ilm. Pressures have been calculated using equilibria in the Sirf system for granulites from the Grenville Province of Ontario and for granulite facies xenoliths from central Mexico. Pressures are consistent with other well-calibrated geobarometers for orthopyroxeneilmenite pairs from two Mexican samples in which oxide textures appear to represent equilibrium. Geologically unreasonable pressures are obtained, however, where oxide textures are complex. Application of data from this study on the equilibrium distribution of iron and magnesium between ilmenite and orthopyroxene suggests that some ilmenite in deep crustal xenoliths is not equilibrated with coexisting pyroxene, while assemblages from exposed granulite terranes have reequilibrated during retrogression. The Sirf equilibria are sensitive to small changes in composition and may be used for determination of activity/composition (a/X) relations of orthopyroxene if an ilmenite model is specified. A symmetric regular solution model has been used for orthopyroxene in conjunction with activity models for ilmenite available from the literature to calculatea/X relations in orthopyroxene of intermediate composition. Data from this study indicate that FeSiO₃?MgSiO₃ orthopyroxene exhibits small, positive deviations from ideality over the range 800–1100°C.     

Global change and the intensification of agriculture in the tropics

Bridging understanding of local environmental change with regional and global patterns of land-use and land-cover change (LUCC) remains a key goal and challenge for our understanding of global environmental change. This meta-analysis attempts to bridge local and regional scales of LUCC by demonstrating the ways in which previously published case studies can be compared and used for a broader regional synthesis in the tropics. In addition to providing results from a meta-analysis, this paper suggests ways to make future case studies more widely comparable.     

Atmospheric forcing of the eastern tropical Pacific: A review

The increase in marine, land surface, atmospheric and satellite data during recent decades has led to an improved understanding of the air–sea interaction processes in the eastern tropical Pacific. This is also thanks to extensive diagnoses from conceptual and coupled ocean–atmosphere numerical models. In this paper, mean fields of atmospheric variables, such as incoming solar radiation, sea level pressure, winds, wind stress curl, precipitation, evaporation, and surface energy fluxes, are derived from global atmospheric data sets in order to examine the dominant features of the low level atmospheric circulations of the region. The seasonal march of the atmospheric circulations is presented to depict the role of radiative forcing on atmospheric perturbations, especially those dominating the atmosphere at low levels.In the tropics, the trade winds constitute an important north–south energy and moisture exchange mechanism (as part of the low level branch of the Hadley circulation), that determines to a large extent the precipitation distribution in the region, i.e., that associated with the Inter-Tropical Convergence Zone (ITCZ). Monsoonal circulations also play an important role in determining the warm season precipitation distribution over the eastern tropical Pacific through a large variety of air–sea–land interaction mechanisms. Westward traveling waves, tropical cyclones, low latitude cold air intrusions, and other synoptic and mesoscale perturbations associated with the ITCZ are also important elements that modulate the annual rainfall cycle. The low-level jets of the Gulf of California, the Intra-Americas Sea (Gulf of Mexico and Caribbean Sea) and Chocó, Colombia are prominent features of the eastern tropical Pacific low-level circulations related to sub-regional and regional scale precipitation patterns. Observations show that the Intra-Americas Low-Level Jet intensity varies with El Niño/Southern Oscillation (ENSO) phases, however its origin and role in the westward propagation and development of disturbances that may hit the eastern tropical Pacific, such as easterly waves and tropical cyclones, are still unclear. Changes in the intensity of the trade winds in the Caribbean Sea and the Gulf of Mexico (associated with eastern tropical Pacific wind jets) exert an important control on precipitation by means of wind–topography interactions. Gaps in the mountains of southern Mexico and Central America allow strong wind jets to pass over the continent imprinting a unique signal in sea surface temperatures and ocean dynamics of the eastern tropical Pacific.The warm pools of the Americas constitute an important source of moisture for the North American Monsoon System. The northeastern tropical Pacific is a region of intense cyclogenetic activity, just west of the coast of Mesoamerica. Over the oceanic regions, large-scale properties of key variables such as precipitation, moisture, surface energy fluxes and wind stress curl are still uncertain, which inhibits a more comprehensive view of the region and stresses the importance of regional field experiments. Progress has been substantial in the understanding of the ocean and atmospheric dynamics of the eastern tropical Pacific, however, recent observational evidence such as that of a shallow meridional circulation cell in that region, in contrast to the classic concept of the Hadley-type deep meridional circulation, suggests that more in situ observations to validate theories are still necessary.This paper is part of a comprehensive review of the oceanography of the eastern tropical Pacific Ocean.