Paleohydrology of arid southeastern maui, hawaiian islands, and its implications for prehistoric human settlement

Arid slopes on the southeastern side of Maui are densely covered with archaeological remains of Hawaiian settlement from the late prehistoric to early postcontact period (ca. A.D. 1500-1860). Permanent habitation sites, agricultural features, and religious structures indicate perennial occupation and farming in a subregion called Kahikinui, yet there is presently no year-round water source. We explore the possibility that postcontact deforestation led to the loss of either (1) perennial channel flow or (2) perennial springs or seeps. To investigate the first possibility, we estimated ancient peak flows on 11 ephemeral channels in Kahikinui using field measurements and paleohydrology. Peak-flow estimates (3-230 m³/s) for a given drainage area are smaller than those for current perennial Maui streams, but are equivalent to gauged peak flows from ephemeral and intermittent streams in the driest regions of Hawai’i and Maui islands. This is consistent with the long-term absence of perennial channel flow in Kahikinui. On the other hand, others have shown that canopy fog-drip in Hawai’i can be greater than rainfall and thus a large part of groundwater recharge. Using isolated live remnants and snags, we estimate the former extent of the forest upstream from archaeological sites. We use rough estimates of the loss of fog-drip recharge caused by deforestation and apply a simple, steady-state hydrologic model to calculate potential groundwater table fall. These order-of-magnitude estimates indicate that groundwater could have fallen by a minimum of several meters, abandoning perennial seeps. This is consistent with archaeological evidence for former perennial seeps, such as stonewalls enclosing potential seeps to protect them. Although longer-term reductions in rainfall cannot be ruled out as a factor, deforestation and loss of fog-drip recharge are obvious and more immediate reasons for a recent loss of perennial water in Kahikinui, Maui.     

The postglacial history of three Picea species in New England, USA

Given the difficulty of separating the three Picea species—P. glauca, P. mariana, and P. rubens (white, black, and red spruce)—in the pollen record, little is known about their unique histories in eastern North America following deglaciation. Here we report the first use of a classification tree analysis (CART) to distinguish pollen grains of these species. It was successfully applied to fossil pollen from eight sites in Maine and one in Massachusetts. We focused on the late glacial/early Holocene (14,000 to 8000 cal yr B.P.) and the late Holocene (1400 cal yr B.P. to present)—the two key periods since deglaciation when Picea has been abundant in the region. The result shows a shift from a Picea forest of P. glauca and P. mariana in the late glacial to a forest of P. rubens and P. mariana in the late Holocene. The small number of P. rubens grains identified from the late glacial/early Holocene samples (<5%) suggests that that species was either absent or rare at most of the sites. The occurrence and distribution of the three species do not reveal any geographic or temporal trend during late glacial time, but the data suggest that they were distributed in local patches on the landscape. The results of this study indicate that the recent population expansion of Picea (1000 to 500 cal yr B.P.) was likely the first time since deglaciation that P. rubens was abundant in the region.     

Channel Bed Mobility Downstream from the Elwha Dams,Washington*

Dams are a major source of fragmentation and degradation of rivers. Although substantial research has been conducted on the environmental impacts of large structures in the United States, smaller dams have received less attention. This study evaluated the impact of two dams of moderate size, the Elwha Dams, on the downstream channel system using field data collection at river cross‐sections. The relationship of average boundary shear stress (τ_o) to critical shear stress (τ_cr) served as the basis for determining channel bed material mobility under the two‐year and ten‐year flood events. The channel had the greatest channel bed mobility at the natural cross‐section upstream from the dams, low bed mobility between the structures, and an increase in channel bed mobility in the low gradient river segment near the mouth of the river. Low bed mobility tended to be associated with a lack of channel system complexity, including reduction or loss of bars and low alluvial terraces and their associated young riparian communities. Although these run‐of‐the‐river dams do not modify streamflow greatly, the loss of sediment from the channel system has had a substantial impact on bed mobility and geomorphic and biotic complexity of the Elwha River.     

Potential Salinity and Temperature Futures for the Chesapeake Bay Using a Statistical Downscaling Spatial Disaggregation Framework

Estuaries are productive and ecologically important ecosystems, incorporating environmental drivers from watersheds, rivers, and the coastal ocean. Climate change has potential to modify the physical properties of estuaries, with impacts on resident organisms. However, projections from general circulation models (GCMs) are generally too coarse to resolve important estuarine processes. Here, we statistically downscaled near-surface air temperature and precipitation projections to the scale of the Chesapeake Bay watershed and estuary. These variables were linked to Susquehanna River streamflow using a water balance model and finally to spatially resolved Chesapeake Bay surface temperature and salinity using statistical model trees. The low computational cost of this approach allowed rapid assessment of projected changes from four GCMs spanning a range of potential futures under a high CO₂ emission scenario, for four different downscaling methods. Choice of GCM contributed strongly to the spread in projections, but choice of downscaling method was also influential in the warmest models. Models projected a ~2–5.5 °C increase in surface water temperatures in the Chesapeake Bay by the end of the century. Projections of salinity were more uncertain and spatially complex. Models showing increases in winter-spring streamflow generated freshening in the Upper Bay and tributaries, while models with decreased streamflow produced salinity increases. Changes to the Chesapeake Bay environment have implications for fish and invertebrate habitats, as well as migration, spawning phenology, recruitment, and occurrence of pathogens. Our results underline a potentially expanded role of statistical downscaling to complement dynamical approaches in assessing climate change impacts in dynamically challenging estuaries.     

