Transient Leakance and Infiltration Characteristics during Lake Bank Filtration

Infiltration capacity of bank filtration systems depends on water extraction and hydraulic resistance of the bed sediments. Lakebed hydraulics may be especially affected by clogging, which is dependent on settlement of fine particles, redox potential, and other factors. In the field, most of these processes are difficult to quantify, and thus, when calculating response to pumping the water flux across the sediment surface is assumed to be linearly dependent on the hydraulic gradient. However, this assumption was not adequate to describe conditions at a bank filtration site located at Lake Tegel, Berlin, Germany. Hence, we first assumed the leakage coefficient (or leakance) is spatially distributed and also temporally variant. Furthermore, observations show that the leakance is considerably higher in shallow than in deeper areas; hence, leakance was assumed to be dependent on the existence and thickness of an unsaturated zone below the lake. The proposed explanation of spatial and temporal variability in leakance involves a hypothesis for redox dependent and reversible biogeochemical clogging, supported by geochemical observations in surface water and ground water. Four leakance approaches are implemented in the ground water flow code MODFLOW2000 and calibrated by inverse modeling using the parameter estimation software PEST. These concepts are evaluated by examining the fit to the hydraulic heads, to infiltration measurements, transport modeling results, and considering the degrees of freedom due to the number of calibration parameters. The leakage concept based on the assumption of the influence of an unsaturated zone on clogging processes best explains the field data.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪—集安深部地质调查重磁电综合地质地球物理解释中.结果表明:该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪—集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

Shallow geology characterization using Rayleigh and Love wave dispersion curves derived from seismic noise array measurements

The local geology and shallow S-wave velocity structure of a site are recognized to be key factors for the increase in the damaging potential of seismic waves. Indeed, seismic amplitudes may be amplified in frequency ranges unfavorable for building stock by the presence of soft sedimentary covers over lying hard bedrock. Hence, microzonation activities, which aim at assessing the site response as accurately as possible, have become a fundamental task for the seismic risk reduction of urbanized areas. Methods based on the measurement of seismic noise, which typically are fast, non-invasive, and low cost, have become a very attractive option in microzonation studies.Using observations derived from seismic noise recordings collected by two-dimensional arrays of seismic stations, we present a novel joint inversion scheme for surface wave curves. In particular, the Love wave, the Rayleigh wave dispersion and the HVSR curves are innovatively combined in a joint inversion procedure carried out following a global search approach (i.e., the Genetic Algorithm).The procedure is tested using a data set of seismic noise recordings collected at the Bevagna (Italy) test-site. The results of the novel inversion scheme are compared with the inversion scheme proposed by Parolai et al. (2005), where only Rayleigh wave dispersion and HVSR curves are used, and with a cross-hole survey.     

Towards improved post‐processing of hydrologic forecast ensembles

Forecast ensembles of hydrological and hydrometeorologial variables are prone to various uncertainties arising from climatology, model structure and parameters, and initial conditions at the forecast date. Post‐processing methods are usually applied to adjust the mean and variance of the ensemble without any knowledge about the uncertainty sources. This study initially addresses the drawbacks of a commonly used statistical technique, quantile mapping (QM), in bias correction of hydrologic forecasts. Then, an auxiliary variable, the failure index (γ), is proposed to estimate the ineffectiveness of the post‐processing method based on the agreement of adjusted forecasts with corresponding observations during an analysis period prior to the forecast date. An alternative post‐processor based on copula functions is then introduced such that marginal distributions of observations and model simulations are combined to create a multivariate joint distribution. A set of 2500 hypothetical forecast ensembles with parametric marginal distributions of simulated and observed variables are post‐processed with both QM and the proposed multivariate post‐processor. Deterministic forecast skills show that the proposed copula‐based post‐processing is more effective than the QM method in improving the forecasts. It is found that the performance of QM is highly correlated with the failure index, unlike the multivariate post‐processor. In probabilistic metrics, the proposed multivariate post‐processor generally outperforms QM. Further evaluation of techniques is conducted for river flow forecast of Sprague River basin in southern Oregon. Results show that the multivariate post‐processor performs better than the QM technique; it reduces the ensemble spread and is a more reliable approach for improving the forecast. Copyright © 2012 John Wiley & Sons, Ltd.     

