Exposure dating of precariously balanced rocks

Precariously balanced rocks (PBRs) are freestanding boulders that are precarious or fragile in the sense that they could be toppled by relatively low-amplitude earthquake ground motion. They are important in paleoseismology because their continued existence limits the amplitude of ground motion experienced at their location during their lifetime. In order to make quantitative use of PBRs for seismic hazard studies, one must determine when they attained their present state of fragility, that is, the point in time when the contact between the rocks and the pedestals on which they rest was exhumed from surrounding soil and the rock became vulnerable to earthquake ground motions. Cosmogenic-nuclide exposure dating can be used for this purpose, but is complicated because nuclide production occurs throughout exhumation of the PBR, so the apparent exposure age of any part of the rock surface exceeds the time that the rock has actually been precariously balanced. Here we describe a method for determining the length of time that a PBR has been fragile by measuring cosmogenic-nuclide concentrations at several locations on the PBR surface, and linking them together with a forward model that accounts for nuclide production before, during, and after exhumation of the PBR. Fitting model to data yields the rate and timing of rock exhumation and thus the length of time the rock has been fragile. We use this method to show that an example PBR in southern California has been fragile for 18.7 ± 2.8 ka.     

A complete and easily accessible means of calculating surface exposure ages or erosion rates from Be and Al measurements

We codify previously published means of calculating exposure ages and erosion rates from ¹⁰Be and ²⁶Al concentrations in rock surfaces, and present a single complete and straightforward method that reflects currently accepted practices and is consistent with existing production rate calibration measurements. It is intended to enable geoscientists, who wish to use cosmogenic-nuclide exposure age or erosion rate measurements in their work to: (a) calculate exposure ages and erosion rates; (b) compare previously published exposure ages or erosion rate measurements on a common basis; (c) evaluate the sensitivity of their results to differences between published production rate scaling schemes. The method is available online at http://hess.ess.washington.edu.     

Investigation of pressure variations over stepped spillways using smooth particle hydrodynamics

In this paper the smoothed particle hydrodynamics, (SPH), technique is used to investigate the pressure distribution on steps located in the non-aerated flow region of a stepped spillway for different discharges typical of skimming flow conditions. The open source code 2D SPHysics has been employed after being validated against the laboratory model studies of flow over broad crested weirs and flow over stepped spillways. The numerical results, in terms of the water surface and velocity profiles at different sections, are in good agreement with the corresponding experimental results. The code is then applied to determine the pressure distribution on the vertical and horizontal step faces. Also, the aspects of the pressure pattern are described and the positions/magnitudes of the maximum and minimum pressure values are presented.     

Parametric investigation of lateral spreading of gently sloping liquefied ground

This paper investigates the main parameters affecting the anticipated maximum surface displacements due to earthquake-induced lateral spreading of mildly sloping ground. The main tool used for this purpose is a numerical methodology employing a bounding surface plasticity model implemented in a finite difference code, which has been thoroughly validated against 16 published centrifuge lateral spreading experiments. This study shows that important problem parameters are the mean ground (surface) acceleration, the duration of strong shaking following the onset of liquefaction, the corrected SPT blowcount, the depth to the sliding plane, the inclination of the ground surface and the fines content of the liquefied soil layers. A new approximate multi-variable relation is proposed for the estimation of ground surface displacements due to lateral spreading in gently sloping ground, which includes the foregoing parameters. The form of the relation builds upon sliding block theory, but its final formulation is based on statistical analysis of the input data and the results from 120 parametric analyses performed with the validated numerical methodology. Comparison of the predictions of the proposed relation for ground surface displacement against pertinent field data (from 256 case histories) and centrifuge test measurements shows satisfactory accuracy. Furthermore, the variation of lateral displacements with depth is explored and distinct displacement patterns are proposed for uniform, 2-layer and 4-layer ground profiles.     

复杂场源形态的海洋可控源电磁三维正演

在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响.     

Intercomparison of 3D pore-scale flow and solute transport simulation methods

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.     

Regional beryllium-10 production rate calibration for late-glacial northeastern North America

The major uncertainty in relating cosmogenic-nuclide exposure ages to ages measured by other dating methods comes from extrapolating nuclide production rates measured at globally scattered calibration sites to the sites of unknown age that are to be dated. This uncertainty can be reduced by locating production rate calibration sites that are similar in location and age to the sites to be dated. We use this strategy to reconcile exposure age and radiocarbon deglaciation chronologies for northeastern North America by compiling ¹⁰Be production rate calibration measurements from independently dated late-glacial and early Holocene ice-marginal landforms in this region. ¹⁰Be production rates measured at these sites are 6–12% lower than predicted by the commonly accepted global ¹⁰Be calibration data set used with any published production rate scaling scheme. In addition, the regional calibration data set shows significantly less internal scatter than the global calibration data set. Thus, this calibration data set can be used to improve both the precision and accuracy of exposure dating of regional late-glacial events. For example, if the global calibration data set is used to calculate exposure ages, the exposure-age deglaciation chronology for central New England is inconsistent with the deglaciation chronology inferred from radiocarbon dating and varve stratigraphy. We show that using the regional data set instead makes the exposure age and radiocarbon chronologies consistent. This increases confidence in correlating exposure ages of ice-marginal landforms in northeastern North America with glacial and climate events dated by other means.     

