Variations of b‐values at the western edge of the Ryukyu Subduction Zone,north‐east Taiwan

Using earthquakes relocated in north‐east Taiwan, we estimated b‐value distribution along a cross‐section located near the Ryukyu slab edge, and four b‐value anomalous areas are evidenced: (1) a high b‐value body lying on top of a low Vp, low Vs and high Vp/Vs sausage‐like body was considered as a region of enhanced partial melt or water supply above which seismicity occurs; (2) beneath the Ilan Plain, an anomalous area characterized by b‐values slightly higher than 1.1 might give evidence to the magma conduits to the Kueishantao Island; (3) above the Ryukyu Wadati‐Benioff zone, at depths ranging from 90 to 110 km, a high b‐value anomaly might correspond to the depth where dehydration occurs in the subducting oceanic plate; and (4) a low b‐value area located within the Ryukyu slab, at depths ranging from 70 to 90 km, might be linked to the compressive mechanisms shown by focal mechanisms and the bending of the subducting plate.     

Direct numerical simulations of flow over two-dimensional and three-dimensional ripples and implication to sediment transport: Steady flow

A well resolved and highly accurate direct numerical simulation (DNS) solver has been developed to understand the implication of hydrodynamics to sediment transport. In the first part of the study we focus on steady flow over two-dimensional and three-dimensional ripples at two Reynolds numbers Re_τ = 180 and 400 (defined by channel half-height and wall-friction velocity) in a channel geometry. The DNS scheme is based on a fourth-order vertical velocity and second-order vertical vorticity formulation, which resolves the difficulties in pressure boundary condition encountered when solving the Navier–Stokes equations. The complex boundary introduced due to the ripples has been imposed in the Cartesian domain using an elegant immersed boundary method. Detailed hydrodynamic analysis has revealed turbulence statistics (in particular, the higher order) and henceforth, the flow structures are sensitive — whether the ripples are two-dimensional or three-dimensional. The importance of fluctuating component of the bottom stress in addition to its mean component; and its significance to sediment transport and ripple migration speed have been investigated.     

Development and validation of a finite element morphological model for shallow water basins

Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.     

Petrology,mineralogy, and oxygen isotope compositions of aluminum‐rich chondrules from CV3 chondrites

Bulk major element composition, petrography, mineralogy, and oxygen isotope compositions of twenty Al‐rich chondrules (ARCs) from five CV3 chondrites (Northwest Africa [NWA] 989, NWA 2086, NWA 2140, NWA 2697, NWA 3118) and the Ningqiang carbonaceous chondrite were studied and compared with those of ferromagnesian chondrules and refractory inclusions. Most ARCs are marginally Al‐richer than ferromagnesian chondrules with bulk Al₂O₃ of 10–15 wt%. ARCs are texturally similar to ferromagnesian chondrules, composed primarily of olivine, pyroxene, plagioclase, spinel, Al‐rich glass, and metallic phases. Minerals in ARCs have intermediate compositions. Low‐Ca pyroxene (Fs_0.6–8.8Wo_0.7–9.3) has much higher Al₂O₃ and TiO₂ contents (up to 12.5 and 2.3 wt%, respectively) than that in ferromagnesian chondrules. High‐Ca pyroxene (Fs_0.3–2.0Wo_33–54) contains less Al₂O₃ and TiO₂ than that in Ca,Al‐rich inclusions (CAIs). Plagioclase (An_77–99Ab_1–23) is much more sodic than that in CAIs. Spinel is enriched in moderately volatile element Cr (up to 6.7 wt%) compared to that in CAIs. Al‐rich enstatite coexists with anorthite and spinel in a glass‐free chondrule, implying that the formation of Al‐enstatite was not due to kinetic reasons but is likely due to the high Al₂O₃/CaO ratio (7.4) of the bulk chondrule. Three ARCs contain relict CAIs. Oxygen isotope compositions of ARCs are also intermediate between those of ferromagnesian chondrules and CAIs. They vary from ?39.4‰ to 13.9‰ in δ¹⁸O and yield a best fit line (slope = 0.88) close to the carbonaceous chondrite anhydrous mineral (CCAM) line. Chondrules with 5–10 wt% bulk Al₂O₃ have a slightly more narrow range in δ¹⁸O (?32.5 to 5.9‰) along the CCAM line. Except for the ARCs with relict phases, however, most ARCs have oxygen isotope compositions (>?20‰ in δ¹⁸O) similar to those of typical ferromagnesian chondrules. ARCs are genetically related to both ferromagnesian chondrules and CAIs, but the relationship between ARCs and ferromagnesian chondrules is closer. Most ARCs were formed during flash heating and rapid cooling processes like normal chondrules, only from chemically evolved precursors. ARCs extremely enriched in Al and those with relict phases could have had a hybrid origin (Krot et al. 2002) which incorporated refractory inclusions as part of the precursors in addition to ferromagnesian materials. The occurrence of melilite in ARCs indicates that melilite‐rich CAIs might be present in the precursor materials of ARCs. The absence of melilite in most ARCs is possibly due to high‐temperature interactions between a chondrule melt and the solar nebula.     

