Modelling converted seismic waveforms in isotropic and anisotropic 1-D gradients: discontinuous versus continuous gradient representations

Over the past decade, there have been numerous receiver function studies directed at imaging the lithosphere-asthenosphere boundary (LAB). Although it is generally accepted that receiver function phases observed in these studies are derived from physical mode conversions at depth within the lithosphere-asthenosphere transition, it is still debatable as to whether these phases are directly indicative of the LAB. This is because interpretation of receiver function LAB signals relies on understanding the elastic characteristics of the Earth??s outer thermal boundary layer. The main issues for receiver function imaging are the sharpness of the elastic material property transition and, more importantly, what specifically are the material gradients. To test the various transition models, a forward modelling approach is required that allows accurate waveform synthetics for a range of discontinuous and continuous gradients in anisotropic, elastic media. We present a derivation of the reflection and transmission response for continuous one-dimensional (1-D) gradients in generally anisotropic elastic media. We evaluate the influence of 1-D isotropic and anisotropic elastic gradients on the seismic waveform by comparing numerical results of models for discontinuous and continuous transitions. The results indicate that discontinuous representations using layers each with uniform parameters and with thicknesses on the order of approximately 1/3 to 1/8 of the dominant seismic wavelength can be used to accurately model P-to-S and S-to-P mode conversions due to continuous transitions of both isotropic and anisotropic elastic properties. From a practical point of view, when comparing synthetic modelling with observation, this constraint can be relaxed further. The presence of signal noise and/or the result of receiver function stacking techniques will likely obscure these subtle waveform e ff ects. Hence this study suggests that accurate synthetic waveforms for LAB transitions can be modelled with discontinuous gradient representations using a reasonable number of discrete transition layers with layer thicknesses no greater than 1/2 to 1/3 the dominant seismic wavelength.     

Stability analysis of the 2007 Chehalis lake landslide based on long-range terrestrial photogrammetry and airborne LiDAR data

On December 4th 2007, a 3-Mm³ landslide occurred along the northwestern shore of Chehalis Lake. The initiation zone is located at the intersection of the main valley slope and the northern sidewall of a prominent gully. The slope failure caused a displacement wave that ran up to 38 m on the opposite shore of the lake. The landslide is temporally associated with a rain-on-snow meteorological event which is thought to have triggered it. This paper describes the Chehalis Lake landslide and presents a comparison of discontinuity orientation datasets obtained using three techniques: field measurements, terrestrial photogrammetric 3D models and an airborne LiDAR digital elevation model to describe the orientation and characteristics of the five discontinuity sets present. The discontinuity orientation data are used to perform kinematic, surface wedge limit equilibrium and three-dimensional distinct element analyses. The kinematic and surface wedge analyses suggest that the location of the slope failure (intersection of the valley slope and a gully wall) has facilitated the development of the unstable rock mass which initiated as a planar sliding failure. Results from the three-dimensional distinct element analyses suggest that the presence, orientation and high persistence of a discontinuity set dipping obliquely to the slope were critical to the development of the landslide and led to a failure mechanism dominated by planar sliding. The three-dimensional distinct element modelling also suggests that the presence of a steeply dipping discontinuity set striking perpendicular to the slope and associated with a fault exerted a significant control on the volume and extent of the failed rock mass but not on the overall stability of the slope.     

Particle fluxes in San Pedro Basin, California: A four-year record of sedimentation and physical forcing

