Ray tomography: Errors and error functions

Tomography is the inversion of boundary projections to reconstruct the internal characteristics of the medium between the source and detector boreholes. Tomography is used to image the structure of geological formations and localized inhomogenieties. This imaging technique may be applied to either seismic or electromagnetic data, typically recorded as transmission measurements between two or more boreholes. Algebraic algorithms are error-driven solutions where the goal is to minimize the error between measured and predicted projections. The purpose of this study is to assess the effect of the ray propagation model, the measurement errors, and the error functions on the resolving ability of algebraic algorithms. The problem under consideration is the identification of a two-dimensional circular anomaly surveyed using crosshole measurements. The results show that: (1) convergence to the position of the circular anomaly in depth between vertical boreholes is significantly better than for convergence in the horizontal direction; (2) error surfaces may not be convex, even in the absence of measurement and model errors; (3) the distribution of information content significantly affects the convexity of averaging error functions; (4) measurement noise and model inaccuracy manifest in increased residuals and in reduced convergence gradients near optimum convergence; (5) the maximum ray error function increases convergence gradients compared with the average error function, and is unaffected by the distribution of information content; however, it has higher probability of local minima. Therefore, inversions based on the minimization of the maximum ray error may be advantageous in crosshole tomography but it requires smooth projections. These results are applicable to both electromagnetic and seismic data for wavelengths significantly smaller than the size of anomalies.     

Spatial interpolation methods for nonstationary plume data

Plume interpolation consists of estimating contaminant concentrations at unsampled locations using the available contaminant data surrounding those locations. The goal of ground water plume interpolation is to maximize the accuracy in estimating the spatial distribution of the contaminant plume given the data limitations associated with sparse monitoring networks with irregular geometries. Beyond data limitations, contaminant plume interpolation is a difficult task because contaminant concentration fields are highly heterogeneous, anisotropic, and nonstationary phenomena. This study provides a comprehensive performance analysis of six interpolation methods for scatter-point concentration data, ranging in complexity from intrinsic kriging based on intrinsic random function theory to a traditional implementation of inverse-distance weighting. High resolution simulation data of perchloroethylene (PCE) contamination in a highly heterogeneous alluvial aquifer were used to generate three test cases, which vary in the size and complexity of their contaminant plumes as well as the number of data available to support interpolation. Overall, the variability of PCE samples and preferential sampling controlled how well each of the interpolation schemes performed. Quantile kriging was the most robust of the interpolation methods, showing the least bias from both of these factors. This study provides guidance to practitioners balancing opposing theoretical perspectives, ease-of-implementation, and effectiveness when choosing a plume interpolation method.     

On coastal flow in the northeast gulf of Alaska near yakutat 1

Abstract

Seasonal maps of geopotential topography and salinity at 10 m were prepared for a region in the northeast Gulf of Alaska near Yakutat. Coastal waters in this region receive considerable runoff in late summer, but the dilution does not result in a greatly intensified coastal flow. Baroclinic coastal flow in the northeast Gulf appears to be weaker than the coastal flow (Kenai Current) to the west along the Kenai Peninsula and in Shelikof Strait.     

The distribution of rare earth elements between monzogranitic melt and the aqueous volatile phase in experimental investigations at 800 °C and 200 MPa

The partitioning of the rare earth elements between a peraluminous monzogranitic melt and a chloride-bearing, sulfur- and carbon dioxide-free, aqueous volatile phase was examined experimentally as a function of chloride and major element concentrations at 800 °C and 200 MPa. The light rare earth elements (e.g. La, Ce) partition into the aqueous volatile phase to a greater extent than the heavy rare earth elements (e.g. Yb, Lu). Distribution of the rare earth elements and the major elements H, Na, K, Ca, and Al between the melt phase (mp) and aqueous volatile phase (aq) is a function of the chlorine concentration in the system, and our data are consistent with the rare earth and major elements occurring as chloride complexes in the aqueous volatile phase. Apparent equilibrium constants for experiments at 800 °C and 200 MPa, K ^′ _REE,Na ^aq/mp , expressed as the ratio of the concentration of a given rare earth element in the aqueous volatile phase to the concentration of the same element in the melt phase, divided by the cubed ratio of sodium in the aqueous volatile phase to the concentration of sodium in the melt phase, decrease systematically with increasing atomic number from K ^′ _La,Na ^aq/mp = 0.41(±0.03) to K ^′ _Lu,Na ^aq/mp =0.11(±0.01), except for Eu. These experimentally derived apparent equilibrium constants for the rare earth elements can be used in a numerical simulation of magmatic volatile exsolution. The simulation gave results consistent with the elemental distribution in the potassic alteration zone of a deep porphyry copper deposit, but higher concentrations of heavy rare earth elements are released into the magmatic aqueous solution than are captured in the secondary mineralization. Received: 1 November 1999 / Accepted: 7 June 2000     

