Response of groundwater chemistry to water deliveries in the lower reaches of Tarim River,Northwest China

In this paper, we analysed the monitored data from nine groundwater-monitoring transects in the lower reaches of Tarim River during the five times of stream water deliveries to the river transect where the stream flow ceased. The results showed that the groundwater depth in the lower reaches of Tarim River rose from −9.30 m before the conveyances to −8.17 and −6.50 m after the first and second conveyances, −5.81 and −6.00 m after the third and fourth the conveyance, and −4.73 m after the fifth. The horizontal extent of groundwater recharge was gradually enlarged along both sides of the channel of conveyance, i.e., from 250 m in width after the first conveyance to 1,050 m away from the channel after the fourth delivery. With the rising groundwater level, the concentrations of major anions Cl⁻, SO₄²⁻ and cations Ca²⁺, Mg²⁺, Na⁺, as well as total dissolved solids (TDS) in groundwater underwent a significant change. The spatial variations in groundwater chemistry indicated that the groundwater chemistry at the transect near Daxihaizi Reservoir changed earlier than that farther from it. In the same transect, the chemical variations were earlier in the monitoring well close to watercourse than that farther away from the stream. In general, the concentration of the major ions and TDS at each monitoring well increased remarkably when the water delivery started, and decreased with the continued water delivery, and then increased once again at the end of the study period. Hence, the whole study period may be divided into three stages: the initial stage, the intermediate stage and the later stage. According to the three stages of groundwater chemistry reaction to water delivery and the relationships between groundwater chemical properties and groundwater depths, we educe that under the situation of water delivery, the optimum groundwater depth in the lower reaches of the Tarim River should be −5 m.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

Depths of Partial Crystallization of H2O-bearing MORB: Phase Equilibria Simulations of Basalts at the MAR near Ascension Island (7 11{degrees}S)

Phase equilibria simulations were performed on naturally quenchedbasaltic glasses to determine crystallization conditions priorto eruption of magmas at the Mid-Atlantic Ridge (MAR) east ofAscension Island (7–11°S). The results indicate thatmid-ocean ridge basalt (MORB) magmas beneath different segmentsof the MAR have crystallized over a wide range of pressures(100–900 MPa). However, each segment seems to have a specificcrystallization history. Nearly isobaric crystallization conditions(100–300 MPa) were obtained for the geochemically enrichedMORB magmas of the central segments, whereas normal (N)-MORBmagmas of the bounding segments are characterized by polybariccrystallization conditions (200–900 MPa). In addition,our results demonstrate close to anhydrous crystallization conditionsof N-MORBs, whereas geochemically enriched MORBs were successfullymodeled in the presence of 0·4–1 wt% H₂O in theparental melts. These estimates are in agreement with direct(Fourier transform IR) measurements of H₂O abundances in basalticglasses and melt inclusions for selected samples. Water contentsdetermined in the parental melts are in the range 0·04–0·09and 0·30–0·55 wt% H₂O for depleted and enrichedMORBs, respectively. Our results are in general agreement (within±200 MPa) with previous approaches used to evaluate pressureestimates in MORB. However, the determination of pre-eruptiveconditions of MORBs, including temperature and water contentin addition to pressure, requires the improvement of magma crystallizationmodels to simulate liquid lines of descent in the presence ofsmall amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

Global Correlations of Ocean Ridge Basalt Chemistry with Axial Depth: a New Perspective

