Groundwater storage variability and annual recharge using well-hydrograph and GRACE satellite data

Most studies using GRACE (Gravity Recovery and Climate Experiment) data for examining water storage anomalies have rich hydrogeological databases. Here, GRACE data are analyzed for southern Mali, Africa, a region with sparse hydrogeological data. GRACE data (2002?C2008) did not overlap with observed groundwater-level data (1982?C2002). Terrestrial water storage from GRACE was corrected for soil moisture using the Global Land Data Assimilation System (GLDAS) model to obtain monthly groundwater storage anomalies and annual net recharge. Historical storage anomalies and net recharge were determined using the water-table fluctuation method for available observation wells. Average annual net recharge averaged 149.1?mm (or 16.4% of annual rainfall) and 149.7?mm (14.8%) from historical water level and GRACE data, respectively. Monthly storage anomaly lows and peaks were observed in May and September, respectively, but have a shift in peak to November using the corrected GRACE data, suggesting that the GLDAS model may poorly predict the timing of soil-water storage in this region. Notwithstanding problems with the GLDAS model, the soil moisture-corrected GRACE data accurately predict the relative timing and magnitude of groundwater-storage changes, suggesting that GRACE data are valuable for identifying long-term regional changes in groundwater storage in areas with sparse hydrogeological data.     

Tectonic relation between northeastern China and the Korean peninsula revealed by interpretation of GRACE satellite gravity data

The major continental blocks in northeastern Asia are the North China block and the South China block, which have collided starting from the Korean peninsula. Geologic and geophysical interpretations reveal a well defined suture zone in northeastern China from Qinling through Dabie to Jiaodong. The discovery of high-pressure metamorphic rocks in the Hongseong area of the Korean peninsula, prominent evidence for the collision zone, indicates extension of the collision zone in northeastern China into the Korean peninsula. Interpretation of the GRACE satellite gravity dataset shows two prominent structural boundaries in the Yellow Sea. One extends from the Jiaodong Belt in eastern China to the Imjingang Belt in the Korean peninsula. The other extends from near Nanjing, eastern China, to Hongseong. Tectonic movement in or near the suture zone may be responsible for seismic activity in the western Korean peninsula and the development of the Yellow Sea sedimentary basin.     

Intracontinental mountain building in Central Asia as inferred from satellite geodetic data

The results of longstanding GPS measurements in the northwestern part of Central Asia are discussed. These results impose certain constraints for modeling of intraplate tectonic processes. In the territory covered by observations, the velocity vectors of recent motions of the Earth’s surface relative to the stable portion of Eurasia decrease northward. The plane field of velocities, which rules out the development of extension zones, indicates the impossibility of the mountain building driven by ascending mantle flows beneath the lithosphere of these regions. The nonuniform spatial distribution of the motions is suggestive of the discrete character of the Earth’s crust and its deformation. The crust is brittle, at least in its upper part, and capable of breaking into blocks. The blocks, which move at different velocities, interact with one another and change their original orientation and position, while experiencing independent deformations. This phenomenon has been exemplified in the Tarim Block and the Tien Shan. Within the limits of the constraints imposed by the GPS measurements, the mechanism of intracontinental mountain building related to the lateral flow of asthenospheric material and to the drag of the overlying lithospheric layers is discussed. This mechanism springs from Argand’s ideas [2, 29] and the plate tectonic concept [10, 23]. The upper-mantle convective flow in the direction of the Indian Plate’s motion was the main cause of the crustal deformation. The detachment of the lithospheric mantle from the Indian Plate approximately 25 Ma ago and its subduction beneath the Himalayas and Tibet, along with simultaneous ascent of the remaining crust and uplift of the Tibetan Plateau, allowed the mantle flow to spread far northward beneath the Asian continent. This process is accompanied by consecutive separation and sinking of the cooling asthenospheric material over the entire area from the Himalayas to Siberia as the subcrustal material cools. As a result, the flow velocity decreases, the roof of the active flow plunges, and the lithosphere becomes thicker. The motion and deformation of the lithospheric layers dragged by deep flow cannot follow the asthenospheric flow strictly, owing to the rigidity of the layers. Therefore, a difference of tangential velocities originates between the flow and the lithosphere, thus giving rise to horizontal shear stresses. These stresses affect the overlying lithospheric layers, including the crustal ones, and bring about their drag and tectonic delamination. Simultaneously, the decreasing velocity in the direction of the mantle flow results in bending of the lithospheric layers that is accompanied by local warping of the crust and its stacking and fragmentation into blocks. The different velocities of block motions lead to their mechanical interactions. This scenario of intracontinental mountain building allows an explanation of the many specific features of tectonic processes and orogeny in within-plate mountainous regions.     

