Rift-related Jurassic basaltic phreatomagmatic volcanism in the central Transantarctic Mountains: precursory stage to flood-basalt effusion

 The Middle Jurassic Kirkpatrick flood basalts and comagmatic Ferrar intrusions in the Transantarctic Mountains represent a major pulse of tholeiitic magmatism related to early stages in the breakup of Gondwana. A record of the volcano-tectonic events leading to formation of this continental flood-basalt province is provided by strata underlying and only slightly predating the Kirkpatrick lavas. In the central Transantarctic Mountains, the lavas rest on widespread (≥7500 km²) tholeiitic pyroclastic deposits of the Prebble Formation. The Prebble Formation is dominated by lahar deposits and is an unusual example of a regionally developed basaltic lahar field. Related, partly fault-controlled pyroclastic intrusions cut underlying strata, and vents are represented by the preserved flanks of two small tephra cones associated with a volcanic neck. Lahar and air-fall deposits typically contain 50–60% accidental lithic fragments and sand grains derived from underlying Triassic – Lower Jurassic strata in the upper part of the Beacon Supergroup. Juvenile basaltic ash and fine lapilli consist of nonvesicular to scoriaceous tachylite, sideromelane, and palagonite, and have characteristics indicating derivation from hydrovolcanic eruptions. The abundance of accidental debris from underlying Beacon strata points to explosive phreatomagmatic interaction of basaltic magma with wet sediment and groundwater, which appears to have occurred in particular where rising magma intersected upper Beacon sand aquifers. Composite clasts in the lahar deposits exhibit complex peperitic textures formed during fine-scale intermixing of basaltic magma with wet sand and record steps in subsurface fuel-coolant interactions leading to explosive eruption. The widespread, sustained phreatomagmatic activity is inferred to have occurred in a groundwater-rich topographic basin linked to an evolving Jurassic rift zone in the Transantarctic Mountains. Coeval basaltic phreatomagmatic deposits of the Mawson and Exposure Hill Formations, which underlie exposures of the Kirkpatrick Basalt up to 1500 km to the north along strike in Victoria Land, appear to represent other parts of a regional, extension-related Middle Jurassic phreatomagmatic province which developed immediately prior to rapid outpouring of the flood basalts. This is consistent with models which assign an important role to lithospheric stretching in the generation of flood-basalt provinces. Received: 28 August 1995 / Accepted: 18 April 1996     

Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

Modeling bed-load transport of coarse sediments in the Great Bay Estuary, New Hampshire

Current, sea level and bed-load transport are investigated in the Lower Piscataqua River section of the Great Bay Estuary, New Hampshire, USA—a well-mixed and geometrically complex system with low freshwater input, having main channel tidal currents ranging between 0.5 and 2 m s⁻¹. Current and sea level forced by the M₂M₄M₆ tides at the estuarine mouth are simulated by a vertically averaged, non-linear, time-stepping finite element model. The hydrodynamic model uses a fixed boundary computational domain and accounts for flooding–drying of tidal flats by making use of a groundwater component. Inertia terms are neglected in comparison with pressure gradient and bottom friction terms, which is consistent with the observed principal dynamic balance for this section of the system. The accuracy of hydrodynamic predictions in the study area is demonstrated by comparison with four tidal elevation stations and two cross-section averaged current measurements. Simulated current is then used to model bed-load transport in the vicinity of a rapidly growing shoal located in the main channel of the lower system. Consisting of coarse sand and gravel, the shoal must be dredged every five to eight years. Two approaches are taken—an Eulerian parametric method in which nodal bed-load flux vectors are averaged over the tidal cycle and a Lagrangian particle tracking approach in which a finite number of sediment particles are released and tracked. Both methods yield pathways and accumulations in agreement with the observed shoal formation and the long-term rate of sediment accumulation in the shoal area.     

On matching the spectrum of Io: Variations in the photometric properties of sulfur-containing mixtures

The problem of comparing laboratory spectra of sulfur-containing binary mixtures with the spectrum of Io is discussed. For the satellite, the observable is the geometric albedo as a function of wavelength, whereas in the laboratory one often measures some other type of albedo. In a previous paper we demonstrated that for pure sulfur the multiplicative factor which converts the laboratory albedos to geometric albedos can be strongly wavelength dependent. The present paper demonstrates that this is also true for binary sulfur-containing mixtures. Furthermore, there is no universal conversion factor applicable to all binary mixtures, nor can the factor be interpolated for a particular mixture from the conversion factors of the two end members. The conversion factor is a function not only of the specific composition of a binary mixture, but of the relative particle size distributions of the two components, and must be measured specifically for each individual sample if a quantitative comparison between a laboratory sample and Io's surface is desired.     

Analysis of stellar occultation data: Effects of photon noise and initial conditions

A new inversion technique for obtaining temperature, pressure, and number density profiles of a planetary atmosphere from an occultation light curve is described. This technique employs an improved boundary condition to begin the numerical inversion and permits the computation of errors in the profiles caused by photon noise in the light curve. We present our assumptions about the atmosphere, optics, and noise and develop the equations for temperature, pressure, and number density and their associated errors. By inverting in equal increments of altitude, Δh, rather than in equal increments of time, Δt, the inversion need not be halted at the first negative point on the light curve as required by previous methods. The importance of the boundary condition is stressed, and a new initial condition is given. Numerical results are presented for the special case of inversion of a noisy isothermal light curve. From these results, simple relations are developed which can be used to predict the noise quality of an occultation. It is found that fractional errors in temperature profiles are comparable to those of pressure and number density profiles. An example of the inversion method is shown. Finally, we discuss the validity of our assumptions. In an appendix we demonstrate that minimum fractional errors in scale height determined from the inversion are comparable to those from an isothermal fit to a noisy isothermal light curve.     

