Effects of fluids on the interaction of granites with limestones: The Notch Peak stock,Utah

In this paper we consider the mechanisms by which the mineralogy and composition of the margins of the Notch Peak granitic stock, Utah, were affected by calcareous argillite country rocks. The contact zone of the granite relative to the unaffected granite is enriched in: K₂O from about 4 to 10 wt.%, Rb from 250 to 510 ppm, Sr from 150 to 790 ppm and Ba from 250 to 2160 ppm. Locally, some of the intrusive rocks at the contact are nearly devoid of quartz and can be classified as syenites. The initial ⁸⁷Sr/⁸⁶Sr ratios range from 0.7069 in the unaffected rocks to 0.7100 in the potassium-enriched samples, approaching the values of the calc-silicate country rocks of about 0.7110.Calculations show that the characteristics of the contact zone near the top of the stock are the result of a number of sequential processes. CO₂-rich fluids derived from the heated calcareous argillites apparently caused a shift in the phase boundaries in the magma, enhancing accumulation of K-feldspar and high-Ca augite at the expense of other phases. The accumulation resulted in the high Ba and Sr concentrations in some samples. However, the high K₂O and Rb concentrations and magmatic ¹⁸O values indicate infiltration of magmatic fluid emanating from the solidifying lower portions of the pluton subsequent to solidification of the cap. The minimum fluid-rock ratios of 4.6 by mass, calculated on the basis of K₂O and Rb concentration shifts, indicate that a substantial amount of the fluid was channeled through this contact zone. The desilication of the rocks forming the syenitic samples at the contact apparently occurred when a chemical potential gradient of silica between the granite and wall-rocks was established as quartz was consumed in the wall-rocks during calc-silicate reactions. The infiltrating magmatic fluid probably acted as a medium for transport of silica across the contact and perhaps exchange of Sr between the country rocks and the intrusion where up to 30% of strontium in the granitic and syenitic samples from the contact zone was derived from the calc-silicates. The syenitic rocks cannot be explained by desilication reactions involving assimilation of the calc-silicates by the granite magma. The results of this study show that fluids interacting with the country rocks need to be considered to explain the effects of country rocks on the composition of the margins of granitic intrusions.     

Generalized thin sheet analysis in magnetotellurics: an extension of Price's analysis

Summary. Price's thin sheet analysis for electromagnetic fields has been extended in order to model the effects of crustal resistivity and conductivity variations on magnetotelluric fields. These extensions allow for a general layered medium below the crust and also account for the vertical resistance of the crust as well as its horizontal conductance. An important parameter emerges from the analysis which determines the distance it takes for the crustal current levels to readjust to changes in the crustal conductance. This adjustment distance is given by the square root of the conductivity thickness product multiplied by the resistivity thickness product. Approximate analytical solutions were developed for two-dimensional geometries in order to demonstrate these effects as well as the modifications produced by finite source wavelengths.     

Petrochemistry of the Fish Cove rhyolite,Keweenaw peninsula,Michigan, U.S.A.

A small (360 × 180 m) rhyolitic intrusive body in the lower portion of the Portage Lake Lava Series of Michigan's Keweenaw peninsula was mapped and sampled in detail. The rhyolite is one of a number of similar bodies which make up less than 1% of the total volume of this thick Late-Precambrian plateau basalt pile. The rock is a low-calcium rhyolite with fine-grained homogeneous texture and sparse phenocrysts of plagioclase and quartz. Analyses of selected trace and major elements for 21 samples taken from the body reveal a chemical zonation consisting of a core zone enriched in K, Rb and Ba, and a border zone relatively poor in these elements. Little areal difference is found with respect to other elements tested (Mn, Sr, Zr, Ca, Ti, and Fe). This apparently primary zonation seems to result from the migration of K, Rb and Ba during crystallization of the shallow intrusive. Though zoned, the trace-element chemistry of the Fish Cove body is distinct from that of eight other rhyolites in the Portage Lake Lava Series, and suggests that fingerprinting by trace elements might be a fruitful method for identifying and correlating the sources of numerous rhyolitic pebbles in conglomerates interbedded with the basaltlava flows of the Portage Lake Series.     

Volume integrals of the products of spherical harmonics and their application to viscous dissipation phenomena in fluids