Style of active intraplate deformation from gravity and seismicity data: the Baikal Rift, Asia

To better understand how active deformation localizes within a continental plate in response to extensional and transtensional tectonics, a combined analysis of high-quality gravity (Bouguer anomaly) and seismicity data is presented consisting of about 35000 earthquakes recorded in the Baikal Rift Zone. This approach allows imaging of deformation patterns from the surface down to the Moho. A comparison is made with heat flow variations in order to assess the importance of lithospheric rheology in the style of extensional deformation. Three different rift sectors can be identified. The southwestern rift sector is characterized by strong gravity and topography contrasts marked by two major crustal faults and diffuse seismicity. Heat flow shows locally elevated values, correlated with recent volcanism and negative seismic P-velocity anomalies. Based on earthquake fault plane solutions and on previous stress field inversions, it is proposed that strain decoupling may occur in this area in response to wrench-compressional stress regime imposed by the India–Asia collision. The central sector is characterized by two major seismic belts; the southernmost one corresponds to a single, steeply dipping fault accommodating oblique extension; in the centre of lake Baikal, a second seismic belt is associated with several dip-slip faults and subcrustal thinning at the rift axis in response to orthogonal extension. The northern rift sector is characterized by a wide, low Bouguer anomaly which corresponds to a broad, high topographic dome and seismic belts and swarms. This topography can be explained by lithospheric buoyancy forces possibly linked to anomalous upper mantle. At a more detailed scale, no clear correlation appears between the surficial fault pattern and the gravity signal. As in other continental rifts, it appears that the lithospheric rheology influences extensional basins morphology. However, in the Baikal rift, the inherited structural fabric combined with stress field variations results in oblique rifting tectonics which seem to control the geometry of southern and northeastern rift basins.     

A proto-monsoonal climate in the late Eocene of Southeast Asia: Evidence from a sedimentary record in central Myanmar

The Burma Terrane has yielded some of the earliest pieces of evidence for monsoonal rainfall in the Bay of Bengal. However, Burmese ecosystems and their potential monsoonal imprint remain poorly studied. This study focuses on the late Eocene Yaw Formation (23° N) in central Myanmar, which was located near the equator (c. 5° N) during the Eocene. We quantitatively assessed the past vegetation, climate, and depositional environments with sporomorph diagrams, bioclimatic analysis, and sequence biostratigraphy. We calculated the palynological diversity and drew inferences with rarefaction analysis by comparing with four other middle to late Eocene tropical palynofloras. Palynological results highlight a high floristic diversity for the palynoflora throughout the section formed by six pollen zones characterized by different vegetation. They indicate that lowland evergreen forests and swamps dominated in the Eocene Burmese deltaic plains while terra firma areas were occupied by seasonal evergreen, seasonally dry, and deciduous forests. This vegetation pattern is typical to what is found around the Bay of Bengal today and supports a monsoon-like climate at the time of the Yaw Formation. Bioclimatic analysis further suggests that in the late Eocene, the Yaw Formation was more seasonal, drier, and cooler compared to modern-day climate at similar near-equatorial latitude. More seasonal and drier conditions can be explained by a well-marked seasonal migration of the Intertropical Convergence Zone (ITCZ), driver of proto-monsoonal rainfall. Cooler temperatures in the late Eocene of central Myanmar may be due to the lack of adequate modern analogues for the Eocene monsoonal climate, while those found at other three Eocene Asian paleobotanical sites (India and South China) may be caused by the effect of canopy evapotranspirational cooling. Our data suggest that paleoenvironmental change including two transgressive–regressive depositional sequences is controlled by global sea level change, which may be driven by climate change and tectonics. The high diversity of the Yaw Formation palynoflora, despite well-marked seasonality, is explained by its crossroads location for plant dispersals between India and Asia.     

