Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

Mineralisation of amethyst-bearing geodes in Ametista do Sul (Brazil) from low-temperature sedimentary brines: evidence from monophase liquid inclusions and stable isotopes

Fluid inclusion studies in combination with hydrogen, oxygen and sulphur isotope data provide novel insights into the genesis of giant amethyst-bearing geodes in Early Cretaceous Paraná continental flood basalts at Amestita do Sul, Brazil. Monophase liquid inclusions in colourless quartz, amethyst, calcite, barite and gypsum were analysed by microthermometry after stimulating bubble nucleation using single femtosecond laser pulses. The salinity of the fluid inclusions was determined from ice-melting temperatures and a combination of prograde and retrograde homogenisation temperatures via the density maximum of the aqueous solutions. Four mineralisation stages are distinguished. In stage I, celadonite, chalcedony and pyrite formed under reducing conditions in a thermally stable environment. Low δ³⁴S_V-CDT values of pyrite (?25 to ?32?‰) suggest biogenic sulphate reduction by organotrophic bacteria. During the subsequent stages II (amethyst, goethite and anhydrite), III (early subhedral calcite) and IV (barite, late subhedral calcite and gypsum), the oxidation state of the fluid changed towards more oxidising conditions and microbial sulphate reduction ceased. Three distinct modes of fluid salinities around 5.3, 3.4 and 0.3 wt% NaCl-equivalent characterise the mineralisation stages II, III and IV, respectively. The salinity of the stage I fluid is unknown due to lack of fluid inclusions. Variation in homogenisation temperatures and in δ¹⁸O values of amethyst show evidence of repeated pulses of ascending hydrothermal fluids of up to 80–90 °C infiltrating a basaltic host rock of less than 45 °C. Colourless quartz and amethyst formed at temperatures between 40 and 80 °C, while the different calcite generations and late gypsum precipitated at temperatures below 45 °C. Calculated oxygen isotope composition of the amethyst-precipitating fluid in combination with δD values of amethyst-hosted fluid inclusions (?59 to ?51?‰) show a significant ¹⁸O-shift from the meteoric water line. This ¹⁸O-shift, high salinities of the fluid inclusions with chloride-sulphate composition, and high δ³⁴S values of anhydrite and barite (7.5 to 9.9?‰) suggest that sedimentary brines from deeper parts of the Guaraní aquifer system must have been responsible for the amethyst mineralisation.     

柴北缘绿梁山地区早古生代弧后盆地型蛇绿岩的年代学、地球化学及大地构造意义

岩石学、地球化学、年代学及Lu-Hf同位素综合研究表明在柴北缘西段绿梁山大平沟地区出露一套弧后盆地型蛇绿岩,岩石类型主要包括变质橄榄岩、变火山岩、变辉长岩及斜长花岗岩。其中变火山岩具有LREE亏损,类似N-MORB的稀土配分模式,同时又具有富集大离子亲石元素,亏损Nb、Ta等高场强元素的岛弧火山岩的地球化学特征,应形成在弧后盆地环境。斜长花岗岩为低钾准铝质花岗岩,具有LREE略微富集,HREE平坦的稀土配分型式,显示强烈Eu正异常,其εHf(t)值介于13.7～15.3之间,为变辉长岩部分熔融的产物,熔融温压条件可能为P=0.8～0.9GPa和T=～800℃。年代学研究结果表明变辉长岩的形成时代为535±2Ma,斜长花岗岩的形成时代为493±3Ma,指示本地区弧后盆地拉张时限至少介于493～535Ma之间,而柴北缘地区古大洋俯冲消减作用应早于535Ma。     

从洋壳俯冲到陆壳俯冲和碰撞：来自羌塘中西部地区榴辉岩和蓝片岩地球化学的证据

