Kinematic characteristics and the influence of reference velocities of phase‐shift‐plus‐interpolation and extended‐split‐step‐Fourier migration methods

The phase‐shift‐plus‐interpolation and extended‐split‐step‐Fourier methods are wavefield‐continuation algorithms for seismic migration imaging. These two methods can be applied to regions with complex geological structures. Based on their unified separable formulas, we show that these two methods have the same kinematic characteristics by using the theory of pseudodifferential operators. Numerical tests on a Marmousi model demonstrate this conclusion. Another important aspect of these two methods is the selection of reference velocities and we explore the influence of the selection of reference velocities by comparing the geometric progression method and the statistical method. We show that the geometric progression method is simple but does not take into account the velocity distribution while the statistical approach is relatively complex but reflects the velocity distribution.     

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.     

Spectral decomposition and de‐noising via time‐frequency and space‐wavenumber reassignment

The reassignment method remaps the energy of each point in a time‐frequency spectrum to a new coordinate that is closer to the actual time‐frequency location. Two applications of the reassignment method are developed in this paper. We first describe time‐frequency reassignment as a tool for spectral decomposition. The reassignment method helps to generate more clear frequency slices of layers and therefore, it facilitates the interpretation of thin layers. The second application is to seismic data de‐noising. Through thresholding in the reassigned domain rather than in the Gabor domain, random noise is more easily attenuated since seismic events are more compactly represented with a relatively larger energy than the noise. A reconstruction process that permits the recovery of seismic data from a reassigned time‐frequency spectrum is developed. Two approaches of the reassignment method are used in this paper, one of which is referred to as the trace by trace time reassignment that is mainly used for seismic spectral decomposition and another that is the spatial reassignment that is mainly used for seismic de‐noising. Synthetic examples and two field data examples are used to test the proposed method. For comparison, the Gabor transform method, inversion‐based method and common deconvolution method are also used in the examples.     

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.     

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.     

Local‐ and field‐scale infiltration into vertically non‐uniform soils with spatially‐variable surface hydraulic conductivities

We studied the problem of local‐ and field‐scale infiltration over a particular class of heterogeneous soils. At the local scale, the soils are described as being vertically non‐uniform, with the saturated hydraulic conductivity continuously decreasing with depth according to a power law function. Analogous to the Green–Ampt model, analytical expressions are first developed for local‐scale infiltration using a sharp front approximation, and model results are compared with numerical solutions of the Richards equation. These results show that saturation does not occur from below in soils with such vertical non‐uniformity, thereby allowing for the use of a sharp front approximation. Because of vertical non‐uniformity, ponding conditions are achieved locally even for rainfall rates less than the surface saturated hydraulic conductivity. Furthermore, infiltration rates asymptotically approach zero at long times. To determine field‐scale infiltration properties, the spatial variability in the surface saturated hydraulic conductivity is represented by a log‐normal random field. Using cumulative infiltration as the independent variable, expressions are developed for the ensemble mean of field‐scale infiltration and the expected time for a given depth of water to infiltrate over the field. Surface horizontal heterogeneity is found to control field‐scale infiltration at small times, whereas local vertical non‐uniformity exerts a strong control at long times. Copyright © 2011 John Wiley & Sons, Ltd.     

Tests on low‐ductility RC frames under high‐ and low‐frequency excitations