半刚性连接钢框架抗震性能的模型试验研究

通过两榀1／2比例的一跨两层半刚性连接钢框架模型的低周反复荷载试验,研究了框架模型的破坏机制、、出饺顺序、滞回特性、耗能能力等抗震性能。从而得出结论：具有腹板双角钢的顶底角钢连接的钢框架由于节点连接刚度较小,因此框架在侧向力作用下变形较大,但同时此连接也改善了梁柱的内力分布,提高了框架的耗能能力,增加了结构抗震性能。     

Development and validation of a metallic haunch seismic retrofit solution for existing under‐designed RC frame buildings

The feasibility and efficiency of a seismic retrofit solution for existing reinforced concrete frame systems, designed before the introduction of modern seismic‐oriented design codes in the mid 1970s, is conceptually presented and experimentally investigated. A diagonal metallic haunch system is introduced at the beam–column connections to protect the joint panel zone from extensive damage and brittle shear mechanisms, while inverting the hierarchy of strength within the beam–column subassemblies and forming a plastic hinge in the beam. A complete step‐by‐step design procedure is suggested for the proposed retrofit strategy to achieve the desired reversal of strength hierarchy. Analytical formulations of the internal force flow at the beam–column‐joint level are derived for the retrofitted joints. The study is particularly focused on exterior beam–column joints, since it is recognized that they are the most vulnerable, due to their lack of a reliable joint shear transfer mechanism. Results from an experimental program carried out to validate the concept and the design procedure are also presented. The program consisted of quasi‐static cyclic tests on four exterior, ? scaled, beam–column joint subassemblies, typical of pre‐1970 construction practice using plain round bars with end‐hooks, with limited joint transverse reinforcement and detailed without capacity design considerations. The first (control specimen) emulated the as‐built connection while the three others incorporated the proposed retrofitted configurations. The experimental results demonstrated the effectiveness of the proposed solution for upgrading non‐seismically designed RC frames and also confirmed the applicability of the proposed design procedure and of the analytical derivations. Copyright © 2006 John Wiley & Sons, Ltd.     

Appraisal of a new 1D weighted joint inversion of ground based and helicopter‐borne electromagnetic data

In order to couple spatial data from frequency‐domain helicopter‐borne electromagnetics with electromagnetic measurements from ground geophysics (transient electromagnetics and radiomagnetotellurics), a common 1D weighted joint inversion algorithm for helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics data has been developed. The depth of investigation of helicopter‐borne electromagnetics data is rather limited compared to time‐domain electromagnetics sounding methods on the ground. In order to improve the accuracy of model parameters of shallow depth as well as of greater depth, the helicopter‐borne electromagnetics, transient electromagnetics, and radiomagnetotellurics measurements can be combined by using a joint inversion methodology. The 1D joint inversion algorithm is tested for synthetic data of helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics. The proposed concept of the joint inversion takes advantage of each method, thus providing the capability to resolve near surface (radiomagnetotellurics) and deeper electrical conductivity structures (transient electromagnetics) in combination with valuable spatial information (helicopter‐borne electromagnetics). Furthermore, the joint inversion has been applied on the field data (helicopter‐borne electromagnetics and transient electromagnetics) measured in the Cuxhaven area, Germany. In order to avoid the lessening of the resolution capacities of one data type, and thus balancing the use of inherent and ideally complementary information content, a parameter reweighting scheme that is based on the exploration depth ranges of the specific methods is proposed. A comparison of the conventional joint inversion algorithm, proposed by Jupp and Vozoff ( 1975 ), and of the newly developed algorithm is presented. The new algorithm employs the weighting on different model parameters differently. It is inferred from the synthetic and field data examples that the weighted joint inversion is more successful in explaining the subsurface than the classical joint inversion approach. In addition to this, the data fittings in weighted joint inversion are also improved.     

