An urban radiation obstruction model

An urban street canyon radiation obstruction model has been developed. The model can describe community structure in terms of the type and dimensions of every building, block, road, park, etc. The need for massive data acquisition in regard to obstruction modeling calls for computerized algorithms, relieving the researcher of the needless tedium of hand calculations and the accompanying high degree of error and labor costs. The model program OBSTRUCT was written in FORTRAN IV for use on the IBM 3033. To facilitate changes or modifications, OBSTRUCT was written in modular form.     

Stepwise use of GFLOW and MODFLOW to determine relative importance of shallow and deep receptors

A stepwise modeling approach is implemented in which a regional one-layer analytic element model is used to simulate the flow system and to furnish boundary conditions for an extracted local three-dimensional model. In this case study the stepwise approach is used to evaluate the fate of recharge in the Menomonee Valley adjacent to Lake Michigan. Two major receptors exist for recharge that flows through contaminated valley fill: the surface water estuary and a tunnel system constructed approximately 75 to 110 m below land surface to store storm runoff. The primary objective of the modeling is to delineate the contributing areas of recharge to each receptor. Of interest is the ability of the one-layer regional model to furnish flux boundary conditions to the local three-dimensional model despite the presence of vertical flow conditions at the boundaries of the local model. Sensitivity analysis suggests that the local model was insensitive to the vertical distribution of the flux. Each step of the modeling approach demonstrates that both receptors play an important role in capturing valley recharge. The pattern of capture of the one-layer model differed in shape from that delineated by the multi-layer local model in the presence of a flow system with pronounced vertical anisotropy and with sinks drawing water from different elevations.     

Simulating ground water-lake interactions: approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAKI problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAKI problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAKI problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.     

Mapping Soils,Vegetation, and Landforms: An Integrative Physical Geography Field Experience*

Students in a graduate seminar at Michigan State University produced a series of detailed vegetation, soils, and landform maps of a 1.5‐square‐mile (3.9 km²) study area in southwest Lower Michigan. The learning outcomes (maps) and skill development objectives (sampling strategies and various GIS applications) of this field‐intensive mapping experience were driven by the assumption that students learn and understand relationships among physical landscape variables better by mapping them than they would in a classroom‐based experience. The group‐based, problem‐solving format was also intended to foster collaboration and camaraderie. The study area lies within a complex, interlobate moraine. Fieldwork involved mapping in groups of two or three, as well as soil and vegetation sampling. Spatial data products assembled and used in the project included topographic maps, a digital elevation model (DEM), aerial photographs, and NRCS (National Resource Conservation Service) soil maps. Most of the soils are dry and sandy, with the main differentiating characteristic being the amount of, and depth to, subsurface clay bands (lamellae) or gravelly zones. The presettlement (early 1830s) vegetation of the area was oak forest, oak savanna, and black oak “barrens.” Upland sites currently support closed forests of white, black, and red oak, with a red maple, dogwood, and sassafras understory. Ecological data suggest that these oak forests will, barring major disturbance, become increasingly dominated by red maple. This group‐based, problem‐solving approach to physical geography education has several advantages over traditional classroom‐based teaching and could also be successfully applied in other, field‐related disciplines.     

Residues of toxaphene decrease in estuarine fish after removal of contaminated sediments

To determine if toxaphene residues in edible fish tissue decreased after removal of contaminated sediments from an estuarine site in 1999, 51 composite samples representing six finfish species were collected in 2001 and analyzed using gas chromatography with electron capture and negative ion mass spectrometric detection. The grand mean total toxaphene residue concentration on a wet weight basis (ΣTOX_wet) was 1,400 ± 3,500 ng g⁻¹ (range: < 18 to 18,000 ng g⁻¹) and was positively correlated with extractable lipid. On a lipid basis, the mean ΣTOX_lip was 26 ± 33 μg g¹, which decreased with increasing distance from the study site. Although benthically-oriented species, such as spot (Leiostomus xanthurus) and striped mullet (Mugil cephalus), exhibited higher mean ΣTOX_wet than those of higher trophic level fish, mean ΣTOX_lip were not significantly different among species. The grand mean ΣTOX for 2001 was 3.8 (wet) and 2.6 (lipid) times less than corresponding preremedial action (1997) concentrations, suggesting that bioavailable toxaphene residues in this system have been reduced. Forage species, such as croaker (Micropogonias undulatus), mullet, and spot, preferentially accumulate toxaphene residues in this system and may serve as vectors of organochlorine contaminants in the estuarine and coastal ocean food web.     