Fluorine in lunar samples: implications concerning lunar fluorapatite

The F contents of a number of Apollo 14 and 15 samples range from less than a ppm for anorthosite rock fragments to ~165 ppm for some soils and breccias. Apollo 15 soils tend to have lower F contents (50–70 ppm) than soils from other sites. In most cases samples were run simultaneously with W-1 in which F was determined to be 216 (±11) ppm.The $\text{F}\text{P}{}_2\text{O}{}_5$ ratio is 0·032 ± 0·005 in soils and rocks. A correlation exists in soils between F, P₂O₅, and that fraction of the Cl which is insoluble in hot water. The $\text{F}\text{Cl}{}_{\mathrm{r}}$ ratio in soils and rocks, though different, requires that the phosphate phase involved be fluorapatite; this is consistent with mineralogical observations. F, like Cl, is correlated with KREEP elements at all sites for which data are available.     

Burial and exhumation history of Pennsylvanian strata, central Appalachian basin: an integrated study

An inferred burial and exhumation history of Pennsylvanian strata in the central Appalachian foreland basin is constrained by integrating palaeothermometers, geochronometers and estimated palaeogeothermal gradients. Vitrinite reflectance data and fluid inclusion homogenization temperatures indicate that burial of Lower and Upper Pennsylvanian strata of the Appalachian Plateau in West Virginia exceeded ～4.4 km during the late Permian and occurred at a rate of ～100 m Myr^?1. Exhumation rates of ～10 m Myr^?1 from the late Permian to the early Cretaceous are constrained using maximum burial conditions and published apatite fission track (AFT) ages. AFT and radiogenic helium ages indicate exhumation rates of ～30–50 m Myr^?1 from the early to late Cretaceous. Radiogenic helium dates and present day sampling depths indicate that exhumation rates from the late Cretaceous to present were ～25 m Myr^?1. Exhumation rates for Upper and Lower Pennsylvanian strata within the Appalachian Plateau are remarkably similar. Early slow exhumation was possibly driven primarily by isostatic rebound associated with Triassic rifting. The later, more rapid exhumation can be attributed to thermal expansion followed by lithospheric flexure related to sediment loading along the passive margin.     

Pito Rift: How a large-offset rift propagates

The Pito Rift area is the site of actively deforming oceanic lithosphere that has been primarily under extension for at least the past million years, based on kinematic reconstructions. The major morphologic features, Pito Deep and Pito Seamount, are aligned toward the Euler pole for relative motion between the Easter and Nazca plates. SeaMARC II side-scan and bathymetry data indicate that there are two general modes of faulting currently active in the Pito Rift area. One is associated with incipient rifting of old (3 Ma) Nazca lithosphere by large NW-SE normal faults, and the other is associated with a broad area of right-lateral transform shear between the Nazca and Easter plates. This transform shear is distributed over a broad region because of the northward growth of the East Rift and parallel tectonic rifting within the Pito Rift area. The majority of the Pito Rift area is composed of preexisting blocks of Nazca plate that are back-tilted away from Pito Deep and strike perpendicular to present and previous relative plate motions. This observation suggests that block-faulting and back-tilting are the primary mechanisms responsible for the distributed lithospheric extension, in agreement with gravity and magnetic analyses (Martinez et al., this issue).The only recent volcanic flows observed in side-scan data are from the Pito Seamount area and to the outside of the outer pseudofault of the East Rift. The significance of the young flows near the outer pseudofault is not understood. We interpret the flows extending northwest from the Pito Seamount as representing a newly formed seafloor spreading axis within the Pito Rift area. Gravity and magnetic analyses (Martinez et al., this issue) together with SeaMARC II bathymetry and side-scan data support this interpretation.Based on the tectonic evolution of the Easter microplate, we propose an evolutionary model for the formation of the Pito Rift area, where new tectonic grabens form immediately west of the previous graben and with slightly more counterclockwise orientation. The duration and history of tectonic activity for each graben are not well constrained.     