The response of low‐ductility reinforced concrete (RC) frames, designed typically for a non‐seismic region, subjected to two frequencies of base excitations is studied. Five half‐scaled, two‐bay, two‐storey, RC frames, each approximately 5 m wide by 3.3 m high, were subjected to both horizontal and/or vertical base excitations with a frequency of 40 Hz as well as a lower frequency of about 4 Hz (close to the fundamental frequency) using a shake table. The imposed acceleration amplitude ranged from 0.2 to 1.2g. The test results showed that the response characteristics of the structures differed under high‐ and low‐frequency excitations. The frames were able to sustain high‐frequency excitations without damage but were inadequate for low‐frequency excitations, even though the frames exhibited some ductility. Linear‐elastic time‐history analysis can predict reasonably well the structural response under high‐frequency excitations. As the frames were not designed for seismic loads, the reinforcement detailing may not have been adequate, based on the crack pattern observed. The effect of vertical excitation can cause significant additional forces in the columns and moment reversals in the beams. The ‘strong‐column, weak‐beam’ approach for lateral load RC frame design is supported by experimental observations. Copyright © 2001 John Wiley & Sons, Ltd.     

Experimental implementation and verification of multi‐degrees‐of‐freedom effective force testing

This paper presents an experimental implementation and verification of multi‐degrees‐of‐freedom effective force testing (MDOF‐EFT). An experimental setup that consists of a two‐degrees‐of‐freedom structural system and two hydraulic actuators at the Johns Hopkins University was utilized in this study. First, experimental system identification was performed to develop compatible analytical models for the multi‐input and multi‐output systems. Dynamics of the control plant, that is, the valve‐to‐force relations, were modeled with a rational polynomial transfer function matrix and delay components. By using the analytical model, a centralized decoupling loop‐shaping force feedback controller was designed such that the forces are uncoupled and the loop transfer functions have desirable dynamic characteristics in the frequency domain. Then, a series of harmonic force and earthquake simulation tests were performed to assess capabilities and limitations of MDOF‐EFT. Experimental results showed that the dynamic forces in the two actuators were accurately controlled to provide tracking while the system was stable and robust for the entire period of the experiment. Furthermore, earthquake simulation tests with increased levels of the reference forces demonstrated the feasibility of MDOF‐EFT with highly nonlinear test structures. Copyright © 2013 John Wiley & Sons, Ltd.     

Triangulation hand‐held laser‐scanning (TriHaLaS) for micro‐ and meso‐habitat surveys in streams

A novel hand‐held laser‐based stream bed survey system is presented. The system facilitates the capture of detailed 3D mapping of shallow (< 0.7 m) riverbed topography in sections approximately 4 m by 2 m. The system includes a stationary reference system, which projects three laser sheets (two at offset angles), within which a hand‐held monitoring pole is moved. The unique configuration of the light sheets intercepts with the pole as it moves within the survey area providing an exact horizontal location. Pole tilt is compensated for by an inertial measurement unit on the pole, and the height above the bed of the pole and submerged scanning laser are monitored relative to the horizontal laser sheet. Verification and application measurements demonstrate high resolution and accuracy in the horizontal (~5 mm) and vertical (~1 mm) direction. The system can be applied at sites where a free view is blocked and other optical through‐water methods fail. It is appropriate for studies on riverbed statistics and dynamics, which necessitate non‐invasive in‐situ surveys. Copyright © 2017 John Wiley & Sons, Ltd.     

Near‐fault ground motions,and the response of elastic and inelastic single‐degree‐of‐freedom (SDOF) systems

In order to investigate the response of structures to near‐fault seismic excitations, the ground motion input should be properly characterized and parameterized in terms of simple, yet accurate and reliable, mathematical models whose input parameters have a clear physical interpretation and scale, to the extent possible, with earthquake magnitude. Such a mathematical model for the representation of the coherent (long‐period) ground motion components has been proposed by the authors in a previous study and is being exploited in this article for the investigation of the elastic and inelastic response of the single‐degree‐of‐freedom (SDOF) system to near‐fault seismic excitations. A parametric analysis of the dynamic response of the SDOF system as a function of the input parameters of the mathematical model is performed to gain insight regarding the near‐fault ground motion characteristics that significantly affect the elastic and inelastic structural performance. A parameter of the mathematical representation of near‐fault motions, referred to as ‘pulse duration’ (T_P), emerges as a key parameter of the problem under investigation. Specifically, T_P is employed to normalize the elastic and inelastic response spectra of actual near‐fault strong ground motion records. Such normalization makes feasible the specification of design spectra and reduction factors appropriate for near‐fault ground motions. The ‘pulse duration’ (T_P) is related to an important parameter of the rupture process referred to as ‘rise time’ (τ) which is controlled by the dimension of the sub‐events that compose the mainshock. Copyright © 2004 John Wiley & Sons, Ltd.     