羌塘中部晚三叠世低温/高压变质带是目前青藏高原内部延伸规模最大的高压变质带,但大量关键高压变质岩石出露地区地球化学资料匮乏,严重制约了对高压变质带原岩建造以及构造演化的全面认识。本文以羌塘中西部地区尚无地球化学资料的果干加年山榴辉岩和红脊山蓝片岩为研究对象,进行了系统的地球化学以及原岩恢复工作。研究表明,果干加年山榴辉岩呈透镜状产于围岩石榴石多硅白云母片岩和少量大理岩中,其原岩为亚碱性玄武岩,具有较低的稀土总量(∑REE=51.19×10-6~59.43×10-6)和轻稀土亏损的特征[(La/Yb)N=0.59~0.70],不具有Nb、Ta、Ti的亏损,与典型的N-MORB特征一致,暗示其原岩可能来源于亏损的地幔源区,形成于洋中脊环境。红脊山地区基性蓝片岩的原岩为碱性玄武岩-亚碱性玄武岩,具有高的TiO2(2.97%~4.14%)和P2O5(0.29%~0.48%)含量,富集轻稀土元素[(La/Yb)N=6.10~11.6]和高场强元素,地球化学特征类似于OIB。但是这些基性蓝片岩与大量的陆源碎屑岩伴生产出,且具有明显的硅铝质上地壳物质混染的特征,与南羌塘地区二叠纪大陆板内基性岩墙的产出特征以及地壳混染特征一致,可能是其俯冲消减的产物。通过本文研究结果并结合区域内已识别出的E-MORB型洋壳和洋岛/海山物质深俯冲的证据,我们认为羌塘中部晚三叠世高压变质带以洋壳物质深俯冲为主,同时亦保留了部分陆壳物质俯冲的证据,暗示大洋向北俯冲消减结束之后,又牵引至少一部分南羌塘北缘陆壳物质经历了随后的俯冲过程。     

黄河流域城市群扩张的时空格局演化及情景模拟——以呼包鄂榆城市群为例

黄河流域城镇扩张对区域景观格局影响显著,城市群人口聚集与增长引发了流域“人-地”矛盾和“空间冲突”等一系列生态环境问题。基于土地利用数据和FLUS模型对2025年和2035年呼包鄂榆城市群城镇化与土地利用时空演化特征进行多情景模拟预测。结果表明：① 1990—2018年呼包鄂榆城市群整体发展水平较低,建设用地面积经历了“平稳增加-缓慢增加-急剧增加”的变化过程,区域总体以草地为主,其占土地总面积的50%以上,其次是未利用土地和耕地,林地和建设用地次之。② 城市群扩张最剧烈地区在空间上主要发生在呼和浩特市、包头市等城市主城区,且扩张模式以外延式扩张为主,扩张来源主要是耕地、草地等生态用地。③ 三种情景模拟发现,2025年和2035年区域土地利用变化的空间结构和特征差异明显。自然发展情景下,城市扩张不受约束,高速增长占据了大量生态用地;加入生态约束条件很好的控制了对草地和林地的占用;经济发展情景下,城市扩张将进一步占据更多的未利用土地和耕地。本研究通过城市群扩张时空格局演化及情景模拟分析,尝试为区域规划、城市空间规划和区域生态空间保护提供多角度、多情景和可选择的政策决策参考。     

Multi-Conjugate Adaptive Optics with laser guide stars: performance in the infrared and visible

The assumption of the Gaussianity of primordial perturbations plays an important role in modern cosmology. The most direct test of this hypothesis consists of testing the Gaussianity of cosmic microwave background (CMB) maps. Counting the pixels with the temperatures in given ranges and thus estimating the one-point probability function of the field is the simplest of all the tests. Other usually more complex tests of Gaussianity generally use a great deal of the information already contained in the probability function. However, the most interesting outcome of such a test would be the signal of non-Gaussianity independent of the probability function. It is shown that the independent information has purely morphological character i.e. it depends on the geometry and topology of the level contours only. As an example we discuss in detail the quadratic model   v = u + α ( u ²-1)  ( u is a Gaussian field with   u¯ =0  and  〈 u ²〉=1  , α is a parameter) that may arise in slow-roll or two-field inflation models. We show that in the limit of small amplitude α the full information about the non-Gaussianity is contained in the probability function. If other tests are performed on this model they simply recycle the same information. A simple procedure allowing us to assess the sensitivity of any statistics to the morphological information is suggested. We provide an analytic estimate of the statistical limit for detecting the quadratic non-Gaussianity α _c as a function of the map size in the ideal situation when the scale of the field is resolved. This estimate is in a good agreement with the results of the Monte Carlo simulations of 256² and 1024² maps. The effect of resolution on the detection quadratic non-Gaussianity is also briefly discussed.     

