Kriging in a global neighborhood

The kriging estimator is usually computed in a moving neighborhood; only the data near the point to be estimated are used. This moving neighborhood approach creates discontinuities in mapping applications. An alternative approach is presented here, whereby all points are estimated using all the available data. To solve the resulting large linear system the kriging estimator is expressed in terms of the inverse of the covariance matrix. The covariance matrix has the advantage of being positive definite and the size of system which can be solved without encountering numerical instability is substantially increased. Because the kriging matrix does not change, the estimator can be written in terms of scalar products, thus avoiding the more time-consuming matrix multiplications of the standard approach. In the particular case of a covariance which is zero for distances greater than a fixed value (the range), the resulting banded structure of the covariance matrix is shown to lead to substantial computational savings in both run time and storage space. In this case the calculation time for the kriging variance is also substantially reduced. The present method is extended to the nonstationary case.     

The petrogenesis of 30–20 Ma basic and intermediate volcanics from the Mogollon-Datil Volcanic Field,New Mexico,USA

In the western USA calcalkaline magmas were generated hundreds of kilometres from the nearest destructive plate margin, and in some areas during regional extension several Ma after the cessation of subduction. The Mogollon-Datil Volcanic Field (MDVF) in southern New Mexico was a centre of active magmatism in the mid- to late-Tertiary, and a detailed field, petrographic and geochemical study has been undertaken to evaluate the relations between extensional tectonics and calcalkaline magmatism in the period 30–20 Ma. The rocks comprise alkalic to high-K calcalkaline lavas, ranging from basalt to high silica andesitc. Most of the basaltic rocks have relatively low HFSE abundances, elevated ⁸⁷Sr/⁸⁶Sr and low ¹⁴³Nd/¹⁴⁴Nd, similar to many Tertiary basalts across the western USA, and they are inferred to have been derived from the continental mantle lithosphere. Two differentiation trends are recognised, with the older magmas having evolved to more calcalkaline compositions by magma mixing between alkalic basaltic andesites and silicic crustal melts, and the younger rocks having undergone 30–40% fractional crystallisation to more alkalic derivatives. The younger basalts also exhibit a shift to relatively higher HSFE abundances, with lower ⁸⁷Sr/⁸⁶Sr and higher ¹⁴³Nd/¹⁴⁴Nd, and these have been modelled as mixtures between an average post-5 Ma Basin and Range basalt and the older MDVF lithosphere-derived basalts. It is argued that the presence of subduction-related geochemical signatures and the development of calcalkaline andesites in the 30–20 Ma lavas from the MDVF are not related to the magmatic effects of Tertiary subduction. Rather, basic magmas were generated by partial melting of the lithospheric mantle which had been modified during a previous subduction event. Since these basalts were generated at the time of maximum extension in the upper crust it is inferred that magma generation was in response to lithospheric extension. The association of the 30–20 Ma calcalkaline andesites with the apparently anorogenic tectonism of late mid-Tertiary extension, is the result of crustal contamination, in that fractionated, mildly alkaline, basaltic andesite magmas were mixed with silicic crustal melts, generating hybrid andesite lavas with calcalkaline affinities.     

Prediction of hydrocarbons in sedimentary basins

To estimate the undiscovered hydrocarbon potential of sedimentary basins, quantitative play assessments specific for each location in a region may be obtained using geostatistical methods combined with the theory of classification of geological objects, a methodology referred to as regionalization. The technique relies on process modeling and measured borehole data as well as probabilistic methods to exploit the relationship between geology (the predictor) and known hydrocarbon productivity (the target) to define prospective stratigraphic intervals within a basin. It is demonstrated in case studies from the oil-producing region of the western Kansas Pennsylvanian Shelf and the gas-bearing Rotliegend sediments of the Northeast German Basin.     

Acidity of twelve northern New England (U.S.A.) lakes in recent centuries

