Monitored Natural Attenuation of Contaminants in the Subsurface: Applications

In recent years there has been increasing interest in the application of passive technologies to reduce or remove contaminants from the subsurface environment including soil and ground water. In most cases, the impetus for this interest lies in a perceived savings compared with more traditional remedial alternatives. In a few cases, the infrastructure at contaminated sites, such as buildings, paved areas, and utilities, makes the use of conventional remedial measures difficult and expensive.
To demonstrate that natural processes are effective in reaching established goals, it is necessary to determine that transformation processes are taking place at a rate that is protective of human health and the environment and that these processes will continue for an acceptable period of time. The basic conditions that must be present to confirm natural attenuation processes arc taking place are discussed along with the behavior of contaminated plumes, monitoring requirements, data analysis, rates of degradation, and mathematical modeling.     

Phase relations between Ca3Fe 2 3 +Si3O12-Fe 3 2 +Fe 2 3 +Si3O12 garnet and CaFeSi2O6-Fe2Si2O6 pyroxene solid solutions

Editorial responsibility: J. Hoefs     

The transition from subduction to continental collision: crustal structure in the North Canterbury region, New Zealand

The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.     

Isotopic evidence for meteoric-hydrothermal alteration of plutonic igneous rocks in the Yakutat Bay and Skagway areas,Alaska

Oxygen and hydrogen isotope ratios were measured on coexisting minerals from quartz diorites and quartz monzonites from a section across the Coast Range batholith in the Skagway area, Alaska, including a variety of outlying plutons west of the batholith in the Yakutat Bay-Mt. St. Elias region (latitudes 59–60°N). The extremely low and variable δ¹⁸Oand δD results indicate widespread meteoric-hydrothermal alternaiton of the Coast Range batholith, and to a lesser extent, of the Yakutat Bay plutons as well. In the Yakutat Bay area, the plutons with K—Ar ages younger than 50 m.y. have widely varying δD values of ?72 to ?148, compared to δD = ?69to?90 for all but one sample in the 50–225 m.y. age grouping (one biotite has δD = ?109). This suggests that the major meteoric-hydrothermal episodes in this area occurred during the Eocene and Miocene. This involved relatively small meteoric water/rock ratios(<0.1), as none of the δ¹⁸O values show any clear-cut evidence of alteration (δ¹⁸O_quartz= 7.4 ?11.8; δ¹⁸O_feldspar= 5.7?10.0). However, in the section across the Coast Range batholith, 85% of the plutonic rocks have very low δD values of ?100 to ?167, and the δ¹⁸O values are extremely variable δ¹⁸O_feldspar= + 10.3to?4.0 and Δ¹⁸O_{quartz-feldspar}= 0.4?10.5. These data indicate that a major portion of the batholith, particularly the quartz monzonite-rich eastern part, but also including many of the quartz diorite plutons as well, interacted with meteoric-hydrothermal convective systems that involved water/rock ratios of about 0.3–1.4. The quart diorite plutons are most depleted in¹⁸O near their northeast contacts against the younger quartz monzonite intrusions. The primary igneous δ¹⁸O values of the quartz diorites were apparently higher than those of the quartz monzonites; they are also unusually high in¹⁸O compared to most other analyzed quartz diorites, suggesting derivation from, exchange with, or assimilation of high-¹⁸O metasediments or altered volcanic rocks. These data and conclusions are very similar to those reached previously on a similar isotopic study of the Coast Range batholith in British Columbia, 700 km to the southeast at latitudes 54–55°N, except that in the Skagway area an even greater proportion of the batholith was apparently depleted in deuterium. This implies that deep (?5km?) circulation of meteoric groundwaters is probably a characteristic of the later stages of emplacement of the Cordilleran batholiths of western North America, suggesting that the eastern sections of these batholiths in particular were emplaced at relatively shallow depths.     

A neodymium,strontium, and oxygen isotopic study of the Cretaceous Samail ophiolite and implications for the petrogenesis and seawater-hydrothermal alteration of oceanic crust

