Lateral lithospheric conductivity patterns,from European synchronous observatory data

Much information on the regional lithospheric structure may come from MTS data acquired by synchronous 2D arrays and processed with regard to the nonlocal response of a laterally inhomogeneous subsurface. We suggest to invert the nonlocal MT responses applying correlation of all surface horizontal and vertical components of the geomagnetic field recorded simultaneously at all stations. The inversion algorithm has been applied to 2004–2005 European observatory data of diurnal S_q variations for first five harmonics and yielded lateral conductivity patterns for different periods. The maps show spatial correlation between conductivity maxima and lithospheric thickness minima and, specifically, highlight the contours of the Pannonian basin, where lithosphere is as thin as ～50 km, from seismic data.     

Quasi-periodic geomagnetic secular variation (from 1985–2005 world observatory data)

We have discovered a 2–4 year periodicity in geomagnetic secular variation (SV) from data of 110 world magnetic observatories. The periodicity in the horizontal component (H) is most prominent and appears to be globally uniform in different regions, on all continents, and in both hemispheres. The quasi-periodic short-wavelength variations show up in the vertical component (Z) as well but locally superpose on long-wavelength regional anomalies. We presume that the short-period fluctuations may be produced by instability of the eccentric dipole (ED) axis proceeding from the analysis of the SV field and optimization modeling of the dipole field with varied ED parameters.     

Effect of inclined conductivity anisotropy on frequency induction and TEM data

The harmonic electromagnetic field of a vertical magnetic dipole above an anisotropic half-space has been simulated using a forward algorithm for layered conductive media with inclined anisotropy. Inclined anisotropy has been found out to change the typical behavior of frequency and transient responses. Qualitative interpretation of FD loop–loop responses of a conducting earth with inclined anisotropy requires taking into account the receiver azimuth dependence of apparent resistivities. In the case of time-domain measurements, this dependence is absent but the apparent resistivities are higher at late times.     

Evaluation of crustal and upper mantle structures using receiver function analysis: ISM broadband observatory data

A three-component broadband seismograph is in operation since January 2007 at the Indian School of Mines (ISM) campus, Dhanbad. We have used the broadband (BB) seismograms of 17 teleseismic events (M ≥ 5.8) recorded by this single BB station during 2008–09 to estimate the crust and upper mantle discontinuities in Dhanbad area which falls in the peninsular India shield. The converted wave technique and the Receiver function analysis are used. A 1-D velocity model has been derived using inversion. The Mohorovicic (Moho) discontinuity (crustal thickness) below the ISM observatory is estimated to be ∼41 km, with an average Poisson ratio of ∼0.28, suggesting that the crust below the Dhanbad area is intermediate to mafic in nature. The single station BB data shed new light to the estimate of crustal thickness beneath the eastern India shield area, which was hitherto elusive. Further, it is observed that the global upper mantle discontinuity at 410 km is delayed by ∼0.6 sec compared to the IASP-91 global model; this may be explained by a slower/hotter upper mantle; while the 660 km discontinuity is within the noise level of data.     

An application of complex demodulation technique to geomagnetic data and conductivity anomaly studies

The method of complex demodulation has been used here to compute the amplitude and phase of a signal present in a geomagnetic series using, Banks’ method. It is found that the results are quite accurate when four or more continuous cycles of the signal of interest are present. The same limitations have also been observed when band-pass filters are used to isolate a signal. The information on phase is always correct. It is concluded that complex demodulation will give correct results for periodic variations like 11-year, 27-day, Sq or pulsations; but will give small values of amplitude for signals like sudden impulses or bays. The latter is not a limitation in conductivity studies where ratios likeZ/H, Z/D orH/D are mostly used in calculations. It has been shown that even with the records of one magnetic storm the cause of anomalies can be accurately identified which otherwise would need a large number of events.     

Measurement of streambed hydraulic conductivity and its anisotropy

 A method is described for the measurement of streambed hydraulic conductivity. Unlike permeameter methods, this method applies straight and l-shaped standpipes directly to streambeds for measurements of in-situ hydraulic conductivity in the vertical (K _v ) and horizontal (K _h ) directions, as well as in other oblique directions (K _s ). This method has advantages in determination of K _v values over grain-size analysis, permeameter tests, or slug test methods. Also unique to this method is that it provides K _s values of a streambed. The measured results can be used to construct a hydraulic conductivity ellipse and to evaluate the anisotropy of streambed sediments. Field examples from the Republican River, Nebraska, demonstrated the usefulness of this method in the determination of streambed hydraulic conductivity and anisotropy along or across a river channel. Results indicate that the K _h is about three to four times larger than K _v , whereas K _s values are larger than K _v but smaller than K _h . Received: 6 March 2000 · Accepted: 18 April 2000     

