Waves, circulation and vertical dependence

Longuet-Higgins and Stewart (J Fluid Mech 13:481–504, 1962; Deep-Sea Res 11:529–562, 1964) and later Phillips (1977) introduced the problem of waves incident on a beach, from deep to shallow water. From the wave energy equation and the vertically integrated continuity equation, they inferred velocities to be Stokes drift plus a return current so that the vertical integral of the combined velocities was nil. As a consequence, it can be shown that velocities of the order of Stokes drift rendered the advective term in the momentum equation negligible resulting in a simple balance between the horizontal gradients of the vertically integrated elevation and wave radiation stress terms; the latter was first derived by Longuet-Higgins and Stewart. Mellor (J Phys Oceanogr 33:1978–1989, 2003a), noting that vertically integrated continuity and momentum equations were not able to deal with three-dimensional numerical or analytical ocean models, derived a vertically dependent theory of wave–circulation interaction. It has since been partially revised and the revisions are reviewed here. The theory is comprised of the conventional, three-dimensional, continuity and momentum equations plus a vertically distributed, wave radiation stress term. When applied to the problem of waves incident on a beach with essentially zero turbulence momentum mixing, velocities are very large and the simple balance between elevation and radiation stress gradients no longer prevails. However, when turbulence mixing is reinstated, the vertically dependent radiation stresses produce vertical velocity gradients which then produce turbulent mixing; as a consequence, velocities are reduced, but are still larger by an order of magnitude compared to Stokes drift. Nevertheless, the velocity reduction is sufficient so that elevation set-down obtained from a balance between elevation gradient and radiation stress gradients is nearly coincident with that obtained by the aforementioned papers. This paper includes four appendices. The first appendix demonstrates the numerical process by which Stokes drift is excluded from the turbulence stress parameterization in the momentum equation. A second appendix determines a bottom slope criterion for the application of linear wave relations to the derivation of the wave radiation stress. The third appendix explores the possibility of generalizing results by non-dimensionalization. The final appendix applies the basic theory to a problem introduced by Bennis and Ardhuin (J Phys Oceanogr 41:2008–2012, 2011).     

An evaluation of 3-D velocity models of the Kanto basin for long-period ground motion simulations

We performed three-dimensional (3-D) finite difference simulations of long-period ground motions (2–10 s) in the Kanto basin using the Japan Seismic Hazard Information Station (J-SHIS 2009), Yamada and Yamanaka (Exploration Geophysics 65(3):139–150, 2012) (YY), and Head Quarter for Earthquake Research Promotion (HERP 2012) velocity models for two intermediate depth (68–80 km) moderate earthquakes (Mw 5.8–5.9), which occurred beneath the Kanto basin. The models primarily differ in the basic data set used in the construction of the velocity models. The J-SHIS and HERP models are the results of integration of mainly geological, geophysical, and earthquake data. On the other hand, the YY model is oriented towards the microtremor-array-observation data. We obtained a goodness of fit between the observed and synthetic data based on three parameters, peak ground velocities (PGVs), smoothed Fourier spectra (FFT), and cross-correlations, using an algorithm proposed by Olsen and Mayhew (Seism Res Lett 81:715–723, 2010). We found that the three models reproduced the PGVs and FFT satisfactorily at most sites. However, the models performed poorly in terms of cross-correlations especially at the basin edges. We found that the synthetics using the YY model overestimate the observed waveforms at several sites located in the areas having V _s 0.3 km/s in the top layer; on the other hand, the J-SHIS and HERP models explain the waveforms better at the sites and perform similarly at most sites. We also found that the J-SHIS and HERP models consist of thick sediments beneath some sites, where the YY model is preferable. Thus, we have concluded that the models require revisions for the reliable prediction of long-period ground motions from future large earthquakes.     

Analysis of STCC eddies using the Okubo–Weiss parameter on model and satellite data

