Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan     

Akinobu Miyakoshi  Takeshi Hayashi  Atsunao Marui  Yasuo Sakura  Shinichi Kawashima  Masafumi Kawai 《International Journal of Earth Sciences》2008,97(2):401-411

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.  相似文献   

A virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations: a new approach for understanding spatial and temporal complexity of magma dynamics     

Diego Perugini  Maurizio Petrelli  Giampiero Poli 《Visual Geosciences》2008,13(1):1-24

In this contribution, we present a virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations with the aim to highlight the power of 3D representations in the understanding of this geological phenomenon. In particular, samples of mixed juveniles from Salina island (Southern Italy) are reconstructed in 3D by serial lapping and digital montage and numerical simulations are performed by using a 3D chaotic dynamical system. Natural and simulated magma mixing structures are visualized by using several multimedia tools including animations and “virtual reality” models. It is shown that magma interaction processes can generate large spatial and temporal compositional heterogeneities in magmatic systems. The same topological structures are observed in both 3D reconstructed rock samples and chaotic numerical simulations, indicating that the mixing of magmas is governed by chaotic dynamics. The use of 3D multimedia models gives the opportunity to penetrate into magma mixing structures and to understand their significance in the context of magma dynamics. Such an approach is very powerful since multimedia tools can strongly capture the attention of the reader bringing him/her into an interactive and memorable geological experience. Electronic supplementary material  enclosed:  相似文献   

山东玲珑金矿床成矿流体地球化学特征   总被引：3，自引：0，他引：3  

王可勇  张春燕  樊岳铭  张晓东 《吉林大学学报(地球科学版)》2008,38(2):194-201

玲珑金矿床第一成矿阶段与含金黄铁矿共生的石英中主要发育4种类型的原生流体包裹体:Ⅰ气液两相,Ⅱ含CO2三相,ⅢCO2,Ⅳ单液相包裹体。流体包裹体成分激光拉曼光谱分析及测温结果显示:①Ⅰa型包体,气液比10%~15%,均一温度为162.7~235.6℃,w(NaCl)(盐度)为4.65%~7.59%,气相平均摩尔分数为:H2O 96.48%,CO22.4%;②Ⅰb型包体,气液比30%~45%,均一温度266.9~349.2℃,w(NaCl)为10.8%~13.4%,气相平均摩尔分数为H2O 69.75%,CO224.74%;③Ⅱ型含CO2包体,CO2相所占比例为20%~90%,其均一温度为193.5~321.6℃,w(NaCl)2.9%~5.3%,CO2相中,H2O的摩尔分数为27.72%,CO2为70.6%。包裹体成分分析及测温结果综合研究认为,玲珑金矿成矿过程中存在大气降水热液与地幔来源流体的混合作用,前者与从流体中分离出的富CO2流体混合,以不同比例被捕获形成Ⅱ型包体;而与分异出CO2后的CO2不饱和地幔流体混合,被捕获形成Ⅰb型包体。两种流体混合导致的含矿热液物化条件变化对金的沉淀成矿具有重要意义。  相似文献   

Seismic hazard of Egypt   总被引：1，自引：0，他引：1  

A. El-Sayed  R. Wahlströ;m  O. Kulhánek 《Natural Hazards》1994,10(3):247-259

Earthquake hazard parameters such as maximum expected magnitude,M _max, annual activity rate,, andb value of the Gutenberg-Richter relation have been evaluated for two regions of Egypt. The applied maximum likelihood method permits the combination of both historical and instrumental data. The catalogue used covers earthquakes with magnitude 3 from the time interval 320–1987. The uncertainties in magnitude estimates and threshold of completeness were taken into account. The hazard parameter determination is performed for two study areas. The first area, Gulf of Suez, has higher seismicity level than the second, all other active zones in Egypt.b-values of 1.2 ± 0.1 and 1.0 ± 0.1 are obtained for the two areas, respectively. The number of annually expected earthquakes with magnitude 3 is much larger in the Gulf of Suez, 39 ± 2 than in the other areas, 6.1 ± 0.5. The maximum expected magnitude is calculated to be 6.5 ± 0.4 for a time span of 209 years for the Gulf of Suez and 6.1 ± 0.3 for a time span of 1667 years for the remaining active areas in Egypt. Respective periods of 10 and 20 years were reported for earthquakes of magnitude 5.0 for the two subareas.  相似文献   

Estimation of covariance matrix from geochemical data with observations below detection limits     

Chang-Jo F. Chung 《Mathematical Geology》1993,25(7):851-865

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.  相似文献   

硅质碎屑砂岩中胶结物来源的新认识     

杨香华 《地质科技情报》1993,12(4):34-39

分析了浅埋藏砂岩的成岩环境,并从岩石学、矿物学、同位素及质量平衡估算方面论证硅质碎屑砂岩中的胶结物主要来源于砂岩浅埋藏期间的地下水,进而分析浅埋藏带内地下水的流动机制及早期胶结作用对晚期成岩作用的影响。  相似文献   

施工顺序对软土地坪复合地基的影响探讨     

魏海涛  蔡智 《岩土工程技术》2021,(1)

场地回填土与搅拌桩施工顺序将直接影响软土地坪复合地基的处理效果,其中沉降控制是影响地坪安全运行的关键。以近海软土地区工业厂房水泥土搅拌桩复合地基工程为例,探讨了复合地基水泥土搅拌桩、回填土两种不同施工顺序的影响因素及相应处理效果;通过数值模拟试验,计算分析了不同填土厚度下两种工况的地基固结总沉降、施工沉降、工后沉降。结果表明,采用“先土后桩”的施工顺序的施工期沉降大于“先桩后土”,但工后沉降大大减小,提高了软土地坪地基处理效果,有利于处理后地坪的安全运行,验证了前述分析结论,为合理设计方案的确定提供了重要的参考依据。  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.  相似文献   

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data     

G. Grünthal  R. Wahlström 《Natural Hazards》2006,38(1-2):159-176

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.  相似文献   

Re-evaluation of focal depths and source mechanisms of selected earthquakes in the Afar depression   总被引：1，自引：0，他引：1  

L. Hagos  H. Shomali  R. Roberts 《Geophysical Journal International》2006,167(1):297-308

We present a stepwise inversion procedure to assess the focal depth and model earthquake source complexity of seven moderate-sized earthquakes  (6.2 > M _w > 5.1)  that occurred in the Afar depression and the surrounding region. The Afar depression is a region of highly extended and intruded lithosphere, and zones of incipient seafloor spreading. A time-domain inversion of full moment tensor was performed to model direct P and SH waves of teleseismic data. Waveform inversion of the selected events estimated focal depths in the range of 17–22 km, deeper than previously published results. This suggests that the brittle–ductile transition zone beneath parts of the Afar depression extends more than 22 km. The effect of near-source velocity structure on the moment tensor elements was also investigated and was found to respond little to the models considered. Synthetic tests indicate that the size of the estimated, non-physical, non-isotropic source component is rather sensitive to incorrect depth estimation. The dominant double couple part of the moment tensor solutions for most of the events indicates that their occurrence is mainly due to shearing. Parameters associated with source directivity (rupture velocity and azimuth) were also investigated. Re-evaluation of the analysed events shows predominantly normal faulting consistent with the relative plate motions in the region.  相似文献