Analysis of fracture induced scattering of microseismic shear-waves

Fractures are pervasive features within the Earth’s crust and have a significant influence on the multi-physical response of the subsurface. The presence of coherent fracture sets often leads to observable seismic scattering enabling seismic techniques to remotely locate and characterise fracture systems. In this study, we confirm the general scale-dependence of seismic scattering and provide new results specific to shear-wave propagation. We do this by generating full waveform synthetics using finite-difference wave simulation within an isotropic background model containing explicit fractures. By considering a suite of fracture models having variable fracture density and fracture size, we examine the widening effect of wavelets due to scattering within a fractured medium by using several different approaches, such as root-mean-square envelope analysis, shear-wave polarisation distortion, differential attenuation analysis and peak frequency shifting. The analysis allows us to assess the scattering behavior of parametrised models in which the propagation direction is either normal or parallel to the fracture surfaces. The quantitative measures show strong observable deviations for fractures size on the order of or greater than the dominant seismic wavelength within the Mie and geometric scattering regime for both propagation normal and parallel to fracture strike. The results suggest that strong scattering is symptomatic of fractures having size on the same order of the probing seismic wave.     

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.     

Modeling patterns of coral bleaching at a remote Central Pacific atoll

A mild bleaching event (9.2% prevalence) at Palmyra Atoll occurred in response to the 2009 ENSO, when mean water temperature reached 29.8-30.1 °C. Prevalence among both abundant and sparse taxa varied with no clear pattern in susceptibility relating to coral morphology. Seven taxon-specific models showed that turbidity exacerbated while prior exposure to higher background temperatures alleviated bleaching, with these predictors explaining an average 16.3% and 11.5% variation in prevalence patterns, respectively. Positive associations occurred between bleaching prevalence and both immediate temperature during the bleaching event (average 8.4% variation explained) and increased sand cover (average 3.7%). Despite these associations, mean unexplained variation in prevalence equalled 59%. Lower bleaching prevalence in areas experiencing higher background temperatures suggests acclimation to temperature stress among several coral genera, while WWII modifications may still be impacting the reefs via shoreline sediment re-distribution and increased turbidity, exacerbating coral bleaching susceptibility during periods of high temperature stress.     

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.     

Proximate environmental drivers of coral communities at Palmyra Atoll: establishing baselines prior to removing a WWII military causeway

A management proposal aims to partly remove a WWII military causeway at Palmyra Atoll to improve lagoon water circulation and alleviate sedimentation stress on the southeast backreef, an area of high coral cover and diversity. This action could result in a shift in sedimentation across reef sites. To provide management advice, we quantified the proximate environmental factors driving scleractinian coral cover and community patterns at Palmyra. The proportion of fine sedimentation was the optimal predictor of coral cover and changes in community structure, explaining 23.7% and 24.7% of the variation between sites, respectively. Scleractinian coral cover was negatively correlated with increases in fine sedimentation. Removing the causeway could negatively affect the Montipora corals that dominate the western reef terrace, as this genus was negatively correlated with levels of fine sedimentation. The tolerance limits of corals, and sediment re-distribution patterns, should be determined prior to complete removal of the causeway.     

Local permit ownership in Alaska salmon fisheries

Changes in ownership of limited entry permits by “local” residents of the region where a fishery occurs may have significant economic and social implications for fishery-dependent regions. This paper examines changes in local permit ownership in Alaska salmon fisheries, for which a long-term decline in rural local permit ownership is an important policy concern. Theoretically, permit markets allocate permits over time to the individuals who are willing to pay the most for them. Any factors that differentially affect what local and non-local residents are willing to pay for permits may affect the equilibrium share of permits held by local residents. For remote rural fisheries in particular, these may include differences between local and non-local residents with respect to access to and costs of financing permits and boats, costs of travel to the fishery, opportunity costs of participation in the fishery, and many other factors. As a fishery increases in profitability, differences between local and non-local residents in access to financing matter more while differences in costs of travel and opportunity costs matter less in the relative ranking of what local and non-local residents are willing to pay for permits. This tends to increase the share of non-local residents among buyers willing to pay the market price for permits, reducing the equilibrium share of permits held by local residents. This leads to a conflict between two important policy goals: increasing fishery profitability and maintaining rural local permit ownership. Consistent with predictions of this theory, the local share of permit ownership in Alaska salmon fisheries is negatively related to permit prices (an indicator of fishery profitability).     

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).     

Response of the shortgrass steppe plant community to fire

Fire is an important driver of ecological pattern and process in grasslands worldwide, although its role in semi-arid systems is less well known. We used published studies and new experimental research to 1) provide a synthesis of existing knowledge of fire in the semi-arid grasslands of the North American Great Plains, and 2) assess the degree of similarity in semi-arid and mesic grassland responses to fire in this region. Based on published studies, burning has neutral to negative effects on aboveground productivity in semi-arid grasslands and variable effects on plant communities. To more rigorously assess fire effects, replicated experimental plots were established in ungrazed shortgrass steppe in northern Colorado and prescribed spring fire was applied in 2006 and 2007, 2006 only, or not at all. Aboveground net primary productivity decreased or remained unchanged with burning. Plant community changes included increases in perennial forbs, decreases in annual grasses and a positive response in annual forbs to the combination of fire and a wet spring in 2007. Combined, these results indicate that post-fire changes in productivity in semi-arid grasslands are neutral to negative, in contrast to positive responses in mesic grasslands, and not strongly negative as previously assumed.     

The Sloan Digital Sky Survey monitor telescope pipeline

The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi‐step process which involves data from three different telescopes: the 1.0‐m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which was used to establish the SDSS standard star network); the SDSS 0.5‐m Photometric Telescope (PT) at the Apache Point Observatory (APO), NewMexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and the SDSS 2.5‐m telescope at APO (which obtains the imaging data for the SDSS proper). In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data from the single‐CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0‐m and the PT). We also describe transformation equations that convert photometry on the USNO‐1.0m u ′g ′r ′i ′z ′ system to photometry the SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the semi‐automated PT factory, which runs MTPIPE in the day‐to‐day standard SDSS operations at Fermilab. Finally, we discuss the use of MTPIPE in current SDSS‐related projects, including the Southern u ′g ′r ′i ′z ′ Standard Star project, the u ′g ′r ′i ′z ′ Open Star Clusters project, and the SDSS extension (SDSS‐II). (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.