Investigating the impact of the Chi‐Chi earthquake on the occurrence of debris flows using artificial neural networks

Debris flows have caused enormous losses of property and human life in Taiwan during the last two decades. An efficient and reliable method for predicting the occurrence of debris flows is required. The major goal of this study is to explore the impact of the Chi‐Chi earthquake on the occurrence of debris flows by applying the artificial neural network (ANN) that takes both hydrological and geomorphologic influences into account. The Chen‐Yu‐Lan River watershed, which is located in central Taiwan, is chosen for evaluating the critical rainfall triggering debris flows. A total of 1151 data sets were collected for calibrating model parameters with two training strategies. Significant differences before and after the earthquake have been found: (1) The size of landslide area is proportioned to the occurrence of debris flows; (2) the amount of critical rainfall required for triggering debris flows has reduced significantly, about half of the original critical rainfall in the study case; and (3) the frequency of the occurrence of debris flows is largely increased. The overall accuracy of model prediction in testing phase has reached 96·5%; moreover, the accuracy of occurrence prediction is largely increased from 24 to 80% as the network trained with data from before the Chi‐Chi earthquake sets and with data from the lumped before and after the earthquake sets. The results demonstrated that the ANN is capable of learning the complex mechanism of debris flows and producing satisfactory predictions. Copyright © 2009 John Wiley & Sons, Ltd.     

Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

Review and new interpretation for the propagation characteristics associated with the 1999 Chi‐Chi earthquake faulting event

The Chi‐Chi earthquake (M_W = 7.6) took place in central western Taiwan in 1999. The earthquake caused reactivation of the Chelungpu Fault and resulted in 100‐km‐long surface ruptures. The fault strikes mostly north–south to NNE–SSW; however, the northern tip of the southern segment of the surface ruptures rotates clockwise to define an east–west trend, then jumps to a shorter NNW‐trending rupture. The largest vertical displacement is recorded in the Shihkang area of the Shihkang–Shangchi Fault Zone, where vertical slips are up to 8–10 m. The Shihkang–Shangchi Fault Zone displays a complex fault pattern as a linkage damage zone between two fault segments with the greatest concentration of faults and fractures. Our new interpretation, based on recent published geometric, kinematic, and geophysical studies on the Chi‐Chi earthquake fault, suggests that the Shihkang–Shangchi Fault Zone is not a simple termination zone, but may be an ‘overstep zone’ or a ‘transfer zone’. Slip analysis along the surface ruptures indicates that they are composed of three fault segments and the amount of slip partly depends on the intersection angle between slip direction and fault strike. Our numerical modeling for the area indicates that Coulomb stress changes are mainly concentrated on tips and bends of the surface ruptures. Slip patterns indicate that the fault propagates toward the northeast. Therefore, this study suggests high potential for future earthquake activity along the unruptured Shangchi segment. Hence, future geohazard studies should focus on the Shangchi segment to evaluate potential earthquakes, determine recurrence intervals, and reduce future earthquake hazards.     

Closed‐form aftershock reliability of damage‐cumulating elastic‐perfectly‐plastic systems

Major earthquakes (i.e., mainshocks) typically trigger a sequence of lower magnitude events clustered both in time and space. Recent advances of seismic hazard analysis stochastically model aftershock occurrence (given the main event) as a nonhomogeneous Poisson process with rate that decays in time as a negative power law. Risk management in the post‐event emergency phase has to deal with this short‐term seismicity. In fact, because the structural systems of interest might have suffered some damage in the mainshock, possibly worsened by damaging aftershocks, the failure risk may be large until the intensity of the sequence reduces or the structure is repaired. At the state‐of‐the‐art, the quantitative assessment of aftershock risk is aimed at building tagging, that is, to regulate occupancy. The study, on the basis of age‐dependent stochastic processes, derived closed‐form approximations for the aftershock reliability of simple nonevolutionary elastic‐perfectly‐plastic damage‐cumulating systems, conditional on different information about the structure. Results show that, in the case hypotheses apply, the developed models may represent a basis for handy tools enabling risk‐informed tagging by stakeholders and decision makers. Copyright © 2013 John Wiley & Sons, Ltd.     

Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration

In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.     

Visco‐acoustic prestack reverse‐time migration based on the time‐space domain adaptive high‐order finite‐difference method

