Seismic vertical array analysis for phase decomposition

We propose a vertical array analysis method that decomposes complex seismograms into body and surface wave time histories by using a velocity structure at the vertical array site. We assume that the vertical array records are the sum of vertically incident plane P and S waves, and laterally incident Love and Rayleigh waves. Each phase at the surface is related to that at a certain depth by the transfer function in the frequency domain; the transfer function is obtained by Haskell's matrix method, assuming a 1-D velocity structure. Decomposed P , S and surface waves at the surface are estimated from the vertical array records and the transfer functions by using a least-squares method in the frequency domain; their time histories are obtained by the inverse Fourier transform. We carried out numerical tests of this method based on synthetic vertical array records consisting of vertically incident plane P and S waves and laterally incident plane Love and Rayleigh waves. Perfect results of the decomposed P , S , Love and Rayleigh waves were obtained for synthetic records without noise. A test of the synthetic records in which a small amount of white noise was added yielded a reasonable result for the decomposed P , S and surface waves. We applied this method to real vertical array records from the Ashigara valley, a moderate-sized sedimentary valley. The array records from two earthquakes occurring at depths of 123 and 148 km near the array (epicentral distance of about 31 km) exhibited long-duration later phases. The analysis showed that duration of the decomposed S waves was a few seconds and that the decomposed surface waves appeared a few seconds after the direct S -wave arrival and had very long duration. This result indicated that the long-duration later phases were generated not by multireflected S waves, but by basin-induced surface waves.     

Relationship between radar reflectivity and the time scale of warm rain formation in a global cloud-resolving model

The connections between radar reflectivity and the time scale of warm rain formation are examined within a global cloud-resolving model. The parameterizations formulae of auto-conversion and accretion processes in the model reveal specific relationships between the time scale for auto-conversion and radar reflectivity of cloud water as well as between the time scale for accretion and radar reflectivity of rain water. The overall time scale for warm rain formation, determined by combined contributions from these processes, is found to relate with total radar reflectivity in the manner that varies with cloud-rain composition between auto-conversion and accretion limits. The global statistics from the model output reveals that the time scale is closely related with the total radar reflectivity, thus suggesting that the radar reflectivity is a gross measure of the warm rain time scale. The relationship developed is applied to both model-simulated and CloudSat-observed radar reflectivities to compare the time scales of warm rain formation between observation and the model. Comparison of the time scale so derived reveals significant differences between the model and CloudSat observations. These differences suggest that the simulated cloud-rain composition in the model is biased toward larger rain water contents and smaller content of cloud water compared to reality due to an accelerated cloud-to-rain water conversion in the model.     

Mineralogy,Fluid Inclusions,and Sulfur Isotopes of the Huanzala Deposits,Peru: Early Skarn and Late Polymetallic Replacement Style Mineralizations

Carbonate‐replacement polymetallic mineralization at the Huanzala deposits (9°51′S, 77°00′W) was conducted in two contrasting stages that occurred in almost the same location. Early‐stage mineralization has a relation with a granodiorite porphyry stock, whereas the late‐stage mineralization is genetically associated with quartz porphyry sills. The early stage involved low to intermediate sulfidation Cu–Zn–(Pb) mineralization associated with metasomatic skarn, and the late stage involved high to intermediate sulfidation Cu–Zn–Pb–(Mn) mineralization associated with hydrothermal alteration characterized by paragonitic sericitization. The orebodies are hosted by steeply dipping (approximately 60°NE) Lower Cretaceous carbonate rocks in a relatively narrow range of approximately 4 km in horizontal extent and less than 1 km in depth. The pathway of the early‐stage brine‐derived fluids (300–>400°C, >33 wt% NaCl equivalent) along a plot of log against 1000/T is best explained by the progressive dual decline of the value and the temperature under rock‐buffering conditions; this decline saw the pathway progress through the stability field of pyrrhotite to reach that of pyrite and promoted a decrease in FeS from 14.5 to 1.6 mol% in the sphalerite. In contrast, an explanation for the pathway of the late‐stage fluids (140–290°C, 3–13 wt% NaCl equivalent) is given by an almost isothermal decline at approximately 270°C, with passing through the stability field of pyrite–bornite to reach that of chalcopyrite, promoting an increase in FeS from 0.1 to 1.6 mol% in the sphalerite, suggesting gas‐buffering conditions. The ore formation pressure records in the fluid inclusions illustrate an approximately 2‐km erosion during the roughly 2‐Myr total lifetime of the hydrothermal system.     

Characteristics of amorphous kerogens fractionated from terrigenous sedimentary rocks