The petrological parameters Na₈ and Fe₈, which are Na₂O andFeO contents in mid-ocean ridge basalt (MORB) melts correctedfor fractionation effects to MgO = 8 wt%, have been widely usedas indicators of the extent and pressure of mantle melting beneathocean ridges. We find that these parameters are unreliable.Fe₈ is used to compute the mantle solidus depth (P_o) and temperature(T_o), and it is the values and range of Fe₈ that have led tothe notion that mantle potential temperature variation of T_P= 250 K is required to explain the global ocean ridge systematics.This interpreted T_P = 250 K range applies to ocean ridges awayfrom ‘hotspots’. We find no convincing evidencethat calculated values for P_o, T_o, and T_P using Fe₈ have anysignificance. We correct for fractionation effect to Mg^# = 0·72,which reveals mostly signals of mantle processes because meltswith Mg^# = 0·72 are in equilibrium with mantle olivineof Fo_89·6 (vs evolved olivine of Fo_{88·1–79·6}in equilibrium with melts of Fe₈). To reveal first-order MORBchemical systematics as a function of ridge axial depth, weaverage out possible effects of spreading rate variation, local-scalemantle source heterogeneity, melting region geometry variation,and dynamic topography on regional and segment scales by usingactual sample depths, regardless of geographical location, withineach of 22 ridge depth intervals of 250 m on a global scale.These depth-interval averages give Fe₇₂ = 7·5–8·5,which would give T_P = 41 K (vs 250 K based on Fe₈) beneathglobal ocean ridges. The lack of Fe₇₂–Si₇₂ and Si₇₂–ridgedepth correlations provides no evidence that MORB melts preservepressure signatures as a function of ridge axial depth. We thusfind no convincing evidence for T_P > 50 K beneath globalocean ridges. The averages have also revealed significantcorrelations of MORB chemistry (e.g. Ti₇₂, Al₇₂, Fe₇₂,Mg₇₂, Ca₇₂, Na₇₂ and Ca₇₂/Al₇₂) with ridge axial depth. Thechemistry–depth correlation points to an intrinsic linkbetween the two. That is, the 5 km global ridge axial reliefand MORB chemistry both result from a common cause: subsolidusmantle compositional variation (vs T_P), which determines themineralogy, lithology and density variations that (1) isostaticallycompensate the 5 km ocean ridge relief and (2) determine thefirst-order MORB compositional variation on a global scale.A progressively more enriched (or less depleted) fertileperidotite source (i.e. high Al₂O₃ and Na₂O, and low CaO/Al₂O₃)beneath deep ridges ensures a greater amount of modal garnet(high Al₂O₃) and higher jadeite/diopside ratios in clinopyroxene(high Na₂O and Al₂O₃, and lower CaO), making a denser mantle,and thus deeper ridges. The dense fertile mantle beneath deepridges retards the rate and restricts the amplitude of the upwelling,reduces the rate and extent of decompression melting, givesway to conductive cooling to a deep level, forces melting tostop at such a deep level, leads to a short melting column,and thus produces less melt and probably a thin magmatic crustrelative to the less dense (more refractory) fertile mantlebeneath shallow ridges. Compositions of primitive MORB meltsresult from the combination of two different, but geneticallyrelated processes: (1) mantle source inheritance and (2) meltingprocess enhancement. The subsolidus mantle compositional variationneeded to explain MORB chemistry and ridge axial depth variationrequires a deep isostatic compensation depth, probably in thetransition zone. Therefore, although ocean ridges are of shalloworigin, their working is largely controlled by deep processesas well as the effect of plate spreading rate variation at shallowlevels. KEY WORDS: mid-ocean ridges; mantle melting; magma differentiation; petrogenesis; MORB chemistry variation; ridge depth variation; global correlations; mantle compositional variation; mantle source density variation; mantle potential temperature variation; isostatic compensation     

Flood hazard delineation combining geomorphological and hydrological methods: an example in the Northern Iberian Peninsula

Flood mapping requires the combination and integration of geomorphological and hydrological-hydraulic methods; however, despite this, there is very little scientific literature that compares and validates both methods. Two types of analysis are addressed in the present article. On the one hand, maps of flood plains have been elaborated using geomorphological evidence and historical flood data in the mountainous area of northwestern Spain, covering an area of more then 232 km² of floodplains. On the other hand, a hydrometeorological model has been developed (Clark semidistributed unit hydrograph) in the Sarria River basin (155 km², NW Spain). This basin is not gauged, hence the model was subjected to a goodness-of-fit test of its parameter (curve number) by means of Monte Carlo simulation. The peak flows obtained by means of the hydrological model were used for hydraulic modeling (one-phase, one-dimensional and steady flow) in a 4 km² urban stretch of the river bed. The delineation of surface areas affected by floods since 1918, as well as those analyzed subsequent to the geomorphological study, reveals a high degree of reliability in the delineation of the flooded areas with frequent recurrence intervals (<50 years). If we compare these flooded surface areas with the estimate obtained by the hydrological-hydraulic method we can see that the latter method overestimates the extent of the surface water by 144% for very frequent recurrence intervals (>10 years) and underestimates it as the recurrence interval increases, by up to 80% less floodplain for exceptional events (>500 years). Finally, a management map is put forth combining the most reliable results available by integrating both methods. Originally presented at the Sixth International Conference on Geomorphology.     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.     

东营凹陷北部沙河街组四段深部储层多重成岩环境及演化模式

主要通过薄片鉴定、扫描电镜观察、包裹体分析、粘土矿物分析、镜质体反射率测试和岩心物性分析等手段,结合构造发育史和有机质热演化史等研究成果,以揭示东营凹陷北带古近系深部碎屑岩储层成岩环境及演化模式为目的展开研究工作。结果表明:东营凹陷北带古近系深部碎屑岩储层存在酸性、碱性和酸性碱性交替等多重成岩环境。酸性成岩环境以碳酸盐矿物溶解、长石溶解蚀变为高岭石并伴生石英次生加大等为标志,碱性成岩环境以石英质颗粒及其次生加大边溶解、长石次生加大和晚期碳酸盐矿物沉淀为标志。研究区古近系深部碎屑岩储层成岩环境由浅至深大致经历了碱性—酸性—酸性碱性交替(局部碱性较强)—碱性—弱碱性的演化过程,并建立了沙河街组四段的成岩演化和储层改造模式。     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.