Mapping groundwater storage variations with GRACE: a case study in Alberta,Canada

The applicability of the Gravity Recovery and Climate Experiment (GRACE) to adequately represent broad-scale patterns of groundwater storage (GWS) variations and observed trends in groundwater-monitoring well levels (GWWL) is examined in the Canadian province of Alberta. GWS variations are derived over Alberta for the period 2002–2014 using the Release 05 (RL05) monthly GRACE gravity models and the Global Land Data Assimilation System (GLDAS) land-surface models. Twelve mean monthly GWS variation maps are generated from the 139 monthly GWS variation grids to characterize the annual GWS variation pattern. These maps show that, overall, GWS increases from February to June, and decreases from July to October, and slightly increases from November to December. For 2002–2014, the GWS showed a positive trend which increases from west to east with a mean value of 12 mm/year over the province. The resulting GWS variations are validated using GWWLs in the province. For the purpose of validation, a GRACE total water storage (TWS)-based correlation criterion is introduced to identify groundwater wells which adequately represent the regional GWS variations. GWWLs at 36 wells were found to correlate with both the GRACE TWS and GWS variations. A factor f is defined to up-scale the GWWL variations at the identified wells to the GRACE-scale GWS variations. It is concluded that the GWS variations can be mapped by GRACE and the GLDAS models in some situations, thus demonstrating the conditions where GWS variations can be detected by GRACE in Alberta.     

Crustal structure under the central High Atlas Mountains (Morocco) from geological and gravity data

Seismic wide angle and receiver function results together with geological data have been used as constraints to build a gravity-based crustal model of the central High Atlas of Morocco. Integration of a newly acquired set of gravity values with public data allowed us to undertake 2–2.5D gravity modelling along two profiles that cross the entire mountain chain. Modelling suggests moderate crustal thickening, and a general state of Airy isostatic undercompensation. Localized thickening appears restricted to the vicinity of a north-dipping crustal-scale thrust fault, that offsets the Moho discontinuity and defines a small crustal root which accounts for the minimum Bouguer gravity anomaly values. Gravity modelling indicates that this root has a northeasterly strike, slightly oblique to the ENE general orientation of the High Atlas belt. A consequence of the obliquity between the High Atlas borders and its internal and deep structure is the lack of correlation between Bouguer gravity anomaly values and topography. Active buckling affecting the crust, a highly elevated asthenosphere, or a combination of both are addressed as side mechanisms that help to maintain the high elevations of the Atlas mountains.     

Formation of diapiric structure in the deformation zone,central Indian Ocean: A model from gravity and seismic reflection data

Analyses of bathymetry, gravity and seismic reflection data of the diffusive plate boundary in the central Indian Ocean reveal a new kind of deformed structure besides the well-reported structures of long-wavelength anticlinal basement rises and high-angle reverse faults. The structure (basement trough) has a length of about 150 km and deepens by up to 1 km from its regional trend (northward dipping). The basement trough includes a rise at its center with a height of about 1.5km. The rise is about 10 km wide with rounded upper surface and bounded by vertical faults. A broad freeair gravity low of about 20 mGal and a local high of 8 mGal in its center are associated with the identified basement trough and rise structure respectively. Seismic results reveal that the horizontal crustal compression prevailing in the diffusive plate boundary might have formed the basement trough possibly in early Pliocene time. Differential loading stresses have been generated from unequal crust/sediment thickness on lower crustal and upper mantle rocks. A thin semi-ductile serpentinite layer existing near the base of the crust that is interpreted to have been formed at mid-ocean ridge and become part of the lithosphere, may have responded to the downward loading stresses generated by the sediments and crustal rocks to inject the serpentinites into the overlying strata to form a classic diapiric structure.     

Comparisons of drought variability between central High Asia and monsoonal Asia: Inferred from tree rings

