Source model of the 11th July 2004 Zhongba earthquake revealed from the joint inversion of InSAR and seismological data

We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain ...     

Sewers as a Source and Sink of Chlorinated‐Solvent Groundwater Contamination,Marine Corps Recruit Depot,Parris Island,South Carolina

Groundwater contamination by tetrachloroethene and its dechlorination products is present in two partially intermingled plumes in the surficial aquifer near a former dry‐cleaning facility at Site 45, Marine Corps Recruit Depot, Parris Island, South Carolina. The northern plume originates from the vicinity of former above‐ground storage tanks. Free‐phase tetrachloroethene from activities in this area entered the groundwater. The southern plume originates at a nearby, new dry‐cleaning facility, but probably was the result of contamination released to the aquifer from a leaking sanitary sewer line from the former dry‐cleaning facility. Discharge of dissolved groundwater contamination is primarily to leaking storm sewers below the water table. The strong influence of sanitary sewers on source distribution and of storm sewers on plume orientation and discharge at this site indicates that groundwater‐contamination investigators should consider the potential influence of sewer systems at their sites.     

Magnitude and rupture duration from high frequency teleseismic P wave with projected landweber deconvolution

Earthquake magnitude and rupture duration are among the most important parameters characterizing an earthquake for the purpose of early tsunami warning. While they can be routinely determined from broadband P waveforms with iterative inver- sion procedures, the inversion procedures may fail when the rupture either lasts longer than the interval between P and later arrivals or requires too much time or human intervention. Little contaminated by later arrivals, high frequency P waves are useful for modeling earthquake source processes, though the envelope waveform is affected by strong scattering in lithosphere. With high frequency envelopes from aftershocks as Empirical Green’s Function (EGF), the coda effects can be removed and more accurate relative source time function (RSTF) of the main shock can be obtained. Assuming that RSTFs cannot be negative, we use the projected Landweber deconvolution method (PLD) to obtain high frequency RSTFs because PLD method has the advantage of non-negativity, causality, and compactness (finite duration). We are able to determine rupture durations of four large earthquakes: the 2004 Sumatra-Andaman earthquake, the 2005 Nias event, the 2006 Java event, and the 2011 Tokuko earthquake. The rupture durations of the Sumatra-Andaman, Nias, and Tohuko events are found to be around 550, 110, and 120 s respectively, consistent with previous studies. The rupture duration of the Java event is about 130 s, supporting that the Java event is a tsunami earthquake. The magnitudes of these earthquakes are found to depend on both the amplitude and the duration of the deconvolved waveforms, and can be approximated by integrating these waveforms.     

Reflection point sideslip and smear in imaging below dipping anisotropic media

Dipping anisotropic clastic strata are ubiquitous in fold and thrust belts. Geological structures below these strata will be mispositioned laterally and vertically on seismic images if we do not properly correct for seismic anisotropy during migration. The magnitude of this lateral mispositioning of a target structure varies with source‐receiver offset, so reflection points will be smeared in the final stacked image. Raytracing demonstrates the lateral‐position and smear phenomena when imaging structures below tilted transversely isotropic media. Analysis of the raytracing results predicts the quantity of lateral‐position error and reflection‐point smear on a seismic image. We created numerical‐model seismic data to show reflection‐point smear on synthetic seismic images and to evaluate the accuracy of the predictions from raytracing.     

Tomographic measurements of pore filling at infiltration fronts

For certain porous media and initial conditions, constant flux infiltrations show a saturation profile which exhibits overshoot. This overshoot is the cause of gravity driven fingering, cannot be described by standard models of unsaturated flow, and is likely controlled by the exact nature of the pore filling at the initial front. Here we report synchrotron X-ray micro-tomography measurements of the porous medium and measure which pores are filled by water and air at the initial wetting front as a function of flux. We find that at high fluxes all the pores are filled with water; for intermediate fluxes, the pores along the edge of the column remain unsaturated; and for low fluxes the pores in the bulk of the experimental column remain unsaturated. This suggests that the unsaturated overshoot conditions observed at higher fluxes are primarily an edge effect of the column. The results can help delineate the correct continuum models that can capture overshoot and gravity driven fingering.     

Simulating Gravity Changes in Topologically Realistic Driven Earthquake Fault Systems: First Results