Direct imaging of extra-solar planets with stationary occultations viewed by a space telescope

The feasibility of detecting planets outside the solar system through imaging at optical wavelengths by a telescope in space is considered. The “black” limb of the moon can be used as an occulting edge to greatly reduce the background light from the planet's star. With this technique and if certain other technical requirements can be realized, a hypothetical Jupiter-Sun system could be detected at a distance of 10 pc. For this system, a signal-to-noise ratio of 9 could be achieved in less than 20 min with a 2.4-m telescope in space. For diffraction limited optics, the required integration time is inversely proportional to the fifth power of the telescope aperture. An orbit for the telescope is described that could achieve a stationary lunar occultation of any star that would last nearly two hours, providing six times more integration time than required by the hypothetical Jupiter-Sun example.     

How to compare the surface of Io to laboratory samples

Since one does not know the photometric functions of various parts of Io, one cannot convert the observed geometric albedo of the satellite to a parameter more directly measurable in the laboratory. One must therefore convert laboratory reflectances to geometric albedos before quantitative comparisons between Io's surface and a laboratory sample are made. This procedure involves determining the wavelength dependence of the sample's photometric function. For substances such as sulfur, whose reflectance varies strongly with wavelength, it is incorrect to assume that the photometric function, and hence the ratio (laboratory reflectance/geometric albedo) is independent of wavelength. To illustrate this point, measurements of the color dependence of this ratio for sulfur are presented for the specific case in which the measured laboratory reflectance is the sample's normal reflectance. In general, unless the laboratory reflectance is precisely the geometric albedo, a wavelength-dependent correction factor must be determined before the laboratory sample can be compared quantitatively with Io's surface.     

The occultation of ϵ Gem by Mars as observed from Agassiz Station

Observations of the April 8, 1976 occultation of ? Gem by Mars made at the Agassiz Station of the Harvard College Observatory have been analyzed to yield temperature profiles of the Martian atmosphere for number densities between 10¹³ and 10¹⁵ cm^?3. Pronounced wavelike structure is evident in both immersion and emersion profiles, with a peak-tto-peak variation of up to 50°K and a vertical scale of 20 km.     

Groundwater control of mangrove surface elevation: Shrink and swell varies with soil depth

We measured monthly soil surface elevation change and determined its relationship to groundwater changes at a mangrove forest site along Shark River, Everglades National Park, Florida. We combined the use of an original design, surface elevation table with new rod-surface elevation tables to separately track changes in the mid zone (0–4 m), the shallow root zone (0–0.35 m), and the full sediment profile (0–6 m) in response to site hydrology (daily river stage and daily groundwater piezometric pressure). We calculated expansion and contraction for each of the four constituent soil zones (surface [accretion and erosion; above 0 m], shallow zone [0–0.35 m], middle zone [0.35–4 m], and bottom zone [4–6]) that comprise the entire soil column. Changes in groundwater pressure correlated strongly, with changes in soil elevation for the entire profile (Adjusted R² = 0.90); this relationship was not proportional to the depth of the soil profile sampled. The change in thickness of the bottom soil zone accounted for the majority (R² = 0.63) of the entire soil profile expansion and contraction. The influence of hydrology on specific soil zones and absolute elevation change must be considered when evaluating the effect of disturbances, sea level rise, and water management decisions on coastal wetland systems.     

Long‐term streamflow trends in Hawai'i and implications for native stream fauna

Climate change has fundamentally altered the water cycle in tropical islands, which is a critical driver of freshwater ecosystems. To examine how changes in streamflow regime have impacted habitat quality for native migratory aquatic species, we present a 50‐year (1967–2016) analysis of hydrologic records in 23 unregulated streams across the five largest Hawaiian Islands. For each stream, flow was separated into direct run‐off and baseflow and high‐ and low‐flow statistics (i.e., Q10 and Q90) with ecologically important hydrologic indices (e.g., frequency of flooding and low flow duration) derived. Using Mann–Kendall tests with a running trend analysis, we determined the persistence of streamflow trends through time. We analysed native stream fauna from ~400 sites, sampled from 1992 to 2007, to assess species richness among islands and streams. Declines in streamflow metrics indicated a general drying across the islands. In particular, significant declines in low flow conditions (baseflows), were experienced in 57% of streams, compared with a significant decline in storm flow conditions for 22% of streams. The running trend analysis indicated that many of the significant downward trends were not persistent through time but were only significant if recent decades (1987–2016) were included, with an average decline in baseflow and run‐off of 10.90% and 8.28% per decade, respectively. Streams that supported higher native species diversity were associated with moderate discharge and baseflow index, short duration of low flows, and negligible downward trends in flow. A significant decline in dry season flows (May–October) has led to an increase in the number of no‐flow days in drier areas, indicating that more streams may become intermittent, which has important implications for mauka to makai (mountain to ocean) hydrological connectivity and management of Hawai'i's native migratory freshwater fauna.