3D finite volume groundwater and heat transport modeling with non-orthogonal grids,using a coordinate transformation method

Many popular groundwater modeling codes are based on the finite differences or finite volume method for orthogonal grids. In cases of complex subsurface geometries this type of grid either leads to coarse geometric representations or to extremely fine meshes. We use a coordinate transformation method (CTM) to circumvent this shortcoming. In computational fluid dynamics (CFD), this method has been applied successfully to the general Navier–Stokes equation. The method is based on tensor analysis and performs a transformation of a curvilinear into a rectangular unit grid, on which a modified formulation of the differential equations is applied. Therefore, it is not necessary to reformulate the code in total. We applied the CTM to an existing three-dimensional code (SHEMAT), a simulator for heat conduction and advection in porous media. The finite volume discretization scheme for the non-orthogonal, structured, hexahedral grid leads to a 19-point stencil and a correspondingly banded system matrix. The implementation is straightforward and it is possible to use some existing routines without modification. The accuracy of the modified code is demonstrated for single phase flow on a two-dimensional analytical solution for flow and heat transport. Additionally, a simple case of potential flow is shown for a two-dimensional grid which is increasingly deformed. The result reveals that the corresponding error increases only slightly. Finally, a thermal free-convection benchmark is discussed. The result shows, that the solution obtained with the new code is in good agreement with the ones obtained by other codes.     

Where now? Reflections on future directions for cosmogenic nuclide research from the CRONUS Projects

The CRONUS-Earth and CRONUS-EU Projects have recently examined the reproducibility of cosmogenic nuclide measurements and the ability of cosmogenic-production models to fit high-quality global calibration data sets. In both cases, although results are adequate for present geochronological needs, they must be improved to meet future demands. We recommend that the cosmogenic-nuclide community embark on a continuing effort to understand the sources of the remaining model-versus-data discrepancy, and that in the meanwhile all cosmogenic ages be benchmarked against reference ages calculated in a consistent fashion. We also recommend that interlaboratory reference materials be routinely analyzed along with samples and blanks, that the results be tracked, and that methods of achieving consistency and precision in sample preparation and analysis be shared within the community. We also recommend additional organization of the community to facilitate the use of intercomparison materials, foster analytical improvements, oversee the evolution of the community age calculator, and synthesize ongoing calibration efforts.     

Seismic behavior of a partially grouted reinforced masonry structure: Shake-table testing and numerical analyses

In regions of low to moderate seismicity in North America, reinforced masonry structures are mostly partially grouted. The behavior of such structures under lateral seismic loads is complicated because of the interaction of the grouted and ungrouted masonry. As revealed in past experimental studies, the performance of partially grouted masonry (PGM) walls under in-plane cyclic lateral loading is inferior to that of fully grouted walls. However, the dynamic behavior of a PGM wall system under severe seismic loads is not well understood. In this study, a full-scale, one-story, PGM building designed for a moderate seismic zone according to current code provisions was tested on a shake table. It was shown that the structure was able to develop an adequate base shear capacity and withstand two earthquake motions that had an effective intensity of two times the maximum considered earthquake with only moderate cracking in mortar joints. However, the structure eventually failed in a brittle manner in a subsequent motion that had a slightly lower effective intensity. A detailed finite element model of the test structure has been developed and validated. The model has been used to understand the distribution of the lateral force resistance among the wall components and to evaluate the shear-strength equation given in the design code. The code equation has been found to be adequate for this structure. Furthermore, a parametric study conducted with the finite element model has shown that the introduction of a continuous bond beam right below a window opening is highly beneficial.     

Extracting palaeoflood data from coarse‐grained Pleistocene river terrace archives: an example from SE Spain

