Development and Applications of Double-difference Seismic Tomography

Double-difference (DD) tomography is a generalization of DD location; it simultaneously solves for the three-dimensional velocity structure and seismic event locations. DD tomography uses a combination of absolute and more accurate differential arrival times and hierarchically determines the velocity structure from larger scale to smaller scale. This method is able to produce more accurate event locations and velocity structure near the source region than standard tomography, which uses only absolute arrival times. We conduct a stability and uncertainty analysis of DD tomography based on a synthetic data set. Currently three versions of the DD tomography algorithms exist: tomoDD, tomoFDD and tomoADD. TomoDD assumes a flat earth model and uses a pseudo-bending ray-tracing algorithm to find rays between events and stations while tomoFDD uses a finite-difference travel-time algorithm and the curvature of the Earth is considered. Both codes are based on a regularly distributed inversion grid, with the former for a local scale and the latter for a regional scale. In contrast, tomoADD adapts the inversion mesh to match with the data distribution based on tetrahedral and Voronoi diagrams. We discuss examples of applying DD tomography to characterize fault zone structure, image high-resolution structure of subduction zones, and determine the velocity structure of volcanoes.     

Planform dynamics of the Lower Mississippi River

This paper presents an analysis of the planform behaviour of the Lower Mississippi River (LMR) using a series of maps and hydrographic surveys covering the period 1765–1975. Data allow analysis at various time and space scales, using fixed and statistically defined reaches, both before and after extensive channel modification. Previous research has interpreted planform change in relation to geomorphological or engineering regime‐type analyses of channel length and width for the LMR as a ‘single system’. The analysis here is broadly consistent with these approaches, but highlights the importance of meander geometry, in the form of the radius of curvature:width ratio. This neglected factor helps resolve paradoxes relating to observed changes in sediment transport and channel stability. When viewed over smaller time and space scales, analysis of dynamics using fixed reach boundaries reveals a downstream trend in the pattern of planform behaviour, which is closely related to the distribution of valley floor deposits, and which also reflects neotectonic influences. Analysis of changes using statistically determined reach boundaries shows that, over shorter time scales, meander trains are continually formed and modified over a period of approximately 120 years. Zones of more‐or‐less dynamic behaviour thus move through the LMR. The research also provides a context for 20th century engineering interventions to the river. These have constrained the magnitude of planform adjustment, but also altered the kind of response that is now possible in relation to changes in discharge and sediment load, and as a consequence of internal feedbacks within the LMR system. Copyright © 2006 John Wiley & Sons, Ltd.     

Imaging the heterogeneous source area of the 2003 M6.4 northern Miyagi earthquake, NE Japan, by double-difference tomography

A shallow M6.4 inland earthquake occurred on 26 July 2003 in the northern part of Miyagi Prefecture, northeastern Japan. This earthquake was a typical inland thrust earthquake, a type that is common in NE Japan. We obtained a detailed seismic velocity structure in the focal area of this earthquake by the double-difference tomography method. Arrival-time data came from temporary seismic stations deployed above the mainshock fault plane. Both the P-wave and S-wave velocities in the hanging wall were lower than those in the footwall. Aftershocks were aligned along a zone where the seismic velocity changes rapidly. This is consistent with the interpretation that the 2003 northern Miyagi earthquake occurred along a fault that acted as a normal fault in the Miocene and has been reactivated as a reverse fault under the present compressional stress regime. The large slip area by the main shock rupture (asperity) corresponds to an area with relatively high P- and S-wave velocities. A zone with low V_p/V_s was detected along the aftershock area. One of the possible causes of this low-V_p/V_s zone is the existence of high-aspect-ratio pores that contain water. Hypocenters of the main shock, largest foreshock, and largest aftershock are also located within the low-V_p/V_s zone.     

A comparison of flow intensities in alluvial rivers: characteristics and implications for modelling flow processes

