Emplacement of the 1907 Mauna Loa basalt flow as derived from precision topography and satellite imaging

An eruption in January of 1907, from the southwest rift zone of Mauna Loa, produced a substantial lava flow field. Satellite images and Differential Global Positioning System (DGPS) survey data, along with observations and photographs from the field, are combined to provide a new perspective on the 1907 eruption. Boundaries of the flow field from the satellite data, combined with field measurements of flow thickness, indicate an area of 25.1 km² and a volume of 86.6 million m³. The eastern lobe of the flow field covers an area of 13.1 km², with a volume of 55.0 million m³, and was emplaced with an average effusion rate of 119 m³/s (at least, for the upper portion of the lobe). Ten DGPS topographic profiles across the eastern lobe aid in distinguishing the characteristics of, and transitions between, the zones identified during the emplacement of the 1984 Mauna Loa flow. Several subdivisions have been built directly on top of or adjacent to the 1907 lava flow. The strong likelihood of future eruptions from the Mauna Loa southwest rift zone makes these housing developments of particular importance for assessments of potential volcanic hazards.     

DyeLIF™: A New Direct‐Push Laser‐Induced Fluorescence Sensor System for Chlorinated Solvent DNAPL and Other Non‐Naturally Fluorescing NAPLs

DyeLIF? is a new version of laser‐induced fluorescence (LIF) for high‐resolution three‐dimensional subsurface mapping of nonaqueous phase liquids (NAPLs) in the subsurface. DyeLIF eliminates the requirement that the NAPL contains native fluorophores (such as those that occur in compounds like polynuclear aromatic hydrocarbons [PAHs]) and can therefore be used to detect chlorinated solvents and other nonfluorescing NAPLs that had previously been undetectable with conventional LIF tools. With DyeLIF, an aqueous solution of water and nontoxic hydrophobic dye is continuously injected ahead of the sapphire detection window while the LIF probe is being advanced in the subsurface. If soil containing NAPL is penetrated, the injected dye solvates into the NAPL within a few milliseconds, creating strong fluorescence that is transmitted via fiber‐optic filaments to aboveground optical sensors. A detailed field evaluation of the novel DyeLIF technology was performed at a contaminated industrial site in Lowell, Massachusetts, USA where chlorinated solvent dense nonaqueous phase liquid (DNAPL) persists below the water table in sandy sediments. Continuously cored boreholes were drilled adjacent to 5 of 30 DyeLIF probes that were advanced at that site. The cores were subsampled in high resolution to generate discrete‐depth soil samples as splits at the same depths where DNAPL was detected in the colocated DyeLIF probes. The cores were analyzed above ground using (1) colorimetric screening using hydrophobic dye tests, (2) laboratory extraction and quantitative chemical analysis, (3) “Benchtop” DyeLIF, and (4) volumetric moisture content. Correlation between DyeLIF and aboveground analyses of the soil cores was excellent: 98% agreement with positive DNAPL detections in samples where DNAPL pore saturations were >0.7% (based on quantitative soil analyses) and the ex situ tests. DyeLIF produced the equivalent of one aboveground colorimetric dye test every 0.2 inch (0.5 cm) of probing. With average daily probing of 395 linear feet (120.4 m), this was the equivalent of 12,039 discrete‐depth colorimetric dye tests/day. Because DyeLIF is an in situ measurement, there are no issues with soil core recovery like there would be for conventional ex situ colorimetric dye tests and 100% characterization of the probed intervals is achieved. Tracking the injection rate and pressure of the dye solution provides simultaneous data regarding relative soil permeability, similar to other direct push (DP) hydraulic profiling tools. Conventional LIF is considered the premier DP tool to identify and map NAPL containing PAHs in the subsurface or confirm its absence. While chlorinated solvent DNAPLs at some field sites contain impurities (e.g., solvated greases or oils) that make them detectable with conventional LIF techniques, at other sites, the DNAPL cannot be detected with conventional LIF. At such sites, the injection of a hydrophobic dye ahead of the sapphire window with the DyeLIF system now makes the LIF technology applicable to the many types of NAPLs that were previously invisible using conventional LIF techniques.     

Natural Stimuli Calibration with Fining Direction Regularization in an Integrated Hydrologic Model

The interaction between surface water and groundwater during flood events is a complex process that has traditionally been described using simplified analytical solutions, or abstracted numerical models. To make the problem tractable, it is common to idealize the flood event, simplify river channel geometry, and ignore bank soil heterogeneity, often resulting in a model that only loosely represents the site, thus limiting its applicability to any specific river cross-section. In this study, we calibrate a site-specific fully-integrated surface and subsurface HydroGeoSphere model using flood events for a cross-section along the South River near Waynesboro, VA. The calibration approach presented in this study demonstrates the incorporation of fining direction regularization with a highly parameterized inversion driven by natural stimuli, to develop several realistic realizations of hydraulic conductivity fields that reflect the depositional history of the system. Specifically, we calibrate a model with 365 unique material zones to multiple flood events recorded in a dense well network while incorporating possible fining sequences consistent with the depositional history of the riverbank. Over 25,000 individual simulations were completed using calibration software and a cloud platform specifically designed for highly parallelized computing environments. The results of this study demonstrate the use of fining direction regularization during model calibration to generate multiple calibrated model realizations that account for the depositional environment of the system.     