Field‐based palaeoflood event reconstruction has the potential to contribute to the development of our understanding of long‐term landscape evolution. However, the reconstruction of past flow event histories (magnitude and frequency) over long‐term (Quaternary) timescales is fraught with difficulties. Here we make a preliminary exploration of some of the practicalities of flood reconstruction from fluvial terrace archives using commonly available sedimentological and geomorphological observations from a field perspective. We utilize Manning and palaeostage indicators to reconstruct historic events that can be used as benchmarks for a lesser used competence based approach, which is applied to coarse‐grained strath terrace deposits. We evaluate the results against gauged records for extreme and catastrophic events that affected the same region in 1973 and 2012. The findings suggest that the competence approach is most effectively applied to terrace deposits if the channel geometry is taken into account when sampling both in cross‐section and in longitudinal section and calibrated against the sedimentology for palaeo‐flow depth. Problems can arise where constrictive channel geometries allow boulder jams to develop, acting as sediment traps for the coarsest material and leading to downstream ‘boulder starvation’. Useful sites to target for palaeoflood reconstruction, therefore, would be upstream of such constrictive reaches where the coarsest transportable bedload has been effectively trapped. Sites to avoid would be downflow, where the deposited material would poorly represent palaeoflood competence. Underestimation from maximum boulder preservation and limited section exposure issues would appear to outweigh possible overestimation concerns related to fluid density and unsteady flow characteristics such as instantaneous acceleration forces. Flood data derived from river terrace deposits suggests that basal terrace geometries and coarse boulder lags common to many terrace sequences are likely the result of extreme flow events which are subsequently filled by lesser magnitude flood events, in this environmental setting. Copyright © 2016 John Wiley & Sons, Ltd.     

Analytical prediction of experimental building pounding

Valuable insights on the problem of seismic pounding have been obtained recently from analytical studies. So far, the proposed analytical models have not been validated experimentally. This paper presents the results of shake table tests of pounding between adjacent three- and eight-storey single-bay steel framed model structures. The pounding response of the frames was measured for various earthquake intensities and initial separations. The experimental results were compared to the predictions resulting from two existing pounding analysis programs. The solution strategy of the first program, SLAM-2, is based on a modal superposition technique. The second program, PC-ANSR, is a non-linear time-step analysis code in which an elastic gap element has been included. Modelling the pounding effect by elastic gap elements in the two programs produced accurate displacement and impact force results. Amplitudes of short acceleration pulses were not well predicted, however, for practical time-step increments. Relative rotations between adjacent floors induced grinding contacts which cannot be captured by uni-axial gap elements.     

Probabilistic seismic loss estimation via endurance time method

Probabilistic Seismic Loss Estimation is a methodology used as a quantitative and explicit expression of the performance of buildings using terms that address the interests of both owners and insurance companies. Applying the ATC 58 approach for seismic loss assessment of buildings requires using Incremental Dynamic Analysis (IDA), which needs hundreds of time-consuming analyses, which in turn hinders its wide application. The Endurance Time Method (ETM) is proposed herein as part of a demand propagation prediction procedure and is shown to be an economical alternative to IDA. Various scenarios were considered to achieve this purpose and their appropriateness has been evaluated using statistical methods. The most precise and efficient scenario was validated through comparison against IDA driven response predictions of 34 code conforming benchmark structures and was proven to be sufficiently precise while offering a great deal of efficiency. The loss values were estimated by replacing IDA with the proposed ETM-based procedure in the ATC 58 procedure and it was found that these values suffer from varying inaccuracies, which were attributed to the discretized nature of damage and loss prediction functions provided by ATC 58.     

Combination of FEM and CMA‐ES algorithm for transmissivity identification in aquifer systems

