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.     

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.     

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.     

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.     

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.     

Catchment area‐based evaluation of the AMC‐dependent SCS‐CN‐based rainfall–runoff models

Using a large set of rainfall–runoff data from 234 watersheds in the USA, a catchment area‐based evaluation of the modified version of the Mishra and Singh (2002a) model was performed. The model is based on the Soil Conservation Service Curve Number (SCS‐CN) methodology and incorporates the antecedent moisture in computation of direct surface runoff. Comparison with the existing SCS‐CN method showed that the modified version performed better than did the existing one on the data of all seven area‐based groups of watersheds ranging from 0·01 to 310·3 km². Copyright © 2005 John Wiley & Sons, Ltd.     

Shaking‐table tests of a full‐scale single‐story masonry veneer wood‐frame structure

This paper presents the findings of shaking‐table experiments conducted to examine the seismic performance of a full‐scale, one‐story, wood‐framed structure with masonry veneer. The structure was designed and constructed in accordance with current U.S. code provisions. The veneer was attached to the wood backing with two kinds of metal anchors, corrugated ties fastened with 8d nails and rigid ties fastened with #8 screws. The tests have shown that the use of nails to fasten veneer anchors to the wood studs is highly unreliable due to the high variation of the nail extraction capacity, which can be influenced by the moisture content of the wood. Other than this, both the wood frame and the masonry veneer performed well under severe ground motions far exceeding a design level earthquake for Seismic Design Category D. Good performance was observed for the rigid veneer ties, which were attached to the wood studs with screws. The results have shown that the veneer walls parallel to the direction of shaking helped to restrain the motion of the wood structure and therefore should not be simply treated as added mass. The detailing of wood roof diaphragms requires special attention in consideration of the out‐of‐plane inertia force of the veneer that can be transmitted through the top plate of the wood‐stud wall to the rim joist. Copyright © 2010 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.     

Image‐domain wavefield tomography with extended common‐image‐point gathers

Waveform inversion is a velocity‐model‐building technique based on full waveforms as the input and seismic wavefields as the information carrier. Conventional waveform inversion is implemented in the data domain. However, similar techniques referred to as image‐domain wavefield tomography can be formulated in the image domain and use a seismic image as the input and seismic wavefields as the information carrier. The objective function for the image‐domain approach is designed to optimize the coherency of reflections in extended common‐image gathers. The function applies a penalty operator to the gathers, thus highlighting image inaccuracies arising from the velocity model error. Minimizing the objective function optimizes the model and improves the image quality. The gradient of the objective function is computed using the adjoint state method in a way similar to that in the analogous data‐domain implementation. We propose an image‐domain velocity‐model building method using extended common‐image‐point space‐ and time‐lag gathers constructed sparsely at reflections in the image. The gathers are effective in reconstructing the velocity model in complex geologic environments and can be used as an economical replacement for conventional common‐image gathers in wave‐equation tomography. A test on the Marmousi model illustrates successful updating of the velocity model using common‐image‐point gathers and resulting improved image quality.     

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.     

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.     

An adaptive free‐surface expression for three‐dimensional finite‐difference frequency‐domain modelling of elastic wave

Finite‐difference frequency‐domain modelling of seismic wave propagation is attractive for its efficient solution of multisource problems, and this is crucial for full‐waveform inversion and seismic imaging, especially in the three‐dimensional seismic problem. However, implementing the free surface in the finite‐difference method is nontrivial. Based on an average medium method and the limit theorem, we present an adaptive free‐surface expression to describe the behaviour of wavefields at the free surface, and no extra work for the free‐surface boundary condition is needed. Essentially, the proposed free‐surface expression is a modification of density and constitutive relation at the free surface. In comparison with a direct difference approximate method of the free‐surface boundary condition, this adaptive free‐surface expression can produce more accurate and stable results for a broad range of Poisson's ratio. In addition, this expression has a good performance in handling the lateral variation of Poisson's ratio adaptively and without instability.     

Poro‐acoustoelasticity of fluid‐saturated rocks

We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data.     

Gabor‐frame‐based Gaussian packet migration

We present a Gaussian packet migration method based on Gabor frame decomposition and asymptotic propagation of Gaussian packets. A Gaussian packet has both Gaussian‐shaped time–frequency localization and space–direction localization. Its evolution can be obtained by ray tracing and dynamic ray tracing. In this paper, we first briefly review the concept of Gaussian packets. After discussing how initial parameters affect the shape of a Gaussian packet, we then propose two Gabor‐frame‐based Gaussian packet decomposition methods that can sparsely and accurately represent seismic data. One method is the dreamlet–Gaussian packet method. Dreamlets are physical wavelets defined on an observation plane and can represent seismic data efficiently in the local time–frequency space–wavenumber domain. After decomposition, dreamlet coefficients can be easily converted to the corresponding Gaussian packet coefficients. The other method is the Gabor‐frame Gaussian beam method. In this method, a local slant stack, which is widely used in Gaussian beam migration, is combined with the Gabor frame decomposition to obtain uniform sampled horizontal slowness for each local frequency. Based on these decomposition methods, we derive a poststack depth migration method through the summation of the backpropagated Gaussian packets and the application of the imaging condition. To demonstrate the Gaussian packet evolution and migration/imaging in complex models, we show several numerical examples. We first use the evolution of a single Gaussian packet in media with different complexities to show the accuracy of Gaussian packet propagation. Then we test the point source responses in smoothed varying velocity models to show the accuracy of Gaussian packet summation. Finally, using poststack synthetic data sets of a four‐layer model and the two‐dimensional SEG/EAGE model, we demonstrate the validity and accuracy of the migration method. Compared with the more accurate but more time‐consuming one‐way wave‐equation‐based migration, such as beamlet migration, the Gaussian packet method proposed in this paper can correctly image the major structures of the complex model, especially in subsalt areas, with much higher efficiency. This shows the application potential of Gaussian packet migration in complicated areas.     

