MANGALA VALLES, MARS: ASSESSMENT OF EARLY STAGES OF FLOODING AND DOWNSTREAM FLOOD EVOLUTION

The Mangala Valles system is an ∼ ∼900 km fluvially carved channel system located southwest of the Tharsis rise and is unique among the martian outflow channels in that it heads at a linear fracture within the crust as opposed to a collapsed region of chaos as is the case with the circum-Chryse channels. Mangala Valles is confined within a broad, north–south trending depression, and begins as a single valley measuring up to 350 km wide that extends northward from a Memnonia Fossae graben, across the southern highlands toward the northern lowlands. Approximately 600 km downstream, this single valley branches into multiple channels, which ultimately lose their expression at the dichotomy boundary. Previous investigations of Mangala Vallis suggested that many of the units mapped interior to the valley were depositional, related to flooding, and that a minimum of two distinct periods of flooding separated by tens to hundreds of millions of years were required to explain the observed geology. We use infrared and visible images from the THermal EMission Imaging System (THEMIS), and topographic data from the Mars Orbiting Laser Altimeter (MOLA), to investigate the nature of the units mapped within Mangala Vallis. We find that the geomorphology of the units, as well as their topographic and geographic distribution, are consistent with most of them originating from a single assemblage of volcanic flow deposits, once continuous with volcanic flows to the south of the Memnonia Fossae source graben. These flows resurfaced the broad, north–south trending depression into which Mangala Vallis formed prior to any fluvial activity. Later flooding scoured and eroded this volcanic assemblage north of the Mangala source graben, resulting in the present distribution of the units within Mangala Vallis. Additionally, our observations suggest that a single period of catastrophic flooding, rather than multiple periods separated by tens to hundreds of millions of years, is consistent with and can plausibly explain the interior geology of Mangala Vallis. Further, we present a new scenario for the source and delivery of water to the Mangala source graben that models flow of groundwater through a sub-cryosphere aquifer and up a fracture that cracks the cryosphere and taps this aquifer. The results of our model indicate that the source graben, locally enlarged to a trough near the head region of Mangala, would have required less than several days to fill up prior to any spill-over of water to the north. Through estimates of the volume of material missing from Mangala (13,000–20,000 km³), and calculation of mean discharge rates through the channel system (∼ ∼5 × 10⁶ m³ s⁻¹), we estimate that the total duration of fluvial activity through the Mangala Valles was 1–3 months.     

Structure and dynamics of the Bok globule B335

CO maps of the Bok globule B335 are presented and used to derive its density profile, mass distribution, and rotational velocity structure. It is found that the cloud is in nearly hydrostatic equilibrium with a density profile that varies roughly as r^?1 in the core and r^?3 in the envelope. The observed rotation is unimportant in the force balance at the present stage of evolution.     

A moment method for radiative transfer in an anisotropically-scattering slab medium with space-dependent albedoω(x)

The moment method is used to solve the radiative transfer problem in an anisotropic scattering plane medium with arbitrary space-dependent albedo(x). The results are compared with those obtained recently by Cengel and Özisik.     

On ‘bimodality of the solar cycle’ and the duration of cycle 21

The period-growth dichotomy of the solar cycle predicts that cycle 21, the present solar cycle, will be of long duration (>133 mo), ending after July 1987. Bimodality of the solar cycle (i.e., cycles being distributed into two groups according to cycle length, based on a comparison to the mean cycle period) is clearly seen in a scatter diagram of descent versus ascent durations. Based on the well-observed cycles 8–20, a linear fit for long-period cycles (being a relatively strong inverse relationship that is significant at the 5% level and having a coefficient of determination r ² 0.66) suggests that cycle 21, having an ascent of 42 mo, will have a descent near 99 mo; thus, cycle duration of about 141 mo is expected. Like cycle 11, cycle 21 occurs on the downward envelope of the sunspot number curve, yet is associated with an upward first difference in amplitude. A comparison of individual cycle, smoothed sunspot number curves for cycles 21 and 11 reveals striking similarity, which suggests that if, indeed, cycle 21 is a long-period cycle, then it too may have an extended tail of sustained, low, smoothed sunspot number, with cycle 22 minimum occurring either in late 1987 or early 1988.     

On the proposed associations of solar neutrino flux with solar particles,cosmic rays,and the solar activity cycle