Seismic velocity estimation from post-critical wide-angle reflections in layered structures

The travel time inversion of wide-angle seismic data is a technique commonly used in the deep seismic sounding. We propose an application of this technique to a smaller scale of a sedimentary layer, where the characteristics of seismic observations changes significantly. Field observations confirmed by synthetic analysis recognize the dominant amplitudes of wide-angle post-critical reflections. A case study is presented in this paper, of a joint interpretation of conventional reflection seismic with reflection imaging, combined with the wide-angle travel time inversion of additional full-spread observations. A joint interpretation results in a precise recognition of the seismic velocity distribution, that is further used for the seismic depth conversion with the uncertainty analysis of the depth of the reflecting horizons. Despite the salt layer in the studied structure this method is able to precisely recognize the seismic velocities of the sub-salt structures.     

Residual drift spectra for RC bridge columns subjected to near-fault earthquakes

Permanent displacement of a bridge column can be directly measured during the inspection after near-fault earthquakes.However,the engineer needs to estimate the expected residual drift at the design stage to determine if the bridge seismic performance is satisfactory.The most direct method to estimate the residual displacement is nonlinear response history analysis,which is time consuming and cumbersome.Alternatively,an attractive but indirect method is generating estimated residual displacement spectra that depend on displacement ductility demand,column period,site conditions,and earthquake characteristics.Given the period and the expected displacement ductility demand for the column,the residual drift response spectra curves can be utilized to estimate the residual drift demand.Residual drift spectra that are applicable to RC bridge columns in different parts of the United States were developed based on nonlinear response history analyses using a comprehensive collection of recorded and synthetic near-fault ground motions and were linked to one-second spectral acceleration(S1)of the AASHTO maps.It was also found that the residual drift ratio is below one percent when S1 is less than 0.6 g.     

Obliquely incident wave propagation across rock joints with virtual wave source method

Due to the presence of joints, waves are greatly attenuated when propagating across rock masses. Zhu et al. (2011) (Normally incident wave propagation across a joint set with virtual wave source method. J. Appl. Geophys.73, 283–288.) studied normally incident wave propagation across a joint set with the virtual wave source method (VWSM). The introduced VWSM has merits in some aspects, especially the capability of separating differently arriving transmitted waves. However, normal wave incidence is only the special case for wave incidence with arbitrary incident angles. Obliquely incident wave propagation across a joint set is more complicated than normally incident wave propagation due to wave transformation at the joints. As a continuation of the previous paper, this work is extended to analytically study obliquely incident wave propagation across joints with VWSM. Complete theoretical reflection and transmission coefficients across single joint described by displacement discontinuity model are derived through plane wave analysis. The superposition of P wave and S wave is for the first time mathematically expressed and studied. The VWSM is verified through comparison with the propagation matrix method. Through extensive parametric studies on wave transmission across single and multiple parallel joints, it is shown that transmitted wave energy is mainly constrained in the transmitted wave of the same type as the incident wave. And with increasing joint stiffness, the transmission coefficients across single joint increases except those whose wave type is different from the incident wave. The amplitude of superposed transmitted wave for P wave incidence increases with incident angle, which is coincident with field observations. Both joint spacing and number of joints have significant effects on transmission coefficients. We find that when joint spacing is sufficiently large, the transmission coefficient is no longer a constant as the normally incident wave propagation case (Zhu et al., 2011). And when joints are very closely spaced, wave attenuation depends little on the number of joints, which is different from the conclusions from equivalent medium method.     

一种高效的基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演策略

基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率.     

Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models

In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information.     

Nonlinear modeling of cyclic response of RC beam–column joints reinforced by plain bars