基于周期项参数模型的GPS流动观测改正方法研究

利用GPS流动观测站200 km范围内站点的连续观测时间序列提高GPS流动观测的定位精度,探讨GPS流动观测结果的可靠性。结果表明：1）流动观测站点速度测量值与连续观测结果差值标准差在改正后平均减小37.97%;2）TAIN、CASH两个站点的同址观测实测值改正后偏离值平均减小38.9%。     

Extreme value modelling of daily areal rainfall over Mediterranean catchments in a changing climate

Heavy rainfall events during the fall season are causing extended damages in Mediterranean catchments. A peaks‐over‐threshold model is developed for the extreme daily areal rainfall occurrence and magnitude in fall over six catchments in Southern France. The main driver of the heavy rainfall events observed in this region is the humidity flux (FHUM) from the Mediterranean Sea. Reanalysis data are used to compute the daily FHUM during the period 1958–2008, to be included as a covariate in the model parameters. Results indicate that the introduction of FHUM as a covariate can improve the modelling of extreme areal precipitation. The seasonal average of FHUM can improve the modelling of the seasonal occurrences of heavy rainfall events, whereas daily FHUM values can improve the modelling of the events magnitudes. In addition, an ensemble of simulations produced by five different general circulation models are considered to compute FHUM in future climate with the emission scenario A1B and hence to evaluate the effect of climate change on the heavy rainfall distribution in the selected catchments. This ensemble of climate models allows the evaluation of the uncertainties in climate projections. By comparison to the reference period 1960–1990, all models project an amplification of the mean seasonal FHUM from the Mediterranean Sea for the projection period 2070–2099, on average by +22%. This increase in FHUM leads to an increase in the number of heavy rainfall events, from an average of 2.55 events during the fall season in present climate to 3.57 events projected for the period 2070–2099. However, the projected changes have limited effects on the magnitude of extreme events, with only a 5% increase in the median of the 100‐year quantiles. Copyright © 2011 John Wiley & Sons, Ltd.     

Floodplain hydrology of the Mekong Delta,Vietnam

The Mekong Delta is one of the largest and most intensively used estuaries in the world. Each year it witnesses widespread flooding which is both the basis of the livelihood for more than 17 million people but also the major hazard. Therefore, a thorough understanding of the hydrologic and hydraulic features is urgently required for various planning purposes. While the general causes and characteristics of the annual floods are understood, the inundation dynamics in the floodplains in Vietnam which are highly controlled by dikes and other control structures have not been investigated in depth. Especially, quantitative analyses are lacking, mainly due to scarce data about the inundation processes in the floodplains. Therefore, a comprehensive monitoring scheme for channel and floodplain inundation was established in a study area in the Plain of Reeds in the northeastern part of the Vietnamese Delta. This in situ data collection was complemented by a series of high‐resolution inundation maps derived from the TerraSAR‐X satellite for the flood seasons 2008 and 2009. Hence, the inundation dynamics in the channels and floodplains, and the interaction between channels and floodplains, could be quantified for the first time. The study identifies the strong human interference which is governed by flood protection levels, cropping patterns and communal water management. In addition, we examine the tidal influence on the inundation in various parts of the Delta, since it is expected that climate change‐induced sea level rise will increase the tidal contribution to floodplain inundation. Copyright © 2011 John Wiley & Sons, Ltd.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.