Evaluation of out‐of‐plane seismic performance of column‐to‐column precast concrete cladding panels in one‐storey industrial buildings

Recent earthquakes in Italy (L'Aquila 2009 and Emilia 2012) highlighted the vulnerability of precast cladding panels, typically associated with a connection system not designed to account for displacement and rotation compatibility between the panels and the supporting structure. Experimental investigations were performed in the past to investigate the in‐plane performance of cladding panels and design recommendations have been made accordingly; however, in the case of out‐of‐plane seismic loads, the load demand is commonly evaluated in the design practice by means of formulations for nonstructural components. This paper summarizes the results obtained from parametric analyses conducted to estimate the out‐of‐plane load demand in column‐to‐column cladding panels typical of one‐storey commercial and industrial buildings. Empirical equations suitable for both new and existing panels are proposed and compared with the design equations given in Eurocode 8 and ASCE 7. The paper also considers the effects of the development of plastic hinges at the column base and of the roof flexibility on the load demand in panel‐to‐column connections. The roof flexibility may generate the torsion of the panels; consequently, an analytical procedure to account for such effects is proposed. Finally, general design recommendations are made.     

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.     

Nested‐Scale Nutrient Flux in a Mixed‐Land‐Use Urbanizing Watershed

To improve quantitative understanding of mixed‐land‐use impacts on nutrient yields, a nested‐scale experimental watershed study design (n = 5) was applied in a 303(d), clean water act impaired urbanizing watershed of the lower Missouri River Basin, USA. From 2010 to 2013, water samples (n = 858 sample days per site) were analysed for total inorganic nitrogen (TIN‐N), nitrite (NO₂–N) nitrate (NO₃–N), ammonia (NH₃–N), and total phosphorus (TP‐P). Annual, seasonal, and monthly flow‐weighted concentrations (FWCs) and nutrient yields were estimated. Mean nutrient concentrations were highest where agricultural land use comprised 58% of the drainage area (NH₃ = 0.111 mg/l; NO₂ = 0.045 mg/l; NO₃ = 0.684 mg/l, TIN = 0.840 mg/l; TP = 0.127 mg/l). Average TP‐P increased by 15% with 20% increased urban land use area. Highly variable annual precipitation was observed during the study with highest nutrient yields during 2010 (record setting wet year) and lowest nutrient yields during 2012 (extreme drought year). Annual TIN‐N and TP‐P yields exceeded 10.3 and 2.04 kg ha^?1 yr^?1 from the agricultural dominated headwaters. Mean annual NH₃–N, NO₂–N, NO₃–N, TIN‐N, and TP‐P yields were 0.742, 0.400, 4.24, 5.38, and 0.979 kg ha^?1 yr^?1, respectively near the watershed outlet. Precipitation accounted for the majority of the explained variance in nutrient yields (R² values from 0.68 to 0.85). Nutrient yields were also dependent on annual precipitation of the preceding year (R² values from 0.87 to 0.91) thus enforcing the great complexity of variable mixed‐land‐use mediated source‐sink nutrient yield relationships. Study results better inform land managers and best management practices designed to mitigate nutrient pollution issues in mixed‐land‐use freshwater ecosystems. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantitative gas saturation estimation by frequency‐dependent amplitude‐versus‐offset analysis

Seismic amplitudes contain important information that can be related to fluid saturation. The amplitude‐versus‐offset analysis of seismic data based on Gassmann's theory and the approximation of the Zoeppritz equations has played a central role in reservoir characterization. However, this standard technique faces a long‐standing problem: its inability to distinguish between partial gas and “fizz‐water” with little gas saturation. In this paper, we studied seismic dispersion and attenuation in partially saturated poroelastic media by using frequency‐dependent rock physics model, through which the frequency‐dependent amplitude‐versus‐offset response is calculated as a function of porosity and water saturation. We propose a cross‐plotting of two attributes derived from the frequency‐dependent amplitude‐versus‐offset response to differentiate partial gas saturation and “fizz‐water” saturation. One of the attributes is a measure of “low frequency”, or Gassmann, of reflectivity, whereas the other is a measure of the “frequency dependence” of reflectivity. This is in contrast to standard amplitude‐versus‐offset attributes, where there is typically no such separation. A pragmatic frequency‐dependent amplitude‐versus‐offset inversion for rock and fluid properties is also established based on Bayesian theorem. A synthetic study is performed to explore the potential of the method to estimate gas saturation and porosity variations. An advantage of our work is that the method is in principle predictive, opening the way to further testing and calibration with field data. We believe that such work should guide and augment more theoretical studies of frequency‐dependent amplitude‐versus‐offset analysis.