Experience of previous earthquakes shows that a considerable portion of buildings reinforced with plain bars sustain relatively large damages especially at the beam–column joints where the damages are mostly caused by either diagonal shear cracks or intersectional cracks caused by bar slippage. While previous works mainly focus on shear failure mode, in this study, the emphasis is placed on slip based cracks as the dominant failure mode. A systematic procedure is introduced to predict the dominant failure mode at the joint which is based on the dimensional properties, reinforcement details, and axial and shear load at the joint. In addition, a relatively simple and efficient nonlinear model is proposed to simulate pre- and post-elastic behavior of the joints which fail under bar slippage mode. In this model, beam and column components are represented by linear elastic elements, dimensions of the joint panel are defined by rigid elements, and effect of slip is taken into account by a nonlinear rotational spring at the end of the beam. The proposed method is validated by experimental results for both internal and external joints .     

GOCE Quick-Look Gravity Solution: Application of the Semianalytic Approach in the Case of Data Gaps and Non-Repeat Orbits

The satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer), the first Core Mission of the Earth Explorer Programme funded by ESA (European Space Agency), is dedicated to the precise modelling of the Earth's gravity field, with its launch planned for 2006. The mathematical models for parameterizing the Earth's gravity field are based on a series expansion into spherical harmonics, yielding a huge number of unknown coefficients. Their computation leads to the solution of very large normal equation systems. An efficient way to handle these equation systems is the so-called semianalytic or lumped coefficients approach, which theoretically requires an uninterrupted, continuous time series of observations, recorded along an exact circular repeat orbit. In this paper the consequences of violating these conditions are analyzed. The effects of an interrupted observation stream onto the estimated spherical harmonic coefficients are demonstrated, and an iterative strategy, which reduces the negative influence depending on the characteristics of the data gaps, is proposed. Additionally, the impact of an imperfectly closing orbit (non-repeat orbit) on the gravity field model is analyzed, and a strategy to minimize the corresponding errors is presented. The applicability of the semianalytic approach also to a joint inversion of satellite-to-satellite tracking data in high-low mode (hl-SST) and satellite gravity gradiometry (SGG) observations is demonstrated, where the analysis of the former component is based on the energy conservation law. Several realistic case studies prove that the semianalytic approach is a feasible tool to generate quick-look gravity solutions, i.e. fast coefficient estimates using only partial data sets. This quick-look analysis shall be able to detect potential distortions of statistical significance (e.g. systematic errors) in the input data, and to give a fast feedback to the GOCE mission control.     

带隅撑胶合木梁-双肢柱框架结构节点抗震性能试验研究

针对胶合木梁-柱节点抗弯能力弱的问题,提出了一种带隅撑的胶合木梁-双肢柱框架结构节点,以提高节点的转动刚度,改善梁柱结构框架的抗侧性能.对3组6个足尺胶合木梁-双肢柱框架结构节点进行了单调和低周反复加载试验以研究节点的抗震性能.结果表明:无隅撑节点类似于铰接,抗弯承载力很小,增设隅撑后显著提高了节点的转动刚度和抗弯承载...     

Scales of columnar jointing in igneous rocks: field measurements and controlling factors

Columnar jointing is a common feature of solidified lavas, sills and dikes, but the factors controlling the characteristic stoutness of columns remain debated, and quantitative field observations are few in number. In this paper, we provide quantitative measurements on sizing of columnar joint sets and our assessment of the principal factors controlling it. We focus on (1) chemistry, as it is the major determinant of the physical (mechanical and thermal) properties of the lava, and (2) geology, as it influences the style of emplacement and lava geometry, setting boundary conditions for the cooling process and the rate of heat loss. In our analysis, we cover lavas with a broad range of chemical compositions (from basanite to phonolite, for six of which we provide new geochemical analyses) and of geological settings. Our field measurements cover 50 columnar jointing sites in three countries. We provide reliable, manually digitized data on the size of individual columns and focus the mathematical analysis on their geometry (23,889 data on side length, of which 17,312 are from full column sections and 3,033 data on cross-sectional area and order of polygonality). The geometrical observations show that the variation in characteristic size of columns between different sites exceeds one order of magnitude (side length ranging from 8 to 338 cm) and that the column-bounding polygons’ average order is less than 6. The network of fractures is found to be longer than required by a minimum-energy hexagonal configuration, indicating a non-equilibrium, geologically quick process. In terms of the development and characteristic sizing of columnar joint sets, our observations suggest that columns are the result of an interplay between the geological setting of emplacement and magma chemistry. When the geological setting constrains the geometry of the emplaced body, it exerts a stronger control on characteristic column stoutness. At unconstrained geometries (e.g. unconfined lava flows), chemistry plays the major role, resulting in stouter columns in felsic lavas and slenderer columns in mafic lavas.     

