Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

Nonmodal Growth on a Sphere at Various Horizontal Scales

Nonmodal growth (NG) and unstable normal mode growth are considered in spherical geometry. Two groups of initial conditions (IC) are studied: "connected" IC (common in Cartesian studies) and "separated" IC (based on observed conditions prior to cyclogenesis). Time series of growth rates are emphasized in conjunction with eigenmode projections. Projections show that early on normal mode growth was much stronger for connected IC and that NG caused negative growth early on of some variables for separated IC. Projections explain why amplitude, kinetic energy (KE), and potential vorticity have more NG than available potential energy (APE). Though varying between ICs and with initial phase shift, NG increases with wavenumber. For middle wavelengths, NG is significant and positive using connected IC but negative or small using separated IC. Total energy and KE growth rates of short waves are very similar during the first 2 days for both ICs. Amplitude time series closely follow KE in all cases studied. APE has less overlap than does KE between the main modes present, so less NG occurs for APE than for KE. In separated IC cases, APE growth rates evolve consistent with emergence of an unstable normal mode and little NG.     

Modeling the Permeability Loss of Metallic Iron Water Filtration Systems

Over the past 30 years the literature has burgeoned with in situ approaches for groundwater remediation. Of the methods currently available, the use of metallic iron (Fe⁰) in permeable reactive barrier (PRB) systems is one of the most commonly applied. Despite such interest, an increasing amount of experimental and field observations have reported inconsistent Fe⁰ barrier operation compared to contemporary theory. In the current work, a critical review of the physical chemistry of aqueous Fe⁰ corrosion in porous media is presented. Subsequent implications for the design of Fe⁰ filtration systems are modeled. The results suggest that: (i) for the pH range of natural waters (>4.5), the high volumetric expansion of Fe⁰ during oxidation and precipitation dictates that Fe⁰ should be mixed with a non‐expansive material; (ii) naturally occurring solute precipitates have a negligible impact on permeability loss compared to Fe⁰ expansive corrosion; and (iii) the proliferation of H₂ metabolizing bacteria may contribute to alleviate permeability loss. As a consequence, it is suggested that more emphasis must be placed on future work with regard to considering the Fe⁰ PRB system as a physical (size‐exclusion) water filter device.     

EARLY HOLOCENE CONDITIONS AT THE WESTERN MARGIN OF CAPE COD BAY

A paleoenvironmental model is proposed for the western shoreline of Cape Cod Bay, Massachusetts, from about 9,000 through 5,500 years B.P. Central to this reconstruction is the integration and correlation of the glacial sediments, radiometric chronology and pollen stratigraphy of an exposed peat deposit located near Nameloc Heights, Plymouth, Massachusetts. Four Carbon-14 dates provide a chronology and suggest an average rate of peat accumulation of about 6 cm/100 years from 8,800 through 7,400 years B.P., decreasing to 3 cm/100 years from 7,400 to 5,500 years B.P. This sequence of radiometric dates provides a temporal framework within which the pollen taxa of the peat deposit can be correlated to the established floristic zones of the southern New England region.     

SHORELINE CHANGE AT THOMPSON ISLAND,BOSTON HARBOR,MASSACHUSETTS, 1938-1977

A photogrammetric and sediment analysis is presented to illustrate the relationship between beach erosion and seacliff recession on Thompson Island, Boston Harbor, Massachusetts. Aerial photographs taken in 1938, 1952, 1963, and 1977 were measured to determine rates of shoreline change around the island. The 39-year average rate of beach erosion is 0.3 m/yr ± 2% with an average rate of cliff recession at 0.2 m/yr ± 2%. Rates of beach erosion between six orientations that reflect principal wave approach to the island were not found to be significantly different (0.05) for the 39-year period. The rank order correlation between beach erosion and cliff recession for these six orientations was moderately well correlated (0.63). To determine possible controls of cliff erosion, 58 sediment samples were collected from glacial cliffs along the shoreline. The textural composition was determined, and then tested with discriminant function analysis. Partial correlation analysis between beach erosion and cliff recession holding a surrogate for sediment size constant improved the rank order from 0.63 to 0.84. The results indicate that coarser-grained cliffs recede at faster rates, but with less erosion occurring on adjacent beaches. Alternatively, finer-grained cliffs recede at slower rates, but with greater erosion occurring on adjacent beaches.     

EVOLUTIONARY CONVERGENCE OR PHENOTYPIC PLASTICITY? RESPONSES TO SUMMER RAIN BY CALIFORNIA CHAPARRAL

An unusual summer rain upon southern California chaparral was followed by rapid water uptake, growth, flowering, and fruiting sequences similar to normal spring growth flushes, owing to genetically fixed traits associated with shrub physiology and reproductive biology. This response to summer rain demonstrates the advantage for shrubs to be flexible in the initiation and duration of growth in an unpredictable environment. Evergreen sclerophyllous chaparral is generally interpreted as a product of evolutionary convergence in which environmental stresses posed by Mediterranean climate limit adaptive evolutionary strategies. Phenological records of chaparral taxa or closely related congeners living in summer rain climates in Arizona, central and southern Baja California, and mainland Mexico, along with fossil evidence, suggest that shrubs evolved in response to a somewhat different Tertiary climate; development of the ecosystem reflects the geographic expansion of similar adaptive forms into the newly developing Mediterranean climate, or migrational convergence. Interpretation of evidence in sprouting behavior and reproductive biology suggests that most evolution in California chaparral derives from selective pressures of burning rather than adaptation to drought.     

DRIFT VOLUME IN THE SOUTHERN PENINSULA OF MICHIGAN—A PRODIGIOUS PLEISTOCENE ENDOWMENT

Michigan's Southern Peninsula is by far the largest drift repository in the Great Lakes Area. Thickness of the glacial sediments averages 85 m and volume, previously unknown, exceeds 9300 km³. The distribution, however, is very uneven, as is the underlying bedrock surface. These variations and known surface characteristics provide a basis for the first qualitative and quantitative identification of distinct drift realms. Even the smallest of the eight realms contains about 35% more drift than the better known Kettle Moraine tract of Wisconsin, which encompasses 50% more area. Comparative analysis within the peninsula shows that Pleistocene glaciation completely transformed the topography to the north, whereas the southern half has palimpsest relationships with subcropping formations. [Key words: drift, drift volume, Michigan, Great Lakes, Pleistocene.]     

BEAVER-POND SEQUENCE MORPHOLOGY AND SEDIMENTATION IN NORTHWESTERN MONTANA

Two adjacent beaver-pond sequences of seven ponds each were examined along the eastern boundary of Glacier National Park, Montana, to determine spatial patterns and amounts of sedimentation, both within individual ponds and along pond sequences. Dams and their associated ponds were categorized into two age groups, young versus old, based on dam morphology and surface vegetation. Old ponds contained significantly more sediment than did young ponds. However, hypothesized patterns of downstream reductions in coarse-grained sediment, both within a pond and along a sequence, were not observed, nor were there any discernible patterns in sediment texture across ponds. [Key words: beaver pond, sedimentation, Rocky Mountains, Montana.]