On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed: Part II-Observations

An analytical expression for the 10 m drag law in terms of the 10 m wind speed at the maximum in the 10 m drag coefficient, and the Charnock constant is presented, which is based on the results obtained from a model of the air-sea interface derived in Bye et al. (2010). This drag law is almost independent of wave age and over the mid-range of wind speeds (5?17 ms^?1) is very similar to the drag law based on observed data presented in Foreman and Emeis (2010). The linear fit of the observed data which incorporates a constant into the traditional definition of the drag coefficient is shown to arise to first-order as a consequence of the momentum exchange across the air-sea boundary layer brought about by wave generation and spray production which are explicitly represented in the theoretical model.     

Extreme storm surge modelling in the North Sea

This study shows that storm surge model performance in the North Sea is mostly unaffected by the application of temporal variations of surface drag due to changes in sea state provided the choice of a suitable constant Charnock parameter in the sea-state-independent case. Including essential meteorological features on smaller scales and minimising interpolation errors by increasing forcing data resolution are shown to be more important for the improvement of model performance particularly at the high tail of the probability distribution. This is found in a modelling study using WAQUA/DCSMv5 by evaluating the influence of a realistic air-sea momentum transfer parameterization and comparing it to the influence of changes in the spatial and temporal resolution of the applied forcing fields in an effort to support the improvement of impact and climate analysis studies. Particular attention is given to the representation of extreme water levels over the past decades based on the example of the Netherlands. For this, WAQUA/DCSMv5 is forced with ERA-Interim reanalysis data. Model results are obtained from a set of different forcing fields, which either (i) include a wave-state-dependent Charnock parameter or (ii) apply a constant Charnock parameter (α _{C h} =?0.032) tuned for young sea states in the North Sea, but differ in their spatial and/or temporal resolution. Increasing forcing field resolution from roughly 79 to 12 km through dynamically downscaling can reduce the modelled low bias, depending on coastal station, by up to 0.25 m for the modelled extreme water levels with a 1-year return period and between 0.1 m and 0.5 m for extreme surge heights.     

Summer surface layer thermal response to surface gravity waves in the Yellow Sea

The Princeton Ocean Model (POM) with generalized coordinate system (POMgcs) is used to study the summer surface-layer thermal response to surface gravity waves in the Yellow Sea (YS). The parameterization schemes of wave breaking developed by Mellor and Blumberg (J Phys Oceanogr 34:693–698, 2004) and Kantha and Clayson (Ocean Model 6:101–124, 2004), respectively, and Stokes production developed by Kantha and Clayson (Ocean Model 6:101–124, 2004) are both included in the Mellor–Yamada turbulence closure model Mellor and Yamada (Rev Geophys 20:851–875, 1982) of POMgcs. Numerical results show that surface gravity waves impact the depth of surface mixed layer of temperature in the YS in summer. The surface mixed layer in the YS cannot be reproduced well and has a visible difference from the observation if the parameterization schemes are not included. A diagnostic analysis of turbulent kinetic energy suggests that both Stokes production and wave breaking play key roles in enhancing the turbulent mixing near the sea surface in the YS. Stokes production seems to have a greater impact throughout the upper mixed layer in the YS in summer than that of wave breaking. In addition, a diagnostic analysis of the momentum balance shows that Coriolis–Stokes forcing has a significant effect on the momentum budget in the upper layer in the YS, and surface gravity waves are able to reduce the velocity of mean flow near the surface and make the mean flow near the surface more homogeneous vertically in the YS.     

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 Elliptical Model for Deformation Due to Groundwater Fluctuations

Historically, surface subsidence as a result of subsurface groundwater fluctuations have produced important and, at times, catastrophic effects, whether natural or anthropogenic. Over the past 30?years, numerical and analytical techniques for the modeling of this surface deformation, based upon elastic and poroelastic theory, have been remarkably successful in predicting the magnitude of that deformation (Le Mouélic and Adragna in Geophys Res Lett 29:1853, 2002). In this work we have extended the formula for a circular-shaped aquifer (Geertsma in J Petroleum Tech 25:734–744, 1973) to a more realistic elliptical shape. We have improved the accuracy of the approximation by making use of the cross terms of the expansion for the elliptic coordinates in terms of the eccentricity, e, and the mean anomaly angle, M, widely used in astronomy. Results of a number of simulations, in terms of e and M developed from the transcendental Kepler equation, are encouraging, giving realistic values for the elliptical approximation of the vertical deformation due to groundwater change. Finally, we have applied the algorithm to modeling of groundwater in southern California.     

The statistics of wave height and crest elevation during the December 2012 storm in the North Sea