Mid Campanian‐Lower Maastrichtian magnetostratigraphy of the James Ross Basin,Antarctica: Chronostratigraphical implications

The James Ross Basin, in the northern Antarctic Peninsula, exposes which is probably the world thickest and most complete Late Cretaceous sedimentary succession of southern high latitudes. Despite its very good exposures and varied and abundant fossil fauna, precise chronological determination of its infill is still lacking. We report results from a magnetostratigraphic study on shelfal sedimentary rocks of the Marambio Group, southeastern James Ross Basin, Antarctica. The succession studied covers a ~1,200 m‐thick stratigraphic interval within the Hamilton Point, Sanctuary Cliffs and Karlsen Cliffs Members of the Snow Hill Island Formation, the Haslum Crag Formation, and the lower López de Bertodano Formation. The basic chronological reference framework is given by ammonite assemblages, which indicate a Late Campanian – Early Maastrichtian age for the studied units. Magnetostratigraphic samples were obtained from five partial sections located on James Ross and Snow Hill islands, the results from which agree partially with this previous biostratigraphical framework. Seven geomagnetic polarity reversals are identified in this work, allowing to identify the Chron C32/C33 boundary in Ammonite Assemblage 8‐1, confirming the Late Campanian age of the Hamilton Point Member. However, the identification of the Chron C32/C31 boundary in Ammonite Assemblage 8‐2 assigns the base of the Sanctuary Cliffs Member to the early Maastrichtian, which differs from the Late Campanian age previously assigned by ammonite biostratigraphy. This magnetostratigraphy spans ~14 Ma of sedimentary succession and together with previous partial magnetostratigraphies on Early‐Mid Campanian and Middle Maastrichtian to Danian columns permits a complete and continuous record of the Late Cretaceous distal deposits of the James Ross Basin. This provides the required chronological resolution to solve the intra‐basin and global correlation problems of the Late Cretaceous in the Southern Hemisphere in general and in the Weddellian province in particular, given by endemism and diachronic extinctions on invertebrate fossils, including ammonites. The new chronostratigraphic scheme allowed us to calculate sediment accumulation rates for almost the entire Late Cretaceous infill of the distal James Ross Basin (the Marambio Group), showing a monotonous accumulation for more than 8 Myr during the upper Campanian and a dramatic increase during the early Maastrichtian, controlled by tectonic and/or eustatic causes.     

The rendition of the Atlantic Warm Pool in the reanalyses

The Atlantic Warm Pool (AWP) region, which is comprised of the Gulf of Mexico, Caribbean Sea and parts of the northwestern tropical Atlantic Ocean, is one of the most poorly observed parts of the global oceans. This study compares three ocean reanalyses, namely the Global Ocean Data Assimilation System of National Centers for Environmental Prediction (NCEP), the Climate Forecast System Reanalysis (CFSR) of NCEP, and the Simple Ocean Data Assimilation (SODA) for its AWP variation. The surface temperature in these ocean reanalyses is also compared with that from the Extended Range SST version 3 and Optimally Interpolated SST version 2 SST analyses. In addition we also compare three atmospheric reanalyses: NCEP-NCAR (R1), NCEP-DOE (R2), and CFSR for the associated atmospheric variability with the AWP. The comparison shows that there are important differences in the climatology of the AWP and its interannual variations. There are considerable differences in the subsurface ocean manifestation of the AWP with SODA (CFSR) showing the least (largest) modulation of the subsurface ocean temperatures. The remote teleconnections with the tropical Indian Ocean are also different across the reanalyses. However, all three oceanic reanalyses consistently show the absence of any teleconnection with the eastern equatorial Pacific Ocean. The influence of the AWP on the tropospheric temperature anomalies last for up to a one season lead and it is found to be relatively weak in R1 reanalyses. A simplified SST anomaly equation initially derived for diagnosing El Niño Southern Oscillation variability is adapted for the AWP variations in this study. The analysis of this equation reveals that the main contribution of the SST variation in the AWP region is from the variability of the net heat flux. All three reanalyses consistently show that the role of the ocean advective terms, including that associated with upwelling in the AWP region, is comparatively much smaller. The covariance of the SST tendency in the AWP with the net heat flux is large, with significant contributions from the variations of the surface shortwave and longwave fluxes.     

An empirical formula for the height of the coastal internal boundary layer

Observations of the height of the daytime coastal internal boundary layer at several sites are used to justify an empirical formula in the Offshore and Coastal Dispersion (OCD) model, which states that the boundary-layer slope is 0.1 in the first 2km from the shoreline, and 0.03 therafter.     

MultiscaleDTM: An open-source R package for multiscale geomorphometric analysis