Analytical Solutions for Free-Hydrocarbon Recovery Using Skimmer and Dual-Pump wells

Accidental release of petroleum hydrocarbons to the subsurface may occur through spills around refineries, leaking pipelines, storage tanks, or other sources. If the spill is large, the hydrocarbon liquids may eventually reach a water table and spread laterally in a pancake-like lens. Hydrocarbons that exist as a separate phase are termed light nonaqueous phase liquids (LNAPLs). The portion of the LNAPL that is mobile, not entrapped as residual saturation, is termed "free product."
This paper presents new analytical solutions for the design of long-term free-product recovery from aquifers with skimmer, single- and dual-pump wells. The solutions are for steady-state flow, based on the assumption of vertical equilibrium, and include the effect of coning of LNAPL, air, and water on flow. The solutions are valid for soils of large hydraulic conductivity where the effect of capillary pressure on coning is small.
The results show how to estimate the maximum rate of inflow of LNAPL for skimmer wells, i.e., wells in which LNAPL is recovered with little or no water production. The paper also shows how to calculate the increase in LNAPL recovery when water is pumped by single- or dual-pump wells. A simple equation is given that can be used to adjust the water rate to avoid smearing of the LNAPL below the water table.     

Impurities and nonstoichiometry in the bulk and on the (101?4) surface of dolomite

Atomistic computer simulation methods have been used to model the nature of nonstoichiometry and impurity defects in the bulk and at the (101?4) surface of dolomite (MgCa(CO₃)₂). Calcium and Mg in the bulk and at the surface have been replaced with divalent Ni, Co, Zn, Fe, Mn, and Cd. The results of these calculations indicate that in the bulk, these impurities will prefer to substitute at the Ca site rather than the Mg site. Ca excess in dolomite is most likely incorporated as basal stacking faults; this nonstoichiometry can influence the site distribution of impurities. The calculated surface segregation energies suggest that of all the impurities studied, only Cd will show a strong preference for the (101?4) surface of dolomite.     

Revisiting “An Exercise in Groundwater Model Calibration and Prediction” After 30 Years: Insights and New Directions

In 1988, an important publication moved model calibration and forecasting beyond case studies and theoretical analysis. It reported on a somewhat idyllic graduate student modeling exercise where many of the system properties were known; the primary forecasts of interest were heads in pumping wells after a river was modified. The model was calibrated using manual trial-and-error approaches where a model's forecast quality was not related to how well it was calibrated. Here, we investigate whether tools widely available today obviate the shortcomings identified 30 years ago. A reconstructed version of the 1988 true model was tested using increasing parameter estimation sophistication. The parameter estimation demonstrated the inverse problem was non-unique because only head data were available for calibration. When a flux observation was included, current parameter estimation approaches were able to overcome all calibration and forecast issues noted in 1988. The best forecasts were obtained from a highly parameterized model that used pilot points for hydraulic conductivity and was constrained with soft knowledge. Like the 1988 results, however, the best calibrated model did not produce the best forecasts due to parameter overfitting. Finally, a computationally frugal linear uncertainty analysis demonstrated that the single-zone model was oversimplified, with only half of the forecasts falling within the calculated uncertainty bounds. Uncertainties from the highly parameterized models had all six forecasts within the calculated uncertainty. The current results outperformed those of the 1988 effort, demonstrating the value of quantitative parameter estimation and uncertainty analysis methods.