Severe drought is a serious natural disaster that frequently strikes East Asia, highlighting the need to understanding its drought regime and the associations with Asian monsoon. Tree-ring-based drought reconstructions provide invaluable paleoclimatic archives for detecting regional and large-scale drought variability and their potential forcings. We herein reviewed many drought reconstructions from central High Asia and monsoonal Asia and compared their similarities and differences, as well as their linkages to Asia monsoon. We compared the decadal-scale variability of six drought reconstructions for the central High Asia, where differing drought variations were found between the western and eastern portions. Seven drought reconstructions were reviewed for monsoonal Asia, from which a difference in drought variability was observed between the northern and southern parts. Therefore, we compared the drought variations of the four sub-regions of western and eastern parts of central High Asia, as well as northern and southern portions of monsoonal Asia. ENSO activity and sea surface temperature of western Pacific and northern Indian Oceans, coupled with Asian monsoon, play an important role in modulating drought variability of much area of the East Asia. An improved denser multi-index tree-ring network of longer length for East Asia is necessary for the establishment of more reliable large-scale drought reconstruction.     

Basin-scale structure control of Carlin-style gold deposits in central Southwestern Guizhou,China: Insights from seismic reflection profiles and gravity data

Carlin-style gold deposits are widely distributed in Southwestern Guizhou, China. Research has dominantly focused on deposit-scale geology and structural control of mineralization. The relationship between the gold deposits and regional structures, in particular their control on the formation and distribution of those deposits is less well understood. Here we use seismic survey and gravity data to reveal the basin-scale structure control of the gold deposits in central Southwestern Guizhou. The seismic reflection profiles show strong reflection horizons including Upper Permian coal measures, and Dongwu and Guangxi unconformities. The seismic data provide important evidence of the basinal structural style and the relationship between fault-related folding and the Carlin-style gold deposits. Regional gravity data also revealed several paleo-uplifts, which appear to control the distribution of Carlin-style gold deposits.The Carlin-style gold deposits are spatially and temporally associated with fault-related folds. Based on their location, we classify the gold deposits into two types, (1) low-angle thrust-controlled deposits; and (2) high-angle reverse fault-related deposits. In the low-angle thrust-controlled deposits, main thrusts are listric, strike NW-SE and dip southwest. These thrusts detach into coal measures, tuffs and volcanic rocks and are bedding-parallel in the vicinity of the Dongwu unconformity. Gold ore bodies are found in the thrusts and characterized by cutting-layer veins near the surface and bedding-parallel veins at depth. The high-angle reverse fault-related deposits are located in the hinge of northeast-southwest-striking fault-related folds, with high-angle reverse faults near the hinge of folds crosscutting the Upper Permian coal measures, volcanic rocks and Dongwu unconformity, and then into the limestone of Middle Permian Maokou Formation. Ore bodies are bedding-parallel in anticline culminations where interlayer fracture zones in Upper Permian coal measures and volcanic rocks are generated by high-angle reverse faulting.We consider that the formation of the gold deposits is part of the tectonic evolution in the area. The two sets of NW-SE- and NE-SW-striking fault-related folds were successively formed during the closure of eastern Tethys Ocean in Late Triassic, and major gold deposits formed in the fault-related folds. The Dongwu unconformity and faults provided conduits for ore fluid flow, and the coal measures and the tuffs of Dachang Member acted as detachment layers and hosts for ore deposition.     

Seafloor spreading in the Weddell Sea from magnetic and gravity data

