地震安评中地震活动性的统计区域与b值

本文讨论我国现行地震安全性评价中,地震活动性研究统计区域的选择及b值计算中存在的问题,提出不一定以地震区带作为统计区域的建议。文中同时指出,累计b值的正确概念应该是无震级上限累计b值,应采用迭代回归计算。文中用实际例子说明作者的建议。     

Impact of artificial subsurface drainage on groundwater travel times and baseflow discharge in an agricultural watershed,Iowa (USA)

Pollutant delivery through artificial subsurface drainage networks to streams is an important transport mechanism, yet the impact of drainage tiles on groundwater hydrology at the watershed scale has not been well documented. In this study, we developed a two‐dimensional, steady‐state groundwater flow model for a representative Iowa agricultural watershed to simulate the impact of tile drainage density and incision depth on groundwater travel times and proportion of baseflow contributed by tile drains. Varying tile drainage density from 0 to 0.0038 m^?1, while maintaining a constant tile incision depth at 1.2 m, resulted in the mean groundwater travel time to decrease exponentially from 40 years to 19 years and increased the tile contribution to baseflow from 0% to an upper bound of 37%. In contrast, varying tile depths from 0.3 to 2.7 m, while maintaining a constant tile drainage density of 0.0038 m^?1, caused mean travel times to decrease linearly from 22 to 18 years and increased the tile contribution to baseflow from 30% to 54% in a near‐linear manner. The decrease in the mean travel time was attributed to decrease in the saturated thickness of the aquifer with increasing drainage density and incision depth. Study results indicate that tile drainage affects fundamental watershed characteristics and should be taken into consideration when evaluating water and nitrate export from agricultural regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Geological interpretation of a deep seismic‐reflection transect across the boundary between the Delamerian and Lachlan Orogens,in the vicinity of the Grampians,western Victoria

A deep seismic‐reflection transect in western Victoria was designed to provide insights into the structural relationship between the Lachlan and the Delamerian Orogens. Three seismic lines were acquired to provide images of the subsurface from west of the Grampians Range to east of the Stawell‐Ararat Fault Zone. The boundary between the Delamerian and Lachlan Orogens is now generally considered to be the Moyston Fault. In the vicinity of the seismic survey, this fault is intruded by a near‐surface granite, but at depth the fault dips to the east, confirming recent field mapping. East of the Moyston Fault, the uppermost crust is very weakly reflective, consisting of short, non‐continuous, west‐dipping reflections. These weak reflections represent rocks of the Lachlan Orogen and are typical of the reflective character seen on other seismic images from elsewhere in the Lachlan Orogen. Within the Lachlan Orogen, the Pleasant Creek Fault is also east dipping and approximately parallel to the Moyston Fault in the plane of the seismic section. Rocks of the Delamerian Orogen in the vicinity of the seismic line occur below surficial cover to the west of the Moyston Fault. Generally, the upper crust is only weakly reflective, but subhorizontal reflections at shallow depths (up to 3 km) represent the Grampians Group. The Escondida Fault appears to stop below the Grampians Group, and has an apparent gentle dip to the east. Farther east, the Golton and Mehuse Faults are also east dipping. The middle to lower crust below the Delamerian Orogen is strongly reflective, with several major antiformal structures in the middle crust. The Moho is a slightly undulating horizon at the base of the highly reflective middle to lower crust at 11–12 s TWT (approximately 35 km depth). Tectonically, the western margin of the Lachlan Orogen has been thrust over the Delamerian Orogen for a distance of at least 25 km, and possibly over 40 km.     

Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.     

Preliminary 3–D geological model of the Kalgoorlie region,Yilgarn Craton,Western Australia,based on deep seismic‐reflection and potential‐field data