It has been proposed that the observed solar neutrino flux exhibits important correlations with solar particles, galactic cosmic rays, and the sunspot cycle, with the latter correlation being opposite in phase and lagging behind the sunspot cycle by about one year. Re-examination of the data-available interval 1971–1981, employing various tests of statistical significance, however, suggests that such a claim is, at present, unwarrantable. For example, on the associations of solar neutrino flux and cosmic-ray flux with the Ap geomagnetic index, neither were found to be statistically significant (at the 95% level of confidence), regardless of the choice of lag (-1, 0, or +1 yr). Presuming linear fits, all correlations with Ap had coefficients of determination (r ², where r is the linear correlation coefficient) less than 16%, meaning that 16% of the variation in the selected test parameters could be explained by the variation in Ap. Similarly, on the associations of solar neutrino flux and cosmic ray flux with sunspot number, only the latter association proved to be of statistical importance. Using the best linear fits, the correlation between yearly averages of solar neutrino flux and sunspot number had r ² 19%, the correlation between weighted moving averages (of order 5) of solar neutrino flux and sunspot number had r ² 45%, and the correlation between cosmic-ray flux and sunspot number had r ² 76%, all correlations being inverse associations. Solar neutrino flux was found not to correlate strongly with cosmic-ray flux, and the Ap geomagnetic index was found not to correlate strongly with sunspot number.     

Accuracy of close binary mass determinations from parallaxes

Absolute masses for W Ursae Majoris and Algol-type close binaries can be determined from their parallax, if observed, and the relative sizes of the stars and their mass ratio, obtained from a light curve solution. An error propagation study compares the typical order of magnitude of the various terms involved, and shows how accurate parallaxes have to be in order to make the procedure work, i.e., making the parallax term not larger than the combined non-parallax terms, and producing reasonably low mass errors. Some comments are made on the possibilities with respect to the HIPPARCOS program.Communication presented at the International Conference on Astrometric Binaries, held on 13–15 June, 1984, at the Remeis-Sternwarte Bamberg, Germany, to commemorate the 200th anniversary of the birth of Friedrich Wilhelm Bessel (1784–1846).     

U–Pb SHRIMP zircon geochronology and T–t–d history of the Kampa Dome, southern Tibet

Structural, petrographic and geochronologic studies of the Kampa Dome provide insights into the tectonothermal evolution of orogenic crust exposed in the North Himalayan gneiss domes of southern Tibet. U–Pb ion microprobe dating of zircons from granite gneiss exposed at the deepest levels within the dome yields concordia ²⁰⁶Pb/²³⁸U age populations of 506 ± 3 Ma and 527 ± 6 Ma, with no evidence of new zircon growth during Himalayan orogenesis. However, the granite contains penetrative deformation fabrics that are also preserved in the overlying Paleozoic strata, implying that the Kampa granite is a Cambrian pluton that was strongly deformed and metamorphosed during Himalayan orogenesis. Zircons from deformed leucogranite sills that cross-cut Paleozoic metasedimentary rocks yield concordant Cambrian ages from oscillatory zoned cores and discordant ages ranging from ca. 491–32 Ma in metamict grains. Since these leucogranites clearly post-date the metasedimentary rocks they intrude, the zircons are interpreted as xenocrysts that are probably derived from the Kampa granite. The Kampa Dome formed via a series of progressive orogenic events including regional ~ N–S contraction and related crustal thickening (D₁), predominately top-to-N ductile shearing and crustal extension (D₂), top-to-N brittle–ductile faulting and related folding on the north limb of the dome, localized top-to-S faulting on the southern limb of the dome, and crustal doming (D₃), and continued N–S contraction, E–W extension and doming (D₄). Structural and geochronologic variability amongst adjacent North Himalayan gneiss domes may reflect changes in the magnitude of crustal exhumation along the North Himalayan antiform, possibly relating to differences in the mid-crustal geometry of the exhuming fault systems.     

Thermal evolution and physics of melt extraction on the ureilite parent body