Studies of sediment cores from 12 acidic lakes in granitic, forested and uninhabited catchments in northern New England, U.S.A. produced diatom-inferred pH (IpH) 5.2 to 5.8 and alkalinity (Ialk) –12 to 31 µeq l^–1, with slowly declining values at some lakes, for one to four centuries prior to logging. Increases of IpH (0.05 to 0.60) and Ialk (5 to 40 µeq l^–1) correlate with logging in the catchments in the early-1800s to early-1900s. Recovery to pre-logging IpH and Ialk correlates with forest succession toward conifers, and is completed in the late-1800s to mid-1900s. Beginning at 1915–1920 (4 lakes), 1930–1950 (4 lakes) and 1965–1970 (4 lakes), IpH and Ialk start decreasing below pre-logging values due to atmospheric acidic deposition, leading to respective total decreases of 0.10 to 0.45 (X=0.25) and 5 to 25 µeq l^–1 (X=15). Inputs of anthropogenic Pb, Zn, V, polycyclic aromatic hydrocarbons, and soot from the atmosphere are first detectable in early- to late-1800s sediment, and Cu in late-1800s sediment, increase rapidly in the late-1800s to mid-1900s, and level off or decrease since the 1960s — partly due to emission controls. Decreasing Ca, Mn, and possibly Zn relative to other metals and normalized to organic content, and increasing flux of Fe to the lakes, indicate soil and water acidification after 1900.This is the fourteenth of a series of papers to be published by this journal which is a contribution of the Paleoecological Investigation of Recent Lake Acidification (PIRLA) project. Drs. D. F. Charles and D. R. Whitehead are guest editors for this series.     

Paleogene island arc collision-related conglomerates, Yarlung–Tsangpo suture zone, Tibet

Coarse clastic rocks in the Liuqu Conglomerate, formed in both terrestrial and subaqueous settings, record a Paleogene phase in the tectonic evolution of Tibet. Facies changes are commonly abrupt with rapid changes in clast types, grain size and stratal patterns. Sediments were derived from the leading (northern) edge of the Indian margin and a Late Jurassic–Cretaceous intraoceanic island arc that lay within Tethys. The coarse clastic sedimentary rocks of the Liuqu conglomerates are extremely proximal, but are locally offset relative to their original source terranes. They record aspects of the history of collision between these terranes and are interpreted to have been deposited in oblique–slip basins that developed along the zone of collision. The absence of clasts derived from terranes to the north of the Yarlung–Tsangpo suture suggests that basins associated with deposition of the Liuqu Conglomerate developed prior to the final collision between India and Asia.     

Salinity tolerance of brown shrimpFarfantepenaeus aztecus as it relates to postlarval and juvenile survival, distribution, and growth in estuaries

The brown shrimp,Farfantepenaeus aztecus, is the major component of the Gulf of Mexico shrimp fishery, and it is critical that we understand its environmental requirements. Brown shrimp spend a large portion of their post-larval (PL) and juvenile life within estuaries distributed along salinity gradients and yet our understanding of the salinity tolerance of various age groups is limited. A series of 48-hr bioassays were conducted in which various ages ofF. aztecus (PL-10, PL-13, PL-15, PL-17, PL-20, and PL-23) were acclimated from a salinity of 26‰ to 1‰, 2‰, 4‰, 8‰, 12‰, and 26‰ in order to determine their tolerance to these salinities. Finally, PL-80.F. aztecus were transferred directly from 25‰ to 2‰, 4‰, and 8‰ waters to study the effects of rapid salinity reductions on juvenile survival. Survival of 10-and 13-day-old PLs was significantly, different from the control (26‰) for all salinities tested. Survival of PL-15 shrimp and older was significantly lower than survival of the controls at 1‰ and 2‰ but similar to the control at all other salinities tested. A 4-wk growth trial was conducted with juvenile shrimp at 2‰, 4‰, 8‰, and 12‰. There was no significant difference in survival among treatments, although shrimp maintained at 8‰ and 12‰ grew significantlymore than shrimp maintained at 2‰ and 4‰. There was no growth difference between shrimp at the two low salinities or between shrimp at the two high salinities. Survival of juveniles transferred directly from 25‰ to various salinities were 100% at 25‰, 94.2% at 8‰, 67.3% at 4‰, and 63.5% at 2‰. These results suggest that PL-13 and younger brown shrimp would have a better chance of survival by delaying entry into estuaries susceptible to rapid salinity declines. The brown shrimp juveniles would, be more densely distributed in areas with salinities greater than 4‰ than in salinities less than 4‰. Although food availability, and bottom type also affect shrimp distribution survival and growth, salinity may also greatly affect the shrimp and its fishery.     