The granite‐greenstone terranes of the Eastern Goldfields Province, Yilgarn Craton, Western Australia, are a major Australian and world gold and nickel source. The Kalgoorlie region, in particular, hosts several world‐class gold deposits. To attempt to understand why these deposits occur where they do, it is important to understand the crustal architecture in the region and how the major mineral systems operate in this architecture. One way to understand these relationships is to develop a detailed 3–D geological model for the region. The best method to map the 3–D geometry of major geological structures is by acquisition and interpretation of seismic‐reflection profiles. To contribute to this aim, a grid of deep seismic‐reflection traverses was acquired in 1999 to examine the 3–D geometry of the region in an area including the Kalgoorlie mineral region and mineral fields to the north and west. This grid was tied to the 1991 regional deep seismic traverse and 1997 high‐resolution seismic profiles in the same region. The grid covers an area measuring approximately 50 km wide by 50 km long and extended to a depth of approximately 50 km (below the base of the crust in this region). The resulting 3–D geological model was further constrained by both surface geological data and geophysical interpretations, with the seismic interpretations themselves also constrained by gravity and magnetic modelling. The 3–D model was used to investigate the geometric relationships between the major faults and shear zones in the area, the relationship between the granite‐greenstone succession and the basement, and the spatial relationships between the greenstones and the granites. Interpretation of the grid of seismic lines and construction of the 3–D geological model confirmed the existence of the detachment surface and led to the recognition that the granite‐greenstone contact usually occurs at a much shallower level than the detachment. Also, west‐dipping faults in the vicinity of the Golden Mile, including the Abattoir Shear through to Boulder‐Lefroy Fault, appear to be more important than previously thought in controlling the structure of that area. An antiformal thrust stack occurs beneath a triangle zone centred on the Golden Mile. The Black Flag Group was deposited in a probable extensional setting, and late extension was also probably more important than previously thought. The granite‐gneiss domes were uplifted by the formation of antiformal thrust stacks at depth beneath them.     

Screening the current soil drainage class map of Flanders (Belgium) for its predictive quality

The predictive quality of the current drainage class map of Flanders was evaluated using data from two monitoring networks: one with good spatial coverage but poor temporal coverage and another with better temporal but poor spatial coverage. We combine both networks to obtain 1678 point predictions for mean highest water (MHW) and mean lowest water (MLW) tables by applying time series modelling and total least squares regression. The resulting MHW and MLW point data set was used to evaluate the currency of the existing map and to identify regional differences. The quality of the current map is moderate, and large differences occur between regions. Especially the Campine region shows large and systematic differences, whereas the southeastern hills and chalk–loam region is relatively accurate. If more weight is given to errors in the wetter drainage classes, about 50% of the area of Flanders would benefit from remapping. Copyright © 2011 John Wiley & Sons, Ltd.     

A generalized conditional intensity measure approach and holistic ground‐motion selection

A generalized conditional intensity measure (GCIM) approach is proposed for use in the holistic selection of ground motions for any form of seismic response analysis. The essence of the method is the construction of the multivariate distribution of any set of ground‐motion intensity measures conditioned on the occurrence of a specific ground‐motion intensity measure (commonly obtained from probabilistic seismic hazard analysis). The approach therefore allows any number of ground‐motion intensity measures identified as important in a particular seismic response problem to be considered. A holistic method of ground‐motion selection is also proposed based on the statistical comparison, for each intensity measure, of the empirical distribution of the ground‐motion suite with the ‘target’ GCIM distribution. A simple procedure to estimate the magnitude of potential bias in the results of seismic response analyses when the ground‐motion suite does not conform to the GCIM distribution is also demonstrated. The combination of these three features of the approach make it entirely holistic in that: any level of complexity in ground‐motion selection for any seismic response analysis can be exercised; users explicitly understand the simplifications made in the selected suite of ground motions; and an approximate estimate of any bias associated with such simplifications is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

Palaeoproterozoic mid‐basin inversion in the northern Mt Isa terrane,Queensland

