A note on frame transformations with applications to geodetic datums

Rigorous equations in compact symbolic matrix notation are introduced to transform coordinates and velocities between ITRF frames and modern GPS-based geocentric geodetic datums. The theory is general but, after neglecting higher than second-order terms, it is shown that the equations revert to the formulation currently applied in most major continental datums. We discuss several examples: the North American Datum of 1983 (NAD83), the European Terrestrial Reference System of 1989 (ETRS89), the Geodetic Datum of Australia of 1994 (GDA94), and the South American Geocentric Reference System (SIRGAS). Electronic Publication     

The permeability of the Elkhorn fault zone,South Park,Colorado

The purposes of this study are to use both field and modeling approaches to characterize the permeability of a fault and to assess the role of the fault on regional ground water flow. The study subject is the Elkhorn fault, a low-angle reverse fault that brings Precambrian crystalline rocks over the sediments of Colorado's South Park Basin. The fault is hypothesized to act as a low-permeability barrier to flow, restricting interaction between the crystalline aquifer and the basin sediments. To test this hypothesis and to better predict the permeability structure of the fault, we synthesized geologic data to create a geologic model of the fault, conducted aquifer tests to estimate the hydrogeologic properties of the fault zone, and used ground water modeling to test the influence of a range of hydraulic properties for the fault zone on ground water flow in the region. Our study suggests that the fault is a low-permeability feature. Estimated heads are best matched to observations by modeling the fault as a 10-foot-thick interval of low-permeability fault gouge. Steady-state flow models show that much of the flow in the study area is topographically driven near land surface. Flow rates decrease with depth in the aquifers. In the footwall, ground water moves updip in the Michigan-San Isabel syncline to discharge in the South Park Basin. In the hanging wall, ground water moves east to a regional ground water divide. Sensitivity analyses indicate that hydraulic heads are most sensitive to changes in hydraulic conductivity and recharge.     

From models to performance assessment: the conceptualization problem

Today, models are ubiquitous tools for ground water analyses. The intent of this paper is to explore philosophically the role of the conceptual model in analysis. Selection of the appropriate conceptual model is an a priori decision by the analyst. Calibration is an integral part of the modeling process. Unfortunately a wrong or incomplete conceptual model can often be adequately calibrated; good calibration of a model does not ensure a correct conceptual model. Petroleum engineers have another term for calibration; they refer to it as history matching. A caveat to the idea of history matching is that we can make a prediction with some confidence equal to the period of the history match. In other words, if we have matched a 10-year history, we can predict for 10 years with reasonable confidence; beyond 10 years the confidence in the prediction diminishes rapidly. The same rule of thumb applies to ground water model analyses. Nuclear waste disposal poses a difficult problem because the time horizon, 1000 years or longer, is well beyond the possibility of the history match (or period of calibration) in the traditional analysis. Nonetheless, numerical models appear to be the tool of choice for analyzing the safety of waste facilities. Models have a well-recognized inherent uncertainty. Performance assessment, the technique for assessing the safety of nuclear waste facilities, involves an ensemble of cascading models. Performance assessment with its ensemble of models multiplies the inherent uncertainty of the single model. The closer we can approach the idea of a long history with which to match the models, even models of nuclear waste facilities, the more confidence we will have in the analysis (and the models, including performance assessment). This thesis argues for prolonged periods of observation (perhaps as long as 300 to 1000 years) before a nuclear waste facility is finally closed.     

Dip and anisotropy effects on flow using a vertically skewed model grid

Darcy flow equations relating vertical and bedding-parallel flow to vertical and bedding-parallel gradient components are derived for a skewed Cartesian grid in a vertical plane, correcting for structural dip given the principal hydraulic conductivities in bedding-parallel and bedding-orthogonal directions. Incorrect-minus-correct flow error results are presented for ranges of structural dip (0 < or = theta < or = 90) and gradient directions (0 < or = phi < or = 360). The equations can be coded into ground water models (e.g., MODFLOW) that can use a skewed Cartesian coordinate system to simulate flow in structural terrain with deformed bedding planes. Models modified with these equations will require input arrays of strike and dip, and a solver that can handle off-diagonal hydraulic conductivity terms.     