It is important to include the viscous effect in seismic numerical modelling and seismic migration due to the ubiquitous viscosity in an actual subsurface medium. Prestack reverse‐time migration (RTM) is currently one of the most accurate methods for seismic imaging. One of the key steps of RTM is wavefield forward and backward extrapolation and how to solve the wave equation fast and accurately is the essence of this process. In this paper, we apply the time‐space domain dispersion‐relation‐based finite‐difference (FD) method for visco‐acoustic wave numerical modelling. Dispersion analysis and numerical modelling results demonstrate that the time‐space domain FD method has great accuracy and can effectively suppress numerical dispersion. Also, we use the time‐space domain FD method to solve the visco‐acoustic wave equation in wavefield extrapolation of RTM and apply the source‐normalized cross‐correlation imaging condition in migration. Improved imaging has been obtained in both synthetic and real data tests. The migration result of the visco‐acoustic wave RTM is clearer and more accurate than that of acoustic wave RTM. In addition, in the process of wavefield forward and backward extrapolation, we adopt adaptive variable‐length spatial operators to compute spatial derivatives to significantly decrease computing costs without reducing the accuracy of the numerical solution.     

Non‐linear‐Maxwell‐element‐type hysteretic control force

A simple non‐linear control law is proposed for reducing structural responses against seismic excitations. This law defines control force dynamics by one differential equation involving a non‐linear term that restrains the control force amplitude. If non‐linearity is neglected, the control force becomes the force in a Maxwell element, so it is called the non‐linear‐Maxwell‐element‐type (NMW) control force. The NMW control force vs. deformation relation plots hysteretic curves. The basic performance of an SDOF model with the NMW control force is examined for various conditions by numerical analyses. Furthermore, the control law is extended to fit an MDOF structural model, and an application example is shown. The computational results show that the NMW control force efficiently reduces structural responses. Copyright © 2000 John Wiley & Sons, Ltd.     

Tests and model calibration of high‐strength steel tubular beam‐to‐column and column‐base composite joints for moment‐resisting structures

Performance‐based engineering (PBE) methodologies allow for the design of more reliable earthquake‐resistant structures. Nonetheless, to implement PBE techniques, accurate finite element models of critical components are needed. With these objectives in mind, initially, we describe an experimental study on the seismic behaviour of both beam‐to‐column (BTC) and column‐base (CB) joints made of high‐strength steel S590 circular columns filled with concrete. These joints belonged to moment‐resisting frames (MRFs) that constituted the lateral‐force‐resisting system of an office building. BTC joints were conceived as rigid and of partial strength, whereas CB joints were designed as rigid and of full strength. Tests on a BTC joint composed of an S275 steel composite beam and high‐strength steel concrete‐filled tubes were carried out. Moreover, two seismic CB joints were tested with stiffeners welded to the base plate and anchor bolts embedded in the concrete foundation as well as where part of a column was embedded in the foundation with no stiffeners. A test programme was carried out with the aim of characterising these joints under monotonic, cyclic and random loads. Experimental results are presented by means of both force–interstory drift ratio and moment–rotation relationships. The outcomes demonstrated the adequacy of these joints to be used for MRFs of medium ductility class located in zones of moderate seismic hazard. Then, a numerical calibration of the whole joint subassemblies was successfully accomplished. Finally, non‐linear time‐history analyses performed on 2D MRFs provided useful information on the seismic behaviour of relevant MRFs. Copyright © 2015 John Wiley & Sons, Ltd.     

A comparison of snowmelt‐derived streamflow from temperature‐index and modified‐temperature‐index snow models

Reliable estimation of the volume and timing of snowmelt runoff is vital for water supply and flood forecasting in snow‐dominated regions. Snowmelt is often simulated using temperature‐index (TI) models due to their applicability in data‐sparse environments. Previous research has shown that a modified‐TI model, which uses a radiation‐derived proxy temperature instead of air temperature as its surrogate for available energy, can produce more accurate snow‐covered area (SCA) maps than a traditional TI model. However, it is unclear whether the improved SCA maps are associated with improved snow water equivalent (SWE) estimation across the watershed or improved snowmelt‐derived streamflow simulation. This paper evaluates whether a modified‐TI model produces better streamflow estimates than a TI model when they are used within a fully distributed hydrologic model. It further evaluates the performance of the two models when they are calibrated using either point SWE measurements or SCA maps. The Senator Beck Basin in Colorado is used as the study site because its surface is largely bedrock, which reduces the role of infiltration and emphasizes the role of the SWE pattern on streamflow generation. Streamflow is simulated using both models for 6 years. The modified‐TI model produces more accurate streamflow estimates (including flow volume and peak flow rate) than the TI model, likely because the modified‐TI model better reproduces the SWE pattern across the watershed. Both models also produce better performance when calibrated with SCA maps instead of point SWE data, likely because the SCA maps better constrain the space‐time pattern of SWE.     

Empirical‐based out‐of‐plane URM infill wall model accounting for the interaction with in‐plane demand

