Damage Survey of Water Supply Systems and Fragility Curve of PVC Water Pipelines in the Chi–Chi Taiwan Earthquake

Substantial damage to water supply systems, including water delivery pipelines, water treatment plants, reservoirs, and water storage tanks, was reported after the 1999 Chi–Chi Taiwan Earthquake. This paper first summarizes the damage survey and then presents the results of seismic fragility analysis for underground pipelines. Construction blueprints of the water delivery pipelines and repair work orders of 11 townships and cities in the disastrous area were digitized into a Geographical Information System (GIS) for analysis and assessment. With the aid of the GIS system, we found that PVC pipes made up 86% of water delivery pipelines while steel, cast iron, ductile iron, PE and others took the rest. Therefore, this paper focuses on the fragility analysis of PVC pipes. Three different methods were applied to derive the fragility relations between the PVC water pipes having nominal diameters (approximately inner diameters) greater than or equal to 65 mm and earthquake intensity parameters such as peak ground acceleration and peak ground velocity. The results were then examined with those of other countries. The discrepancy between our results and the empirical equation used by HAZUS, an earthquake loss estimation software developed by the Federal Emergency Management Agency was not significant.     

Structural damage detection using the optimal weights of the approximating artificial neural networks

This work presents a novel neural network‐based approach to detect structural damage. The proposed approach comprises two steps. The first step, system identification, involves using neural system identification networks (NSINs) to identify the undamaged and damaged states of a structural system. The partial derivatives of the outputs with respect to the inputs of the NSIN, which identifies the system in a certain undamaged or damaged state, have a negligible variation with different system errors. This loosely defined unique property enables these partial derivatives to quantitatively indicate system damage from the model parameters. The second step, structural damage detection, involves using the neural damage detection network (NDDN) to detect the location and extent of the structural damage. The input to the NDDN is taken as the aforementioned partial derivatives of NSIN, and the output of the NDDN identifies the damage level for each member in the structure. Moreover, SDOF and MDOF examples are presented to demonstrate the feasibility of using the proposed method for damage detection of linear structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Factorial and indicator kriging methods using a geographic information system to delineate spatial variation and pollution sources of soil heavy metals

This study uses factorial and indicator kriging methods with a geographic information system (GIS), to overlay and delineate the variation and pollution sources of soil heavy metals in fields in the north of Changhua County, Taiwan. Accordingly, the spatial variation and pollution probability patterns of soil cadmium, chromium and copper are estimated to support further environmental monitoring, remediation and planning. The overlay maps of industrial plants, irrigation channels and kriging estimates, constitute a GIS. This reveals that the short- and long-range variations and soil pollution by Cd, Cr and Cu strongly corresponded to the locations of irrigation systems and industrial plants in the area of study. The maps also verify most locations of suspected pollution sources from the GIS database. Moreover, the probability spatial patterns estimated by indicator kriging, based on the pollution thresholds of Cd, Cr and Cu, were also highly correlated with irrigation channels and industrial plants. The GIS, with factorial and indicator kriging under different combinations of variations and pollution probabilities of soil Cd, Cr and Cu, provide alternative scenarios for delineating the areas that contain high variation and high pollution probability of these heavy metals, to enable decision makers to identify pollution sources, and to further investigate and remedy soil pollution.     

Identification of soil degradation during earthquake excitations by Bayesian inference

A Bayesian inference approach is introduced to identify soil degradation behaviours at four downhole array sites. The approach of inference is based on a parametric time‐varying infinite impulse response filter model. The approach is shown to be adaptive to the changes of filter parameters and noise amplitudes. Four sites, including the Lotung (Taiwan), Chiba (Japan), Garner Valley (California), and Treasure Island (California) sites with downhole seismic arrays are analysed. Our results show two major types of soil degradation behaviour: the well‐known strain‐dependent softening, and reduction in stiffness that is not instantaneously recoverable. It is also found that both types of soil degradation are more pronounced in sandy soils than in clayey soils. The mechanism for the second type of soil degradation is not yet clear to the authors and suggested to be further studied. Copyright © 2003 John Wiley & Sons, Ltd.     

Earth's rotation variations: a wavelet analysis

Historical to up‐to‐date data of the minute variations in the solid Earth's rotation are subjected to a comprehensive time‐frequency wavelet analysis. The length‐of‐day for the period 1962–2012 confirms the presence of a prominent, robust 6‐year periodicity and reveals an anomalously strong 18.6‐year tidal oscillation as well as a ~13‐year quasi‐periodic signal. The polar‐motion excitation for the period 1900–2012 validates the existence of the ~26‐year Markowitz wobble and shows an ~8‐year quasi‐periodic signal, but no appreciable 18.6‐year periodicity. Although it is known that exchanges of angular momentum with the geophysical fluids are responsible for the rotational variations of the solid Earth, the exact physical mechanisms involved on interannual‐to‐decadal timescales are still far from understood.     