Weathering of sulfidic shale and copper mine waste: secondary minerals and metal cycling in Great Smoky Mountains National Park,Tennessee, and North Carolina,USA

Metal cycling via physical and chemical weathering of discrete sources (copper mines) and regional (non-point) sources (sulfide-rich shale) is evaluated by examining the mineralogy and chemistry of weathering products in Great Smoky Mountains National Park, Tennessee, and North Carolina, USA. The elements in copper mine waste, secondary minerals, stream sediments, and waters that are most likely to have negative impacts on aquatic ecosystems are aluminum, copper, zinc, and arsenic because these elements locally exceed toxicity guidelines for surface waters or for stream sediments. Acid-mine drainage has not developed in streams draining inactive copper mines. Acid-rock drainage and chemical weathering processes that accompany debris flows or human disturbances of sulfidic rocks are comparable to processes that develop acid-mine drainage elsewhere. Despite the high rainfall in the mountain range, sheltered areas and intermittent dry spells provide local venues for development of secondary weathering products that can impact aquatic ecosystems.Electronic Supplementary Material  Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.

Episodic colonization of an intertidal mudflat by native cordgrass (<Emphasis Type="Italic">Spartina foliosa</Emphasis>) at Tijuana Estuary

Following heavy winter storms and sedimentation in 1993,Spartina foliosa (Pacific cordgrass) clones established on a 6.5-ha mudflat in Tijuana Estuary, with over 80 new clones counted by 1997. El Niño Southern Oscillation (ENSO) storms in 1993 apparently facilitated the habitat conversion through river flooding, which caused a temporary reduction in soil salinity and delivered large volumes of sediment. Extreme sedimentation likely raised mudflat elevations enough to allowSpartina establishment. We hypothesized that clones, once established, increased sedimentation in a positive feedback loop leading to accelerated habitat conversion. We collected data on elevation,Spartina expansion, and sediment accretion in two consecutive years (1998–1999). The elevation range of the mudflats in 1998 (0.5–0.83 m NGVD) was within the elevation range ofSpartina at this site (0.39–0.83 m NGVD), indicating that remaining mudflats are at elevations suitable for further expansion. Sediment accretion ranged from 4.0–12.7 cm between 1997–1998 (ENSO conditions), but was close to long-term averages (? cm yr^?1) in 1998–1999 (nonflood year) indicating how susceptible Tijuana Estuary is to sedimentation from episodic storms. Although accretion rates were similar withinSpartina clones and on bare mudflats over the ENSO winter, clones were typically dome-shaped suggesting higher sediment retention rates within clones. The radial expansion rates of clones (1.31±0.25 m in 1998; 1.12±0.07 m in 1999) approximated the maxima reported for this species and were not related to clone size or vigor. Conditions on the mudflat appear ideal forSpartina growth, masking differences that might otherwise be observed. Given the likelihood of sedimentation-driven habitat conversion in southern California and other Mediterranean-type estuaries, management efforts are needed to address sedimentation issues on a watershed scale. We recommend that wetland restoration projects in southern California include large areas of intertidal mudflat, both to maintain habitat for shorebird feeding and to allow colonization by salt marsh vegetation.     

Instructions for Authors

Instructions for Authors

Instructions for Authors     

The Northern North Atlantic: A Changing Environment, edited by P. Schäfer, W. Ritzrau, M. Schläter and J. Thiede (Book Review)

Volume Contents

Volume contents     

Stable Isotope Ratios: Hurricane Olivia

