地震波的场方程矩阵和能量的正定二次型及其意义

场方程是能够表述半空间地震波场的整体特征.波动方程、速度方程和能量方程.通过分析可知每个场方程都具有各自的“场方程矩阵”.能量方程能够对所有场方程矩阵进行综合和贯通,给出了能量方程以“弹性矩阵”为核心的普适性表达形式.最后,运用矩阵的正定二次型理论阐述了“能量矩阵与弹性矩阵”之间一致的对称性和正定性.能量矩阵蕴含的动态力的平衡关系、速度的时间_空间分布和能量的传播及变化的物理意义,能够从能量矩阵的正定二次型特性表述出来.本文研究分析问题的方法完全适用于复杂介质模型,相关的认识和结论可以拓展到均匀黏弹性各向同性介质、均匀弹性各向异性介质、均匀黏弹性各向异性介质以及比奥饱和流体介质.     

小波变换应用于日长变化的研究进展

随着空间测地技术的发展,日长变化的观测资料已经达到了空前高的程度,迫切需要合适的数学工具和方法从这些高精度的观测资料中提取有价值的信息.近些年,小波变换被应用于天文学和地球科学的资料分析中,它通过伸缩和平移等运算对时间信号进行多尺度的细化分析,从中提取有效的信息.本文主要介绍了小波变换应用于日长变化的研究进展,并初步探讨该领域以后的研究方向.     

塔里木河流域径流变化趋势及其对气候变化的响应

This paper has studied the change of streamflow and the impact of climatic variability conditions on regional hydrological cycle in the headwater of the Tarim River Basin. This study investigates possible causes of observed trends in streamflow in an environment which is highly variable in terms of atmospheric conditions, and where snow and ice melt play an important role in the natural hydrological regime. The discharge trends of three head streams have a significant increase trend from 1957 to 2002 with the Mann-Kendall test. Complex time-frequency distributions in the streamflow regime are demonstrated especially by Morlet wavelet analysis over 40 years. The purpose is to ascertain the nature of climatic factors spatial and temporal distribution, involved the use of EOF （Empirical Orthogonal Function） to compare the dominant temperature, precipitation and evaporation patterns from normally climatic records over the Tarim＇s headwater basin. It shows that the first principal component was dominated since the 1990s for temperature and precipitation, which identifies the significant ascending trend of spatial and temporal pattern characteristics under the condition of the global warming. An exponential correlation is highlighted between surface air temperature and mean river discharge monthly, so the regional runoff increases by 10%-16% when surface air temperature rises by 1 ℃. Results suggest that headwater basins are the most vulnerable environments from the point of view of climate change, because their watershed properties promote runoff feeding by glacier and snow melt water and their fundamental vulnerability to temperature changes affects rainfall, snowfall, and glacier and ice melt.     

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.     

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.     

Stage level,volume and time‐frequency information content of Lake Tana using stochastic and wavelet analysis methods

Lake Tana is the largest fresh water body situated in the north‐western highlands of Ethiopia. In addition to its ecological services, it serves for local transport, electric power generation, fishing, recreational purposes, and source of dry season irrigation water supply. Evidence shows that the lake has dried at least once at about 15,000–17,000 before present owing to a combination of high evaporation and low precipitation events. Past attempts to understand and simulate historical fluctuation of Lake Tana based on simplistic water balance approach of inflow, outflow, and storage have failed to capture well‐known events of drawdown and rise of the lake that have happened in the last 44 years. This study tested different stochastic methods of lake level and volume simulation for supporting Lake Tana operational planning decision support. Three stochastic methods (perturbations approach, Monte Carlo methods, and wavelet analysis) were employed for lake level and volume simulation, and the results were compared with the stage level measurements. Forty‐four years of daily, monthly, and mean annual lake level data have shown a Gaussian variation with goodness of fit at 0.01 significant levels of the Kolmogorov–Smirnov test. The stochastic simulations predicted the lake stage level of the 1972, 1984, and 2002/2003 historical droughts 99% of the time. The information content (frequency) of fluctuation of Lake Tana for various periods was resolved using Wigner's Time‐Frequency Decomposition method. The wavelet analysis agreed with the perturbations and Monte Carlo simulations resolving the time (1970s, 1980s, and 2000s) in which low frequency and high spectral power fluctuation has occurred. The Monte Carlo method has shown its superiority for risk analysis over perturbation and deterministic method whereas wavelet analysis reconstructed historical record of lake stage level at daily and monthly time scales. Copyright © 2012 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.     

Neuro-fuzzy models employing wavelet analysis for suspended sediment concentration prediction in rivers

Abstract

The study of sediment load is important for its implications to the environment and water resources engineering. Four models were considered in the study of suspended sediment concentration prediction: artificial neural networks (ANNs), neuro-fuzzy model (NF), conjunction of wavelet analysis and neuro-fuzzy (WNF) model, and the conventional sediment rating curve (SRC) method. Using data from a US Geological Survey gauging station, the suspended sediment concentration predicted by the WNF model was in satisfactory agreement with the measured data. Also the proposed WNF model generated reasonable predictions for the extreme values. The cumulative suspended sediment load estimated by this model was much higher than that predicted by the other models, and is close to the observed data. However, in the current modelling, the ANN, NF and SRC models underestimated sediment load. The WNF model was successful in reproducing the hysteresis phenomenon, but the SRC method was not able to model this behaviour. In general, the results showed that the NF model performed better than the ANN and SRC models.

Citation Mirbagheri, S. A., Nourani, V., Rajaee, T. & Alikhani, A. (2010) Neuro-fuzzy models employing wavelet analysis for suspended sediment concentration prediction in rivers. Hydrol. Sci. J. 55(7), 1175–1189.     

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.     

An answer to an important controversy and the need for caution when using simple models to predict inelastic earthquake response of buildings with torsion

This paper presents evidence that the extension of conclusions based on the widely used simplified, one story, eccentric systems of the shear‐beam type, to actual, nonsymmetric buildings and consequent assessments of the pertinent code provisions, can be quite erroneous, unless special care is taken to match the basic properties of the simplified models to those of the real buildings. The evidence comes from comparisons of results obtained using three variants of simplified models, with results from the inelastic dynamic response of three‐ and five‐story eccentric buildings computed with detailed MDOF systems, where the members are idealized with the well‐known plastic hinge model. In addition, a convincing answer is provided on a pertinent hanging controversy: For frame‐type buildings, designed in accordance with the dynamic provisions of modern codes (such as EC8 or IBC2000), which allow reduced shears at the stiff edge due to torsion, the frames at the flexible sides are the critical elements in terms of ductility demands. Copyright © 2009 John Wiley & Sons, Ltd.