In this paper, we propose a coupling of a finite element model with a metaheuristic optimization algorithm for solving the inverse problem in groundwater flow (Darcy's equations). This coupling performed in 2 phases is based on the combination of 2 codes: This is the HySubF‐FEM code (hydrodynamic of subsurface flow by finite element method) used for the first phase allowing the calculation of the flow and the CMA‐ES code (covariance matrix adaptation evolution strategy) adopted in the second phase for the optimization process. The combination of these 2 codes was implemented to identify the transmissivity field of groundwater by knowing the hydraulic head in some point of the studied domain. The integrated optimization algorithm HySubF‐FEM/CMA‐ES has been validated successfully on a schematic case offering an analytical solution. As realistic application, the integrated optimization algorithm HySubF‐FEM/CMA‐ES was applied to a complex groundwater in the north of France to identify the transmissivity field. This application does not use zonation techniques but solves an optimization problem at each internal node of the mesh. The obtained results are considered excellent with high accuracy and fully consistent with the hydrogeological characteristics of the studied aquifer.However, the various numerical simulations performed in this paper have shown that the CMA‐ES algorithm is time‐consuming. Finally, the paper concludes that the proposed algorithm can be considered as an efficient tool for solving inverse problems in groundwater flow.     

Modelling fluxes of water and sediment between Venice Lagoon and the sea

To describe the exchange of water and sediment through the Venice Lagoon inlets a 3-D hydrodynamic and sediment transport model has been developed and applied to a domain comprising Venice Lagoon and a part of the Adriatic Sea. The model has been validated for both current velocities and suspended particle concentration against direct observations and from observations empirically derived fluxes from upward-looking acoustic Doppler current profiler probes installed inside each inlet. The model provides estimates of the suspended sediment transport in the lower 3 m of the water column that is not detected by acoustic Doppler current profiler sensors. The bedload model prediction has been validated against measured sand transport rates collected by sand traps deployed in the Lido and Chioggia inlets. Results indicate that, in the Lido inlet, 87% of the total load is in suspension, while the rest moves as bedload.     

Constitutive model for soil amplification of ground shaking: Parameter calibration,comparisons, validation

A one-dimensional constitutive model, developed for the nonlinear ground response analysis of layered soil deposits, is calibrated and validated experimentally in this paper. The small number of parameters renders the model easily implementable, yet quite flexible in effectively reproducing almost any type of experimentally observed hysteretic soil behavior. In particular, the model generates realistic shear modulus and damping curves as functions of shear strain, as well as stress–strain hysteresis loops. The model is calibrated against three sets of widely-used published shear modulus and damping (G : γ and ξ : γ) curves and a library of parameter values is assembled to facilitate its use. The model, along with a developed explicit finite-difference code, NL-DYAS, for analyzing the wave propagation in layered hysteretic soil deposits, is tested against established constitutive models and numerical tools such as Cyclic1D [12] and SHAKE [42], and validated against experimental data from two centrifuge tests. Emphasis is given on the proper assessment of the V_s profile in the centrifuge tests, on the role of soil nonlinearity, and on comparisons of two inelastic codes (NL-DYAS and Cyclic1D) with equivalent linear (SHAKE) analysis.     

Simulating seasonal to multi-decadal variation in lake thermal response to meteorological forcing using the UCLAKE 1-dimensional model code

Lake temperature responses to climate forcing are of interest on account of the important linkages between water temperature and ecosystem processes. This paper describes a new 1-dimensional (1D) numerical model code and its application to investigations of multi-scale linkages between the vertical temperature structure and meteorological forcing. UCLAKE is implemented as highly portable open-source software, based on computationally efficient algorithms, and able to resolve sub-daily (e.g., hourly) dynamics while retaining the efficiency to simulate multi-decadal time scales.A UCLAKE model is calibrated and validated against thermistor profile time series for a small upland lake in North Wales, UK. Some of the challenges in 1D model calibration are explored and a sensitivity analysis reveals a dependence of optimal parameter set values on water column depth and time. An exploratory 52-year hindcast simulation demonstrates the computational efficiency of UCLAKE for multi-decadal studies of trends in lake temperature that vary with depth. A supplementary application of UCLAKE to Windermere, in the English Lake District, demonstrates its performance for larger and deeper lakes.     

Direct numerical simulation of pore-scale flow in a bead pack: Comparison with magnetic resonance imaging observations

A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier–Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental uncertainty in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods. © 2012 Elsevier Science. All rights reserved.     

Total electron content monitoring using triple frequency GNSS data: A three-step approach

Global Navigation Satellite Systems (GNSS) are very useful tools to study the ionosphere. Nevertheless, the precision of the usual dual frequency total electron content (TEC) monitoring technique is affected by code delays (hardware and multipath), and is therefore limited. This paper introduces a TEC monitoring technique based on triple frequency GPS and Galileo measurements. The three steps of this technique are validated on triple frequency simulated data. In fact, as it is not affected by code delays, the precision of the reconstructed TEC is improved in regards with the dual frequency technique.     

A new boundary element solution to evaluate the geometric effects of the canyon site on the displacement response spectrum

This paper presents a step-by-step procedure using the three-dimensional boundary element approach to study the behavior of semi-circular canyons under seismic shear waves. The boundary element code TDASC allows utilization for various canyon geometries, evaluation of concurrent seismic waves and calculation of the ground motions on canyons due to an excitation at any arbitrary point of the incident field. Considering the widening ratio of the canyon(including prismatic, semi-prismatic and non-prismatic canyons), wave characteristics(wavelength, dimensionless period, direction) and maximum amplification pattern, the solution was applied to carry out a series of parametric studies. It was shown that canyon form can significantly affect the displacement amplification, especially at the points located on its edges. By increasing the wave dimensionless frequency(η 1), the amplification pattern becomes more complex. On the basis of the results from a variety of considered cases, a new expression has been presented for the limiting wavelength beyond which the widening of the canyon will not have a major effect on the displacement amplification. To verify the reliability of the proposed approach, the obtained results, expressed in terms of displacement amplitude, were compared with those from the available published literature and a reasonably good agreement was observed.