The North Pacific Subtropical Counter Current (STCC) is a weak zonal current comprising of a weak eastward flow near the surface (with speeds of less than 0.1 m/s and a thickness of approximately 50–100 m) and westward flow (the North Equatorial Current) beneath. Previous studies (e.g., Qiu J Phys Oceanogr 29: 2471–2486, 1999) have shown that the STCC is baroclinically unstable. Therefore, despite its weak mean speeds, nonlinear STCC eddies with diameters ~300 km or larger and rotational speeds exceeding the eddy propagation speeds develop (Samelson J Phys Oceanogr 27: 2645–2662, 1997; Chelton et al. Prog Oceanogr 91: 167–216, 2011). In this study, the authors present numerical experiments to describe and explain the instability and eddy-generation processes of the STCC and the seasonal variation. Emphasis is on finite-amplitude eddies which are analyzed based on the parameter of Okubo (Deep-Sea Res 17: 445–454, 1970) and Weiss (Physica D 48: 273–294, 1991). The temperature and salinity distribution in March and April offer the favorable condition for eddies to grow, while September and October are unfavorable seasons for the generation of eddies. STCC is maintained not only by subsurface front but also by the sea surface temperature (SST) front. The seasonal variation of the vertical shear is dominated by the seasonal surface STCC velocity. The SST front enhances the instability and lead to the faster growth of STCC eddies in winter and spring. The near-surface processes are therefore crucial for the STCC system.     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $\text{ V }_\mathrm{S30}$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $\text{ M }_\mathrm{w}$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $\text{ R }_\mathrm{epi}$ , and hypocentral distance, $\text{ R }_\mathrm{hyp}$ ) and finite-fault (distance to the surface projection of the rupture, $\text{ R }_\mathrm{JB}$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.     

Ground motion simulations for the 23 October 2011 Van,Eastern Turkey earthquake using stochastic finite fault approach

The 23 October 2011 Van (Mw 7.1) earthquake that occurred in Eastern Turkey resulted in heavy damage particularly in the city of Van and town of Ercis. This paper presents ground motion simulations of Van earthquake by using stochastic finite fault method (EXSIM, Motazedian and Atkinson in Bull Seismol Soc Am 95:995–1010, 2005; Boore in Bull Seismol Soc Am 99:3202–3216, 2009) that provides a simple and effective tool to generate high frequency strong motion. The input parameters related to source, path, and site effects are calibrated on the basis of minimizing the error functions between simulations and observations both in time and frequency domain. Validated model parameters are used to produce synthetics in regional extent with the aim of understanding the level and distribution of the ground shaking particularly in the near fault region where no recordings are available within the 40 km of the epicenter. This paper evaluates the effect of two different slip models on ground motion intensity measures over the area of interest and addresses the variability in the near fault region associated with the source effect. The synthetics are compared with the corresponding estimations of ground motion prediction equations by Boore and Atkinson (Earthq Spectra 24:99–138, 2008), Akkar and Bommer (Seismol Res Lett 81:195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100:2978–2995, 2010). Our results indicate that despite the limitation of the method for incorporating the directivity effect and inadequate representation of the soil conditions at the individual stations, a satisfactory match between synthetics and observations are obtained both in time and frequency domain. Spatial distributions of the synthetics in regional level also show reasonable correlation with ground motion prediction equations and damage observations.     

Real-time risk assessment in seismic early warning and rapid response: a feasibility study in Bishkek (Kyrgyzstan)

Earthquake early warning systems (EEWS) are considered to be an effective, pragmatic, and viable tool for seismic risk reduction in cities. While standard EEWS approaches focus on the real-time estimation of an earthquake’s location and magnitude, innovative developments in EEWS include the capacity for the rapid assessment of damage. Clearly, for all public authorities that are engaged in coordinating emergency activities during and soon after earthquakes, real-time information about the potential damage distribution within a city is invaluable. In this work, we present a first attempt to design an early warning and rapid response procedure for real-time risk assessment. In particular, the procedure uses typical real-time information (i.e., P-wave arrival times and early waveforms) derived from a regional seismic network for locating and evaluating the size of an earthquake, information which in turn is exploited for extracting a risk map representing the potential distribution of damage from a dataset of predicted scenarios compiled for the target city. A feasibility study of the procedure is presented for the city of Bishkek, the capital of Kyrgyzstan, which is surrounded by the Kyrgyz seismic network by mimicking the ground motion associated with two historical events that occurred close to Bishkek, namely the 1911 Kemin (M?=?8.2; ±0.2) and the 1885 Belovodsk (M?=?6.9; ±0.5) earthquakes. Various methodologies from previous studies were considered when planning the implementation of the early warning and rapid response procedure for real-time risk assessment: the Satriano et al. (Bull Seismol Soc Am 98(3):1482–1494, 2008) approach to real-time earthquake location; the Caprio et al. (Geophys Res Lett 38:L02301, 2011) approach for estimating moment magnitude in real time; the EXSIM method for ground motion simulation (Motazedian and Atkinson, Bull Seismol Soc Am 95:995–1010, 2005); the Sokolov (Earthquake Spectra 161: 679–694, 2002) approach for estimating intensity from Fourier amplitude spectra; and the Tyagunov et al. (Nat Hazard Earth Syst Sci 6:573–586, 2006) approach for risk computation. Innovatively, all these methods are jointly applied to assess in real time the seismic risk of a particular target site, namely the city of Bishkek. Finally, the site amplification and vulnerability datasets considered in the proposed methodology are taken from previous studies, i.e., Parolai et al. (Bull Seismol Soc Am, 2010) and Bindi et al. (Soil Dyn Earthq Eng, 2011), respectively.     

Attenuation Tomography Based on Strong Motion Data: Case Study of Central Honshu Region, Japan

Three-dimensional frequency dependent S-wave quality factor (Q_β(f)) value for the central Honshu region of Japan has been determined in this paper using an algorithm based on inversion of strong motion data. The method of inversion for determination of three-dimensional attenuation coefficients is proposed by Hashida and Shimazaki (J Phys Earth. 32, 299–316, 1984) and has been used and modified by Joshi (Curr Sci. 90, 581–585, 2006; Nat Hazards. 43, 129–146, 2007) and Joshi et al. (J. Seismol. 14, 247–272, 2010). Twenty-one earthquakes digitally recorded on strong motion stations of Kik-net network have been used in this work. The magnitude of these earthquake ranges from 3.1 to 4.2 and depth ranging from 5 to 20 km, respectively. The borehole data having high signal to noise ratio and minimum site effect is used in the present work. The attenuation structure is determined by dividing the entire area into twenty-five three-dimensional blocks of uniform thickness having different frequency-dependent shear wave quality factor. Shear wave quality factor values have been determined at frequencies of 2.5, 7.0 and 10 Hz from record in a rectangular grid defined by 35.4°N to 36.4°N and 137.2°E to 138.2°E. The obtained attenuation structure is compared with the available geological features in the region and comparison shows that the obtained structure is capable of resolving important tectonic features present in the area. The proposed attenuation structure is compared with the probabilistic seismic hazard map of the region and shows that it bears some remarkable similarity in the patterns seen in seismic hazard map.     

Evaluation of nearshore wave models in steep reef environments

To provide coastal engineers and scientists with a quantitative evaluation of nearshore numerical wave models in reef environments, we review and compare three commonly used models with detailed laboratory observations. These models are the following: (1) SWASH (Simulating WAves till SHore) (Zijlema et al. 2011), a phase-resolving nonlinear shallow-water wave model with added nonhydrostatic terms; (2) SWAN (Simulating WAve Nearshore) (Booij et al. 1999), a phase-averaged spectral wave model; and (3) XBeach (Roelvink et al. 2009), a coupled phase-averaged spectral wave model (applied to modeling sea-swell waves) and a nonlinear shallow-water model (applied to modeling infragravity waves). A quantitative assessment was made of each model’s ability to predict sea-swell (SS) wave height, infragravity (IG) wave height, wave spectra, and wave setup ( \( \overline{\eta} \) ) at five locations across the laboratory fringing reef profile of Demirbilek et al. (2007). Simulations were performed with the “recommended” empirical coefficients as documented for each model, and then the key wave-breaking parameter for each model (α in SWASH and γ in both SWAN and XBeach) was optimized to most accurately reproduce the observations. SWASH, SWAN, and XBeach were found to be capable of predicting SS wave height variations across the steep fringing reef profile with reasonable accuracy using the default coefficients. Nevertheless, tuning of the key wave-breaking parameter improved the accuracy of each model’s predictions. SWASH and XBeach were also able to predict IG wave height and spectral transformation. Although SWAN was capable of modeling the SS wave height, in its current form, it was not capable of modeling the spectral transformation into lower frequencies, as evident in the underprediction of the low-frequency waves.     

An analytical method for the evaluation of the in-plan irregularity of non-regularly asymmetric buildings

This paper describes a new method for the evaluation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ of uncoupled torsional to lateral frequencies in real multi-storey buildings. The above-mentioned parameters greatly affect the lateral-to-torsional coupling of the response of asymmetric systems and thus are of paramount importance in the assessment of the in-plan irregularity of buildings. The proposed method, which is a generalization of that suggested by Calderoni et al. (Earthq Spectra 18(2):219–231, 2002), allows the calculation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ from the structural response to arbitrary distributions of forces and torsional couples. The effectiveness of the method is validated on some regularly and non-regularly asymmetric buildings characterised by different in-plan irregularity. The analyses demonstrate that the results of the method are rigorous in the case of regularly asymmetric systems and only slightly depend upon the heightwise distribution of the forces in the case of non-regularly asymmetric systems. Finally, the values of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ resulting from the proposed method are compared to those obtained by means of the procedure suggested by Makarios and Anastassiadis in (Struct Des Tall Spec Build 7(1):33–55, 1998a; Struct Des Tall Spec Build 7(1):57–71, 1998b) .     

On the selection of GMPEs for Vrancea subcrustal seismic source

The Vrancea subcrustal seismic source is characterized by large magnitude ( $M_{W} \ge 7$ ) intermediate-depth earthquakes that occur two or three times during a century on average. In this study several procedures are used to grade four candidate ground motion prediction equations proposed for Vrancea source in the SHARE project. In the work of Delavaud et al. (J Seismol 16(3):451–473, 2012) four ground motion prediction models developed for subduction zones (Zhao et al. in Bull Seism Soc Am 96(3):898–913, 2006; Atkinson and Boore in Bull Seism Soc Am 93(4):1703–1729, 2003; Youngs et al. in Seism Res Lett 68(1):58–73, 1997; Lin and Lee in Bull Seism Soc Am 98(1):220–240, 2008) are suggested as suitable for Vrancea subcrustal seismic source. The paper presents the appropriateness analysis of the four suggested ground motion prediction equations done using a dataset of 109 triaxial accelerograms recorded during seven Vrancea seismic events with moment magnitude $M_{W}$ between 5.4 and 7.4, occurred in the past 35 years. The strong ground motions were recorded in Romania, as well as in Bulgaria, Republic of Moldova and Serbia. Based on the ground motion dataset several goodness-of-fit measures are used in order to quantify how well the selected models match with the recorded data. The compatibility of the four ground motion prediction models with respect to magnitude scaling and distance scaling implied by strong ground motion dataset is investigated as well. The analyses show that the Youngs et al. (Seism Res Lett 68(1):58–73, 1997) and Zhao et al. (Bull Seism Soc Am 96(3):898–913, 2006) ground motion prediction models have a better fit with the data and can be candidate models for Probabilistic Seismic Hazard Assessment.     

Scaling Transition in Earthquake Sources: A Possible Link Between Seismic and Laboratory Measurements

We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M _W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M _W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M _W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M _W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M _W ~ 5.5. A representative fault slip for M _W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M _W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.     

Compatible ground-motion prediction equations for damping scaling factors and vertical-to-horizontal spectral amplitude ratios for the broader Europe region

In a companion article Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4, 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) present a new ground-motion prediction equation (GMPE) for estimating 5 %-damped horizontal pseudo-acceleration spectral (PSA) ordinates for shallow active crustal regions in Europe and the Middle East. This study provides a supplementary viscous damping model to modify 5 %-damped horizontal spectral ordinates of Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) for damping ratios ranging from 1 to 50 %. The paper also presents another damping model for scaling 5 %-damped vertical spectral ordinates that can be estimated from the vertical-to-horizontal (V/H) spectral ratio GMPE that is also developed within the context of this study. For consistency in engineering applications, the horizontal and vertical damping models cover the same damping ratios as noted above. The article concludes by introducing period-dependent correlation coefficients to compute horizontal and vertical conditional mean spectra (Baker in J Struct Eng 137:322–331, 2011). The applicability range of the presented models is the same as of the horizontal GMPE proposed by Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b): as for spectral periods $0.01 \hbox { s}\le \,\hbox {T}\le \,4\hbox { s}$ as well as PGA and PGV for V/H model; and in terms of seismological estimator parameters $4\le \hbox {M}_\mathrm{w} \le 8, \hbox { R} \le 200 \hbox { km}, 150\hbox { m/s}\le \hbox { V}_\mathrm{S30}\le $ 1,200 m/s, for reverse, normal and strike-slip faults. The source-to-site distance measures that can be used in the computations are epicentral $(\hbox {R}_\mathrm{epi})$ , hypocentral $(\hbox {R}_\mathrm{hyp})$ and Joyner–Boore $(\hbox {R}_\mathrm{JB})$ distances. The implementation of the proposed GMPEs will facilitate site-specific adjustments of the spectral amplitudes predicted from probabilistic seismic hazard assessment in Europe and the Middle East region. They can also help expressing the site-specific design ground motion in several formats. The consistency of the proposed models together with the Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) GMPE may be advantageous for future modifications in the ground-motion definition in Eurocode 8 (CEN in Eurocode 8, Design of structures for earthquake resistance—part 1: general rules, seismic actions and rules for buildings. European Standard NF EN 1998-1, Brussels, 2004).     