The role of masonry infills in the seismic behavior of reinforced concrete buildings has been widely studied in terms of their strength and stiffness contribution in the in‐plane (IP) direction, while fewer studies have been carried out on their response and modeling in the out‐of‐plane (OOP) direction. In this paper, the state of the art in code and literature provisions regarding infills' OOP capacity and seismic demand is presented, together with a review of the experimental tests that have been carried out to investigate infills' OOP behavior and the effects of IP‐OOP interaction. This review aims to collect an experimental database that is used to evaluate the effectiveness of literature and code provisions and to propose a semiempirical approach both for predicting infills' OOP strength, stiffness, and displacement capacity and for modeling the effects of IP displacement demand on OOP behavior and vice versa. Then, the state of the art on modeling of infills' OOP behavior and IP‐OOP interaction is presented together with a new macro model based on the proposed formulations and conceived to represent the IP and OOP behavior by taking into account the mutual interaction effects. Finally, the proposed model is used for an example application on two case‐study buildings, showing the effects of taking into account or neglecting the IP‐OOP interaction phenomena.     

Handling dip discrimination phenomenon in common‐reflection‐surface stack via combination of output‐imaging‐scheme and migration/demigration

In the application of a conventional common‐reflection‐surface (CRS) stack, it is well‐known that only one optimum stacking operator is determined for each zero‐offset sample to be simulated. As a result, the conflicting dip situations are not taken into account and only the most prominent event contributes to any a particular stack sample. In this paper, we name this phenomenon caused by conflicting dip problems as ‘dip discrimination phenomenon’. This phenomenon is not welcome because it not only leads to the loss of weak reflections and tips of diffractions in the final zero‐offset‐CRS stacked section but also to a deteriorated quality in subsequent migration. The common‐reflection‐surface stack with the output imaging scheme (CRS‐OIS) is a novel technique to implement a CRS stack based on a unified Kirchhoff imaging approach. As far as dealing with conflicting dip problems is concerned, the CRS‐OIS is a better option than a conventional CRS stack. However, we think the CRS‐OIS can do more in this aspect. In this paper, we propose a workflow to handle the dip discrimination phenomenon based on a cascaded implementation of prestack time migration, CRS‐OIS and prestack time demigration. Firstly, a common offset prestack time migration is implemented. Then, a CRS‐OIS is applied to the time‐migrated common offset gather. Afterwards, a prestack time demigration is performed to reconstruct each unmigrated common offset gather with its reflections being greatly enhanced and diffractions being well preserved. Compared with existing techniques dealing with conflicting dip problems, the technique presented in this paper preserves most of the diffractions and accounts for reflections from all possible dips properly. More importantly, both the post‐stacked data set and prestacked data set can be of much better quality after the implementation of the presented scheme. It serves as a promising alternative to other techniques except that it cannot provide the typical CRS wavefield attributes. The numerical tests on a synthetic Marmousi data set and a real 2D marine data set demonstrated its effectiveness and robustness.     

Large‐scale real‐time hybrid simulation of a three‐story steel frame building with magneto‐rheological dampers

A series of large‐scale real‐time hybrid simulations (RTHSs) are conducted on a 0.6‐scale 3‐story steel frame building with magneto‐rheological (MR) dampers. The lateral force resisting system of the prototype building for the study consists of moment resisting frames and damped brace frames (DBFs). The experimental substructure for the RTHS is the DBF with the MR dampers, whereas the remaining structural components of the building including the moment resisting frame and gravity frames are modeled via a nonlinear analytical substructure. Performing RTHS with an experimental substructure that consists of the complete DBF enables the effects of member and connection component deformations on system and damper performance to be accurately accounted for. Data from these tests enable numerical simulation models to be calibrated, provide an understanding and validation of the in‐situ performance of MR dampers, and a means of experimentally validating performance‐based seismic design procedures for real structures. The details of the RTHS procedure are given, including the test setup, the integration algorithm, and actuator control. The results from a series of RTHS are presented that includes actuator control, damper behavior, and the structural response for different MR control laws. The use of the MR dampers is experimentally demonstrated to reduce the response of the structure to strong ground motions. Comparisons of the RTHS results are made with numerical simulations. Based on the results of the study, it is concluded that RTHS can be conducted on realistic structural systems with dampers to enable advancements in resilient earthquake resistant design to be achieved. Copyright © 2014 John Wiley & Sons, Ltd.     

Least‐squares reverse‐time migration in a matrix‐based formulation

