Towards a nonperturbation transport theory in hetereogeneous aquifers

A large body of existing theories of flow and contaminant transport in aquifers ignore the presence of recharge, eliminate the boundary conditions, neglect transient conditions in groundwater flow, conceive hydraulic gradients as linear, and require parameter variability to be stationary and Gaussian. The most outstanding and difficult to justify assumption is the subjective small size of the stochastic terms (i.e., small perturbation methods), which usually is forced by considering the logarithm of the hydraulic conductivity. Several problems in flow and contaminant subsurface hydrology, such as the enhanced dispersion parameters with plume size or time after injection, remain to be observed in the light of a stochastic theory that allows a more realistic consideration of physical and hydrologic properties. In this article, an attempt is made to reformulate a contaminant transport equation (the variable dispersion equation, VDE) with transport parameters in terms of regional hydrologic and aquifer hydraulic properties, such as recharge rate, spatially random transmissivity, hydraulic gradient, aquifer thickness, and soil porosity. Subsequently, a general analytic procedure, the method of decomposition, is used to derive a solution to the VDE. This procedure does not require small perturbation, logarithmic transformations, or specific probability law assumptions. Comparison tests with existing theoretical and field results are given. The tests illustrate the enhanced dispersion and shifting concentration effects produced by the variable dispersion equation. Finally a generalization of the method to nonstationary dispersion in three-dimensional domains is proposed.     

Trace Element and Isotope Geochemistry of the Volcanic Rocks of Bequia, Grenadine Islands, Lesser Antilles Arc: a Study of Subduction Enrichment and Intra-crustal Contamination

Pliocene volcanics on the island of Bequia comprise two interbeddedsuites of basalts and andesites. The isotopically homogeneoussuite (IHS) has a limited range of Sr—Nd—Pb isotopes(⁸⁷Sr/⁸⁶Sr 0.7040–0.7046, ¹⁴³ Nd/¹⁴⁴ Nd 0.5130 and ²⁰⁶Pb/²⁰⁴Pb 19.36–19.51), and mantle-like ¹⁸O values (5.5in clinopyroxene). The isotopically diverse suite (IDS) is characterizedby much wider ranges of radiogenic isotopes (⁸⁷ Sr/⁸⁶Sr 0.7048–0.7077,¹⁴³ Nd/¹⁴⁴ Nd 0.5128–0.5123 and ²⁰⁶ Pb/²⁰⁴ Pb 19.7–20.2),in which all of the Sr and Pb ratios are higher and Nd ratiosare lower than those of the IHS. The IDS is also characterizedby high ¹⁸ O values, up to 7.6 in clinopyroxene. The Sr andPb isotope ratios are too high, and the Nd isotope ratios aretoo low in the IDS for any of these lavas to be derived fromunmodified depleted mantle. Both suites are petrologically very similar and their majorelement compositions and phenocryst contents suggest that theywere formed largely by fractional crystallization of a hydroustholeiitic melt at pressures <3 kbar. The isotopic ratiosand enrichments in large ion lithophile elements (LILE), andto some extent light rare earth elements (LREE), as comparedwith mid-ocean ridge basalts (MORB), of the IHS lavas suggestthat they were derived from a depleted mantle source which hadbeen re-enriched by the addition of 1–4% of a subductioncomponent. This component probably comprised a mixture of dehydrationfluids, and perhaps minor siliceous melts, released from subductingsediments and mafic crust. The extreme isotopic ranges, largeenrichments in incompatible elements, more fractionated LREEpatterns and higher ¹⁸ O values of the IDS lavas are interpretedas resulting from 10–55% assimilation—fractionalcrystallization of sediments, derived from the Guyana Shield,which are present in the arc crust, by IHS type melts. KEY WORDS: trace elements; radiogenic isotopes; arc lavas; Lesser Antilles ^*Corresponding author.     

Total ozone variations in the tropical belt: An application for quality of ground based measurements