Geology and structural history of the southwest Precordillera margin, northern Mendoza Province, Argentina

Rocks and structures in the southwest Precordillera terrane, located in western Argentina, constrain the Paleozoic distribution of continents and the development of the western margin of Gondwana. Detailed mapping of an area in the southwest Precordillera allowed identification of several pre-Carboniferous rock units formed in distinct tectonic environments and were later tectonically juxtaposed. The pre-Carboniferous rock units comprise carbonate metasiltstone, metasandstone, massive diabase, and quartzo-feldspathic gneiss intruded by ultramafic rocks and layered gabbro. Preliminary structural analysis indicates that the present distribution of units is due to two contractional deformation episodes, an east-directed Devonian ductile event and a west-directed Tertiary brittle event. The metasedimentary rocks, which form the structural base of the area and are part of the western Precordilleran passive margin sequence, were juxtaposed along minor ductile shear zones early in the ductile event. Their contact was then folded during continued ductile deformation; at this time the ultramafic/layered gabbro complex and the massive diabase were emplaced over the metasedimentary units along narrow ductile shear zones. Brittle deformation, associated with the Andean orogeny, involved open folding, thrust faulting, and reactivation of some ductile features.     

Comparison of Electrostatic and Non‐Electrostatic Models for U(VI) Sorption on Aquifer Sediments

A non‐electrostatic generalized composite surface complexation model (SCM) was developed for U(VI) sorption on contaminated F‐Area sediments from the U.S. Department of Energy Savannah River Site, South Carolina. The objective of this study was to test if a simpler, semi‐empirical, non‐electrostatic U(VI) sorption model (NEM) could achieve the same predictive performance as a SCM with electrostatic correction terms in describing U(VI) plume evolution and long‐term mobility. One‐dimensional reactive transport simulations considering key hydrodynamic processes, Al and Fe minerals, as well as H⁺ and U surface complexation, with and without electrostatic correction terms, were conducted. The NEM was first calibrated with laboratory batch H⁺ and U(VI) sorption data on F‐Area sediments, and then the surface area of the NEM was adjusted to match field observations of dissolved U(VI). Modeling results indicate that the calibrated NEM was able to perform as well as the previously developed electrostatic model in predicting the long‐term evolution of H⁺ and U(VI) at the site, given the variability of field‐site data. The electrostatic and NEM models yield somewhat different results for the time period when basin discharge was active; however, it is not clear which modeling approach may be better to model this early time period because groundwater quality data during this period were not available. A key finding of this study is that the applicability of NEM (and thus robustness of its predictions) to the field system evolves with time and is strongly dependent on the pH range that was used to develop the model.     

Estimation of soil shear modulus softening during strong ground shaking using ground surface and downhole acceleration recordings

An algorithm for estimating the average shear modulus for a soil site containing a downhole accelerometer array is described. For distant or weak earthquakes, the usual procedure for estimating shear modulus is to perform time‐series correlation of two downhole records. The vertical distance between instruments divided by the peak correlation lag time gives the average shear wave velocity. The shear modulus follows easily. This method is not applicable for stronger earthquakes where non‐linear softening effects lead to progressively slower shear wave velocities. The method presented here overcomes the softening effect by compressing the time scale of the upper acceleration record. Time compression is accomplished in such a way that the peak correlation of the two records is maximized. The algorithm steps through the records, maximizing the correlation peak by adjusting the time scale within an active time interval. The resulting compressed upper record can be interpreted as the ground motion that would have occurred had softening not taken place. The summation of the various time scale adjustments shows both the amount of softening and the time at which it occurred. Copyright © 2000 John Wiley & Sons, Ltd.