The equations governing the conversion of kinetic energy into heat in moving viscous media are formulated as volume integrals of products of spherical harmonics. Although the formulation of the fundamental equations is classical, difficulties in the integration of certain products of generalized spherical harmonics over a sphere have permitted heretofore the treatment of only two cases. The closed, form evaluation of eight fundamental types of definite integrals of the product of spherical harmonics, some of them new, or at least missing in the literature, makes possible for the first time the evaluation of these volume integrals in closed form for arbitrary order and index. Explicit details are given for the rates of energy dissipation produced by viscous motions characterized by spheroidal as well as toroidal symmetry.     

Future land cover and climate may drive decreases in snow wind-scour and transpiration,increasing streamflow at a Colorado,USA headwater catchment

Understanding how land cover change will impact water resources in snow-dominated regions is of critical importance as these locations produce disproportionate runoff relative to their land area. We coupled a land cover evolution model with a spatially explicit, physics-based, watershed process model to simulate land cover change and its impact on the water balance in a 5.0 km² headwater catchment spanning the alpine–subalpine transition on the Colorado Front Range. We simulated two potential futures both with greater air temperature (+4°C/century) and more precipitation (+15%/century, MP) or less precipitation (−15%/century, LP) from 2000 to 2100. Forest cover in the catchment increased from 72% in 2000 to 84% and 83% in 2050 and to 95% and 92% in 2100 for MP and LP, respectively. Surprisingly, increases in forest cover led to mean increases in annual streamflow production of 12 mm (6%) and 2 mm (1%) for MP and LP in 2050 with an annual control streamflow of 208 mm. In 2100, mean streamflow production increased by 91 mm (44%) and 61 mm (29%) for MP and LP. This result counters previous work as runoff production increased with forested area due to decreases in snow wind-scour and increases in drifting leeward of vegetation, highlighting the need to better understand the impacts of forest expansion on the spatial pattern of snow scour, deposition and catchment effective precipitation. Identifying the hydrologic response of mountainous areas to climate warming induced land cover change is critically important due to the potential water resources impacts on downstream regions.     

An error analysis of profile flux,stability, and roughness length measurements made in the marine atmospheric surface layer

An analytical error analysis of profile-derived fluxes of heat, moisture, and momentum, along with stability and roughness length, is performed using the accuracies of the constituent temperature, humidity, and wind speed measurements. Five experiments, representing more than two thirds of the existing marine profile data presently contained in the literature, are compared. Much of the profile data examined was used to develop the transfer coefficients presently employed by a large number of competing bulk aerodynamic flux schemes. Depending upon the experiment, typical profile-method measurement errors were found to range from 15 to 35% for a sensible heat flux of ± 10 W m^-2, from 15 to 105% for a latent heat flux of ± 100 W m^-2, from 10 to 40% for a stress of 0.05 N m^-2, from 15 to 60% for a Monin-Obukhov stability of ± 0.05, and from 25 to 100% for a roughness length of 2 × 10^-4 m. Smaller magnitude flux values were found to contain typical errors as large as 100% for sensible heat flux, 300% for latent heat flux, and 60% for stress.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

Dual extraction ofR-mode andQ-mode factor solutions

It is mathematically possible to extract both R-mode and Q-mode factors simultaneously (RQ-mode factor analysis)by invoking the Eckhart-Young theorem. The resulting factors will be expressed in measures determined by the form of the scalings that have been applied to the original data matrix. Unless the measures for both solutions are meaningful for the problem at hand, the factor results may be misleading or uninterpretable. Correspondence analysis uses a symmetrical scaling of both rows and columns to achieve measures of proportional similarity between objects and variables. In the literature, the resulting similarity is a χ ² distance appropriate for analysis of enumerated data, the original application of correspondence analysis. Justification for the use of this measure with interval or ratio data is unconvincing, but a minor modification of the scaling procedure yields the profile similarity, which is an appropriate measure. Symmetrical scaling of rows and columns is unnecessary for RQ-mode factor analysis. If the data are scaled so the minor product W'Wis the correlation matrix, the major product WW'is expressed in the Euclidean distances between objects. Therefore, RQ-mode factor analysis can be performed so that the Rmode is a principal components solution and the Qmode is a principal coordinates solution. For applications where the magnitudes of differences are important, this approach will yield more interpretable results than will correspondence analysis.     

On Christoffel roots for nondetached slowness surfaces

The only restriction on the values of the elasticity parameters is the stability condition. Within this condition, we examine the Christoffel equation for nondetached qP slowness surfaces in transversely isotropic media. If the qP slowness surface is detached, each root of the solubility condition corresponds to a distinct smooth wavefront. If the qP slowness surface is nondetached, the roots are elliptical but do not correspond to distinct wavefronts; also, the qP and qSV slowness surfaces are not smooth.