A preliminary attempt to fractionate amorphous kerogens from terrigenous bulk kerogen by a benzene-water two phase partition method under acidic condition was made. Microscopic observation revealed that amorphous kerogens and structured kerogens were fractionated effectively by this method. Characteristics of the amorphous and structured kerogens fractionated by this method were examined by some chemical analyses and compared with those of the bulk kerogen and humic acid isolated from the same rock sample (Haizume Formation, Pleistocene, Japan). The elemental and infrared (IR) analyses showed that the amorphous kerogen fraction had the highest atomic$\text{H}\text{C}$ ratio and the lowest atomic $\text{N}\text{C}$ ratio and was the richest in aliphatic structures and carbonyl and carboxyl functional groups. Quantities of fatty acids from the saponification products of each geopolymer were in agreement with the results of elemental and IR analyses. Distribution of the fatty acids was suggestive that more animal lipids participate in the formation of amorphous kerogens because of the abundance of relatively lower molecular weight fatty acids (such as C₁₆ and C₁₈ acids) in saponification products of amorphous kerogens. On the other hand, although the amorphous kerogen fraction tends to be rich in aliphatic structures compared with bulk kerogen of the same rock samples, van Krevelen plots of elemental compositions of kerogens from the core samples (Nishiyama Oil Field, Tertiary, Japan) reveal that the amorphous kerogen fraction is not necessarily characterized by markedly high atomic $\text{H}\text{C}$ ratio. This was attributed to the oxic environment of deposition and the abundance of biodegraded terrestrial amorphous organic matter in the amorphous kerogen fraction used in this work.     

Intra-oceanic island arc origin for Iratsu eclogites of the Sanbagawa belt,central Shikoku,southwest Japan

New geochemical and Sr–Nd isotopic data for the Iratsu eclogite and surrounding metamorphic rocks of the Sanbagawa belt, Japan, show that, while the protoliths of the metamorphic rocks formed in a variety of tectonic settings, the Iratsu body represents a deeply subducted and accreted island arc. The igneous protoliths of eclogites and garnet amphibolites were probably generated from a mantle source that had components of both a depleted mantle modified by slab-released fluid (as seen in a negative Nb anomaly) and an enriched mantle, similar to that of ocean island basalts (OIB). Fractional crystallization modeling indicates that the protoliths of some garnet clinopyroxenites from the Iratsu body are cumulates from a basaltic magma that crystallized under high O₂ and H₂O fugacities in the middle to lower crust. The source characteristics and crystallization conditions suggest that the protoliths of the Iratsu rocks formed in an oceanic island arc. Quartz eclogites from the marginal zone of the Iratsu body have geochemical signatures similar to turbidites from the Izu–Bonin island arc (as seen in a negative Nb anomaly and a concave REE pattern). The protoliths might be volcaniclastic turbidites that formed in a setting proximal to the oceanic island arc. Geochemical and isotopic signatures of the surrounding mafic schists are similar to normal (N-) and enriched (E-) mid-ocean-ridge basalt (MORB), and distinct from the rocks from the Iratsu body. The protoliths of the mafic schists likely formed in a plume-influenced mid-ocean ridge or back-arc basin. Pelitic schists from the surrounding rocks and pelitic gneisses from the marginal zone of the Iratsu body have evolved, continental geochemical signatures (as seen in a negative ε_Nd(t) value (~?5)), consistent with their origin as continent-derived trench-fill turbidites.     

沙虫干呈味及功能性成分研究

对沙虫干抽提物成分的组成进行的研究结果表明 :游离氨基酸含量 w为 10 .3% ,占抽提物含氮物质的 77.9%。抽提物组成中 ,呈甘味的 Gly(3.2 % )、Ala(2 .5% )和呈鲜味的 IMP(0 .2 7% )、Glu(0 .2 5% )、琥珀酸 (0 .35% )等成分同沙虫干独特的风味有关 ,Na 、Cl-等无机离子也同其呈味相关。此外 ,牛磺酸含量 w高达 3.2 % ,是富含牛磺酸的一种海洋生物 ;精氨酸的含量 w亦较高 ,为 1.2 %。     

Attenuation of compressional waves in peridotite measured as a function of temperature at 200 MPa

