Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Stratigraphic evidence for multiple drainings of glacial Lake Missoula along the Clark Fork River, Montana, USA

Glacial Lake Missoula, a source of Channeled Scabland flood waters, inundated valleys of northwest Montana to altitudes of ∼ 1265 m and to depths of  >600 m, as evidenced by shorelines and silty lacustrine deposits. This study describes previously unrecognized catastrophic lake-drainage deposits that lie stratigraphically beneath the glacial-lake silts. The unconsolidated gravelly flood alluvium contains imbricated boulder-sized clasts, cross-stratified gravel with slip-face heights of 2-> 35 m, and 70- to 100-m-high gravel bars which all indicate a high-energy, high-volume alluvial environment. Gravel bars and high scablands were formed by catastrophic draining of one or possibly more early, high lake stands (1200-1265 m). Most glacial-lake silt, such as the Ninemile section, was deposited stratigraphically above the earlier deposits, represents a lower lake stand(s) (1050-1150 m), and was not deposited in lake(s) responsible for the highest discharge events. The glaciolacustrine silt-covered benches are incised by relict networks of valleys formed during the drainage of the last glacial lake. Significant erosion associated with the last lake draining was confined to the inner Clark Fork River canyon.     

Observations Of Nocturnal Drainage Flow In A Shallow Gully

The formation of cold air drainage flows in a shallow gully is studied during CASES-99 (Cooperative Atmosphere-Surface Exchange Study). Fast and slow response wind and temperature measurements were obtained on an instrumented 10-m tower located in the gully and from a network of thermistors and two-dimensional sonic anemometers, situated across the gully. Gully flow formed on clear nights even with significant synoptic flow. Large variations in surface temperature developed within an hour after sunset and in situ cooling was the dominant factor in wind sheltered locations. The depth of the drainage flow and the height of the down-gully wind speed maximum were found to be largest when the external wind speed above the gully flow is less than 2 m s^-1. The shallow drainage current is restricted to a depth of a few metres, and is deepest when the stratification is stronger and the external flow is weaker. During the night the drainage flow breaks down, sometimes on several occasions, due to intermittent turbulence and downward fluxes of heat and momentum. The near surface temperature may increase by 6 ° C in less than 30 min due to the vertical convergence of downward heat flux. The mixing events are related to acceleration of the flow above the gully flow and decreased Richardson number. These warming events also lead to warming of the near surface soil and reduction of the upward soil heat flux. To examine the relative importance of different physical mechanisms that could contribute to the rapid warming, and to characterize the turbulence generated during the intermittent turbulent periods, the sensible heat budget is analyzed and the behaviour of different turbulent parameters is discussed.     

Anasazi (Pre-Columbian Native-American) Migrations During The Middle-12Th and Late-13th Centuries – Were they Drought Induced?

Severe droughts in the middle-12th and late-13th centuries appear to have affected Anasazi (pre-Columbian Native American) populations. During the first drought most of the great houses in the central San Juan Basin were vacated; the second drought resulted in the abandonment of the Four Corners region. During the first drought, villages may not have been completely abandoned. The multi-year drought periods probably were characterized by reductions in both winter and summer precipitation. Maize is dependent on winter precipitation for its germination and initial growth and on summer (monsoonal) precipitation for its continued growth. Reductions in precipitation are hypothesized to have resulted in low yields of maize, the dietary staple of the Anasazi. A comparison of historic climate data and tree-ring-based reconstructions of precipitation in the Four Corners region with tree-ring-based reconstructions of the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) indicate that severe and persistent drought in the Four Corners region occurs when the PDO is negative and the AMO is positive. Historic climate data from the greater San Juan Basin indicate that a negative PDO is characterized by reductions in both water-year and summer precipitation, reinforcing the concept that at least some multi-year droughts involved weakening of the summer monsoon with attendant decreases in the yields of maize.     

Remote estimation of surficial seafloor properties through the application Angular Range Analysis to multibeam sonar data