Joint Perturbation in Geoacoustic Emission,Radon, Thoron,and Atmospheric Electric Field Based on Observations in Kamchatka

Izvestiya, Physics of the Solid Earth - Based on the integration and analysis of the published theoretical results and field observations, the possibility is considered of a joint perturbation in...     

Large-scale atmospheric patterns associated with mesoscale features leading to extreme precipitation events in Northwestern Italy

The synoptic and large-scale atmospheric conditions for heavy rainfall events in Northwestern Italy are diagnosed through the joint analysis of surface precipitation gauges and reanalysis atmospheric fields. Quantiles of local surface gauge precipitation observations are used to estimate the much larger-scale composite maps (conditional mean fields) of vertically integrated moisture flux, low-level winds, sea-level pressure, and 500 hPa height across the Atlantic and European domains. Remarkably, coarse-resolution reanalysis data show distinct synoptic conditions for heavy precipitation in localized regions that are below the resolution of the reanalysis. In this paper the key attributes of the new approach that is based on the joint analysis of gridded reanalysis and station data are presented. The application of the methodology are used to establish supporting evidence for hydrometeorological processes that lead to extreme precipitation across Northwest Italy. The results confirm the role of large-scale flow features that interact with regional topography in producing localized extreme precipitation. Whereas previous studies were based on a few case studies (modeled or observational), in this study the approach to producing a large ensemble of cases and composite statistics are introduced.     

感应-横向组合测井联合反演及多解性研究

电阻率测井反演问题解的局部收敛性和非唯一性增加了测井解释的难度，影响了实际的应用，把测量原理迥异的感应型和电流型测井结果联合起来反演，是减小这些影响的有效方法，利用经测井理论和反演技术，建立了适用感应－横向组合测井响应的联合演算法，研究了单一的横向测井反演的感应－横向联合反演的局部收敛性、非唯一性及解的偏差，结果表明，联合反演更有利于同时确定地层电阻率、侵入半径和侵入带电阻率，当初值地较大范围内变化时，反演结果均能在真值或真值附近收敛，改善了反演结果的局部收敛性，增加了反演结果的可信度，避免了差异较大乃至相互矛盾的解，现场实例证实了这一结论。     

Turbulent viscosity optimized by data assimilation

As an alternative approach to classical turbulence modelling using a first or second order closure, the data assimilation method of optimal control is applied to estimate a time and space-dependent turbulent viscosity in a three-dimensional oceanic circulation model. The optimal control method, described for a 3-D primitive equation model, involves the minimization of a cost function that quantifies the discrepancies between the simulations and the observations. An iterative algorithm is obtained via the adjoint model resolution. In a first experiment, a k ± L model is used to simulate the one-dimensional development of inertial oscillations resulting from a wind stress at the sea surface and with the presence of a halocline. These results are used as synthetic observations to be assimilated. The turbulent viscosity is then recovered without the k + L closure, even with sparse and noisy observations. The problems of controllability and of the dimensions of the control are then discussed. A second experiment consists of a two-dimensional schematic simulation. A 2-D turbulent viscosity field is estimated from data on the initial and final states of a coastal upwelling event.