Summary The study of the regime of ozone variations in the huge tropical belt (25° S to 25° N), which are, in general, very small and zonally nearly symmetric, permits to establish a statistical model for estimating the ozone deviations using Total Ozone Mapping Spectrometer (TOMS) data. The equatorial stratospheric winds at 25 and 50hPa and the solar flux at 10.7 cm are used as major predictors and the linear trend was also estimated. The 10m/sec stratospheric wind change is related to1.2% ozone change at the equator, to practically no change in the 8–15° belts and up to 1.4% change with opposite phase over the tropics in spring but nearly zero change in fall. The solar cycle related amplitude is about 1.4% per 100 units of 10.7 cm solar flux. The ozone trends are negative: not significant over the equator and about –2% per decade (significant at 95% level) over the tropics. The latter could have been enforced by the 2 to 4% lower ozone values during 1991–1993, part of which might be related to the effects of the Mt. Pinatubo eruption, but might also be due to the strong QBO. The estimated deviations are verified versus reliable observations and the very good agreement permits applying the model for quantitative quality control of the reported ozone data from previous years. The standard deviation of the difference between observed ozone deviations and those estimated from the model is only 0.9–1.6% for yearly mean, that means instruments used for total ozone observations in the tropical belt should have systematic error of less than 1%. Cases when the discrepancies between the model and reported observations at a given station exceed 2–3% for time interval of 2 or more years should be verified.With 17 Figures     

Photometric investigation of the evolved contact binary system DN Aur

We present here the photometric light curve analyses of the eclipsing binary star DN Aur. The CCD photometry, performed at the Behlen observatory using the 0.76 m automated telescope gave 646 individual data points inV andR bandpass filters. From this data we have determined a new epoch and an orbital period of 0.6168891 days. The published spectral classification is F3.The Wilson-Devinney model was used to derive the photometric solutions. DN Aur is a W UMa type contact binary system. The mass ratio, (q=m ₂/m ₁0.210, where star 2 is eclipsed at the primary minimum) suggests that the system has A-type configuration. The computed light curve has a third light of about 22 percent and a total eclipse in the secondary minimum. A solution with a cool spot on the secondary component is also found. We recommend spectroscopic study of DN Aur even though the light curve analysis show it to be a single line spectroscopic system. Generally contact systems of spectral type F3 have periods ranging from 0.25 to 0.5 days. The longer period of DN Aur suggests that it is an evolved contact system with case A mass transfer.     

The contact binary system YY Eri

Relatively new and unanalysed photometric data-sets of the contact binary system YY Eri are presented. The light curves have been analysed using information limit optimization techniques, and the Binary Maker program of Bradstreet (1992).Comparison of the results allows insight into determinacy questions affecting the W UMa type of light curve, which relates to our general understanding of contact binaries. A limiting contact configuration cannot be ruled out on the basis of empirical data-analysis alone.Very careful analysis of the photospheric flux distribution over the surface of the Sun may help establish appropriate values of the gravity brightening parameter for cool dwarf stars. At present, however, empirical photometric information on contact binaries is non-discriminatory: i.e. alternative cosmogonies can find alternative support from available evidence. The implication is then for more and better observational data to allow better real independent parameter determination.     

Photochemistry and Vertical Transport in Io's Atmosphere and Ionosphere

We present an updated model for the photochemistry of Io's atmosphere and ionosphere and use this model to investigate the sensitivity of the chemical structure to vertical transport rates. SO₂is assumed to be the dominant atmospheric gas, with minor molecular sodium species such as Na₂S or Na₂O released by sputtering or venting from the surface. Photochemical products include SO, O₂, S, O, Na, NaO, NaS, and Na₂. We consider both “thick” and “thin” SO₂atmospheres that encompass the range allowed by recent HST and millimeter-wave observations, and evaluate the possibility that O₂and/or SO may be significant minor dayside constituents and therefore likely dominant nightside gases. The fast reaction between S and O₂limits the column abundance of O₂to ∼10⁴less than that calculated by Kumar (J. Geophys. Res.87, 1677–1684, 1982; 89(A9), 7399–7406, 1984) for a pure sulfur/oxygen atmosphere. If a significant source of NaO₂or Na₂O were supplied by the surface and mixed rapidly upward, then oxygen liberated in the chemical reactions which also liberate free Na would provide an additional source of O₂. Fast eddy mixing will enhance the transport of molecular sodium species to the exobase, in addition to increasing the vertical transport rate of ions. Ions produced in the atmosphere will be accelerated by the reduced corotation electric field penetrating the atmosphere. These ions experience collisions with the neutral gas, leading to enhanced vertical ion diffusion. The dominant ion, Na⁺, is lost primarily by charge exchange with Na₂O and/or Na₂S in the lower atmosphere and by diffusion through the ionopause in the upper atmosphere. The atmospheric column abundance of SO, O₂, and the upper atmosphere escape rates of Na, S, O, and molecular sodium species are all strong functions of the eddy mixing rate. Most atmospheric escape, including that of molecular sodium species, probably occurs from the low density “background” SO₂atmosphere, while a localized high density “volcanic” SO₂atmosphere can yield an ionosphere consistent with that detected by the Pioneer 10 spacecraft.     