This paper describes least‐squares reverse‐time migration. The method provides the exact adjoint operator pair for solving the linear inverse problem, thereby enhancing the convergence of gradient‐based iterative linear inversion methods. In this formulation, modified source wavelets are used to correct the source signature imprint in the predicted data. Moreover, a roughness constraint is applied to stabilise the inversion and reduce high‐wavenumber artefacts. It is also shown that least‐squares migration implicitly applies a deconvolution imaging condition. Three numerical experiments illustrate that this method is able to produce seismic reflectivity images with higher resolution, more accurate amplitudes, and fewer artefacts than conventional reverse‐time migration. The methodology is currently feasible in 2‐D and can naturally be extended to 3‐D when computational resources become more powerful.     

Simulation of masonry out‐of‐plane failure modes by multi‐body dynamics

The seismic assessment of the local failure modes in existing masonry buildings is currently based on the identification of the so‐called local mechanisms, often associated with the out‐of‐plane wall behavior, whose stability is evaluated by static force‐based approaches and, more recently, by some displacement‐based proposals. Local mechanisms consist of kinematic chains of masonry portions, often regarded as rigid bodies, with geometric nonlinearity and concentrated nonlinearity in predefined contact regions (unilateral no‐tension behavior, possible sliding with friction). In this work, the dynamic behavior of local mechanisms is simulated through multi‐body dynamics, to obtain the nonlinear response with efficient time history analyses that directly take into account the characteristics of the ground motion. The amplification/filtering effects of the structure are considered within the input motion. The proposed approach is validated with experimental results of two full‐scale shaking‐table tests on stone masonry buildings: a sacco‐stone masonry façade tested at Laboratório Nacional de Engenharia Civil and a two‐storey double‐leaf masonry building tested at European Centre for Training and Research in Earthquake Engineering (EUCENTRE). Copyright © 2015 John Wiley & Sons, Ltd.     

Seismic self‐centering steel beam‐to‐column moment connections using bolted friction devices

This paper first presents the force–deformation relationship of a post‐tensioned (PT) steel beam‐to‐column connection constructed with bolted web friction devices (FDs). This paper then describes the test program conducted in the National Center for Research on Earthquake Engineering, Taiwan, on four bolted FDs and four full‐scale PT beam‐to‐column moment connection subassemblies using the FDs. Tests confirm that (1) the hysteretic behavior of four bolted FDs is very stable, (2) the friction coefficient between the steel plate and the brass shim is about 0.34, (3) the proposed force–deformation relationships reasonably predict the experimental responses of the PT connections under cyclically increasing deformations up to a beam peak rotation of 0.05 rad, and (4) the decompression moments do not degrade as beam cyclic deformations increase. Copyright © 2007 John Wiley & Sons, Ltd.     

Simplified macro‐model for infill masonry walls considering the out‐of‐plane behaviour

One of the main challenges in earthquake risk mitigation is the assessment of existing buildings not designed according to modern codes and the development of effective techniques to strengthen these structures. Particular attention should be given to RC frame structures with masonry infill panels, as demonstrated by their poor performance in recent earthquakes in Europe. Understanding the seismic behaviour of masonry‐infilled RC frames presents one of the most difficult problems in structural engineering. Analytical tools to evaluate infill–frame interaction and the failure mechanisms need to be further studied. This research intends to develop a simplified macro‐model that takes into account the out‐of‐plane behaviour of the infill panels and the corresponding in‐plane and out‐of‐plane interaction when subjected to seismic loadings. Finally, a vulnerability assessment of an RC building will be performed in order to evaluate the influence of the out‐of‐plane consideration in the building response. Copyright © 2015 John Wiley & Sons, Ltd.     

Multipathing in three‐parameter common‐image gathers from reverse‐time migration

Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region.     

Solving the complex near‐surface problem using 3D data‐driven near‐surface layer replacement

In many land seismic situations, the complex seismic wave propagation effects in the near‐surface area, due to its unconsolidated character, deteriorate the image quality. Although several methods have been proposed to address this problem, the negative impact of 3D complex near‐surface structures is still unsolved to a large extent. This paper presents a complete 3D data‐driven solution for the near‐surface problem based on 3D one‐way traveltime operators, which extends our previous attempts that were limited to a 2D situation. Our solution is composed of four steps: 1) seismic wave propagation from the surface to a suitable datum reflector is described by parametrized one‐way propagation operators, with all the parameters estimated by a new genetic algorithm, the self‐adjustable input genetic algorithm, in an automatic and purely data‐driven way; 2) surface‐consistent residual static corrections are estimated to accommodate the fast variations in the near‐surface area; 3) a replacement velocity model based on the traveltime operators in the good data area (without the near‐surface problem) is estimated; 4) data interpolation and surface layer replacement based on the estimated traveltime operators and the replacement velocity model are carried out in an interweaved manner in order to both remove the near‐surface imprints in the original data and keep the valuable geological information above the datum. Our method is demonstrated on a subset of a 3D field data set from the Middle East yielding encouraging results.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.