A technique has been developed to determine attenuation in rocks at high temperature using a gas-media, high-pressure apparatus. A pulse transmission technique and a spectral ratio method are used to study compressional seismic properties of rocks. Seismic waves are transmitted to and from the sample through buffer rods of mullite. The effect of seismic wave reflections within the sample assembly are cancelled out by taking ratios of the spectra measured at different temperatures. In order to obtain good signal-to-noise ratio for resolving the attenuation at high pressure and temperature, special care is taken in the sample assembly and the ultrasonic coupling between the sample, buffer rods and transducers. A very tight connection of the sample-buffer rod-transducer is essential for obtaining high frequency signals (>300 kHz) at high temperature. A small mass is attached to each outside end of the transducer to drive low frequency signals (<250 kHz) into the sample. Before attenuation measurements, the sample and the buffer rods are tightly compacted in a platinum tube at high pressure and room temperature to ensure pressure seal of the sample assembly. The frequency range of measurement covers 50 to 450 kHz for the sample. Attenuation is very small in the buffer rod compared to the sample for the entire temperature range of the study. Because of the small attenuation, a wide frequency band of 50 kHz to 3.2 MHz can be covered for investigating the attenuation in the buffer rod. The technique has been used to measure attenuation at high confining pressure, and temperatures including sub- and hyper-solidus of upper mantle rocks. Therefore, effects of partial melting on attenuation can be studied.The method is applied to the attenuation measurement in a peridotite as a function of temperature to 1225°C at 200 MPa confining pressure. At high temperature, signal amplitude decays more rapidly at high frequency than at low frequency, from which attenuation (andQ) can be determined using a spectral ratio method. No frequency dependence ofQ is resolved for both the sample and the buffer rod over the entire temperature and frequency ranges of the measurement. The results show thatQ decreases rapidly with increasing temperature even in the temperature range below the solidus of peridotites. Such temperature sensitivity ofQ is probably more useful to probe thermal structure in the upper mantle than that of conductivity at temperatures below the solidus. The results in this study are compared with available seismic velocity, electrical conductivity and solidus data for peridotites, suggesting that there is no discontinuous change in both mechanical and electrical properties of peridotites at the solidus temperature. Even at hypersolidus temperatures, it appears that velocity drops and conductivity increases continuously (not abruptly) with increasing melt fraction. This implies that mechanical and electrical properties of the upper mantle will gradually change at the boundary where the geotherm crosses the solidus.     

Matrix method evaluating an internal radiation field in a plane-parallel atmosphere

Suppose that the atmosphere consists of homogeneous sublayers whose scattering and transmission matrices are known, then it is shown in this paper how to determine the intensity and polarization parameters of internal radiation field by the matrix method.     

Dust Trails of 8P/Tuttle and the Unusual Outbursts of the Ursid Shower

We calculate the position of dust trails from comet 8P/Tuttle, in an effort to explain unusual Ursid meteor shower outbursts that were seen when the comet was near aphelion. Comet 8P/Tuttle is a Halley-type comet in a 13.6-year orbit, passing just outside of Earth's orbit. We find that the meteoroids tend to be trapped in the 12:14 mean motion resonance with Jupiter, while the comet librates in a slightly shorter period orbit around the 13:15 resonance. It takes 6 centuries to decrease the perihelion of the meteoroid orbits enough to intersect Earth's orbit, during which time the meteoroids and comet separate in mean anomaly by 6 years, thus explaining the 6-year lag between the comet's return and Ursid outbursts. The resonances also prevent dispersion along the comet orbit and limit viewing to only one year in each return. We identified past dust trail encounters with dust trails from 1392 (Dec. 1945) and 1378 (Dec. 1986) and predicted another outburst on 2000 December 22 at around 7:29 and 8:35 UT, respectively, from dust trails dating to the 1405 and 1392 returns. This event was observed from California using video and photographic techniques. At the same time, five Global-MS-Net stations in Finland, Japan, and Belgium counted meteors using forward meteor scatter. The outburst peaked at 8:06±07 UT, December 22, at zenith hourly rate ∼90 per hour, and the Ursid rates were above half peak intensity during 4.2 h. We find that most Ursid orbits do scatter around the anticipated positions, confirming the link with comet 8P/Tuttle and the epoch of ejection. The 1405 and 1392 dust trails appear to have contributed similar amounts to the activity profile. Some orbits provide a hint of much older debris being present as well. This work is the strongest evidence yet for the relevance of mean motion resonances in Halley-type comet dust trail evolution.     

Experimental study of the alteration process of labradorite in acid hydrothermal solutions

Labradorite was altered artificially by HC1 solution ranging from M = 1 to M = 0.003 at 245 and 230°C. The products of alteration were examined by X-ray diffraction, electron microscopy, electron diffraction, infrared spectroscopy and the electron microprobe and the solution was analyzed chemically.Amorphous silica only was formed in solutions with M_HCl = 1 and M_HCl = 0.3. In a solution with M_HCl = 0.2, amorphous silica was initially formed, later dissolved and replaced by kaolinite. A mixture of microcrystalline boehmite and amorphous aluminosilicate was formed, altering to kaolinite in solutions with M_HCl = 0.1 and 0.3. Small amounts of kaolinite were initially formed but the alteration soon stopped in solution with M_HCl = 0.003. Relationships between the alteration processes and pH of the solutions can be roughly explained by using solubility diagrams assuming the congruent dissolution of labradorite and precipitation of the products in partial equilibrium. However, these assumptions are not valid with strongly acid solutions.The rate of dissolution of labradorite is controlled not only by its surface area, but also by the diffusion of matter through the layer of alteration products.