A model of an exploding,radiating star in general relativity

In a previous paper Adams, Cary and Cohen (1994) presented a model of a supernova. In that paper the equations of General Relativity describing the evolution of a spherically symmetric, radiating star were solved analytically. The evolution of the star was determined by the application of boundary conditions at the center and at the edge. Due to lmitations in the presupernova model, only the very slow inward motion of an unstable, degenerate core could be considered. The solution was also limited by the need to exclude a runaway term, one that increased exponentially with time. Without the exclusion of the runaway, the luminosity would have increased without bound and the mass would have become negative.This paper presents a completely analytic solution to the equations of General Relativity describing the evolution of a Type II supernova. Professor S.E. Woosley kindly gave us data on the physical variables of a 12M ₀ presupernova star. In our model the core collapses within 1 s, leaving a 1.3M ₀ remnant. Shortly afterward 10.6M ₀ is ejected to infinity, and 0.17M ₀ is radiated away in the form of neutrinos. The distance of the edge from the center increases proportionally to the two-thirds power of the time. The luminosity decreases proportionally to the inverse four-thirds power.Although the runaway solution was modified by the exploding rather than a static envelope, it must still be excluded by adjusting initial conditions. Its character is changed from an exponential to a very large power (55) of time. The removal of a degree of freedom by this exclusion leads to physically non-sensical results such as negative luminosity. The inclusion of a term describing motion of the mantle due to neutrino interactions provides the additional degree of freedom necessary for physically reasonable results.     

The X-ray coronae of solitary late-type stars

For solar cycles 20 and 21, the longitudinal distribution of the D, G, and H-type solar flares which are related to the final phases of active region evolution, have been analysed for the northern and the southern hemispheres separately. One active zone has been found for D, G, and H-type flares, and one more active zone has been found for the H-type flares of the northern hemisphere for cycle 20. Two active zones have been found for the D and H-type flares of the northern hemisphere for cycle 21. Southern-hemisphere flares are concentrated in two active zones for cycle 20. The active zone in the northern hemisphere, which rotates with a synodic period of about 26.73 days, produced 30% of the examined D-type flares during cycle 20 and persisted in the same position during the two solar cycles, 20 and 21. The active zone in the southern hemisphere rotated with a synodic period of about 27.99 days. Only the active zone producing D-type flares persisted in the same position during the two solar cycles.     

Propagation of MHD body and surface waves in magnetically structured regions of the solar atmosphere

The fact that magnetically structured regions exist in the solar atmosphere has been known for a number of years. It has been suggested that different kinds of magnetohydrodynamic (MHD) waves can be efficiently damped in these regions and that the dissipated wave energy may be responsible for the observed enhancement in radiative losses. From a theoretical point of view, an important task would be to investigate the propagation and dissipation of MHD waves in these highly structured regions of the solar atmosphere. In this paper, we study the behavior of MHD body and surface waves in a medium with either a single or double (slab) magnetic interface by use of a nonlinear, two-dimensional, time-dependent, ideal MHD numerical model constructed on the basis of a Lagrangean grid and semi-implicit scheme. The processes of wave confinement and wave energy leakage are discussed in detail. It is shown that the obtained results depend strongly on the type of perturbations imposed on the interface or slab and on the plasma parameter, . The relevance of the obtained results to the heating problem of the upper parts of the solar atmosphere is also discussed.     

White-light enhancements and small-scale chromospheric activities in an active region

A microflare or a group of Ellerman bombs was found to be associated with several points of white-light enhancements. These points had similar sizes as facular points (d 0.3 arc sec). Temporal evolution of these activities is described. Origins of these activities are discussed to be deeply seated excess heating in chromospheric and photospheric levels.