Oil spill fishery impact assessment modeling: The fisheries recruitment problem

A model to assess the impact of oil spill on fisheries, consisting of an oil spill fates model, a continental shelf hydrodynamics model, an ichthyoplankton transport and fates model, and a fish population model, has been applied to the Georges Bank-Gulf of Maine region to estimate the impact of oil spills on several important commercial fisheries. The model addresses direct impacts of oil on a fishery through hydrocarbon-induced egg and larval mortality. This early life stage mortality is estimated by dynamically mapping the spatial intersection of the surface and subsurface oil concentrations resulting from the spill with the developing eggs and larvae. Ichthyoplankton entering an area with hydrocarbon concentrations in excess of a specified threshold are assumed lost. Model output is given in terms of differential catch, comparing the non-impacted and the hydrocarbon impacted fisheries. Difficulties in establishing stock-recruit relationships, and the inability to predict first year survival even one year ahead make the quantification of absolute catch losses impossible. Output of the model system discussed here is therefore limited to relative rather than absolute catch losses.The paper is organized to demonstrate first the importance of the recruitment question to impact estimation, second that a modeling methodology is necessary to evaluate impacts given the magnitude of unexplained observed recruitment variability, and third a stochastic solution to the problem which places impact estimates in the context of a probability distribution. Lastly, the model system is applied to the problem of attaining better early life history mortality estimates, to ultimately improve impact estimation capabilities.     

River flooding, storm resuspension, and event stratigraphy on the northern California shelf: observations compared with simulations

Coast-hugging surface flood plumes occur on the inner shelf of northern California during the winter season, generating dense, near-bottom suspensions which may attain fluid mud concentrations as particles settle. The period of storm-heightened waves may continue into the flood period, leading to gravity-driven seaward displacement of the bottom suspension; or the wave regime may ameliorate, leaving the suspension to consolidate as a short-lived, inner-shelf flood bed. Such beds tend to be resuspended within days or weeks by subsequent storm events that may recreate the original high concentrations. The sediment is thus dispersed seaward by gravity flows, to be deposited as a muddy flood bed on the central shelf. The locus of deposition of these “high-concentration regimes” is a function of the relative intensities of river discharge and storm wave height. Greater discharge piles thicker storm beds nearer shore, while intense wave regimes allow deposition of the fluid mud further seaward. During events with high values of both parameters, large amount of fluid mud may bypass over the shelf edge. In contrast, “low-concentration regimes” occur during storm periods when there has been no recent flood deposition on the inner shelf. The shelf floor is better consolidated than in the previous case, and the resulting suspended sediment concentrations are lower. As a consequence, low-concentration regimes are winnowing and bypassing regimes, and the beds deposited are thinner and sandier. Algorithms describing deposition by high and low-concentration regimes have been embedded in a probabilistic model. A simulation of a 400-year sequence of beds deposited by winter storms and floods suggests that on the Eel shelf, the Holocene transgressive systems tract consists of back-stepping, seaward-fining event beds, whose timelines (bedding planes) dip more gently than do their gradational facies boundaries. At these longer time scales, flood beds dominate over storm beds.     

Using 2MASS to detect main sequence companions to subdwarf B stars

We explore the binary fraction of subdwarf B (sdB) stars by using the Two Micron All Sky Survey to search for main sequence companions. We have convolved Kurucz models with appropriate filter bandpasses to examine how various combinations separate single sdB stars from sdB+MS binaries. We notice that the 2MASS magnitude limits greatly increase the fraction of inferred sdB+MS binaries and examine a magnitude-limited sample appropriate for single sdB stars.