Moored sediment traps were deployed from January 2004 through December 2007 at depths of 550 and 800 m in San Pedro Basin (SPB), CA (33°33.0′N, 118°26.5′W). Additionally, floating sediment traps were deployed at 100 and 200 m for periods of 12-24 h during spring 2005, fall 2007, and spring 2008. Average annual fluxes of mass, particulate organic carbon (POC), ??¹³C_org, particulate organic nitrogen (PON), ??¹⁵N-PON, biogenic silica (bSiO₂), calcium carbonate (CaCO₃), and detrital material (non-biogenic) were coupled with climate records and used to examine sedimentation patterns, vertical flux variability, and organic matter sources to this coastal region. Annual average flux values were determined by binning data by month and averaging the monthly averages. The average annual fluxes to 550 m were 516±42 mg/m² d for mass (sdom of the monthly averages, n=117), 3.18±0.26 mmol C/m² d for POC (n=111), 0.70±0.05 mmol/m² d for CaCO₃ (n=110), 1.31±0.21 mmol/m² d for bSiO₂ (n=115), and 0.35±0.03 mmol/m² d for PON (n=111). Fluxes to 800 and to 550 m were similar, within 10%. Annual average values of ??¹³C_org at 550 m were −21.8±0.2‰ (n=108), and ??¹⁵N averages were 8.9±0.2‰ (n=95). The timing of both high and low flux particle collection was synchronous between the two traps. Given the frequency of trap cup rotation (4-11 days), this argues for particle settling rates ≥83 m/d for both high and low flux periods. The moored traps were deployed over one of the wettest (2004-2005, 74.6 cm rainfall) and driest (2006-2007, 6.6 cm) rain years on record. There was poor correlation (Pearson's correlation coefficient, 95% confidence interval) of detrital mass flux with: C_org/N ratio (r=0.10, p=0.16); ??¹⁵N (r=−0.19, p=0.02); and rainfall (r=0.5, p=0.43), suggesting that runoff does not immediately cause increases in particle fluxes 15 km offshore. ??¹³C_org values suggest that most POC falling to the basin floor is marine derived. Coherence between satellite-derived chlorophyll a records from the trap location (±9 km² resolution) and SST data indicates that productivity and export occurs within a few days of upwelling and both of these parameters are reasonable predictors of POC export, with a time lag of a few days to 2 weeks (with no time lag—SeaWiFS chlorophyll a and POC flux, r=0.25, p=0.0014; chlorophyll a and bSiO₂ flux, r=0.28, p=0.0002).     

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application

Models for water transfer in the crop–soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop–soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop–soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration.     

The GPS Toolbox

The GPS Toolbox is dedicated to highlighting algorithms utilized by GPS engineers and scientists. If you have an interesting algorithm you would like to share with our readers or if you have a topic you would like to see covered in a future column, contact us at gps-toolbox@ngs.noaa.gov. To comment on the algorithms presented here, or to leave a request for an algorithm you may be looking for, visit our Web site (http://www.ngs.noaa.gov/gps-toolbox). ? 2000 John Wiley & Sons, Inc.     

Tracing the formation history of intermediate-age star clusters in the Small Magellanic Cloud

Colour–magnitude diagrams (CMDs) are presented for the first time for 10 star clusters projected on to the Small Magellanic Cloud (SMC). The photometry was carried out in the Washington system C and T ₁ filters allowing the determination of ages by means of the magnitude difference between the red giant clump and the main-sequence turnoff (MSTO), and metallicities from the red giant branch (RGB) locus. The clusters all have ages in the range 1.5–4 Gyr and metallicities between  −1.3 < [Fe/H] < −0.6  , with respective errors of ∼0.5 Gyr and 0.3 dex. This increases substantially the sample of intermediate-age clusters in the SMC with well-derived parameters. We combine our results with those for other clusters in the literature to derive as large and homogeneous a data base as possible (totalling 26 clusters) in order to study global effects. We find evidence for two peaks in the age distribution of SMC clusters, at ∼6.5 and 2.5 Gyr, in good agreement with previous hints involving smaller samples. The most recent peak occurs at a time that corresponds to a very close encounter between the Large Magellanic Cloud (LMC) and the SMC according to the recent dynamical models of Bekki et al. that they used to explain the enhancement of LMC clusters with this age. It appears cluster formation may have been similarly stimulated in the SMC by this encounter as well. We also find very good agreement between cluster ages and metallicities and the prediction from a bursting model from Pagel and Tautvaišienė with a burst that occurred 3 Gyr ago. These two lines of evidence together favour a bursting cluster formation history as opposed to a continuous one for the SMC.     

Identifying the source of a strong fullerene envelope arising from laser desorption mass spectrometric analysis of meteoritic insoluble organic matter

Insoluble organic matter (IOM) has been obtained from two carbonaceous chondrite meteorites and subjected to analysis by laser desorption mass spectrometry (LDMS) using standard operating conditions that were optimized for fullerene detection (3-6 μJ pulses at 337 nm focused to a spot size of approximately 100 μm in diameter). The preparation process yields no free C₆₀ in the IOM, and other experiments suggest that this material does not contain appreciable amounts of fullerenes. Nevertheless, a pronounced high-mass envelope is observed in LDMS, extending from 720 amu to about 4000 amu, with peaks spaced apart every 24 amu (corresponding to the gain or loss of C₂ units). We attribute this high-mass envelope to the existence of various fullerene molecules. The present work demonstrates that these fullerene molecules are created by the laser desorption laser ionization process under typical laser conditions used for studying free fullerenes in organic solvent extracts of natural samples (toluene and 1,2,4-trichlorobenzene). The implications of this false positive detection of fullerene molecules on the reports of fullerenes in other meteoritic samples have been investigated by introducing IOM into typical fullerene extraction procedures and examining the LDMS results. We found that IOM is capable of producing false positive signals in these experiments. The effect of ambient laboratory contamination producing fullerene signals is also described. It is found that extensive centrifugation of the meteoritic extracts is able to reduce the observed fullerene envelope, which points to an association of this envelope with IOM particulates that have passed through the filtering steps. We suggest the exercise of extreme caution in interpreting fullerene data from LDMS experiments.