An example of solute spreading in nonstationary, bounded geological formations

We analyze the impact of a linear trend in the mean log-conductivity on the transport of a conservative tracer in a bounded domain. The effects of such a linear trend on solute transport were analyzed in depth for unbounded domains (Rubin and Seong, Water Resour Res 30(11):2901–2911, 1994; Indelman and Rubin, Water Resour Res 31(5):1257–1265, 1995; Water Resour Res 32(5):1257–1265, 1996), whereas studies concerning this special case of medium nonstationarity in finite domains usually focus on head or flow statistics (Guadagnini et al., Stoch Environ Res Risk Assess, 17:394–407, 2003). In this study both ensemble and effective plume moments are provided for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis involving a steady velocity field spatially nonstationary is developed by using the stochastic finite element method. Results show that ensemble moments are affected by increasing trends both parallel and normal to the mean flow direction, but the impact on effective plume moments is very different. A parallel trend does not seem to influence the effective second moments; while a normal trend, although modifies the transverse effective moment only weakly, strongly increases the longitudinal one, especially for large initial sizes of the source. Furthermore, the increase of the particle displacement variance produced by a parallel trend in the finite domain disagrees with the results obtained in an unbounded domain, due to the boundary conditions here considered making both head and velocity moments nonstationary and nonsymmetric.     