河床水力传导度及其各向异性的测定

采用直接测定法观测了黑河中游段河床水力传导度及其各向异性,结果表明:河床水力传导度不仅存在较强的各向异性,而且存在空间尺度上的变异性.河床中心位置在垂直、水平和θ=30°方向的平均水力传导度分别为0.45、22.49和1.71 m/d,河床边分别为5.95、29.69、16.80 m/d.在同一测点,水力传导度随着与水平方向的夹角增大呈幂函数曲线下降.试验结果表明河床边是河水的主要渗漏区,并且以侧渗为主.     

Estimating hydraulic conductivity from specific capacity for Hawaii aquifers,USA

Site-specific relationships between specific capacity and hydraulic parameters (transmissivity and hydraulic conductivity) were investigated for volcanic rocks in Maui, Hawaii, USA. Details about well construction were commonly ignored in previous studies. To improve on such efforts, specific-capacity values were normalized by the open interval of the well. Correcting specific capacity for turbulent head losses using step-drawdown tests and including aquifer penetration length improved the correlation between specific capacity and hydraulic conductivity and reduced uncertainty in the prediction of hydraulic parameters. The relationships provide estimates of aquifer parameters with correlation coefficients between 0.81 and 0.99. The relationships for Maui can probably be extended to other Hawaii islands, given the similarity of aquifer formations and a reasonable fit to step-drawdown data from Oahu. Hydraulic conductivity was estimated from 1,257 specific-capacity values in the Hawaii’s well database. Hydraulic-conductivity estimates for dike-free volcanic rocks are consistent on different islands. For all islands, the estimates range from 3 to 8,200 m/d, with a geometric-mean and median value of 272 and 291 m/d, respectively. A geostatistical approach was applied to Maui and Oahu to generate island-wide hydraulic-conductivity maps to facilitate groundwater management efforts.     

Upper mantle anisotropy inferred from shear wave splitting beneath the Eastern Indian Shield region

We estimate the shear wave splitting parameters vis-à-vis the thicknesses of the continental lithosphere beneath the two permanent seismic broadband stations located at Dhanbad (DHN) and Bokaro (BOKR) in the Eastern Indian Shield region. Broadband seismic data of 146 and 131 teleseismic earthquake events recorded at DHN and BOKR stations during 2007–2014 were analyzed for the present measurements. The study is carried out using rotation-correlation and transverse component minimization methods. We retain our “Good”, “Fair” and “Null” measurements, and estimate the splitting parameters using 13 “Good” results for DHN and 10 “Good” results for BOKR stations. The average splitting parameters (ϕ, δt) for DHN and BOKR stations are found to be 50.76°±5.46° and 0.82 ± 0.2 s and 56.30°±5.07° and 0.95 ± 0.17 s, and the estimated average thicknesses of the anisotropic layers beneath these two stations are ∼ 94 and ∼109 km, respectively. The measured deviation of azimuth of the fast axis direction (ϕ) from the absolute motion of the Indian plate ranges from ∼8° to 14°. The measured deviation of azimuth of the fast axis direction (ϕ) from the absolute motion of the Indian plate ranges from ∼8° to 14°. The eastward deviation of the fast axis azimuths from absolute plate motion direction is interpreted to be caused by induced outflow from the asthenosphere. Further, the delay time found in the present analysis is close to the global average for continental shield areas, and also coherent with other studies for Indian shield regions. The five “Null” results and the lower delay time of ∼0.5–0.6 s might be indicating multilayer anisotropy existing in the mantle lithosphere beneath the study area.     

Genesis of diamonds in the lower mantle

The “forbidden” assemblage (ferropericlase + enstatite) as inclusions in diamonds has been taken as evidence to imply that these inclusions and their host diamonds formed initially in the lower mantle. Magnesite is probably the only stable carbonate at depths greater than ∼220 km. Like dehydration reactions, the reaction boundary for the decarbonation of magnesite has a positive dT/dP slope at lower pressures, which becomes negative at higher pressures, if no other phase intervenes. This reaction boundary probably intersects the geotherm between ∼900 and ∼1100 km, below which magnesite decomposes into an assemblage periclase + diamond + oxygen. Thus, ferropericlase is the most likely inclusion in diamond formed in the lower mantle. The high frequency of sole occurrence of ferropericlase in diamonds from Sao Luiz, Brazil seems to substantiate the present speculation. Received: 8 June 1998 / Accepted: 28 September 1998     