The Mellish Park Syncline is located in the northern part of the Mt Isa terrane. It has an axial trace that transects the remnants of the unconformity‐bounded Palaeoproterozoic Leichhardt and Isa Superbasins. The syncline is separated into a lower and upper component based upon variation in fold geometry across the basin‐bounding unconformity. The lower syncline, in the Leichhardt Superbasin, is tight and has an inclined west‐dipping axial plane. The upper syncline, in the Isa Superbasin, is open and upright. The geometry of the lower syncline is a consequence of a period of shortening and basin inversion which post‐dated the Leichhardt Rift Event (ca 1780–1740 Ma) and pre‐dated the Mt Isa Rift Event (ca 1710–1655 Ma), forming an open and upright north‐oriented syncline. Subsequent southeast tilting and half‐graben development during the Mt Isa Rift Event resulted in the lower syncline being tilted into its inclined geometry. Sequences of the Isa Superbasin were then deposited onto the eroded syncline. The geometry of the upper syncline reflects regional east‐west shortening during the Isan Orogeny (ca 1590–1500 Ma). The position of the upper syncline was largely controlled by the pre‐existing lower syncline. At this time the lower syncline was reactivated and tightened by flexural slip folding.     

Multifractal behaviour of long‐term karstic discharge fluctuations

Karstic watersheds are highly complex hydrogeological systems that are characterized by a multiscale behaviour corresponding to the different pathways of water in these systems. The main issue of karstic spring discharge fluctuations consists in the presence and the identification of characteristic time scales in the discharge time series. To identify and characterize these dynamics, we acquired, for many years at the outlet of two karstic watersheds in South of France, discharge data at 3‐mn, 30‐mn and daily sampling rate. These hydrological records constitute to our knowledge the longest uninterrupted discharge time series available at these sampling rates. The analysis of the hydrological records at different levels of detail leads to a natural scale analysis of these time series in a multifractal framework. From a universal class of multifractal models based on cascade multiplicative processes, the time series first highlights two cut‐off scales around 1 and 16 h that correspond to distinct responses of the aquifer drainage system. Then we provide estimates of the multifractal parameters α and C₁ and the moment of divergence q_D corresponding to the behaviour of karstic systems. These results constitute the first estimates of the multifractal characteristics of karstic spingflows based on 10 years of high‐resolution discharge time series and should lead to several improvements in rainfall‐karstic springflow simulation models. Copyright © 2012 John Wiley & Sons, Ltd.     

Crustal complexity in the Lachlan Orogen revealed from teleseismic receiver functions

There is an ongoing debate about the tectonic evolution of southeast Australia, particularly about the causes and nature of its accretion to a much older Precambrian core to the west. Seismic imaging of the crust can provide useful clues to address this issue. Seismic tomography imaging is a powerful tool often employed to map elastic properties of the Earth's lithosphere, but in most cases does not constrain well the depth of discontinuities such as the Mohorovi?i? (Moho). In this study, an alternative imaging technique known as receiver function (RF) has been employed for seismic stations near Canberra in the Lachlan Orogen to investigate: (i) the shear-wave-velocity profile in the crust and uppermost mantle, (ii) variations in the Moho depth beneath the Lachlan Orogen, and (iii) the nature of the transition between the crust and mantle. A number of styles of RF analyses were conducted: H-K stacking to obtain the best compressional–shear velocity (V _P /V _S) ratio and crustal thickness; nonlinear inversion for the shear-wave-velocity structure and inversion of the observed variations in RFs with back-azimuth to investigate potential dipping of the crustal layers and anisotropy. The thick crust (up to 48 km) and the mostly intermediate nature of the crust?mantle transition in the Lachlan Orogen could be due to the presence of underplating at the base of the crust, and possibly to the existing thick piles of Ordovician mafic rocks present in the mid and lower crust. Results from numerical modelling of RFs at three seismic stations (CAN, CNB and YNG) suggest that the observed variations with back-azimuth could be related to a complex structure beneath these stations with the likelihood of both a dipping Moho and crustal anisotropy. Our analysis reveals crustal thickening to the west beneath CAN station which could be due to slab convergence. The crustal thickening may also be related to the broad Macquarie volcanic arc, which is rooted to the Moho. The crustal anisotropy may arise from a strong N–S structural trend in the eastern Lachlan Orogen and to the preferred crystallographic orientation of seismically anisotropic minerals in the lower and middle crust related to the paleo-Pacific plate convergence.