The variation of the backscatter strength with the angle of incidence is an intrinsic property of the seafloor, which can be used in methods for acoustic seafloor characterization. Although multibeam sonars acquire backscatter over a wide range of incidence angles, the angular information is normally neglected during standard backscatter processing and mosaicking. An approach called Angular Range Analysis has been developed to preserve the backscatter angular information, and use it for remote estimation of seafloor properties. Angular Range Analysis starts with the beam-by-beam time-series of acoustic backscatter provided by the multibeam sonar and then corrects the backscatter for seafloor slope, beam pattern, time varying and angle varying gains, and area of insonification. Subsequently a series of parameters are calculated from the stacking of consecutive time series over a spatial scale that approximates half of the swath width. Based on these calculated parameters and the inversion of an acoustic backscatter model, we estimate the acoustic impedance and the roughness of the insonified area on the seafloor. In the process of this inversion, the behavior of the model parameters is constrained by established inter-property relationships. The approach has been tested using a 300 kHz Simrad EM3000 multibeam sonar in Little Bay, NH. Impedance estimates are compared to in situ measurements of sound speed. The comparison shows a very good correlation, indicating the potential of this approach for robust seafloor characterization.     

Coming to the table: Early stakeholder engagement in marine spatial planning

From 2009 to 2011, marine spatial planning (MSP) rapidly gained visibility in the United States as a promising ocean management tool. A few small-scale planning efforts were completed in state waters, and the Obama Administration proposed a framework for large-scale regional MSP throughout the U.S. Exclusive Economic Zone. During that same time period, the authors engaged a variety of U.S ocean stakeholders in a series of dialogs with several goals: to share information about what MSP is or could be, to hear stakeholder views and concerns about MSP, and to foster better understanding between those who depend on ocean resources for their livelihood and ocean conservation advocates. The stakeholder meetings were supplemented with several rounds of in-depth interviews and a survey. Despite some predictable areas of conflict, project participants agreed on a number of issues related to stakeholder engagement in MSP: all felt strongly that government planners need to engage outsiders earlier, more often, more meaningfully, and through an open and transparent process. Equally important, the project affirmed the value of bringing unlike parties together at the earliest opportunity to learn, talk, and listen to others with whom they rarely engage.     

Regolith depth growth on an icy body orbiting Saturn and evolution of bidirectional reflectance due to surface composition changes

Using a Markov chain model, we consider the regolith growth on a small body in orbit around Saturn, subject to meteoritic bombardment, and assuming all impact ejecta are re-collected. We calculate the growth of regolith and the fractional pollution, assuming an initial pure ice body and amorphous carbon as a pollutant. We extend the meteorite flux of Cuzzi and Estrada (Cuzzi, J., Estrada, P. [1998]. Icarus 132, 1-35) to larger sizes to consider the effect of disruption of the moonlet on other moonlets in the ensemble. This is a relatively small effect, completely negligible for moonlets of 1 m radius. For the given impact model, fractional pollution reaches 22% for 1 m bodies, but only 3% for 10 m bodies, 1.7% for 20 m bodies, and 1% for 30 m bodies after 4 byr. By considering an ensemble of moonlets, which have identical cross-sections for releasing and capturing ejecta, this analysis can be extended to a model of particles in Saturn’s rings, where the calculated spectra can be compared to observed ring spectra. The measured spectral reflectance of Saturn’s rings from Cassini observations therefore constrains the size and age of the ring particles. The comparison between 1 m, 10 m, 20 m, and 30 m particles confirms that for larger ring mass, the current rings would be less polluted; for the largest particles, we expect negligible changes in the UV spectrum after 4 byr of meteoritic bombardment. We consider two end members for mixing of the meteoritic material: areal and intimate. Given the uncertainties in the actual mixing of the meteoritic infall and in its composition (as a worst case, we assume the meteoritic material is 100% amorphous carbon, intimately mixed) initially pure ice 30 m ring particles would darken after 4 byr of exposure by 15%.