A Tailored Computation of the Mean Dynamic Topography for a Consistent Integration into Ocean Circulation Models

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography—the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99–110, 2012; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie 2012). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.     

Focal Mechanism Solutions of the 2008 Wenchuan earthquake sequence from P-wave polarities and SH/P amplitude ratios: new results and implications

The 2008 Wenchuan earthquake, a major intraplate earthquake with M _w 7.9, occurred on the slowly deforming Longmenshan fault. To better understand the causes of this devastating earthquake, we need knowledge of the regional stress field and the underlying geodynamic processes. Here, we determine focal mechanism solutions (FMSs) of the 2008 Wenchuan earthquake sequence (WES) using both P-wave first-motion polarity data and SH/P amplitude ratio (AR) data. As P-wave polarities are more reliable information, they are given priority over SH/P AR, the latter of which are used only when the former has loose constraint on the FMSs. We collect data from three categories: (1) permanent stations deployed by the China Earthquake Administration (CEA); (2) the Western Sichuan Passive Seismic Array (WSPSA) deployed by Institute of Geology, CEA; (3) global stations from Incorporated Research Institutions for Seismology. Finally, 129 events with magnitude over M _s 4.0 in the 2008 WES are identified to have well-constrained FMSs. Among them, 83 are well constrained by P-wave polarities only as shown by Cai et al. (Earthq Sci 24(1):115–125, 2011), and the rest of which are newly constrained by incorporating SH/P AR. Based on the spatial distribution and FMSs of the WES, we draw following conclusions: (1) the principle compressional directions of most FMSs of the WES are subhorizontal, generally in agreement with the conclusion given by Cai et al. (2011) but with a few modifications that the compressional directions are WNW–ESE around Wenchuan and ENE–WSW around Qingchuan, respectively. The subhorizontal compressional direction along the Longmenshan fault from SW to NE seems to have a left-lateral rotation, which agrees well with regional stress field inverted by former researchers (e.g., Xu et al., Acta Seismol Sin 30(5), 1987; Acta Geophys Sin 32(6), 1989; Cui et al., Seismol Geol 27(2):234–242, 2005); (2) the FMSs of the events not only reflected the regional stress state of the Longmenshan region, but also were obviously controlled by the faults to some extent, which was pointed out by Cai et al. (2011) and Yi et al. (Chin J Geophys 55(4):1213–1227, 2012); (3) while the 2008 Wenchuan earthquake and some of its strong aftershocks released most of the elastic energy accumulated on the Longmenshan fault, some other aftershocks seem to occur just for releasing the elastic energy promptly created by the 2008 Wenchuan earthquake and some of its strong aftershocks. (4) Our results further suggest that the Longmenshan fault from Wenchuan to Beichuan was nearly fully destroyed by the 2008 Wenchuan earthquake and accordingly propose that there is less probability for great earthquakes in the middle part of the Longmenshan fault in the near future, although there might be a barrier to the southwest of Wenchuan and it is needed to pay some attention on it in the near future.     