Currently, GPS and InSAR measurements are used to monitor deformation produced by slip on earthquake faults. It has been suggested that another method to accomplish many of the same objectives would be through satellite-based gravity measurements. The Gravity Recovery and Climate Experiment (GRACE) mission has shown that it is possible to make detailed gravity measurements from space for climate dynamics and other purposes. To build the groundwork for a more advanced satellite-based gravity survey, we must estimate the level of accuracy needed for precise estimation of fault slip in earthquakes. We turn to numerical simulations of earthquake fault systems and use these to estimate gravity changes. The current generation of Virtual California (VC) simulates faults of any orientation, dip, and rake. In this work, we discuss these computations and the implications they have for accuracies needed for a dedicated gravity monitoring mission. Preliminary results are in agreement with previous results calculated from an older and simpler version of VC. Computed gravity changes are in the range of tens of μGal over distances up to a few hundred kilometers, near the detection threshold for GRACE.     

Oscillating pulsar polar gaps

An analytical model for oscillating pair creation above the pulsar polar cap is presented in which the parallel electric field is treated as a large amplitude, superluminal, electrostatic wave. An exact formalism for such wave is derived in one dimension and applied to both the low-density regime in which the pair plasma density is much lower than the corotating charge density and the high-density regime in which the pair plasma density is much higher than the corotating charge density. In the low-density regime, which is relevant during the phase leading to a pair cascade, a parallel electric field develops resulting in a rapid acceleration of particles. The rapid acceleration leads to bursts of pair production and the system switches to the oscillatory phase, corresponding to the high-density regime, in which pairs oscillate with net drift motion in the direction of wave propagation. Oscillating pairs lead to a current that oscillates with large amplitude about the Goldreich–Julian current. The drift motion can be highly relativistic if the phase speed of large amplitude waves is moderately higher than the speed of light. Thus, the model predicts a relativistic outflow of pairs, a feature that is required for avoiding overheating of the pulsar polar cap and is also needed for the pulsar wind.     

Petrogenesis of the late Cretaceous K‐rich volcanic rocks from the Central Pontide orogenic belt,North Turkey

Subduction‐related volcanic rocks are widespread in the Central Pontides of Turkey, and represented by the Hamsaros volcanic succession in the Sinop area to the north. The volcanic rocks display high‐K calc‐alkaline, shoshonitic and ultra‐K affinities. ⁴⁰Ar/³⁹Ar age data indicate that the rocks occurred during the Late Cretaceous (ca 82 Ma), and the volcanic suites were coeval. Primitive mantle‐normalized trace element patterns of all the lavas are characterized by strong enrichments in large ion lithophile elements (LILE) (Rb, Ba, K, and Sr), Th, U, Pb, and light rare earth elements (LREE; La, Ce) and prominent negative Nb, Ta, and Ti anomalies, all typical of subduction‐related lavas. There is a systematic increase in the enrichment of incompatible trace elements from the high‐K calc‐alkaline lavas through the shoshonitic to the ultra‐K lavas. In addition, the shoshonitic and ultra‐K lavas have significantly higher ⁸⁷Sr/⁸⁶Sr (0.70666–0.70834) and lower ¹⁴³Nd/¹⁴⁴Nd (0.51227–0.51236) initial ratios than coexisting high‐K calc‐alkaline lavas (⁸⁷Sr/⁸⁶Sr 0.70576–0.70613, ¹⁴³Nd/¹⁴⁴Nd 0.51245–0.51253). Geochemical and isotopic data show that the shoshonitic and ultra‐K rocks cannot be derived from the high‐K calc‐alkaline suite by any shallow level differentiation process, and point to a derivation from distinct mantle sources. The shoshonitic and ultra‐K rocks were derived from metasomatic veins related to melting of recycled subducted sediments, but the high‐K calc‐alkaline rocks from a lithospheric source metasomatized by fluids from subduction zone.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

Detecting active comets in the SDSS

Using a sample of serendipitously discovered active comets in the Sloan Digital Sky Survey (SDSS), we develop well-controlled selection criteria for greatly increasing the efficiency of comet identification in the SDSS catalogs. After follow-up visual inspection of images to reject remaining false positives, the total sample of SDSS comets presented here contains 19 objects, roughly one comet per 10 million other SDSS objects. The good understanding of selection effects allows a study of the population statistics, and we estimate the apparent magnitude distribution to r∼18, the ecliptic latitude distribution, and the comet distribution in SDSS color space. The most surprising results are the extremely narrow range of colors for comets in our sample (e.g. root-mean-square scatter of only ∼0.06 mag for the g-r color), and the similarity of comet colors to those of jovian Trojans. We discuss the relevance of our results for upcoming deep multi-epoch optical surveys such as the Dark Energy Survey, Pan-STARRS, and the Large Synoptic Survey Telescope (LSST), and estimate that LSST may produce a sample of about 10,000 comets over its 10-year lifetime.