A re-compilation of magnetic data in the Weddell Sea is presented and compared with the gravity field recently derived from retracked satellite altimetry. The previously informally named ‘Anomaly-T,’ an east–west trending linear positive magnetic and gravity anomaly lying at about 69°S, forms the southern boundary of the well-known Weddell Sea gravity herringbone. North of Anomaly-T, three major E–W linear magnetic lows are shown, and identified with anomalies c12r, c21–29(r) and c33r. On the basis of these, and following work by recent investigators, isochrons c13, c18, c20, c21, c30, c33 and c34 are identified and extended into the western Weddell Sea. Similarly, a linear magnetic low lying along the spine of the herringbone is shown and provisionally dated at 93–96 Ma. Anomaly-T is tentatively dated to be M5n, in agreement with recent tectonic models.Although current tectonic models are generally in good agreement to the north of T, to the south interpretations differ. Some plate tectonic models have only proposed essentially north–south spreading in the region, whilst others have suggested that a period of predominantly east–west motion (relative to present Antarctic geographic coordinates) occurred during the mid-Mesozoic spreading between East and West Gondwana. We identify an area immediately to the south of T which appears to be the southerly extent of N–S spreading in the herringbone. Following recent work, the extreme southerly extent of the N–S directed spreading of the herringbone is provisionally dated M9r/M10. In the oldest Weddell Sea, immediately to the north and east of the Antarctic shelf, we see subtle features in both the magnetic and gravity data that are consistent with predominantly N–S spreading in the Weddell Sea during the earliest opening of East and West Gondwana. In between, however, in a small region extending approximately from about 50 km south of T to about 70°S and from approximately 40° to 53°W, the magnetic and gravity data appear to suggest well-correlated linear marine magnetic anomalies (possible isochrons) perpendicular to T, bounded and offset by less well-defined steps and linear lows in the gravity (possible fracture zones). These magnetic and gravity data southwest of T suggest that the crust here may record an E–W spreading episode between the two-plate system of East and West Gondwana prior to the initiation of the three-plate spreading system of South America, Africa and Antarctica. The E–W spreading record to the east of about 35°W would then appear to have been cut off at about M10 time during the establishment of N–S three-plate spreading along the South American–Antarctic Ridge and then subducted under the Scotia Ridge.     

Topography of the Moho undulations in France from gravity data: their age and origin

The complete gravity data set from France and part of the neighboring countries has been analyzed to compute the topography of the Moho undulations. This work is based on an improved filtering technique and an appropriate assumed density contrast between the crust and the upper mantle. Comparison with deep seismic refraction data reveals that this relief map expresses the continuity and geometry of the Moho undulations better than the sparsely distributed seismic refraction data in France. This gravity Moho map, though may not give absolute depths at places, provides a far better correlation with surface geology than the result from other geophysical techniques. Four domains have been recognized: (a) the Alpine domain where all the Moho undulations are concentric with the Alps; (b) the Armorican domain in which all the undulations are north-west/south-east oriented; (c) the Pyrenean domain, in which the undulations are parallel with the Mountain chain; and (d) the Massif Central Domain which does not show clear structural orientation because of the influence of the strong heat flow located at the lower crust/upper mantle interface. Study of the topography and of the superficial structures associated with these undulations reveals that the undulations delineated in the Alpine Domain result from the Tertiary compression which shaped the Alps. The Armorican Domain was first created during the Lower to Middle Cretaceous opening of the Bay of Biscay. It is now slightly affected by the Tertiary to Quaternary closure of this Bay. The Pyrenean Domain was successively shaped by the Lower Cretaceous oblique opening of the Bay of Biscay and by the Upper Cretaceous to Eocene northward displacement of Spain. Comparison between the Moho undulations map and the stress map of France reveals that most of the undulations are perpendicular to the actual shortening directions. This observation suggests that the Mesozoic, Cenozoic and Quaternary stress directions were roughly the same. Massif Central is characterized by the convergence of these three sets of undulations. Its Post-Oligocene uplift was probably the result of the converging stresses recognized in the three surrounding domains. When the Moho undulations and the topography are compared, two types of periodic crustal instabilities can be recognized. One corresponds to the buckling of the crust developed under compression, the other to boudinage which was associated with extension. Both phenomena show a typical wavelength of 200–250 km which is in agreement with the results of the actual physical and numerical modeling currently available.     

Geological applications potential of DEM,ETM+, and gravity data in arid and semi-arid regions with special reference to south Western Desert,Egypt

The study area comprises the southeastern part of the Western Desert of Egypt. The aims of the present study are to delineate the major surface and subsurface structures and to evaluate tectonic framework of the study area. Surface studies of the present work are based mainly on data from two different types of remote sensing systems: 15 cloud-free Landsat Enhanced Thematic Mapper Plus (ETM+) images and Shuttle Radar Topography Mission (SRTM) DEM data. The remote sensing data were used to recognize and extract geological lineaments by calculating and interpreting DEM derivatives, including shaded relief maps, slope maps, and traverse profiles. A Bouguer anomaly map (scale 1:500,000) was used to delineate a detailed picture of the subsurface structure in the study area. The analysis of surface data, including SRTM DEM and ETM+, was combined with subsurface data, including Bouguer anomaly map, and gave new insight into the tectonic and structure patterns of the study area.     