Assessing rockfall susceptibility in steep and overhanging slopes using three-dimensional analysis of failure mechanisms

Rockfalls strongly influence the evolution of steep rocky landscapes and represent a significant hazard in mountainous areas. Defining the most probable future rockfall source areas is of primary importance for both geomorphological investigations and hazard assessment. Thus, a need exists to understand which areas of a steep cliff are more likely to be affected by a rockfall. An important analytical gap exists between regional rockfall susceptibility studies and block-specific geomechanical calculations. Here we present methods for quantifying rockfall susceptibility at the cliff scale, which is suitable for sub-regional hazard assessment (hundreds to thousands of square meters). Our methods use three-dimensional point clouds acquired by terrestrial laser scanning to quantify the fracture patterns and compute failure mechanisms for planar, wedge, and toppling failures on vertical and overhanging rock walls. As a part of this work, we developed a rockfall susceptibility index for each type of failure mechanism according to the interaction between the discontinuities and the local cliff orientation. The susceptibility for slope parallel exfoliation-type failures, which are generally hard to identify, is partly captured by planar and toppling susceptibility indexes. We tested the methods for detecting the most susceptible rockfall source areas on two famously steep landscapes, Yosemite Valley (California, USA) and the Drus in the Mont-Blanc massif (France). Our rockfall susceptibility models show good correspondence with active rockfall sources. The methods offer new tools for investigating rockfall hazard and improving our understanding of rockfall processes.     

Microearthquake seismicity and active tectonics of northwestern Greece

We carried out a microearthquake survey lasting for six weeks in northwest Greece using 18 portable seismograph stations to examine a region in which normal and thrust faulting have been reported in close proximity to one another. With this array we located 148 events and determined fault plane solutions for eight events using only rays radiated upwards. The seismicity of the region is diffuse with events extending to depths of nearly 30 km, and there is a minimum in activity near a depth of 15 km. The fault plane solutions exhibit a wide spectrum of fault types and orientations and are not consistent with simple zones of shortening or extension. Neither tractions applied to the edges or bottom of the region nor deviatoric stresses that compensate for lateral variations in crustal thickness can account for the variety of fault plane solutions. We think that the complicated behavior is a manifestation of inhomogeneous deformation due, at least in part, to a pre-existing complicated juxtaposition of structures and formations.     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

The Effect of Vegetation Productivity on Millet Prices in the Informal Markets of Mali, Burkina Faso and Niger

Systematic evaluation of food security throughout the Sahel has been attempted for nearly two decades. Food security analyses have used both food prices to determine the ability of the population to access food, and satellite-derived vegetation indices that measure vegetation production to establish how much food is available each year. The relationship between these two food security indicators is explored here using correspondence analysis and through the use of Markov chain models. Two sources of quantitative data were used: 8 km normalized difference vegetation index (NDVI) data from the Advanced Very High Resolution Radiometers (AVHRR) carried on the NOAA series of satellites, and monthly millet prices from 445 markets in Mali, Niger and Burkina Faso. The results show that the growing season vegetation production is related to the price of millet at the annual and the seasonal time scales. If the growing season was characterized by erratic, sparse rainfall, it resulted in higher prices, and well-distributed, abundant rainfall resulted in lower prices. The correspondence between vegetation production and millet prices is used to produce maps of millet prices for West Africa.     

Oxidation of CO on surface hematite in high CO2 atmospheres

We propose a mechanism for the oxidation of gaseous CO into CO₂ occurring on the surface mineral hematite (Fe₂O₃(s)) in hot, CO₂-rich planetary atmospheres, such as Venus. This mechanism is likely to constitute an important source of tropospheric CO₂ on Venus and could at least partly address the CO₂ stability problem in Venus’ stratosphere, since our results suggest that atmospheric CO₂ is produced from CO oxidation via surface hematite at a rate of 0.4 petagrammes (Pg) CO₂ per (Earth) year on Venus which is about 45% of the mass loss of CO₂ via photolysis in the Venusian stratosphere. We also investigated CO oxidation via the hematite mechanism for a range of planetary scenarios and found that modern Earth and Mars are probably too cold for the mechanism to be important because the rate-limiting step, involving CO(g) reacting onto the hematite surface, proceeds much slower at lower temperatures. The mechanism may feature on extrasolar planets such as Gliese 581c or CoRoT-7b assuming they can maintain solid surface hematite which, e.g. starts to melt above about 1200 K. The mechanism may also be important for hot Hadean-type environments and for the emerging class of hot Super-Earths with planetary surface temperatures between about 600 and 900 K.