This paper compares flow intensity data obtained with different instruments from a variety of fluvial environments. It examines associations between the root-mean-square of longitudinal velocity fluctuations (flow intensity), local mean velocity, relative depth, and boundary resistance. Results indicate systematic differences in the behaviour of flow intensity which scale with respect to position in the boundary layer (deep sand-bedded rivers), boundary grain resistance (shallow river environments with coarse beds), and possibly form resistance (shallower sand-bedded rivers). Preliminary approaches to prediction and modelling of variations in flow intensity are suggested based upon linear regression relationships. Intensity values are also compared with theoretical and empirical limits to the use of Taylor's substitution, which allows time and frequency properties of a single-point velocity time series to be used to yield a flow length scale. In general, limits are exceeded in all environments for near-boundary flow measurements, but are met for y/d > 0·3 in most cases in sand-bed rivers, and for y/d > 0·4 in some gravel-bed environments. © 1998 John Wiley & Sons, Ltd.     

A complete dynamical ozone budget measured in the tropical marine boundary layer during PASE

The Pacific Atmospheric Sulfur Experiment (PASE) was a field mission that took place aboard the NCAR C-130 airborne laboratory over the equatorial Pacific Ocean near Christmas Island (Kirimati, Republic of Kiribati) during August?CSeptember, 2007. Eddy covariance measurements of the ozone fluxes at various altitudes above the ocean surface, along with simultaneous mapping of the horizontal gradients provided a unique opportunity to observe all of the dynamical components of the ozone budget in this remote marine environment. The results of six daytime and two sunrise flights indicate that vertical transport into the marine boundary layer from above and horizontal advection by the tradewinds are both important source terms, while photochemical destruction consisting of 82% photolysis (leading to OH production), 11% reaction with HO₂, and 7% reaction with OH provides the main sink. The overall photochemical lifetime of ozone in the marine boundary layer was found to be 6.5 days. Ocean uptake of ozone was observed to be fairly slow (mean deposition velocity of 0.024?±?0.014 cm s^?1) accounting for a diurnally averaged loss rate that was ??30% as large as the net photochemical destruction. From the measurement of deposition velocity an ozone reactivity of ??50 s^?1 in seawater is inferred. Due to the unprecedented measurement accuracy of the dynamical budget terms, unobserved photochemistry was able to be deduced, leading to the conclusion that 3.9?±?3.0 ppt (parts per trillion by volume) of NO is present on average in the daytime tropical marine boundary layer, broadly consistent with several previous studies in similar environments. It is estimated, however, that each ppt of BrO hypothetically present would counter each ppt of NO above the requisite 3.9 ppt needed for budget closure. The long-term budget of ozone is further analyzed in the buffer layer, between the boundary layer and free troposphere, and used to derive an entrainment velocity across the trade wind inversion of 0.51 ± 0.38 cm s^?1.     

Comparison of alternative representations of hydraulic-conductivity anisotropy in folded fractured-sedimentary rock: modeling groundwater flow in the Shenandoah Valley (USA)

A numerical representation that explicitly represents the generalized three-dimensional anisotropy of folded fractured-sedimentary rocks in a groundwater model best reproduces the salient features of the flow system in the Shenandoah Valley, USA. This conclusion results from a comparison of four alternative representations of anisotropy in which the hydraulic-conductivity tensor represents the bedrock structure as (model A) anisotropic with variable strikes and dips, (model B) horizontally anisotropic with a uniform strike, (model C) horizontally anisotropic with variable strikes, and (model D) isotropic. Simulations using the US Geological Survey groundwater flow and transport model SUTRA are based on a representation of hydraulic conductivity that conforms to bedding planes in a three-dimensional structural model of the valley that duplicates the pattern of folded sedimentary rocks. In the most general representation, (model A), the directions of maximum and medium hydraulic conductivity conform to the strike and dip of bedding, respectively, while the minimum hydraulic-conductivity direction is perpendicular to bedding. Model A produced a physically realistic flow system that reflects the underlying bedrock structure, with a flow field that is significantly different from those produced by the other three models.