Holocene environments and climate in the Mongolian Altai reconstructed from the Hoton-Nur pollen and diatom records: a step towards better understanding climate dynamics in Central Asia

This study presents the results of the palynological and diatom analyses of the sediment core recovered in Hoton-Nur Lake (48°37′18″N, 88°20′45″E, 2083 m) in 2004. Quantitative reconstruction of the Holocene vegetation and climate dynamics in the semiarid Mongolian Altai suggests that boreal woodland replaced the primarily open landscape of northwestern Mongolia at about 10 kyr BP (1 kyr = 1000 cal yr) in response to a noticeable increase in precipitation from 200–250 mm/yr to 450–550 mm/yr. A decline of the forest vegetation and a return to a predominance of open vegetation types occurred after 5 kyr BP when precipitation sums decreased to 250–300 mm/yr. Prior to 11.5 kyr BP diatom concentrations are relatively low and the lake is dominated by planktonic Cyclotella and small Fragilariaceae, suggesting the existence of a relatively deep and oligotrophic/mesotrophic lake. The great abundance of Staurosirella pinnata from the beginning of the record until 10.7 kyr BP might imply intensified erosion processes in the catchment and this is fully consistent with the presence of scarce and dry vegetation and the generally arid climate during this period. From about 10.7 kyr BP, more planktonic diatom taxa appeared and increased in abundance, indicating that the lake became more productive as diatom concentration increased. This change correlates well with the development of boreal woodland in the catchment. Decrease in precipitation and changes in the vegetation towards steppe are reflected by the rapid increase in Aulacoseira distans from about 5 kyr BP. The Holocene pollen and diatom records do not indicate soil and vegetation cover disturbances by the anthropogenic activities, implying that the main transformations of the regional vegetation occurred as a result of the natural climate change. Our reconstruction is in agreement with the paleomonsoon records from China, demonstrating an abrupt strengthening of the summer monsoon at 12 kyr BP and an associated increase in precipitation and in lake levels between 11 and 8 kyr BP, followed by the stepwise attenuation of the monsoon circulation and climate aridization towards the modern level. The records from the neighboring areas of Kazakhstan and Russia, situated west and north of Hoton-Nur, demonstrate spatially and temporally different Holocene vegetation and climate histories, indicating that the Altai Mountains as a climate boundary are of pivotal importance for the Holocene environmental and, possibly, habitation history of Central Asia.     

Crustal properties from S-wave and gravity data along a seismic refraction profile in Romania

VRANCEA'99 is a seismic refraction line that was carried out in 1999 to investigate the deep structure and physical properties of the upper lithosphere of the southeastern Carpathians and its foreland. It runs from the city of Bacau to the Danube River, traversing the Vrancea epicentral area of strong intermediate-depth seismicity and the city of Bucharest.

Interpretation of P-wave arrivals led to a velocity model that displays a multi-layered crust with velocities increasing with depth. The range of P wave velocities in the sedimentary cover increases from N to S and a structuring of the autochthonous basement of the Moesian Platform is observed. The crystalline crust displays thickness variations, but at the same time the lateral velocity structure along the seismic line remains almost constant. An intra-crustal boundary separates an upper crust from the lower crust. Within the upper mantle a low velocity zone is detected at a depth of about 55-km.

The interpretation of observable S-waves resulted in a velocity model that shows the same multi-layered crust, with S-velocities increasing similarly with depth as the P-waves. The corresponding Poisson's ratio is highly variable throughout the crust and ranges from 0.20–0.35 for the sedimentary cover to 0.22–0.25 for the crystalline crust. The interpretation of the V_p, V_s and Poisson's ratio in petrological terms suggests a large variety of rocks from sand and clay to sandstone, limestone and dolomite within the sedimentary cover. Within the crystalline crust the most probably rock types are granite, granodiorite, granite–gneiss and/or felsic amphibolite–gneiss in the upper part and gneiss and /or amphibolite in the lower part.

Based on the 2-D seismic velocity model, a density model is developed. Density values are assigned to each layer in agreement with the P-wave velocity model and with values accepted for the geological units in the area. After several iterations a good fit between the computed and observed Bouguer anomalies was obtained along the seismic line.     

Evidence for Triassic salt domes in the Tunisian Atlas from gravity and geological data

