The solubility of quartz in aqueous sodium chloride solution at 350°C and 180 to 500 bars

The solubility of quartz in 2, 3, and 4 molal NaCl was measured at 350°C and pressures ranging from 180 to 500 bars. The molal solubility in each of the salt solutions is greater than that in pure water throughout the measured pressure range, with the ratio of solubility in NaCl solution to solubility in pure water decreasing as pressure is increased. The measured solubilities are significantly higher than solubilities calculated using a simple model in which the water activity in NaCl solutions decreases either in proportion to decreasing vapor pressure of the solution as salinity is increased or in proportion to decreasing mole fraction of water in the solvent.     

MELTWATER FLUXES AT AN ARCTIC FOREST-TUNDRA SITE

Models of surface energy balance and snow metamorphism are utilized to predict the energy and meltwater fluxes at an Arctic site in the forest–tundra transition zone of north-western Canada. The surface energy balance during the melt period is modelled using an hourly bulk aerodynamic approach. Once a snowcover becomes patchy, advection from the bare patches to the snow-covered areas results in a large spatial variation in basin snowmelt. In order to illustrate the importance of small-scale, horizontal advection, a simple parameterization scheme using sensible heat fluxes from snow free areas was tested. This scheme estimates the maximum horizontal advection of sensible heat from the bare patches to the snow-covered areas. Calculated melt was routed through the measured snowcover in each landscape type using a variable flow path, meltwater percolation model. This allowed the determination of the spatial variability in the timing and magnitude of meltwater release for runoff. Model results indicate that the initial release of meltwater first occurred on the shallow upland tundra sites, but meltwater release did not occur until nearly two weeks later on the deep drift snowcovers. During these early periods of melt, not all meltwater is available for runoff. Instead, there is a period when some snowpacks are only partially contributing to runoff, and the spatial variation of runoff contribution corresponds to landscape type. Comparisons of melt with and without advection suggests that advection is an important process controlling the timing of basin snowmelt.     

Sources of dissolved REE in mountainous streams draining granitic rocks, Sierras Pampeanas (Córdoba, Argentina)

Stream waters draining granitic terrains from the highest part (850 to 2200 m a.s.l.) of Sierras Pampeanas (Córdoba, Argentina, ∼32°S, ∼65°W) were sampled in order to define sources and distribution of dissolved rare earth elements (REE), and to describe the geochemical processes that govern their mobility. The contribution of the regional granite to the dissolved REE pool in stream water is limited due to the physical conditions predominating in the area (i.e., steep slopes and semiarid climate). Therefore, precipitation is considered a seasonally significant source controlling REE concentration in stream water. Dissolved REE concentrations are inversely correlated with monthly precipitation and rainfall frequency. During the rainy season (i.e., the austral summer) REE concentrations in stream water are lower than during the dry season (i.e., austral winter). Such low concentrations reflect the balance between the REE input from precipitation and their removal by adsorption. In contrast, during the dry season, the longer residence time of water within fractures and colluvium determines an increased REE concentration in the base flow. Lower pH values also contribute to raise REE concentration through desorption from mineral surfaces.     

On the Occurrence of Zellia from the Maping Limestone of Chengkung, Central Yunnan

INTRODUCTIONThe specimens described in the present paper were collected in 1939 by the late Mr. T. Y. Hsu and Mr. C. S. Pien from the Maping limestone exposed E of Tuikoshan, Chengkung district, Central Yunnan. All specimens belong to Zellia, a subgenus of Pseudoschwagerina; three forms may be referred to known subspecies while two are new.     

Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Lateral continuity of basement seismic reflections in 15 Ma ultrafast-spreading crust at ODP Site 1256

The Ocean Drilling Program (ODP) initiated drilling at Site 1256D in the Guatemala Basin, about 1,000 km off the East Pacific Rise to penetrate plutonic rocks, anticipated to be relatively shallow in this region, formed at an ultra-fast spreading rate. IODP Expedition E312 successfully drilled into gabbros at ~1,150 m in basement. Multi-channel seismic traces show weak laterally coherent sub-basement reflections at borehole depths. Synthetic reflectivity seismograms were computed using a Ricker wavelet and impedance profiles from borehole sonic logs. These seismograms show significant sub-basement amplitude peaks. A zero-offset vertical seismic profile, shot on E312, was processed to investigate the authenticity of these reflections and their relationship to borehole geology. A dual scheme of the median filtering and F–K dip filtering was used. Tests with synthetic seismograms indicate the approach is effective at reasonable SNR levels. Downgoing energy is clearly identified but negligible upgoing energy is visible over random noise. These results indicate that lava flows and igneous contacts in upper ocean crust have significant topography on lateral scales less than the Fresnel Zone (~300 m) due to igneous and tectonic processes.     

Sedimentological analysis of two drill cores through the crater moat‐filling breccia,Flynn Creek impact structure,Tennessee

Sedimentological (line‐logging) analysis of two drill cores, FC77‐3 and FC67‐3, situated, respectively, in the northwestern and southeastern quadrants of the Flynn Creek impact structure's crater‐moat area reveals that the ~27 m thick crater moat‐filling breccia consists of three subequal parts. These parts, which were deposited during early modification stage of this marine‐target impact structure, are distinguished on the basis of vertical trends in sorting, grain size, and counts of clasts per meter in comparison with other well‐known marine‐target impact structures, namely Lockne, Tvären, and Chesapeake Bay. The lower part is interpreted to represent mainly slump deposits, and the middle part is interpreted to represent a stage intermediate between slump and marine resurge, that is, a traction flow driven by overriding suspension flow. The upper part (size graded, and relatively well sorted and fine grained) is interpreted to represent marine resurge flow only. The upper part is capped by a relatively thin and relatively fine‐grained calcarenite to calcisiltite deposit.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.     

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.