Micromorphological Analysis of Coastal Sediments from Willapa Bay, Washington, USA: A Technique for Analysing Inferred Tsunami Deposits

Tsunami deposits are provisionally distinguished in the field on the basis of anomalous sand horizons, fining-up and fining-landward, coupled with organic-rich, fragmented `backwash' sediments. In this paper, micromorphological features of a sediment sequence previously interpreted as being of tsunami origin are described. These characteristics are shown to be consistent with the macro-scale features used elsewhere, but show additional details not seen in standard stratigraphies, including possible evidence for individual waves, possibly wave-magnitude progression, organic fragment alignment and intraclast microstructures. Although replication and more complete studies are needed, this analysis confirms the identification of a tsunami in Willapa Bay in ca.1700 AD, while demonstrating a widely applicable technique for confirming or refuting possible tsunami deposits.     

Lithium and lithium isotope profiles through the upper oceanic crust: a study of seawater-basalt exchange at ODP Sites 504B and 896A

Ocean Drilling Program (ODP) Hole 504B near the Costa Rica Rift is the deepest hole drilled in the ocean crust, penetrating a volcanic section, a transition zone and a sheeted dike complex. The distribution of Li and its isotopes through this 1.8-km section of oceanic crust reflects the varying conditions of seawater alteration with depth. The upper volcanic rocks, altered at low temperatures, are enriched in Li (5.6-27.3 ppm) and have heavier isotopic compositions (δ⁷Li=6.6-20.8‰) relative to fresh mid-ocean ridge basalt (MORB) due to uptake of seawater Li into alteration clays. The Li content and isotopic compositions of the deeper volcanic rocks are similar to MORB, reflecting restricted seawater circulation in this section. The transition zone is a region of mixing of seawater with upwelling hydrothermal fluids and sulfide mineralization. Li enrichment in this zone is accompanied by relatively light isotopic compositions (−0.8-2.1‰) which signify influence of basalt-derived Li during mineralization and alteration. Li decreases with depth to 0.6 ppm in the sheeted dike complex as a result of increasing hydrothermal extraction in the high-temperature reaction zone. Rocks in the dike complex have variable isotopic values that range from −1.7 to 7.9‰, depending on the extent of hydrothermal recrystallization and off-axis low-temperature alteration. Hydrothermally altered rocks are isotopically light because ⁶Li is preferentially retained in greenschist and amphibolite facies minerals. The δ⁷Li values of the highly altered rocks of the dike complex are complementary to those of high-temperature mid-ocean ridge vent fluids and compatible to equilibrium control by the alteration mineral assemblage. The inventory of Li in basement rocks permits a reevaluation of the role of oceanic crust in the budget of Li in the ocean. On balance, the upper 1.8 km of oceanic crusts remains a sink for oceanic Li. The observations at 504B and an estimated flux from the underlying 0.5 km of gabbro suggest that the global hydrothermal flux is at most 8×10⁹ mol/yr, compatible with geophysical thermal models. This work defines the distribution of Li and its isotopes in the upper ocean crust and provides a basis to interpret the contribution of subducted lithosphere to arc magmas and cycling of crustal material in the deep mantle.     

Isotope variations in a Sierra Nevada snowpack and their relation to meltwater

Isotopic variations in melting snow are poorly understood. We made weekly measurements at the Central Sierra Snow Laboratory, California, of snow temperature, density, water equivalent and liquid water volume to examine how physical changes within the snowpack govern meltwater δ¹⁸O. Snowpack samples were extracted at 0.1 m intervals from ground level to the top of the snowpack profile between December 1991 and April 1992. Approximately 800 mm of precipitation fell during the study period with δ¹⁸O values between −21.35 and −4.25‰. Corresponding snowpack δ¹⁸O ranged from −22.25 to −6.25‰. The coefficient of variation of δ¹⁸O in snowpack levels decreased from −0.37 to −0.07 from winter to spring, indicating isotopic snowpack homogenization. Meltwater δ¹⁸O ranged from −15.30 to −8.05‰, with variations of up to 2.95‰ observed within a single snowmelt episode, highlighting the need for frequent sampling. Early snowmelt originated in the lower snowpack with higher δ¹⁸O through ground heat flux and rainfall. After the snowpack became isothermal, infiltrating snowmelt displaced the higher δ¹⁸O liquid in the lower snowpack through a piston flow process. Fractionation analysis using a two-component mixing model on the isothermal snowpack indicated that δ¹⁸O in the initial and final half of major snowmelt was 1.30‰ lower and 1.45‰ higher, respectively, than the value from simple mixing. Mean snowpack δ¹⁸O on individual profiling days showed a steady increase from −15.15 to −12.05‰ due to removal of lower δ¹⁸O snowmelt and addition of higher δ¹⁸O rainfall. Results suggest that direct sampling of snowmelt and snow cores should be undertaken to quantify tracer input compositions adequately. The snowmelt sequence also suggests that regimes of early lower δ¹⁸O and later higher δ¹⁸O melt may be modeled and used in catchment tracing studies.     

Laboratory Study of Air Sparging of TCE-Contaminated Saturated Soils and Ground Water

Air sparging has proven to be an effective remediation technique for treating saturated soils and ground water contaminated by volatile organic compounds (VOCs). Since little is known about the system variables and mass transfer mechanisms important to air sparging, several researchers have recently performed laboratory investigations to study such issues. This paper presents the results of column experiments performed to investigate the behavior of dense nonaqueous phase liquids (DNAPFs). specifically trichloroethylene (TCE), during air sparging. The specific objectives of the study were (1) to compare the removal of dissolved TCE with the removal of dissolved light nonaqueous phase liquids (LNAPLs). such as benzene or toluene; (2) to determine the effect of injected air-flow rate on dissolved TCE removal; (3) to determine the effect of initial dissolved TCE concentration on removal efficiency; and (4) to determine the differences in removal between dissolved and pure-phase TCE. The test results showed that (1) the removal of dissolved TCE was similar to that of dissolved LNAPL: (2) increased air-injection rates led to increased TCE removal at lower ranges of air injection, but further increases at higher ranges of air injection did not increase the rate of removal, indicating a threshold removal rate had been reached; (3) increased initial concentration of dissolved TCE resulted in similar rates of removal: and (4) the removal of pure-phase TCE was difficult using a low air-injection rate, but higher air-injection rates led to easier removal.     

How do streamflow generation mechanisms affect watershed hypsometry?

Several studies have shown that the dominant streamflow generation mechanism in a river basin can leave distinct geomorphological signatures in basin topography. In particular, it has been suggested previously that basins generated by groundwater discharge tend to have a larger hypsometric integral than surface runoff basins because fluvial erosion is more focused in the valleys where groundwater discharge tends to occur. In this analysis, we aim to clarify this relationship by developing an alternative method to quantify the effects of streamflow generation mechanisms on basin hypsometry and by using a numerical model that can generate streamflow by different processes to evaluate the sensitivity of the results to the hydrological and geomorphological properties of the basin. The model results suggest that the hypsometric characteristics that are usually associated with groundwater discharge basins, such as a larger hypsometric integral, occur primarily when drainage networks are still advancing in the watershed. During later stages of development, an additional factor such as lithological controls or a distinct geomorphological process would be needed to preserve these features. The model results also show that the hypsometric effects are stronger when the parameters of the fluvial erosion process promote the influence of small discharge rates. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of flooded rice paddy on remotely sensed evapotranspiration in the Krishna River basin,India

Evapotranspiration (ET) is one of the major water exchange processes between the earth's surface and the atmosphere. ET is a combined process of evaporation from open water bodies, bare soil and plant surfaces, and transpiration from vegetation. Remote sensing-based ET models have been developed to estimate spatially distributed ET over large regions, however, many of them reportedly underestimate ET over semi-arid regions (Jamshidi et al., Journal of Hydrometeorology, 2019, 20, 947–964). In this work, we show that underestimation of ET can occur due to the open water evaporation from flooded rice paddies ignored in the existing ET models. To address the gap in ET estimation, we have developed a novel approach that accounts for the missing ET component over flooded rice paddies. Our method improved ET estimates by a modified Penman-Monteith algorithm that considered the fraction of open water evaporation from flooded rice paddies. Daily ET was calculated using ground based meteorological data and the MODIS satellite data over the Krishna River Basin. Seasonal and annual ET values over the Krishna Basin were compared with two different ET algorithms. ET estimates from these two models were also compared for different crop combinations. Results were validated with flux tower-based measurements from other studies. We have identified a 17 mm/year difference in average annual ET over the Krishna River Basin with this new ET algorithm. This is very critical in basin scale water balance analysis and water productivity studies.     

Water mining from the deep critical zone by apple trees growing on loess

There have been significant recent advances in understanding the ecohydrology of deep soil. However, the links between root development and water usage in the deep critical zone remains poorly understood. To clarify the interaction between water use and root development in deep soil, we investigated soil water and root profiles beyond maximum rooting depth in five apple orchards planted on farmland with stand ages of 8, 11, 15, 18, and 22 years in a subhumid region on the Chinese Loess Plateau. Apple trees rooted progressively deeper for water with increasing stand age and reached 23.2 ± 0.8 m for the 22‐year‐old trees. Soil water deficit in deep soil increased with tree age and was 1,530 ± 43 mm for a stand age of 22 years. Measured root deepening rate was far great than the reported pore water velocity, which demonstrated that trees are mining resident old water. The deficits are not replenished during the life‐span of the orchard, showing a one‐way mining of the critical zone water. The one‐way root water mining may have changed the fine root profile from an exponential pattern in the 8‐year‐old orchard to a relative uniform distribution in older orchards. Our findings enhance our understanding of water‐root interaction in deep soil and reveal the unintended consequences of critical zone dewatering during the lifespan of apple trees.     

Use of long-period surface waves for rapid determination of earthquake-source parameters

A method for rapid retrieval of earthquake-source parameters from long-period surface waves is developed. With this method, the fault geometry and seismic moment can be determined immediately after the surface wave records have been retrieved. Hence, it may be utilized for warning of tsunamis in real time. The surface wave spectra are inverted to produce either a seismic moment tensor (linear) or a fault model (nonlinear). The method has been tested by using the IDA (International Deployment of Accelerographs) records. With these records the method works well for the events larger than M_s = 6, and is useful for investigating the nature of slow earthquakes.For events deeper than 30 km, all of the five moment tensor elements can be determined. For very shallow events (d ? 30 km) the inversion becomes ill-conditioned and two of the five source moment tensor elements become unresolvable. This difficulty is circumvented by a two-step inversion. In the first step, the unresolvable elements are constrained to be zero to yield a first approximation. In the second step, additional geological and geophysical data are incorporated to improve the first approximation. The effect of the source finiteness is also included.