Grid resolution study of ground water flow and transport

Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution.     

Error analyses of resonant orbits for geodesy

An error analysis of resonant orbits for geodesy indicates that attempts to use resonance to recover high order geopotential coefficients may be seriously hampered by errors in the geopotential. This effect, plus the very high correlations (up to .99) of the resonant coefficients with each other and the orbital period in single satellite solutions, makesindividual resonant orbits of limited value for geodesy. Multiple-satellite, single-plane solutions are only a slight improvement over the single satellite case. Accurate determination of high order coefficients from low altitude resonant satellites requires multiple orbit planes and small drift-periods to reduce correlations and effects of errors of non-resonant geopotential terms. Also, the effects of gravity model errors on low-altitude resonant satellites make the use of tracking arcs exceeding two to three weeks of doubtful validity. Because high-altitude resonant orbits are less affected by non-resonant terms in the geopotential, much longer tracking arcs can be used for them.     

On solar system nomenclature

Arguments are presented for naming topographic features on other solar system objects after human beings other than astronomers; and to institute a more consistent scheme for Jovian satellite nomenclature.     

Influence of calibration methodology on ground water flow predictions

We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).     

Relict silicate inclusions in extraterrestrial chromite and their use in the classification of fossil chondritic material

Chromite is the only common meteoritic mineral surviving long-term exposure on Earth, however, the present study of relict chromite from numerous Ordovician (470 Ma) fossil meteorites and micrometeorites from Sweden, reveals that when encapsulated in chromite, other minerals can survive for hundreds of millions of years maintaining their primary composition. The most common minerals identified, in the form of small (<1-10 μm) anhedral inclusions, are olivine and pyroxene. In addition, sporadic merrillite and plagioclase were found.Analyses of recent meteorites, holding both inclusions in chromite and corresponding matrix minerals, show that for olivine and pyroxene inclusions, sub-solidus re-equilibration between inclusion and host chromite during entrapment has led to an increase in chromium in the former. In the case of olivine, the re-equilibration has also affected the fayalite (Fa) content, lowering it with an average of 14% in inclusions. For Ca-poor pyroxene the ferrosilite (Fs) content is more or less identical in inclusions and matrix. By these studies an analogue to the commonly applied classification system for ordinary chondritic matrix, based on Fa in olivine and Fs in Ca-poor pyroxene, can be established also for inclusions in chromite. All olivine and Ca-poor pyroxene inclusions (>1.5 μm) in chromite from the Ordovician fossil chondritic material plot within the L-chondrite field, which is in accordance with previous classifications. The concordance in classification together with the fact that inclusions are relatively common makes them an accurate and useful tool in the classification of extraterrestrial material that lacks matrix silicates, such as fossil meteorites and sediment-dispersed chromite grains originating primarily from decomposed micrometeorites but also from larger impacts.     

Glacial and environmental changes over the last 60 000 years in the Polar Ural Mountains,Arctic Russia,inferred from a high‐resolution lake record and other observations from adjacent areas

Our knowledge about the glaciation history in the Russian Arctic has to a large extent been based on geomorphological mapping supplemented by studies of short stratigraphical sequences found in exposed sections. Here we present new geochronological data from the Polar Ural Mountains along with a high‐resolution sediment record from Bolshoye Shchuchye, the largest and deepest lake in the mountain range. Seismic profiles show that the lake contains a 160‐m‐thick sequence of unconsolidated lacustrine sediments. A well‐dated 24‐m‐long core from the southern end of the lake spans the last 24 cal. ka. From downward extrapolation of sedimentation rates we estimate that sedimentation started about 50–60 ka ago, most likely just after a large glacier had eroded older sediments from the basin. Terrestrial cosmogenic nuclide (TCN) exposure dating (¹⁰Be) of boulders and Optically Stimulated Luminescence (OSL) dating of sediments indicate that this part of the Ural Mountains was last covered by a coherent ice‐field complex during Marine Isotope Stage (MIS) 4. A regrowth of the glaciers took place during a late stage of MIS 3, but the central valleys remained ice free until the present. The presence of small‐ and medium‐sized glaciers during MIS 2 is reflected by a sequence of glacial varves and a high sedimentation rate in the lake basin and likewise from ¹⁰Be dating of glacial boulders. The maximum extent of the mountain glaciers during MIS 2 was attained prior to 24 cal. ka BP. Some small present‐day glaciers, which are now disappearing completely due to climate warming, were only slightly larger during the Last Glacial Maximum (LGM) as compared to AD 1953. A marked decrease in sedimentation rate around 18–17 cal. ka BP indicates that the glaciers then became smaller and probably disappeared altogether around 15–14 cal. ka BP.