Ergodicity and Earthquake Catalogs: Forecast Testing and Resulting Implications

Recently the equilibrium property of ergodicity was identified in an earthquake fault system (Tiampo et al., Phys. Rev. Lett. 91, 238501, 2003; Phys. Rev. E 75, 066107, 2007). Ergodicity in this context not only requires that the system is stationary for these networks at the applicable spatial and temporal scales, but also implies that they are in a state of metastable equilibrium, one in which the ensemble averages can be substituted for temporal averages when studying their behavior in space and time. In this work we show that this property can be used to identify those regions of parameter space which are stationary when applied to the seismicity of two naturally-occurring earthquake fault networks. We apply this measure to one particular seismicity-based forecasting tool, the Pattern Informatics index (Tiampo et al., Europhys. Lett. 60, 481–487, 2002; Rundle et al., Proc. National Acad. Sci., U.S.A., Suppl. 1, 99, 2463, 2002), in order to test the hypothesis that the identification of ergodic regions can be used to improve and optimize forecasts that rely on historic seismicity catalogs. We also apply the same measure to synthetic catalogs in order to better understand the physical process that affects this accuracy. We show that, in particular, ergodic regions defined by magnitude and time period provide more reliable forecasts of future events in both natural and synthetic catalogs, and that these improvements can be directly related to specific features or properties of the catalogs that impact the behavior of their spatial and temporal statistics.     