混合型缓冲材料砌块渗透性及其各向异性研究

缓冲材料砌块具备良好的防渗性能,是保证高放废物处置库封闭性的关键。缓冲材料砌块经单轴加载挤压成型,块体尺寸大,测试与评价其渗透性能的均匀性和各向异性,对缓冲屏障设计具有重要价值。研究制备不同干密度的混合型缓冲材料砌块,定向切割获得砌块内部不同位置、不同方向的小试样,开展刚性壁和柔性壁渗透试验,综合评价了混合型缓冲砌块不同方向的水力传导系数。研究结果表明：砌块干密度越大,砌块的水力传导系数越低。所有砌块的水力传导系数均处于10?10 cm/s数量级,沿砌块z、?、r方向的水力传导系数基本相等,表明砌块渗透性具有良好的均匀性和各向同性,且满足高放缓冲材料防渗要求。两种试验方法获得的渗透结果差异较小,柔性壁渗透条件下砌块产生10%～20%的体积膨胀,导致其水力传导系数略大于刚性壁下的平行试样。引入CT试验,分析了渗透前后砌块的微观特性变化。试验证实,在饱和渗透之后,以CT值显示的砌块密度分布趋于均匀,砌块内部存在的微裂隙愈合消失。     

Structure of the Siberian upper mantle from superlong seismic profile data

The paper discusses the mantle structure along superlong seismic profiles in Russia examined using the method of homogenous functions. Two-dimensional heterogeneous sections of the upper mantle were calculated from travel-time curves to a depth of 500–600 km with allowance for the Earth’s curvature without using any a priory information. The presentation of sections as surfaces with a shaded relief combined with velocity contours allowed discerning the principal interfaces in the lithosphere and in the upper mantle, the internal structure of layers, and local heterogeneities of different shapes (convective cells and slabs) in the sections.     

Fluorine contents of some hydrous minerals derived from upper mantle and lower crust

Fluorine contents have been determined in about forty samples of amphibole, mica and apatite in alkali basalt and kimberlite and their incorporated xenoliths. They show a wide variation ranging from 15,000 to 100 ppm, corresponding to about 40 to 0.2 per cent substitution of F for OH in hydroxyl site of hydrous minerals. Fluorine abundances in these minerals reflect those of their host magmas or rocks; Itinome-gata xenoliths are the lowest and South African kimberlites and their xenoliths are the highest. F/OH and also. D/H (Kuroda et al. 1975) ratios in coexisting phlogopite-potassic richterite from peridotite and mica nodules are thoughts to have formed under no simple equlibrium conditions.     

Thermal and redox state of the subcontinental lithospheric mantle of NE Spain from thermobarometric data on mantle xenoliths

Mantle xenoliths in within-plate Cenozoic alkaline mafic lavas from NE Spain are used to assess the local subcontinental lithospheric mantle geotherm and the influence of melting and metasomatism on its oxidation state. The xenoliths are mainly anhydrous spinel lherzolites and harzburgites and gradations between, with minor pyroxenites. Most types show protogranular textures, but transitional protogranular–porphyroclastic and equigranular lherzolites also exist. Different thermometers used in the estimates provide higher subsolidus equilibrium temperatures for harzburgites (1,062 ± 29°C) than for lherzolites (972 ± 89°C), although there is overlap; the lowest temperatures correspond to porphyroclastic lherzolites, whereas pyroxenites give the highest temperatures (up to 1,257°C). Maximum pressures for subsolidus equilibrium of peridotites are at 2.0–1.8 GPa. Later they followed adiabatic decompression and harzburgites registered lower pressures (1.02 ± 0.19 GPa) than lherzolites (1.41 ± 0.27 GPa). One pyroxenite gives values consistent with the spinel lherzolite field (1.08 GPa). The shallowest barometric data are in agreement with the highest local conductive geotherms, which implies that the lithosphere–asthenosphere boundary is at 70–60 km minimum depth. Higher equilibrium temperatures for the harzburgites could be explained by the existence of mafic magma bodies or dykes at the lower crust–mantle boundary. Paleo-fO₂ conditions during partial melting as inferred from the covariation between V and MgO concentrations are mainly between QFM−1 and QFM−2 in log units. However, most thermobarometric fO₂ estimates are between QFM−1 and QFM+1, suggesting oxidation caused by later metasomatism during uplift and cooling.     