Detailed gravity data were analyzed to constrain two controversial geological models of evaporitic structures within the Triassic diapiric zone (Triassic massifs of Jebel Debadib and Ben Gasseur) of the northern Tunisian Atlas. Based on surface observations, two geological models have been used to explain the origin of the Triassic evaporitic bodies: (1) salt dome/diapiric structure or (2) a “salt glacier”. The gravity analysis included the construction of a complete Bouguer gravity anomaly map, horizontal gravity gradient (HGG) map and two and a half-dimensional (2.5D) forward models. The complete Bouguer gravity anomaly map shows a prominent negative anomaly over the Triassic evaporite outcrops. The HGG map showed the location of the lateral density changes along northeast structural trends caused by Triassic/Cretaceous lithological differences. The modeling of the complete Bouguer gravity anomaly data favored the diapiric structure as the origin of the evaporitic bodies. The final gravity model constructed over Jebel Debadib indicates that the Triassic evaporitic bodies are thick and deeply rooted involving a dome/diapiric structure and that the Triassic material has pulled upward the younger sediment cover by halokinesis. Taking in account kinematic models and the regional tectonic events affecting the northern margin of Africa, the above diapirs formed during the reactive to active to passive stages of continental margin evolution with development of sinks. Otherwise, this study shows that modeling of detailed gravity data adds useful constraints on the evolution of salt structures that may have an important impact on petroleum exploration models.     

Alkali loss and retention in an iron-rich peralkaline phonolite dyke from the Gardar province,south Greenland

An unusually iron-rich, strongly peralkaline dyke occurring near the Ilimaussaq intrusion, south Greenland, is described. The dyke has developed two facies characterised respectively by the mafic assemblages pyroxene+amphibole+aemigmatite, and pyroxene+magnetite. The differences can be explained as an effect of loss of up to 4% Na₁O from the magnetite-bearing facies (relative to the magnetite-free facies). The effect is in accordance with the phase relations at 1 atm in the synthetic system Na₂OAl₂O₂Fe₂O₂SiO₂. Continued alkali loss and magnetic precipitation will change the differentiation trend from agpaitic with iron enrichment towards miaskitic with iron depletion. It is suggested that the extreme degree of iron enrichment in the dyke resulted from the coincidence of low oxygen fugacity and strong peralkalinity tending to stabilise both ferrous and ferric iron in the liquid.     

Hf isotope compositions and HREE variations in off-craton garnet and spinel peridotite xenoliths from central Asia

Garnet-facies continental mantle is poorly understood because the vast majority of mantle xenoliths in continental basalts are spinel peridotite. Peridotite xenoliths from Vitim (southern Siberia) and Mongolia provide some of the best samples of garnet and garnet-spinel facies off-craton lithospheric mantle. Garnets in those fertile to moderately depleted lherzolites show a surprisingly broad range of HREE abundances, which poorly correlate with modal and major oxide compositions. Some garnets are zoned and have Lu-rich cores. We argue that these features indicate HREE redistribution after the partial melting, possibly related to spinel-garnet phase transition on isobaric cooling. Most peridotites from Vitim have depleted to ultra-depleted Hf isotope compositions (calculated from mineral analyses: ε_Hf(0) = +17 to +45). HREE-rich garnets have the most radiogenic ε_Hf values and plot above the mantle Hf-Nd isotope array while xenoliths with normal HREE abundances usually fall within or near the depleted end of the MORB field. Model Hf isotope ages for the normal peridotites indicate an origin by ancient partial melt extraction from primitive mantle, most likely in the Proterozoic. By contrast, an HREE-rich peridotite yields a Phanerozoic model age, possibly reflecting overprinting of the ancient partial melting record with that related to a recent enrichment in Lu. Clinopyroxene-garnet Lu-Hf isochron ages (31-84 Ma) are higher than the likely eruption age of the host volcanic rocks (∼16 Ma). Garnet-controlled HREE migration during spinel-garnet and garnet-spinel phase transitions may be one explanation for extremely radiogenic ¹⁷⁶Hf/¹⁷⁷Hf reported for some mantle peridotites; it may also contribute to Hf isotope variations in sub-lithospheric source regions of mantle-derived magmas.