Layer-stripping reverse-time migration1

We present a layer-stripping method of migration for irregularly layered media in which first-order velocity discontinuities separate regions of constant or smoothly varying velocity. We use the reverse-time method to migrate seismic data layer by layer, from the surface downwards. As part of the migration of a given layer, the bottom boundary of the layer is defined based on power in the migrated signal, and a seismic section is collected along it. This new section serves as the boundary condition for migration in the next layer. This procedure is repeated for each layer, with the final image formed from the individual layer images. Layer-stripping migration consists of three steps: (1) layer definition, (2) wavefield extrapolation and imaging, and (3) boundary determination. The migration scheme when used with reverse-time extrapolation is similar to datuming with an imaging condition. The reverse-time method uses an explicit fourth-order time, tenth-order space, finite-difference approximation to the scalar wave equation. The advantages of layer-stripping reverse-time migration are: (1) it preserves the benefits of the reverse-time method by handling strong velocity contrasts between layers and steeply dipping structures; (2) it reduces computer memory and saves computation time in high-velocity layers, and (3) it allows interpretational control of the image. Post-stack layer-stripping reverse-time migration is illustrated with a synthetic CMP data example. Prestack migration is illustrated with a synthetic data set and with a marine seismic reflection profile across the Santa Maria Basin and the Hosgri Fault in central California.     

Chronology and petrogenesis of young achondrites,Shergotty, Zagami,and ALHA77005: late magmatism on a geologically active planet

A study was undertaken to determine the chronology, petrogenesis and relationships among the shergottites, Shergotty and Zagami and the unique achondrite ALHA77005. These meteorites are the product of a variety of complex processes.Petrogenesis: Chondrite-normalized abundance patterns of Shergotty and Zagami are very similar and show pronounced depletions of both the light REE (La-Nd) and heavy REE (Dy-Lu) relative to Sm-Gd. These characteristic depletions are even more pronounced for ALHA77005. The light REE depletion is qualitatively consistent with the presence of cumulus pyroxene and/or olivine in these meteorites, but trace element models show that the parental magmas of all three meteorites were probably also light REE depleted. Both trace element model calculations and combined Rb-Sr and Sm-Nd isotopic systematics show that the meteorites could not have been co-magmatic nor can ALHA77005 be representative of the source material of the shergottites. Light REE depletion of the parental magmas also implies light REE depletion of the source material. The Sm-Nd systematics of the shergottites require a time-averaged sub-chondritic (light REE enriched) Sm-Nd ratio since 4.6 AE ago. The Sm-Nd systematics of ALHA77005 permit a time-averaged super-chondritic (light REE depleted) Sm/Nd ratio if its crystallization age is less than T_ICE = 0.72 AE.Chronology. Rb-Sr internal isochrons for all three meteorites and a Sm-Nd internal isochron for Zagami are concordant at ~ 180 Myr. ³⁹Ar-⁴⁰Ar plateau ages of Shergotty and Zagami maskelynite are ~250–260 Myr. These ages apparently reflect resetting of these isotopic systems by shock metamorphism which converted the feldspar to maskelynite. The concordance of these ages suggests a single shock event during which the meteorites were in close physical proximity. The time of this event is most precisely given by the Rb-Sr age of 180 ± 4 Myr for Zagami.The crystallization ages of the meteorites were not precisely determined. Extreme upper limits are determined by Sm-Nd model ages relative to an eucrite initial ${}^{143}\text{Nd}{}^{144}\text{Nd = 0.505835}$ at 4.6 AE ago. These model ages for Shergotty, Zagami and ALHA77005 are 3600, 3500 and 2850 Myr, respectively. The Sm-Nd whole rock age of 1340 ± 60 Myr for the three meteorites gives the crystallization age if the Sm/Nd ratios of the precursor materials were always the same. We consider this 1340 Myr age as a “best estimate” upper limit. “Best estimate” lower limits for Shergotty and Zagami are taken from the average ³⁹Ar-⁴⁰Ar ages of 1200 and 900 Myr of pyroxene separates. The average ³⁹Ar-⁴⁰Ar age of a whole rock sample of ALHA77005 was 1600 Myr and can be partitioned between a low temperature (feldspar) phase and a high temperature (olivine + pyroxene + inclusions) “phase”. The average apparent ³⁹Ar-⁴⁰Ar age of the low temperature phase is ~1050 Myr, which is chosen as the “best estimate” lower limit to the age. The crystallization ages of Shergotty, Zagami and ALHA77005 probably lie within the ranges of 1200–1300, 900–1300 and 1000–1300 Myr, respectively. The Rb-Sr whole rock age of 4400 ± 400 Myr and single-stage BABI model ages of ~4800–5100 Myr are interpreted as reflecting differentiation of the parent body at ~4600 Myr ago.The complex geochemical and isotopic evolution recorded by these meteorites suggests a geologically active parent body capable of sustaining melting at two or more epochs in its history.