This paper concerns new field measurements of wave height and crest elevation probability distributions as measured in the North Sea during a storm in December 2012. The water surface elevation was recorded by Saab WaveRadar REX instruments mounted on eight fixed-jacket platforms in addition to a Datawell Directional Waverider buoy. The storm generated an easterly sea state which peaked well in excess of the 100-year wave height for that direction in the region. Furthermore, 19 freak waves occurred during the storm according to the definition as reported by Haver (2000). The present study demonstrates that the significant steepness and spectral bandwidth during the storm remain almost constant. Consequently, there is little change in the commonly applied design wave height and crest elevation probability distributions throughout the storm. Whilst the bulk of the recorded data was in good agreement with the theoretical distributions, it was demonstrated that when the wind speed was larger than 25 m/s, the measured crest elevation lies above the second-order Forristall distribution.     

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.     

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.     

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.     

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.     

Coupled wave and surge modelling for the eastern Irish Sea and implications for model wind-stress

We revisit the surge of November 1977, a storm event which caused damage on the Sefton coast in NW England. A hindcast has been made with a coupled surge-tide-wave model, to investigate whether a wave-dependent surface drag is necessary for accurate surge prediction, and also if this can be represented by an optimised Charnock parameter. The Proudman Oceanographic Laboratory Coastal Modelling System-Wave Model (POLCOMS-WAM) has been used to model combined tides, surges, waves and wave-current interaction in the Irish Sea on a 1.85 km grid. This period has been previously thoroughly studied, e.g. Jones and Davies [Jones, J.E., Davies, A.M., 1998. Storm surge computations for the Irish Sea using a three-dimensional numerical model including wave-current interaction. Continental Shelf Research 18(2), 201–251] and we build upon this previous work to validate the POLCOMS-WAM model to test the accuracy of surge elevation predictions in the study area. A one-way nested approach has been set up from larger scale models to the Irish Sea model. It was demonstrated that (as expected) swell from the North Atlantic does not have a significant impact in the eastern Irish Sea. To capture the external surge generated outside of the Irish Sea a (1/9° by 1/6°) model extending beyond the continental shelf edge was run using the POLCOMS model for tide and surge.     

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.     

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.     

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.     

Locations and magnitudes of earthquakes in Central Asia from seismic intensity data

We apply the Bakun and Wentworth (Bull Seism Soc Am 87:1502–1521, 1997) method to determine the location and magnitude of earthquakes occurred in Central Asia using MSK-64 intensity assignments. The attenuation model previously derived and validated by Bindi et al. (Geophys J Int, 2013) is used to analyse 21 earthquakes that occurred over the period 1885–1964, and the estimated locations and magnitudes are compared to values available in literature. Bootstrap analyses are performed to estimate the confidence intervals of the intensity magnitudes, as well as to quantify the location uncertainty. The analyses of seven significant earthquakes for the hazard assessment are presented in detail, including three large historical earthquakes that struck the northern Tien-Shan between the end of the nineteenth and the beginning of the twentieth centuries: the 1887, M 7.3 Verny, the 1889, M 8.3 Chilik and the 1911, M 8.2 Kemin earthquakes. Regarding the 1911, Kemin earthquake the magnitude values estimated from intensity data are lower (i.e. MI_LH?=?7.8 and MI_W?=?7.6 considering surface wave and moment magnitude, respectively) than the value M?=?8.2 listed in the considered catalog. These values are more in agreement with the value M _S?=?7.8 revised by Abe and Noguchi (Phys Earth Planet In, 33:1–11, 1983b) for the surface wave magnitude. For the Kemin earthquake, the distribution of the bootstrap solutions for the intensity centre reveal two minima, indicating that the distribution of intensity assignments do not constrain a unique solution. This is in agreement with the complex source rupture history of the Kemin earthquake, which involved several fault segments with different strike orientations, dipping angles and focal mechanisms (e.g. Delvaux et al. in Russ Geol Geophys 42:1167–1177, 2001; Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). Two possible locations for the intensity centre are obtained. The first is located on the easternmost sub-faults (i.e. the Aksu and Chon-Aksu segments), where most of the seismic moment was released (Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). The second location is located on the westernmost sub-faults (i.e. the Dzhil'-Aryk segment), close to the intensity centre location obtained for the 1938, M 6.9 Chu-Kemin earthquake (MI_LH?=?6.9 and MI_W?=?6.8).     

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).     

Pre-Observational Evolution of Surface Temperature in Romania as Inferred from Borehole Temperature Measurements

Temperature data from nine boreholes in the Carpathian orogen in Romania were used to obtain information on the ground surface temperature history (GSTH) in the last 250?years. The temperature measurements were taken with a thermistor probe (sensitivity in the 10 mK range) using the stop-and-go technique, at 10 m intervals, in the depth range of 20–580?m. The least squares inverse modelling approach of Tarantola and Valette (J Geophys 50:159–170, 1982) was used to infer the GSTH. Long-term air temperature records available from the Romanian weather station network were used as a comparison term for the first 100–150?years of the GSTH, and as a forcing function in a POM-SAT model that combines borehole temperature profiles (BTPs) and meteorological time series (surface air temperature, SAT) to produce information on the so-called pre-observational mean (POM). Results from a global circulation model for the Romanian area are incorporated in the discussion as well.     

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.     

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.