Bidirectional substructure pseudo‐dynamic tests and analysis of a full‐scale two‐story buckling‐restrained braced frame

A two‐story buckling‐restrained brace (BRB) frame was tested under bidirectional in‐plane and out‐of‐plane loading to evaluate the BRB stability and gusset plate design. The test comprised pseudo‐dynamic loadings using the 1999 Chi‐Chi earthquake scaled to the 50%, 10%, and 2% probability of exceedance in 50 years and a cyclic regime of increasing amplitudes of up to 3.0% story drift ratio (SDR). The specimen had a unique configuration where the beams were connected to the columns through shear tabs welded to the column flanges and bolted to the beam webs. Stable hysteretic behavior with only minor cracking at the gusset‐to‐column welds was observed under the pseudo‐dynamic tests, with maximum in‐plane and out‐of‐plane SDRs of 2.24% and 1.47% respectively. Stable behavior continued into the cyclic test where fracture of the gusset‐to‐column welds occurred in the first cycle to simultaneous bidirectional SDR of 3.0%. The observed BRB stability is consistent with a methodology developed for BRB frames under simultaneous in‐plane and out‐of‐plane drifts. The specimen behavior was studied using a finite element model. It was shown that gusset plates are subjected to a combination of BRB force and frame action demands, with the latter increasing the gusset‐to‐beam and gusset‐to‐column interface demands by an average of 69% and 83% respectively. Consistent with the test results, failure at the gusset‐to‐column interfaces is computed when frame action demands are included, thus confirming that not considering frame action demands may results in unconservative gusset plate designs. Copyright © 2016 John Wiley & Sons, Ltd.     

A Rosenbrock‐W method for real‐time dynamic substructuring and pseudo‐dynamic testing

A variant of the Rosenbrock‐W integration method is proposed for real‐time dynamic substructuring and pseudo‐dynamic testing. In this variant, an approximation of the Jacobian matrix that accounts for the properties of both the physical and numerical substructures is used throughout the analysis process. Only an initial estimate of the stiffness and damping properties of the physical components is required. It is demonstrated that the method is unconditionally stable provided that specific conditions are fulfilled and that the order accuracy can be maintained in the nonlinear regime without involving any matrix inversion while testing. The method also features controllable numerical energy dissipation characteristics and explicit expression of the target displacement and velocity vectors. The stability and accuracy of the proposed integration scheme are examined in the paper. The method has also been verified through hybrid testing performed of SDOF and MDOF structures with linear and highly nonlinear physical substructures. The results are compared with those obtained from the operator splitting method. An approach based on the modal decomposition principle is presented to predict the potential effect of experimental errors on the overall response during testing. Copyright © 2009 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.     

Full waveform inversion of SH‐ and Love‐wave data in near‐surface prospecting

We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls.     

Seismic zoning for initial‐ and total‐cost minimization

We consider two design criteria to study seismic zoning. In the first, codes require that structures be designed for some specified values. Zoning is then optimal when it minimizes the expected present value of the initial costs of all structures to be built in the region being zoned. In the second criterion, it is designed so that the present value of the total cost is minimized, including initial and maintenance costs as well as losses due to damage and failure. We will call these criteria zoning for the initial‐ and total‐cost minimization, respectively. It is shown that under certain conditions, the boundaries coincide with isoparametric curves and the problem may be solved in one dimension. We also deal with problems not reducible to a single dimension. Different methods are proposed to solve the various kinds of problems. The work ends with some illustrative examples. Copyright © 2000 John Wiley & Sons, Ltd.     