The oxygen and hydrogen isotopic compositions of rains from HurricaneOlivia (1994) in the eastern Pacific were measured. The rains werecollected on 24 and 25 September during airplane flights conducted at anelevation of 3 km. Hurricane Olivia peaked in intensity to a category-4storm between the two dates. Isotope ratios of rains from HurricaneOlivia were markedly lower ( ¹⁸O = –13.9to –28.8) than that of rain collected from a thunderstormat an elevation of 2.3 km outside the influence of Olivia (¹⁸O = –3.8). A distinct decrease in isotoperatios from the first day to the next ( ¹⁸O =–18.4 to –21.9) in Hurricane Olivia wasattributed to decreased updraft velocities and outflow aloft. Thisshifted the isotopic water mass balance so that fewer hydrometeors werelifted and more ice descended to flight level. A decrease in the averagedeuterium excess from the first day to the next (d = 15.5 to 7.1)was attributed to an increase in the relative humidity of the watervapor `source' area. We hypothesize that the `source' region for therain was in the boundary layer near the storm center and that becausethe hurricane was at peak intensity prior to the second day the relative humidity was higher.     

Oxygen isotope study of the Long Valley magma system, California: isotope thermometry and convection in large silicic magma bodies

Products of voluminous pyroclastic eruptions with eruptive draw-down of several kilometers provide a snap-shot view of batholith-scale magma chambers, and quench pre-eruptive isotopic fractionations (i.e., temperatures) between minerals. We report analyses of oxygen isotope ratio in individual quartz phenocrysts and concentrates of magnetite, pyroxene, and zircon from individual pumice clasts of ignimbrite and fall units of caldera-forming 0.76 Ma Bishop Tuff (BT), pre-caldera Glass Mountain (2.1-0.78 Ma), and post-caldera rhyolites (0.65-0.04 Ma) to characterize the long-lived, batholith-scale magma chamber beneath Long Valley Caldera in California. Values of '¹⁸O show a subtle 1‰ decrease from the oldest Glass Mountain lavas to the youngest post-caldera rhyolites. Older Glass Mountain lavas exhibit larger (~1‰) variability of '¹⁸O(quartz). The youngest domes of Glass Mountain are similar to BT in '¹⁸O(quartz) values and reflect convective homogenization during formation of BT magma chamber surrounded by extremely heterogeneous country rocks (ranging from 2 to +29‰). Oxygen isotope thermometry of BT confirms a temperature gradient between "Late" (815 °C) and "Early" (715 °C) BT. The '¹⁸O(quartz) values of "Early" and "Late" BT are +8.33 and 8.21‰, consistent with a constant '¹⁸O(melt)=7.8ǂ.1‰ and 100 °C temperature difference. Zircon-melt saturation equilibria gives a similar temperature range. Values of '¹⁸O(quartz) for different stratigraphic units of BT, and in pumice clasts ranging in pre-eruptive depths from 6 to 11 km (based on melt inclusions), and document vertical and lateral homogeneity of '¹⁸O(melt). Worldwide, five other large-volume rhyolites, Lava Creek, Lower Bandelier, Fish Canyon, Cerro Galan, and Toba, exhibit equal '¹⁸O(melt) values of earlier and later erupted portions in each of the these climactic caldera-forming eruptions. We interpret the large-scale '¹⁸O homogeneity of BT and other large magma chambers as evidence of their longevity (>10⁵ years) and convection. However, remaining isotopic zoning in some quartz phenocrysts, trace element gradients in feldspars, and quartz and zircon crystal size distributions are more consistent with far shorter timescales (10²-10⁴ years). We propose a sidewall-crystallization model that promotes convective homogenization, roofward accumulation of more evolved and stagnant, volatile-rich liquid, and develops compositional and temperature gradients in pre-climactic magma chamber. Crystal + melt + gas bubbles mush near chamber walls of variable '¹⁸O gets periodically remobilized in response to chamber refill by new hotter magmas. One such episode of chamber refill by high-Ti, Sr, Ba, Zr, and volatile-richer magma happened 10³-10⁴ years prior to the 0.76-Ma caldera collapse that caused magma mixing at the base, mush thawing near the roof and walls, and downward settling of phenocrysts into this hybrid melt.