Reminiscing my sixty year pursuit of the physics of the Sun and the Galaxy

Reminiscing begins with childhood and passes on to student days through graduate school and the first real contact with research.Then early academic positions and stumbling efforts to pursue my ideas.The first significant progress came as a research associate with Prof.W.M.Elsasser at the University of Utah,beginning with an introduction to magnetohydrodynamics and the generation of the geomagnetic field through induction in the liquid metal core of Earth.A move to the University of Chicago to work with Prof.J.A.Simpson,on the implications of cosmic ray variations and interplanetary magnetic fields,led to the theory of coronal expansion and the solar wind and then to exploring the dynamical effects of cosmic rays on the galactic magnetic field.Spontaneous current sheets and intrinsic rapid reconnection in interlaced magnetic field line topologies were the next big project,leading up to retirement.Finally,it is a pleasure to recall my many associates,whose fresh thinking helped stimulate the daily research activities.     

Depositional environment and OSL chronology of the Homeb silt deposits, Kuiseb River, Namibia

Previous studies suggest that the Homeb silts of the Kuiseb valley, Namibia (i) accumulated in a dune-dammed lake, (ii) are end-point deposits, (iii) represent an aggrading river bed, and (iv) are slackwater deposits. Thus, they have been used alternatively as evidence of past drier conditions or past wetter conditions. Lithostratigraphic analysis of two sediment sequences at Homeb indicates sedimentation by aggradation of the Kuiseb River triggered by a transition from an arid to humid climate. OSL ages for the sequences were obtained by the SAR protocol on aliquots of 9.6-mm and 4.0-mm diameter and on single grains. Four-millimeter aliquot minimum ages closely approximate the single-grain minimum ages and are younger than 9.6-mm aliquot minimum and central ages. Based on these results, the small-aliquot (4-mm) approach appears to provide ages comparable to those obtained by the more laborious and time-consuming single-grain method. Minimum ages indicate rapid deposition of the Homeb Silts in at least two episodes centered at 15 ka and 6 ka during climate transitions from arid to humid. Flash floods eroded the valley fills during slightly more arid conditions.     

A simple method to model the reduced environment of lake bottom sapropel formation

The Kambala and Barchin brackish lakes(Baraba steppe,southern West Siberia) contain an organic-rich sapropel layer that was formed in oxygen-depleted waters.We measured the bulk sediment elemental composition,the water chemistry and determined the mineralogical composition and predominant biota species(Diatoms and Cyanobacteria in phytoplankton community respectively) in the lakes.The result indicates that the first lake has a siliceous type of sapropel and the second a carbonaceous one.A computer thermodynamic model was developed for chemical interaction in water-bottom sediment systems of the Kambala and Barchin Lakes.The surface sodium bicarbonate waters are supersaturated with respect to calcite,magnesite(or low Mg-calcite),quartz and chlorite with minor strontianite,apatite and goethite(pH8.9-9.3,Eh 0.3 V).Nevertheless,it is shown that during sapropel deposition,deep silt waters should be anoxic(Eh0 V).The virtual component CH_2O has been used to create an anoxic environment suitable for pyrite formation due to the biotic community impact and abiotic reduction.Thermodynamic calculation has shown that silt water is not necessarily euxinic(anoxic and sulfidic).Depending on Eh,sulfate sulfur can dominate in solution,causing the formation of gypsum together with pyrite.An attempt was made to find a reason for solution supersaturation with respect to Ca and Mg ions due to their complexation with humic acids.     

Evaluation of an alternative method for numerically modeling nonhydrostatic flows over irregular terrain

Modeling nonhydrostatic atmospheric flow requires the solution of the vertical equation of motion and a prognostic or diagnostic equation for pressure. If the nonhydrostatic components of the flow are relatively small, they can be approximated and incorporated into a purely hydrostatic model, which usually is conceptually simpler and computationally more efficient. A method to do this for a linear model of local thermally-induced circulations is further developed and adapted to a non-linear numerical model of the neutral atmospheric boundary layer. A hydrostatic model and the quasi-nonhydrostatic version were used to simulate neutral flow over simple terrain features. One set of observations taken over a simple change in roughness and another set taken over a change in both roughness and terrain were simulated by both models to assess the capabilities of the quasi-nonhydrostatic technique.It is found that (as expected) the pressure deviation from the hydrostatic state is negligible for the roughness change, but it is an important aspect of neutral flow over terrain. Thus, for flow encountering a simple roughness change, the hydrostatic approximation is good, even for small horizontal scales. However, the quasi-nonhydrostatic model qualitatively produces the features in the observations for flow over a terrain change that the hydrostatic model cannot produce.Journal Paper No. J-12737 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa. Project No. 2779.     

On facies belts and facies mosaics: Holocene isolated platforms,South China Sea

Although the general criteria for recognition and environmental interpretation of different carbonate facies are well‐established, a predictive understanding of the areal extent and spatial patterning of facies bodies and why they might organize into facies belts or facies mosaics is poorly constrained. To explore patterns and process dynamics of facies on isolated carbonate platforms, quantitative analysis of thematic maps derived from remote sensing images of 27 Holocene atolls of the Paracel and Spratly chains in the South China Sea explores variability within and among platforms. On these systems, most annular shelf‐margin reefs are less than 500 m wide on both chains; inboard of the reefs, reef sand aprons range up to 500 m (Spratlys) and 1000 m (Paracels) wide. Around individual platforms, Spratly Chain sand apron widths are wider to the north‐west, whereas apron widths in the Paracel Chain are more symmetrical; collectively, data indicate log‐normal width‐exceedance probability distributions. Platform‐interior patch reefs include area‐exceedance probability distributions and gap size distributions (lacunarity) consistent within chains, but distinct between the chains. To understand the processes underlying distinct distributions, simulations explored distinct growth scenarios. Results suggest that differences may represent distinct process classes: proportional growth processes with multiplicative random effects (reef sand aprons – belts), versus non‐linear, size‐proportional growth of randomly aged and distributed elements (patch reefs – mosaics). The probabilistically distinct sizes and spatial patterns of geomorphic elements within these general process classes are interpreted to represent ‘variations on themes’ related to the different impacts of tropical storms, winter cold fronts and circulation in each chain. The results highlight fundamentally different growth patterns impacting the sizes and distribution of facies belts and mosaics on isolated carbonate platforms. Because these types of bodies ultimately construct stratigraphy, the themes could be applied to understand and predict variability in the architecture of subsurface reservoir analogues.     

Use of archaeology to date liquefaction features and seismic events in the New Madrid seismic zone,Central United States

Prehistoric earthquake-induced liquefaction features occur in association with Native American occupation horizons in the New Madrid seismic zone. Age control of these liquefaction features, including sand-blow deposits, sand-blow craters, and sand dikes, can be accomplished by extensive sampling and flotation processing of datable materials as well as archaeobotanical analysis of associated archaeological horizons and pits. This approach increases both the amount of carbon for radiocarbon dating and the precision dating of artifact assemblages. Using this approach, we dated liquefaction features at four sites northwest of Blytheville, Arkansas, and found that at least one significant earthquake occurred in the New Madrid seismic zone between A.D. 1180 and 1400, probably about A.D. 1300 ± 100 yr. In addition, we found three buried sand blows that formed between 3340 B.C. and A.D. 780. In this region where very large to great earthquakes appear to be closely timed, archaeology is helping to develop a paleoearthquake chronology for the New Madrid seismic zone. © 1996 John Wiley & Sons, Inc.