Electrical resistivity tomography study of Taal volcano hydrothermal system, Philippines

Taal volcano (311?m in altitude) is located in The Philippines (14°N, 121°E) and since 1572 has erupted 33 times, causing more than 2,000 casualties during the most violent eruptions. In March 2010, the shallow structures in areas where present-day surface activity takes place were investigated by DC resistivity surveys. Electrical resistivity tomography (ERT) lines were performed above the two identified hydrothermal areas located on the northern flank of the volcano and in the Main Crater, respectively. Due to rough topography, deep valleys, and dense vegetation, most measurements were collected using a remote method based on a laboratory-made equipment. This allowed retrieval of information down to a depth of 250?m. ERTs results detail the outlines of the two geothermal fields defined by previous self-potential, CO₂ soil degassing, ground temperature, and magnetic mapping (Harada et al. Japan Acad Sci 81:261–266, 2005; Zlotnicki et al. Bull Volcanol 71:29–49, 2009a, Phys Chem Earth 34:294–408, 2009b). Hydrothermal fluids originate mainly from inside the northern part of the Main Crater at a depth greater than the bottom of the Crater Lake, and flow upward to the ground surface. Furthermore, water from the Main Crater Lake infiltrates inside the surrounding geological formations. The hydrothermal fluids, outlined by gas releases and high temperatures, cross the crater rim and interact with the northern geothermal field located outside the Main Crater.     