Obtaining the steady‐state drawdown solutions of constant‐head and constant‐flux tests

The solutions of constant‐head and constant‐flux tests are commonly used to predict the temporal or spatial drawdown distribution or to determine aquifer parameters. Theis and Thiem equations, for instance, are well‐known transient and steady‐state drawdown solutions, respectively, of the constant‐flux test. It is known that the Theis equation is not applicable to the case where the aquifer has a finite boundary or the pumping time tends to infinity. On the other hand, the Thiem equation does not apply to the case where the aquifer boundary is infinite. However, the issue of obtaining the Thiem equation from the transient drawdown solution has not previously been addressed. In this paper, the drawdown solutions for constant‐head and constant‐flux tests conducted in finite or infinite confined aquifers with or without consideration of the effect of the well radius are examined comprehensively. Mathematical verification and physical interpretation of the solutions to these two tests converging or not converging to the Thiem equation are presented. The result shows that there are some finite‐domain solutions for these two tests that can converge to the Thiem equation when the time becomes infinitely large. In addition, the time criteria to give a good approximation to the finite‐domain solution by the infinite‐domain solution and the Thiem equation are investigated and presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Automatic calibration of HEC‐HMS using single‐objective and multi‐objective PSO algorithms

This study presents single‐objective and multi‐objective particle swarm optimization (PSO) algorithms for automatic calibration of Hydrologic Engineering Center‐ Hydrologic Modeling Systems rainfall‐runoff model of Tamar Sub‐basin of Gorganroud River Basin in north of Iran. Three flood events were used for calibration and one for verification. Four performance criteria (objective functions) were considered in multi‐objective calibration where different combinations of objective functions were examined. For comparison purposes, a fuzzy set‐based approach was used to determine the best compromise solutions from the Pareto fronts obtained by multi‐objective PSO. The candidate parameter sets determined from different single‐objective and multi‐objective calibration scenarios were tested against the fourth event in the verification stage, where the initial abstraction parameters were recalibrated. A step‐by‐step screening procedure was used in this stage while evaluating and comparing the candidate parameter sets, which resulted in a few promising sets that performed well with respect to at least three of four performance criteria. The promising sets were all from the multi‐objective calibration scenarios which revealed the outperformance of the multi‐objective calibration on the single‐objective one. However, the results indicated that an increase of the number of objective functions did not necessarily lead to a better performance as the results of bi‐objective function calibration with a proper combination of objective functions performed as satisfactorily as those of triple‐objective function calibration. This is important because handling multi‐objective optimization with an increased number of objective functions is challenging especially from a computational point of view. Copyright © 2012 John Wiley & Sons, Ltd.     

Post‐Synthesis of Ti‐SBA‐15 in Supercritical CO2‐Ethanol Solution

Ti‐SBA‐15 materials with Ti incorporated into the framework of SBA‐15 and controllable Ti contents were successfully prepared via a post‐treatment route in supercritical CO₂‐ethanol solution, followed by calcination. The resultant Ti‐SBA‐15 materials were characterized by means of different techniques including X‐ray diffraction, X‐ray photoelectron spectroscopy, transmission electron microscopy, IR analysis, and N₂ sorption analysis. It was demonstrated that the resultant materials retained a structure similar to that of the parent SBA‐15, and Ti was incorporated into the framework of SBA‐15. At high Ti content, i.e, Si/Ti = 5:1, a TiO₂ phase formed and was coated onto the inner surface of SBA‐15 in addition to the incorporation of the Ti in the framework. The BET surface areas of the Ti‐SBA‐15 samples decreased with increasing Ti content. The presence of small amounts of H₂O in the reaction medium may have resulted in some TiO₂ nanoparticles being uniformly distributed in the pores of the SBA‐15 accompanying the incorporation of Ti in the SBA‐15 framework.