Finite element formulation and algorithms for unsaturated soils. Part I: Theory

This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd.     

Boltysh,another end‐Cretaceous impact

Abstract— The Chixculub impact occurred at the Cretaceous/Tertiary (K/T) boundary, and although several other Late Cretaceous and Paleogene impact craters have, at times, been linked with the K/T boundary, isotope geochronology has demonstrated that all have significantly different ages. The currently accepted age of the 24 km diameter Boltysh crater, a K‐Ar whole‐rock age, places it in the Coniacian at 88 ± 3 Ma. However, comprehensive Ar‐Ar dating of a range of melt samples yields a mean age of 65.17 ± 0.64 Ma, within errors of the K/T boundary. Several of the fresh samples exhibit signs of excess argon but this seems to be concentrated in rapidly crystallized glass‐rich samples. The Ar‐Ar age confirms an earlier fission track measurement and thus two dating techniques have yielded an age within errors of the K/T boundary for this crater. Crucially, although the ages of Boltysh and Chixculub are within errors, they may not have formed synchronously. Craters of 24 km diameter occur much more commonly than impacts of Chixculub dimensions, but their proximity does raise the important question of how many impacts there might have been close to the K/T boundary.     

FOUR NEW IRON METEORITE FINDS

Four new irons are described; Buenaventura (IIIB) from Chihuahua, Mexico: mass 114 kg; Denver City (anomalous) from Texas, USA: mass 26.1 kg; Kinsella (IIIB) from Alberta, Canada: mass 3.7 kg; and Tacoma (IA) from Washington, USA: mass 17 g. Denver City is unique, i.e., not related to any other known iron. Tacoma is the smallest iron meteorite recorded. All were purchased for the UCLA collection following a publicity drive for new meteorites     

Transboundary impacts on regional ground water modeling in Texas

Recent legislation required regional grassroots water resources planning across the entire state of Texas. The Texas Water Development Board (TWDB), the state's primary water resource planning agency, divided the state into 16 planning regions. Each planning group developed plans to manage both ground water and surface water sources and to meet future demands of various combinations of domestic, agricultural, municipal, and industrial water consumers. This presentation describes the challenges in developing a ground water model for the Llano Estacado Regional Water Planning Group (LERWPG), whose region includes 21 counties in the Southern High Plains of Texas. While surface water is supplied to several cities in this region, the vast majority of the regional water use comes from the High Plains aquifer system, often locally referred to as the Ogallala Aquifer. Over 95% of the ground water demand is for irrigated agriculture. The LERWPG had to predict the impact of future TWDB-projected water demands, as provided by the TWDB, on the aquifer for the period 2000 to 2050. If detrimental impacts were noted, alternative management strategies must be proposed. While much effort was spent on evaluating the current status of the ground water reserves, an appropriate numerical model of the aquifer system was necessary to demonstrate future impacts of the predicted withdrawals as well as the effects of the alternative strategies. The modeling effort was completed in the summer of 2000. This presentation concentrates on the political, scientific, and nontechnical issues in this planning process that complicated the modeling effort. Uncertainties in data, most significantly in distribution and intensity of recharge and withdrawals, significantly impacted the calibration and predictive modeling efforts. Four predictive scenarios, including baseline projections, recurrence of the drought of record, precipitation enhancement, and reduced irrigation demand, were simulated to identify counties at risk of low final ground water storage volume or low levels of satisfied demand by 2050.     

Constraints on spatial variability in soft-sediment communities affected by contamination from an Antarctic waste disposal site

A small-scale (<500 m length) transect-based survey was conducted in December 1998 to examine the spatial distribution of soft-sediment communities and of concentrations of heavy metals and hydrocarbons in sediments in Brown Bay, adjacent to an abandoned waste dump, at Casey Station, Antarctica. Samples were taken along three transects at increasing distances (nine stations) from the shore-line waste dump. A gradient of contamination was detected, but concentrations of contaminants were very variable with "hotspots" or high levels of contaminants at some stations. Multivariate analysis revealed that the distribution of soft-sediment communities was distinctly different between the inner, middle and outer stations. Abundances of most taxa were very variable with few patterns apparent, but some fauna displayed an abundance gradient from the inner to the outer part of the bay. Many taxa had maximum abundances at outer stations and minimum at inner stations. Multivariate correlations between environmental variables and soft-sediment communities indicated that combinations of some metals (Cd, Cu, Sn, Pb) and grain size (mainly finer fractions, fine sands and coarse silts) were the variables that best "matched" the community patterns within Brown Bay. This study indicated that there were significant correlations between the presence of contaminants and the distribution and composition of soft-sediment communities over very small spatial scales.     

Rapid expansion of a cylindrical cavity in a rate-type soil

An analytical solution for rapid (undrained) expansion of a long cylindrical cavity from initially zero radius is presented. A simple rate-type soil model possessing the salient features of the modified Cambridge critical state model is used. Application of the solution to determination of the stress field surrounding a pile immediately after driving is discussed.     

Implications of recent numerical calculations for the fission theory of the origin of the Moon

Recent numerical simulations of fission instabilities in rotating fluid masses give results which are consistent with some requirements of the fission hypothesis for the origin of the moon.     

Processes affecting deposition of sediment in a small, morphologically complex lake

Devil Lake is morphologically complex as a result of Pleistocene glacial erosion of the Frontenac Axis of the Canadian Shield. In order to assess the processes causing highly variable sedimentation in the lake, we monitored currents, suspended sediment and temperature in the lake before and during autumn overturn in 2002. Strong summer thermal stratification (stability number to 0.11 s^–1 declining with the approach of overturn) was insufficient to prevent a dynamic response in the hypolimnion to wind forcing. Superimposed on a gradual increase in suspended sediment concentration in the last weeks of stratification from less than 2 g/l to about 30 g/l were shorter-term rises lasting up to several days. Associated with these events was an increase in particle size of the sediment from a mode of 40–50 to 150–200 m ascribed to flocculation from primary particles. These events culminated in rapid (<1 h) clearing of the water associated with strong, sustained winds over the lake, especially from the southwest. After overturn, the events were more frequent, and flocculation was unable to develop as well in the more vigorous circulation. However, currents in the hypolimnion occurred throughout the period before, during and after overturn with speed related to wind speed, but direction largely independent of wind direction. The results represent an approach to understanding the nature of sedimentary processes and thus to strengthening the use of sedimentary records as proxy in environmental and paleoenvironmental assessment.     

Separable spatial modeling of spillovers and disturbances

The single spatial parameter in the spatial autoregressive model affects both the estimation of spillovers and the estimation of spatial disturbances. Consequently, the spatial autoregressive model has the undesirable property that if the degree of spatial dependence in the disturbances differs from that in the spillovers, neither may be estimated correctly. We show theoretically that the dependence structure for the spillovers and disturbances can differ and conduct a Monte Carlo experiment that verifies these findings. In contrast, estimates from a simple separable model show little bias in all the scenarios. We also show differences between the spatial autoregressive model and the separable model on five empirical examples.