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.     

Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to e_max are found to produce final states containing too many bodies with masses which are too small when e_max < 0.10, while simulations with e_max > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The e_max = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets.     

210Pb balance in Long Island Sound

A material balance is constructed for excess ²¹⁰Pb (relative to ²²⁶Ra) as a test of the retentivity of Long Island Sound for a reactive heavy metal. Excess ²¹⁰Pb is supplied to Long Island Sound chiefly by direct atmospheric deposition [1 ± 0.2(dis·min^?1)cm^?2·yr^?1]. Rivers supply less than 20% of the atmospheric flux, and other inputs, from open ocean waters, ²²⁶Ra decay, groundwater seepage, and sewage discharge, appear to be negligible. The total input of excess ²¹⁰Pb represents approximately the flux required to maintain the inventory of excess ²¹⁰Pb measured in sediment cores from central Long Island Sound; that is, excess ²¹⁰Pb is lost from Long Island Sound chiefly by radioactive decay. The retention of excess ²¹⁰Pb within Long Island Sound is achieved in two steps: a rapid removal of soluble ²¹⁰Pb onto suspended particles and the ongoing entrapment of particles in the basin by the residual bottom-water influx from the east.     

Stability of frosts in the solar system

Calculations on the stability of various frosts (against evaporation) for solar system objects in circular and elliptical orbits are made. It is found that the stability of these frosts is dependent on the rate of rotation of the object, the latitude of the area on the object being considered, and the eccentricity of the orbit as well as its mean distance from the Sun. These factors greatly influence the amount of solar radiation incident and reradiated from a given area on the object. The likelihood of finding these frosts on the surfaces of objects and the lifetimes of objects composed of these frosts is discussed.     

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.     

The geochemistry and provenance of Apollo 16 mafic glasses

The regolith of the Apollo 16 lunar landing site is composed mainly of feldspathic lithologies but mafic lithologies are also present. A large proportion of the mafic material occurs as glass. We determined the major element composition of 280 mafic glasses (>10 wt% FeO) from six different Apollo 16 soil samples. A small proportion (5%) of the glasses are of volcanic origin with picritic compositions. Most, however, are of impact origin. Approximately half of the mafic impact glasses are of basaltic composition and half are of noritic composition with high concentrations of incompatible elements. A small fraction have compositions consistent with impact mixtures of mare material and material of the feldspathic highlands. On the basis of major-element chemistry, we identified six mafic glass groups: VLT picritic glass, low-Ti basaltic glass, high-Ti basaltic glass, high-Al basaltic glass, KREEPy glass, and basaltic-andesite glass. These glass groups encompass 60% of the total mafic glasses studied. Trace-element analyses by secondary ion mass spectroscopy for representative examples of each glass group (31 total analyses) support the major-element classifications and groupings. The lack of basaltic glass in Apollo 16 ancient regolith breccias, which provide snapshots of the Apollo 16 soil just after the infall of Imbrium ejecta, leads us to infer that most (if not all) of the basaltic glass was emplaced as ejecta from small- or moderate-sized impacts into the maria surrounding the Apollo 16 site after the Imbrium impact. The high-Ti basaltic glasses likely represent a new type of basalt from Mare Tranquillitatis, whereas the low-Ti and high-Al basaltic glasses possibly represent the composition of the basalts in Mare Nectaris. Both the low-Ti and high-Al basaltic glasses are enriched in light-REEs, which hints at the presence of a KREEP-bearing source region beneath Mare Nectaris. The basaltic andesite glasses have compositions that are siliceous, ferroan, alkali-rich, and moderately titaniferous; they are unlike any previously recognized lunar lithology or glass group. Their likely provenance is within the Procellarum KREEP Terrane, but they are not found within the Apollo 16 ancient regolith breccias and therefore were likely deposited at the Apollo 16 site post-Imbrium. The basaltic-andesite glasses are the most ferroan variety of KREEP yet discovered.     

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.     

应用等效纬度-海拔模型进行地温及多年冻土制图

This research presents a method for permafrost mapping in discontinuous permafrost regions based on equivalent latitude/elevation concept in interior Alaska. In winter months, study site has a strong temperature inversion in air up to 700 m elevation. Air temperature data and the effects of slope, aspect and elevation were used to create an equivalent latitude/elevation model. This model was well correlated with mean annual surface temperature (0.79). In this watershed, the thawing index (It≈1 400 ℃*days) at the ground surface and snow depth do not vary greatly from south facing to north facing slopes. The primary controlled factor that determines the mean annual surface temperature was the winter surface temperature. The permafrost stability is effectively controlled by the freezing index. We determined 37.5% of Caribou-Poker Creeks Research Watershed has unstable or thawing permafrost. At least 2.1% of the permafrost in this watershed may have disappeared in the last 90 years due to climate warming. This method makes it possible to evaluate the permafrost stability in the present, past and future.     

Flexible spectral methods for the generation of random fields with power-law semivariograms

Random field generators serve as a tool to model heterogeneous media for applications in hydrocarbon recovery and groundwater flow. Random fields with a power-law variogram structure, also termed fractional Brownian motion (fBm) fields, are of interest to study scale dependent heterogeneity effects on one-phase and two-phase flow. We show that such fields generated by the spectral method and the Inverse Fast Fourier Transform (IFFT) have an incorrect variogram structure and variance. To illustrate this we derive the prefactor of the fBm spectral density function, which is required to generate the fBm fields. We propose a new method to generate fBm fields that introduces weighting functions into the spectral method. It leads to a flexible and efficient algorithm. The flexibility permits an optimal choice of summation points (that is points in frequency space at which the weighting function is calculated) specific for the autocovariance structure of the field. As an illustration of the method, comparisons between estimated and expected statistics of fields with an exponential variogram and of fBm fields are presented. For power-law semivariograms, the proposed spectral method with a cylindrical distribution of the summation points gives optimal results.