Petrogenesis of the Northwest Africa 4898 high‐Al mare basalt

Northwest Africa (NWA) 4898 is the only low‐Ti, high‐Al basaltic lunar meteorite yet recognized. It predominantly consists of pyroxene (53.8 vol%) and plagioclase (38.6 vol%). Pyroxene has a wide range of compositions (En_12–62Fs_25–62Wo_11–36), which display a continuous trend from Mg‐rich cores toward Ca‐rich mantles and then to Fe‐rich rims. Plagioclase has relatively restricted compositions (An_87–96Or_0–1Ab_4–13), and was transformed to maskelynite. The REE zoning of all silicate minerals was not significantly modified by shock metamorphism and weathering. Relatively large (up to 1 mm) olivine phenocrysts have homogenous inner parts with Fo ~74 and sharply decrease to 64 within the thin out rims (~30 μm in width). Four types of inclusions with a variety of textures and modal mineralogy were identified in olivine phenocrysts. The contrasting morphologies of these inclusions and the chemical zoning of olivine phenocrysts suggest NWA 4898 underwent at least two stages of crystallization. The aluminous chromite in NWA 4898 reveals that its high alumina character was inherited from the parental magma, rather than by fractional crystallization. The mineral chemistry and major element compositions of NWA 4898 are different from those of 12038 and Luna 16 basalts, but resemble those of Apollo 14 high‐Al basalts. However, the trace element compositions demonstrate that NWA 4898 and Apollo 14 high‐Al basalts could not have been derived from the same mantle source. REE compositions of its parental magma indicate that NWA 4898 probably originated from a unique depleted mantle source that has not been sampled yet. Unlike Apollo 14 high‐Al basalts, which assimilated KREEPy materials during their formation, NWA 4898 could have formed by closed‐system fractional crystallization.     

Observation of saturnian stream particles in the interplanetary space

In January 2004 the dust instrument on the Cassini spacecraft detected the first high-velocity grain expelled from Saturn - a so-called stream particle. Prior to Cassini’s arrival at Saturn in July 2004 the instrument registered 801 faint impacts, whose impact signals showed the characteristic features of a high-velocity impact by a tiny grain. The impact rates as well as the directionality of the stream particles clearly correlate with the sector structure of the interplanetary magnetic field (IMF). The Cosmic Dust Analyser (CDA) registered stream particles dominantly during periods when the IMF direction was tangential to the solar wind flow and in the prograde direction. This finding provides clear evidence for a continuous outflow of tiny dust grains with similar properties from the saturnian system. Within the compressed part of co-rotating interaction regions (CIRs) of the IMF, characterized by enhanced magnetic field strength and compressed solar wind plasma, CDA observed impact bursts of faster stream particles. We find that the bursts result from the stream particles being sped up inside the compressed CIR regions. Our analysis of the stream-particle dynamics inside rarefaction regions of the IMF implies that saturnian stream particles have sizes between 2 and 9 nm and exit the saturnian systems closely aligned with the planet’s ring plane with speeds in excess of 70 km s⁻¹.     

Long-term polar motion prediction using normal time–frequency transform

This paper presents normal time–frequency transform (NTFT) application in harmonic/quasi-harmonic signal prediction. Particularly, we use the normal wavelet transform (a special NTFT) to make long-term polar motion prediction. Instantaneous frequency, phase and amplitude of Chandler wobble, prograde and retrograde annual wobbles of Earth’s polar motion are analyzed via the NTFT. Results show that the three main wobbles can be treated as quasi-harmonic processes. Current instantaneous harmonic information of the three wobbles can be acquired by the NTFT that has a kernel function constructed with a normal half-window function. Based on this information, we make the polar motion predictions with lead times of 1 year and 5 years. Results show that our prediction skills are very good with long lead time. An abnormality in the predictions occurs during the second half of 2005 and first half of 2006. Finally, we provide the future (starting from 2013) polar motion predictions with 1- and 5-year leads. These predictions will be used to verify the effectiveness of the method proposed in this paper.     

Simulation of probability distributions commonly used in hydrological frequency analysis

Random variable simulation has been applied to many applications in hydrological modelling, flood risk analysis, environmental impact assessment, etc. However, computer codes for simulation of distributions commonly used in hydrological frequency analysis are not available in most software libraries. This paper presents a frequency‐factor‐based method for random number generation of five distributions (normal, log–normal, extreme‐value type I, Pearson type III and log‐Pearson type III) commonly used in hydrological frequency analysis. The proposed method is shown to produce random numbers of desired distributions through three means of validation: (1) graphical comparison of cumulative distribution functions (CDFs) and empirical CDFs derived from generated data; (2) properties of estimated parameters; (3) type I error of goodness‐of‐fit test. An advantage of the method is that it does not require CDF inversion, and frequency factors of the five commonly used distributions involves only the standard normal deviate. Copyright © 2006 John Wiley & Sons, Ltd.