We develop a physical model of the thermal history of the ureilite parent body (UPB) that numerically tracks the history of its heating, hydration, dehydration, partial melting and smelting as a function of its formation time and the initial values of its composition, formation temperature and water ice content. Petrologic and chemical data from the main group (non-polymict) ureilite meteorites, which sample the interior of the UPB between depths corresponding to pressures in the range 3-10 MPa, are used to constrain the model. We find that to achieve the ∼30% melting inferred for ureilites from all sampled depths, the UPB must have had a radius between ∼80 and ∼130 km and must have accreted about 0.55 Ma after CAI formation. Melting began in the body at ∼1 Ma after CAI, and the time at which 30% melting was reached varied with depth in the asteroid but was always between ∼4.5 and ∼5.8 Ma after CAI. The total rate at which melt was produced in the UPB varied from more than 100 m³ s⁻¹ in the very early stages of melting at ∼1 Ma after CAI to ∼5 m³ s⁻¹ between 2 and 3 Ma after CAI, decreasing to extremely small values as the end of melting was approached beyond ∼5 Ma. Although the initial period of high melt production occupied only a short time around 1 Ma after CAI, it corresponded to ∼half (16%) of total silicate melting, and all strictly basaltic (i.e. plagioclase-saturated) melts must have been produced during this period.A very efficient melt transport network, consisting of a hierarchy of veins and larger pathways (dikes), developed quickly at the start of melting, ensuring rapid (timescales of months) transport of any single parcel of melt to shallow levels, thus ensuring that chemical interaction between melts and the rocks through which they subsequently passed was negligible. Volatile (mainly carbon monoxide) production due to smelting began at the start of silicate melting in the shallowest parts of the UPB and at later times at greater depths. Except at the very start and very end of melting, the volatile content of the melts produced was always high - generally between 15 and 35 mass % - and most of the melt produced was erupted at the surface of the UPB with speeds well in excess of the escape velocity and was lost into space. However, we show that 30% melting at the 3 MPa pressure level was only possible if ∼15% of the total melt produced in the asteroid was retained as a small number (∼5) of very extensive, sill-like intrusions centered at a depth of ∼7 km below the surface, near the base of the ∼8 km thick outer crust of the asteroid that was maintained at temperatures below the basalt solidus by conductive heat loss to the surface. The horizontal extents of these sills occupied about 75% of the surface area of the UPB, and the sills acted as buffers between the steady supply of melt from depth and the intermittent explosive eruption of the melt into space. We infer that samples from these intrusions are preserved as the rare feldspathic (loosely basaltic) clasts in polymict ureilites, and show that the cooling histories of the sills are consistent with these clasts reaching isotopic closure at ∼5 Ma after CAI, as given by ²⁶Al-²⁶Mg, ⁵³Mn-⁵³Cr and Pb-Pb age dates.     

Defect properties of albite

Knowledge of the defect properties of Lunar and Mercurian minerals has recently become important, with the advent of models which attempt to explain the formation of the thin exosphere of these celestial bodies. Here, we have calculated the formation energies of sodium and oxygen vacancies in the mineral albite (NaAlSi₃O₈), as well as the Schottky defect energy for the removal of a Na₂O unit. We have employed both the supercell and Mott–Littleton approaches, using Kohn–Sham density functional theory and classical interatomic potential methods. As well as reporting the defect energies and structures, we comment upon the relative merits of the methods used.     

Short-term and interannual variability of redox-sensitive chemical parameters in hypoxic/anoxic bottom waters of the Chesapeake Bay

A combination of CTD casts, discrete bottle sampling and in situ voltammetric microelectrode profiling was used to examine changing redox conditions in the water column at a single station south of the Bay Bridge in the upper Chesapeake Bay in late July/early August, 2002–2005. Short-term (2–4 h) fluctuations in the oxic/suboxic/anoxic interface were documented using in situ voltammetric solid-state electrodes. Profiles of dissolved oxygen and sulfide revealed tidally-driven vertical fluctuations of several meters in the depth and thickness of the suboxic zone. Bottom water concentrations of sulfide, Mn²⁺ and Fe²⁺ also varied over the tidal cycle by approximately an order of magnitude. These data indicate that redox species concentrations at this site varied more due to physical processes than biogeochemical processes. Based on analysis of ADCP data, tidal currents at this station were strongly polarized, with the principal axis of tidal currents aligned with the mainstem channel. Together with the chemical data, the ADCP analysis suggests tidal flushing of anoxic bottom waters with suboxic water from north of the site. The present study is thus unique because while most previous studies have focused on processes across relatively stable redox interfaces, our data clearly demonstrate the influence of rapidly changing physical mixing processes on water column redox chemistry.Also noted during the study were interannual differences in maximum bottom water concentrations of sulfide, Mn²⁺ and Fe²⁺. In 2003, for example, heavy spring rains resulted in severe hypoxia/anoxia in June and early July. While reported storm-induced mixing in late July/early August 2003 partially alleviated the low-oxygen conditions, bottom water concentrations of sulfide, Mn²⁺ and Fe²⁺ were still much higher than in the previous year. The latter implies that the response time of the microbial community inhabiting the suboxic/anoxic bottom waters to changing redox conditions is slow compared to the time scale of episodic mixing events. Bottom water concentrations of the redox-sensitive chemical species should thus be useful as a tracer to infer prior hypoxic/anoxic conditions not apparent from ambient oxygen levels at the time of sampling.