Effects of Location Errors in Pattern Informatics

The effect of location errors in the performance of seismicity-based forecasting methods was studied here using one particular binary forecast technique, the Pattern Informatics (PI) technique (Rundle et al., Proc Nat Acad Sci USA 99, 2514–2521, 2002; Tiampo et al., Pure Appl Geophys 159, 2429–2467, 2002). The Southern Californian dataset was used to generate a series of perturbed catalogs by adding different levels of noise to epicenter locations. The PI technique was applied to these perturbed datasets to perform retrospective forecasts that were evaluated by means of skill scores, commonly used in atmospheric sciences. These results were then compared to the effectiveness obtained from the original dataset. Isolated instances of decline of the PI performance were observed due to the nature of the skill scores themselves, but no clear trend of degradation was identified. Dependence on the total number of events in a catalog also was studied, with no systematic degradation in the performance of the PI for catalogs with events in the cases studied. These results suggest that the stability of the PI method is due to the invariance of the clustering patterns identified by the TM metric (Thirumalai and Mountain, Phys Rev A 39, 3563–3573, 1989) when applied to seismicity.     

Indian summer monsoon rainfall prediction using artificial neural network

Forecasting the monsoon temporally is a major scientific issue in the field of monsoon meteorology. The ensemble of statistics and mathematics has increased the accuracy of forecasting of Indian summer monsoon rainfall (ISMR) up to some extent. But due to the nonlinear nature of ISMR, its forecasting accuracy is still below the satisfactory level. Mathematical and statistical models require complex computing power. Therefore, many researchers have paid attention to apply artificial neural network in ISMR forecasting. In this study, we have used feed-forward back-propagation neural network algorithm for ISMR forecasting. Based on this algorithm, we have proposed the five neural network architectures designated as BP1, BP2, $\ldots, $ … , BP5 using three layers of neurons (one input layer, one hidden layer and one output layer). Detail architecture of the neural networks is provided in this article. Time series data set of ISMR is obtained from Pathasarathy et al. (Theor Appl Climatol 49:217–224 1994) (1871–1994) and IITM (http://www.tropmet.res.in/, 2012) (1995–2010) for the period 1871–2010, for the months of June, July, August and September individually, and for the monsoon season (sum of June, July, August and September). The data set is trained and tested separately for each of the neural network architecture, viz., BP1–BP5. The forecasted results obtained for the training and testing data are then compared with existing model. Results clearly exhibit superiority of our model over the considered existing model. The seasonal rainfall values over India for next 5 years have also been predicted.