Geomorphic and sedimentologic evidence for the separation of Lake Superior from Lake Michigan and Huron

A common break was recognized in four Lake Superior strandplain sequences using geomorphic and sedimentologic characteristics. Strandplains were divided into lakeward and landward sets of beach ridges using aerial photographs and topographic surveys to identify similar surficial features and core data to identify similar subsurface features. Cross-strandplain, elevation-trend changes from a lowering towards the lake in the landward set of beach ridges to a rise or reduction of slope towards the lake in the lakeward set of beach ridges indicates that the break is associated with an outlet change for Lake Superior. Correlation of this break between study sites and age model results for the strandplain sequences suggest that the outlet change occurred sometime after about 2,400 calendar years ago (after the Algoma phase). Age model results from one site (Grand Traverse Bay) suggest an alternate age closer to about 1,200 calendar years ago but age models need to be investigated further. The landward part of the strandplain was deposited when water levels were common in all three upper Great Lakes basins (Superior, Huron, and Michigan) and drained through the Port Huron/Sarnia outlet. The lakeward part was deposited after the Sault outlet started to help regulate water levels in the Lake Superior basin. The landward beach ridges are commonly better defined and continuous across the embayments, more numerous, larger in relief, wider, have greater vegetation density, and intervening swales contain more standing water and peat than the lakeward set. Changes in drainage patterns, foreshore sediment thickness and grain size help in identifying the break between sets in the strandplain sequences. Investigation of these breaks may help identify possible gaps in the record or missing ridges in strandplain sequences that may not be apparent when viewing age distributions and may justify the need for multiple age and glacial isostatic adjustment models. This is the third in a series of ten papers published in this special issue of Journal of Paleolimnology. These papers were presented at the 47th Annual Meeting of the International Association for Great Lakes Research (2004), held at the University of Waterloo, Waterloo, Ontario, Canada. P.F. Karrow and C.F.M. Lewis were guest editors of this special issue. The U.S. Government's right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Inverse methods for estimating primary input signals from time-averaged isotope profiles

Mammalian teeth are invaluable archives of ancient seasonality because they record along their growth axes an isotopic record of temporal change in environment, plant diet, and animal behavior. A major problem with the intra-tooth method is that intra-tooth isotope profiles can be extremely time-averaged compared to the actual pattern of isotopic variation experienced by the animal during tooth formation. This time-averaging is a result of the temporal and spatial characteristics of amelogenesis (tooth enamel formation), and also results from laboratory sampling. This paper develops and evaluates an inverse method for reconstructing original input signals from time-averaged intra-tooth isotope profiles. The method requires that the temporal and spatial patterns of amelogenesis are known for the specific tooth and uses a minimum length solution of the linear system Am = d, where d is the measured isotopic profile, A is a matrix describing temporal and spatial averaging during amelogenesis and sampling, and m is the input vector that is sought. Accuracy is dependent on several factors, including the total measurement error and the isotopic structure of the measured profile. The method is shown to accurately reconstruct known input signals for synthetic tooth enamel profiles and the known input signal for a rabbit that underwent controlled dietary changes. Application to carbon isotope profiles of modern hippopotamus canines reveals detailed dietary histories that are not apparent from the measured data alone. Inverse methods show promise as an effective means of dealing with the time-averaging problem in studies of intra-tooth isotopic variation.     

Assessment of critical path analyses of the relationship between permeability and electrical conductivity of pore networks

Critical path analysis (CPA) is a method for estimating macroscopic transport coefficients of heterogeneous materials that are highly disordered at the micro-scale. Developed originally to model conduction in semiconductors, numerous researchers have noted that CPA might also have relevance to flow and transport processes in porous media. However, the results of several numerical investigations of critical path analysis on pore network models raise questions about the applicability of CPA to porous media. Among other things, these studies found that (i) in well-connected 3D networks, CPA predictions were inaccurate and became worse when heterogeneity was increased; and (ii) CPA could not fully explain the transport properties of 2D networks. To better understand the applicability of CPA to porous media, we made numerical computations of permeability and electrical conductivity on 2D and 3D networks with differing pore-size distributions and geometries. A new CPA model for the relationship between the permeability and electrical conductivity was found to be in good agreement with numerical data, and to be a significant improvement over a classical CPA model. In sufficiently disordered 3D networks, the new CPA prediction was within ±20% of the true value, and was nearly optimal in terms of minimizing the squared prediction errors across differing network configurations. The agreement of CPA predictions with 2D network computations was similarly good, although 2D networks are in general not well-suited for evaluating CPA. Numerical transport coefficients derived for regular 3D networks of slit-shaped pores were found to be in better agreement with experimental data from rock samples than were coefficients derived for networks of cylindrical pores.     

Hypoxia in Manila Bay, Philippines during the northeast monsoon

Herein we present results from one of the first extensive bay-wide oceanographic surveys of Manila Bay, wherein 31 stations were sampled during the northeast monsoon (cold and dry season). A band of hypoxic bottom water (dissolved oxygen<2.8 mg/L) spanned the midsection of the bay from east to west. Bottom nitrate concentrations (5.7-16.8 μM; avg. 11.1 μM) and total organic carbon values in sediments (1.7-3.1%; avg. 2.4%) were high in the midsection, which coincided with the band of hypoxic bottom water. Physical processes and site-specific accumulation of organic material likely lead to hypoxic conditions in Manila Bay, even during the northeast monsoon period when the water column is relatively well mixed. The results of this study complement the previously reported widespread hypoxia that occurs during the rainy season. Thus, hypoxia may be pervasive in the bay throughout the year, although it varies in intensity and spatial extent.     