松辽盆地南部岩石圈地幔变形特征及各向异性分析——来自橄榄岩包体的约束

SKS测量结果显示松辽盆地快波方向分布较为复杂, 由于缺乏深部岩石变形资料约束, 制约了地震波各向异性成因解释。本文通过对松辽盆地南部双辽地区橄榄岩包体进行详细的岩石学、地球化学、显微结构、各向异性研究, 结果显示双辽地区橄榄岩包体的平衡温度为893~1152℃, 来源于岩石圈地幔。橄榄岩中橄榄石的晶格优选方位(CPO)类型主要为A型、D型和AG型, 其中, AG型和D型CPO可能形成于西太平洋板块向欧亚大陆俯冲回撤导致的岩石圈变形, AG型橄榄石CPO也可能形成于熔体存在下的橄榄石变形。基于CPO计算获得的橄榄岩包体全岩V_P各向异性(AV_P)为4.79%~11.80%, 最大剪切波各向异性(AV_Smax)为3.13%~7.93%。结合地球物理测量结果, 推断松辽盆地南部复杂的SKS各向异性的主要贡献可能来源于面理陡倾或直立的岩石圈地幔。

     

干旱半干旱区土壤含盐量和电导率高光谱估算

以新疆南疆地区渭干河-库车河三角洲绿洲不同盐渍化程度的土壤为研究对象, 将土壤样品的含盐量与电导率数据和使用ASD FieldSpec Pro FR光谱仪测得土壤样品高光谱数据作为本研究原始数据. 首先, 对原始土壤光谱反射率曲线进行Savitzky-Golay滤波以消除光谱曲线噪声可能引起的误差, 对得到的光谱曲线进行对数、 倒数等15种光谱变换; 其次, 对土壤样品含盐量与电导率关系进行分析, 并通过两者与土壤光谱反射率不同光谱变换形式的相关性比较分析, 遴选出一阶微分、 对数倒数一阶微分、 连续统去除和连续统去除一阶微分四种相关性较好的变换形式; 最后, 以此变换形式对含盐量和电导率进行建模, 并对二者进行了高光谱估算精度的比较. 结果表明: 渭干河-库车河三角洲绿洲的含盐量与电导率的相关性较高, 达到0.9975; 电导率与四种不同光谱变换形式之间的相关性要优于含盐量, 特别在一些土壤盐渍化的敏感波段尤为突出; 无论是含盐量还是电导率, 多元逐步回归模型均优于一元线性模型, 且以电导率建立的两种回归模型的决定系数均高于含盐量. 研究表明, 土壤电导率对高光谱信息的反应比含盐量更敏感, 同时, 以电导率建模的估算精度比含盐量更高. 因此, 以电导率替代含盐量进行土壤盐渍化高光谱估算研究是一种精度更高、 速度更快的方法, 可为提高土壤盐渍化高光谱估算提供理论依据.     

回归分析法在日变数据推算中的应用

在海洋磁测数据日变改正过程中,常会因各种因素限制,导致日变数据缺失,这时就需要结合多台站数据推算工区位置处的日变数据。通过借鉴前人研究成果及考察相关分析结果,可以近似认为台站间日变数据具有一定的线性关系,并引入统计学中常用的回归分析法拟合出某一位置日变数据的推算公式,进而得到推算数据。选取南海区域实测数据验证可知该方法具有一定的可靠性,可以解决日变数据缺失等问题,从而提升海洋磁测数据的精度,是一种改进的磁测数据日变改正方法。     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Dissociation of CAS phase in the uppermost lower mantle

The high-pressure stability limit of calcium aluminosilicate (CAS) phase has been examined in its end-member CaAl₄Si₂O₁₁ composition at 18–39 GPa and 1,670–2,300 K in a laser-heated diamond-anvil cell (LHDAC). The in-situ synchrotron X-ray diffraction measurements revealed that the CAS phase decomposes into three-phase assemblage of cubic Al-bearing CaSiO₃ perovskite, Al₂O₃ corundum, and SiO₂ stishovite above 30 GPa and 2,000 K with a positive pressure–temperature slope. Present results have important implications for the subsolidus mineral assemblage of subducted sediment and the melting phase relation of basalt in the lower mantle.