Digital terrain models (DTMs) are datasets containing altitude values above or below a reference level, such as a reference ellipsoid or a tidal datum over geographic space, often in the form of a regularly gridded raster. They can be used to calculate terrain attributes that describe the shape and characteristics of topographic surfaces. Calculating these terrain attributes often requires multiple software packages that can be expensive and specialized. We have created a free, open-source R package, MultiscaleDTM , that allows for the calculation of members from each of the five major thematic groups of terrain attributes: slope, aspect, curvature, relative position, and roughness, from a regularly gridded DTM. Furthermore, these attributes can be calculated at multiple spatial scales of analysis, a key feature that is missing from many other packages. Here, we demonstrate the functionality of the package and provide a simulation exploring the relationship between slope and roughness. When roughness measures do not account for slope, these attributes exhibit a strong positive correlation. To minimize this correlation, we propose a new roughness measure called adjusted standard deviation. In most scenarios tested, this measure produced the lowest rank correlation with slope out of all the roughness measures tested. Lastly, the simulation shows that some existing roughness measures from the literature that are supposed to be independent of slope can actually exhibit a strong inverse relationship with the slope in some cases.     

Distributed hydrological modelling using weather radar in gauged and ungauged basins

Distributed hydrological modelling using space–time estimates of rainfall from weather radar provides a natural approach to area-wide flood forecasting and warning at any location, whether gauged or ungauged. However, radar estimates of rainfall may lack consistent, quantitative accuracy. Also, the formulation of hydrological models in distributed form may be problematic due to process complexity and scaling issues. Here, the aim is to first explore ways of improving radar rainfall accuracy through combination with raingauge network data via integrated multiquadric methods. When the resulting gridded rainfall estimates are employed as input to hydrological models, the simulated river flows show marked improvements when compared to using radar data alone. Secondly, simple forms of physical–conceptual distributed hydrological model are considered, capable of exploiting spatial datasets on topography and, where necessary, land-cover, soil and geology properties. The simplest Grid-to-Grid model uses only digital terrain data to delineate flow pathways and to control runoff production, the latter by invoking a probability-distributed relation linking terrain slope to soil absorption capacity. Model performance is assessed over nested river basins in northwest England, employing a lumped model as a reference. When the distributed model is used with the gridded radar-based rainfall estimators, it shows particular benefits for forecasting at ungauged locations.     

A sediment record of trace metal loadings in the Upper Mississippi River

Sediment cores from Lake Pepin, a natural lake on the Upper Mississippi River, reveal the historical trends in trace metal use and discharge in the watershed. Lead-210 dated concentration profiles of trace metals (Ag, Cd, Cr, Cu, Hg, Pb, V, Zn) in sediment cores from throughout the lake generally showed low and stable concentrations prior to settlement (circa 1830), peak concentrations between 1940 and 1975, and substantial decreases thereafter. Whole-lake sediment accumulation rates increased greatly over the period of record, from 79,000 metric tons year⁻¹ prior to 1830, to 876,000 metric tons year⁻¹ during the 1990s. Whole-lake accumulation rates of most trace metals peaked in the 1960s but decreased sharply after that. Sediment and trace metal accumulation rates decreased in the downstream direction, and approximately two-thirds of the sediment and trace metal mass accumulated in Lake Pepin since 1800 was deposited in the upper 30% (by area) of the lake. The dramatic declines in trace metal concentrations and accumulation rates in Lake Pepin sediments since 1970 coincide with increased pollution control and prevention efforts throughout the watershed, including the implementation of secondary treatment at a large municipal wastewater treatment plant upstream. This is one of eight papers dedicated to the “Recent Environmental History of the Upper Mississippi River” published in this special issue of the Journal of Paleolimnology. D.R. Engstrom served as guest editor of the special issue.     

Influence of rock mass strength on the erosion rate of alpine cliffs

Collapse of cliff faces by rockfall is a primary mode of bedrock erosion in alpine environments and exerts a first‐order control on the morphologic development of these landscapes. In this work we investigate the influence of rock mass strength on the retreat rate of alpine cliffs. To quantify rockwall competence we employed the Slope Mass Rating (SMR) geomechanical strength index, a metric that combines numerous factors contributing to the strength of a rock mass. The magnitude of cliff retreat was calculated by estimating the volume of talus at the toe of each rockwall and projecting that material back on to the cliff face, while accounting for the loss of production area as talus buries the base of the wall. Selecting sites within basins swept clean by advancing Last Glacial Maximum (LGM) glaciers allowed us to estimate the time period over which talus accumulation occurred (i.e. the production time). Dividing the magnitude of normal cliff retreat by the production time, we calculated recession rates for each site. Our study area included a portion of the Sierra Nevada between Yosemite National Park and Lake Tahoe. Rockwall recession rates determined for 40 alpine cliffs in this region range from 0·02 to 1·22 mm/year, with an average value of 0·28 mm/year. We found good correlation between rockwall recession rate and SMR which is best characterized by an exponential decrease in erosion rate with increasing rock mass strength. Analysis of the individual components of the SMR reveals that joint orientation (with respect to the cliff face) is the most important parameter affecting the rockwall erosion rate. The complete SMR score, however, best synthesizes the lithologic variables that contribute to the strength and erodibility of these rock slopes. Our data reveal no strong independent correlations between rockwall retreat rate and topographic attributes such as elevation, aspect, or slope angle. Copyright © 2009 John Wiley & Sons, Ltd.