Hydromechanical behavior during constant-rate pumping tests in fractured gneiss

Deformation and fluid pressure during constant-rate pumping tests were characterized by sensitivity analyses using a theoretical model of a deformable fracture, and by conducting and analyzing field tests at a site underlain by fractured biotite gneiss in western South Carolina, USA. The sensitivity analysis indicates that displacement is a hysteretic function of fluid pressure during constant-rate pumping tests, and the signals are affected by properties of the fracture (normal stiffness, aperture, and heterogeneities in the fracture plane) and enveloping rock (elastic modulus and permeability). The field tests used a removable borehole extensometer to measure axial displacements in the pumping well. The field results are generally similar to simulations, and the hydraulic and mechanical properties obtained by inversion of the pumping test results are essentially the same as results from hydromechanical slug tests. The observed displacements early in the pumping tests are less than, whereas late in the test they are greater than predicted values assuming a uniform fracture. This difference can be explained as a consequence of preferential flow within the fracture. These results indicate that displacements in a wellbore are sensitive to details of fracture connectivity and preferential flow that are difficult to detect through the pressure signal alone.     

The remarkable occurrence of large rainfall-induced debris flows at two different locations on July 12, 2008, Southern Sierra Nevada, CA, USA

On July 12, 2008, two convective cells about 155 km apart produced a brief period of intense rainfall triggering large debris flows in the southern Sierra Nevada. The northernmost cell was centered over Oak Creek Canyon, an east-flowing drainage, and its tributaries near Independence, CA, USA. About 5:00 p.m., debris flows passed down the South Fork and North Fork of Oak Creek to merge into a large single feature whose passage affected the historic Mt. Whitney Fish hatchery and blocked California State Highway 395. At about the same time, the southernmost cell was largely centered over Erskine Creek, a main tributary of the west-flowing Kern River. Debris flows issued from several branches to coalesce into a large debris flow that passed along Erskine Creek, through the town of Lake Isabella, CA, USA and into the Kern River. It was observed reaching Lake Isabella about 6:30 p.m. Both debris flows caused significant disruption and damage to local communities.     

A coronal magnetic field model with horizontal volume and sheet currents

When globally mapping the observed photospheric magnetic field into the corona, the interaction of the solar wind and magnetic field has been treated either by imposing source surface boundary conditions that tacitly require volume currents outside the source surface (Schatten, Wilcox, and Ness, 1969) or by limiting the interaction to thin current sheets between oppositely directed field regions (Wolfson, 1985). Yet observations and numerical MHD calculations suggest the presence of non-force-free volume currents throughout the corona as well as thin current sheets in the neighborhoods of the interfaces between closed and open field lines or between oppositely directed open field lines surrounding coronal helmet-streamer structures. This work presents a model including both horizontal volume currents and streamer sheet currents. The present model builds on the magnetostatic equilibria developed by Bogdan and Low (1986) and the current-sheet modeling technique developed by Schatten (1971). The calculation uses synoptic charts of the line-of-sight component of the photospheric magnetic field measured at the Wilcox Solar Observatory. Comparison of an MHD model with the calculated model results for the case of a dipole field and comparison of eclipse observations with calculations for CR 1647 (near solar minimum) show that this horizontal current-current-sheet model reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces coronal helmet structures better than the current-sheet model.     

Eclipsing Binary Stars in Nearby Galaxies

We briefly discuss the survey programme we are conducting to detect eclipsing binaries in local group galaxies. Some lightcurves from studies of M31, IC 1613 and NGC 6822 are presented along with details of future work.     

Studies of multiple stellar systems – IV. The triple-lined spectroscopic system Gliese 644

We present a radial velocity study of the triple-lined system Gliese 644 and derive spectroscopic elements for the inner and outer orbits with periods of 2.965 5 and 627 d. We also utilize old visual data, as well as modern speckle and adaptive optics observations, to derive a new astrometric solution for the outer orbit. These two orbits together allow us to derive masses for each of the three components in the system: M _A=0.410±0.028 (6.9 per cent), M _Ba=0.336±0.016 (4.7 per cent), and M _Bb=0.304±0.014 (4.7 per cent) M_⊙. We suggest that the relative inclination of the two orbits is very small. Our individual masses and spectroscopic light ratios for the three M stars in the Gliese 644 system provide three points for the mass–luminosity relation near the bottom of the main sequence, where the relation is poorly determined. These three points agree well with theoretical models for solar metallicity and an age of 5 Gyr. Our radial velocities for Gliese 643 and vB 8, two common proper motion companions of Gliese 644, support the interpretation that all five M stars are moving together in a physically bound group. We discuss possible scenarios for the formation and evolution of this configuration, such as the formation of all five stars in a sequence of fragmentation events leading directly to the hierarchical configuration now observed, versus formation in a small N cluster with subsequent dynamical evolution into the present hierarchical configuration.