In the Samail ophiolite,¹⁴⁷Sm-¹⁴³Nd,⁸⁷Rb-⁸⁷Sr, and¹⁸O/¹⁶O isotopic systems have been used to distinguish between sea-floor hydrothermal alteration and primary magmatic isotopic variations. The Rb-Sr and¹⁸O/¹⁶O isotopic systems clearly exhibit sensitivity to hydrothermal interactions with seawater while the Sm-Nd system appears essentially undisturbed. Internal isochrons have been determined by the¹⁴⁷Sm-¹⁴³Nd method using coexisting plagioclase and pyroxene and give crystallization ages of 130 ± 12m.y. from Ibra and 100 ± 20 m.y. from Wadi Fizh. These ages are interpreted as the time of formation of the Samail oceanic crust and are older than the inferred emplacement age of 65–85 m.y. The initial¹⁴³Nd/¹⁴⁴Nd ratios for a tectonized harzburgite, cumulate gabbros, plagiogranite, sheeted dikes and a basalt have a limited range in ε_Nd of from 7.5 to 8.6 for all lithologies, demonstrating a clear oceanic affinity and supporting earlier interpretations based on geologic observations and geochemistry. The⁸⁷Sr/⁸⁶Sr initial ratios on the same rocks have an extremely large range of from 0.70296 to 0.70650 (ε_Sr = ?19.7 to +30.5) and the δ¹⁸O values vary from 2.6 to 12.7. These large variations are clearly consistent with hydrothermal interaction of seawater with the oceanic crust. A model is presented for the closed system exchange of Sr and O, that in principle illustrates how the Sr isotopic data may be utilized to estimate the water/rock ratio and subsequently used to evaluate the temperature of equilibration between the water and silicates from the¹⁸O/¹⁶O water-rock fractionation.     

Fallout radionuclides in the Pacific Ocean: Vertical and horizontal distributions,largely from GEOSECS stations

From GEOSECS stations, largely, the 1974 distributions of Pu and of¹³⁷Cs are described in the Pacific Ocean north of about 20°S latitude. Changes in some of these distributions are described from 1978 cruises by the authors.The Pacific exhibited, everywhere, a shallow subsurface layer of Pu-rich water with its concentration maximum at about 465 m in 1974; over a large portion of the central North Pacific a second layer of Pu-labelled water, less concentrated than the shallow layer, lay just above the bottom. Similar features were not observed in the case of¹³⁷Cs.The inventories of both Pu and¹³⁷Cs in the water column at most 1974 stations are substantially greater than those to be expected from world-wide fallout alone; these inventory excesses appear to be attributable to close-in fallout, but only if the ratio Pu/¹³⁷Cs in this source was much higher than in world-wide fallout. The North Pacific mean ratio of the inventories is 2.2 times that observed in world-wide fallout.Resolubilization of Pu both from sinking particles and from sediments explains peculiarities of its depth distributions.There is little evidence for tracer movement by sliding downward along density surfaces;¹³⁷Cs appears to have moved to depth by downmixing at the edge of the Kuroshio, and then moved horizontally and upward alongσ_t contours. The shallow Pu-rich layer shows no coordination with density, salinity or O₂ isopleths. The deep Pu-rich layer is restricted to a narrow range of O₂ concentrations that confirm its origin in the Aleutian Trench and rapid spread southward and laterally. Near-bottom circulation processes have been much more active than here-to-fore described.     

A mathematical model to predict the inelastic response of a steel frame: Formulation of the model

The purpose of this research is to use data from experiments to formulate a mathematical model that will predict the non-linear response of a single-storey steel frame to an earthquake input. The process used in this formulation is system identification. The form of the model is a second-order non-linear differential equation with linear viscous damping and Ramberg—Osgood type hysteresis. The damping coefficient and the three parameters in the hysteretic model are to be established. An integral weighted mean squared error function is used to evaluate the [goodness of fit] between the model's response and the structure's response when both are subjected to the same excitation. The function includes errors in displacement and acceleration and is integrated from zero to a time T, which may be the full duration of the recorded response or only a portion of it. The parameters are adjusted using a modified Gauss-Newton method until the error function is minimized. The computer program incorporating these steps in the system identification process is verified with simulated data. Results given in the paper show that in every case the program converges in few iterations to the assigned set of parameters.     

A mathematical model to predict the inelastic response of a steel frame: Establishment of parameters from shaking table experiments

The purpose of this research is to use data from experiments to formulate a mathematical model that will predict the non-linear response of a single-storey steel frame to an earthquake input. The process used in this formulation is system identification. In experiments performed on a shaking table, the frame was subjected to two earthquake motions at several intensities. In each case the frame underwent severe inelastic deformation. A computer program which incorporates the concepts of system identification makes use of the recorded data to establish four parameters in a non-linear mathematical model. When different amounts of data are used in the program, parameter sets are established which give the best model response for that amount of test data. The resulting sets of parameters reflect the way in which the properties of the structure change during the excitation. However, when the full durations of the different excitations are used, the sets of parameters are almost identical. For each of these sets of parameters, the correlation of the computed accelerations with the measured is excellent, and the shape of the computed displacement response compares very well with the measured response, although the permanent offset of the displacements is not computed exactly. Suggestions are given on how to overcome this deficiency in the mathematical model.