The Wabar impact craters, Saudi Arabia, revisited

The very young Wabar craters formed by impact of an iron meteorite and are known to the scientific community since 1933. We describe field observations made during a visit to the Wabar impact site, provide analytical data on the material collected, and combine these data with poorly known information discovered during the recovery of the largest meteorites. During our visit in March 2008, only two craters (Philby-B and 11 m) were visible; Philby-A was completely covered by sand. Mapping of the ejecta field showed that the outcrops are strongly changing over time. Combining information from different visitors with our own and satellite images, we estimate that the large seif dunes over the impact site migrate by approximately 1.0–2.0 m yr␣1 southward. Shock lithification took place even at the smallest, 11 m crater, but planar fractures (PFs) and undecorated planar deformation features (PDFs), as well as coesite and stishovite, have only been found in shock-lithified material from the two larger craters. Shock-lithified dune sand material shows perfectly preserved sedimentary structures including cross-bedding and animal burrows as well as postimpact structures such as open fractures perpendicular to the bedding, slickensides, and radiating striation resembling shatter cones. The composition of all impact melt glasses can be explained as mixtures of aeolian sand and iron meteorite. We observed a partial decoupling of Fe and Ni in the black impact glass, probably due to partitioning of Ni into unoxidized metal droplets. The absence of a Ca-enriched component demonstrates that the craters did not penetrate the bedrock below the sand sheet, which has an estimated thickness of 20–30 m.

Core formation and mantle differentiation on Mars

Geochemical investigation of Martian meteorites (SNC meteorites) yields important constraints on the chemical and geodynamical evolution of Mars. These samples may not be representative of the whole of Mars; however, they provide constraints on the early differentiation processes on Mars. The bulk composition of Martian samples implies the presence of a metallic core that formed concurrently as the planet accreted. The strong depletion of highly siderophile elements in the Martian mantle is only possible if Mars had a large scale magma ocean early in its history allowing efficient separation of a metallic melt from molten silicate. The solidification of the magma ocean created chemical heterogeneities whose ancient origin is manifested in the heterogeneous 142Nd and 182W abundances observed in different meteorite groups derived from Mars. The isotope anomalies measured in SNC meteorites imply major chemical fractionation within the Martian mantle during the life time of the short-lived isotopes 146Sm and 182Hf. The Hf-W data are consistent with very rapid accretion of Mars within a few million years or, alternatively, a more protracted accretion history involving several large impacts and incomplete metal-silicate equilibration during core formation. In contrast to Earth early-formed chemical heterogeneities are still preserved on Mars, albeit slightly modified by mixing processes. The preservation of such ancient chemical differences is only possible if Mars did not undergo efficient whole mantle convection or vigorous plate tectonic style processes after the first few tens of millions of years of its history.

Petrogenesis of rift-related tephrites, phonolites and trachytes (Central European Volcanic Province, Rhon, FRG): Constraints from Sr, Nd, Pb and O isotopes

The volcanic rocks of the Rhön area (Central European Volcanic Province, Germany) belong to a moderately alkali basaltic suite that is associated with minor tephriphonolites, phonotephrites, tephrites, phonolites and trachytes. Based on isotope sytematics (87Sr/86Sr: 0.7033–0.7042; 143Nd/144Nd: 0.51279–0.51287; 206Pb/204Pb: 19.1–19.5), the inferred parental magmas formed by variable degrees of partial melting of a common asthenospheric mantle source (EAR: European Asthenospheric Reservoir of Cebriá and Wilson, 1995). Tephrites, tephriphonolites, phonotephrites, phonolites and trachytes show depletions and enrichments in some trace elements (Sr, Ba, Nb, Zr, Y) indicating that they were generated by broadly similar differentiation processes that were dominated by fractionation of olivine, clinopyroxene, amphibole, apatite and titaniferous magnetite ± plagioclase ± alkalifeldspar. The fractionated samples seem to have evolved by two distinct processes. One is characterized by pure fractional crystallization indicated by increasing Nb (and other incompatible trace element) concentrations at virtually constant 143Nd/144Nd ~ 0.51280 and 87Sr/86Sr ~ 0.7035. The other process involved an assimilation–fractional crystallization (AFC) process where moderate assimilation to crystallization rates produced evolved magmas characterized by higher Nb concentrations at slightly lower 143Nd/144Nd down to 0.51275. Literature data for some of the evolved rocks show more variable 87Sr/86Sr ranging from 0.7037 to 0.7089 at constant 143Nd/144Nd ~ 0.51280. These features may result from assimilation of upper crustal rocks by highly differentiated low-Sr (< 100 ppm Sr) lavas. However, based on the displacement of the differentiated rocks from this study towards lower 143Nd/144Nd ratios and modeled AFC processes in 143Nd/144Nd vs. 87Sr/86Sr and 207Pb/204Pb vs. 143Nd/144Nd space assimilation of lower crustal rocks seems more likely. The view that assimilation of lower crustal rocks played a role is confirmed by high-precision double-spike Pb isotope data that reveal higher 207Pb/204Pb ratios (15.62–15.63) in the differentiated rocks than in the primitive basanites (15.58–15.61). This is compatible with incorporation of radiogenic Pb from lower crustal xenoliths (207Pb/204Pb: 15.63–15.69) into the melt. However, 206Pb/204Pb ratios are similar for the differentiated rocks (19.13–19.35) and the primitive basanites (19.12–19.55) implying that assimilation involved an ancient crustal end member with a higher U/Pb ratio than the mantle source of the basanites. In addition, alteration-corrected δ18O values of the differentiated rocks range from c. 5 to 7‰ which is the same range as observed in the primitive alkaline rocks. This study confirms previous interpretations that highlighted the role of AFC processes in the evolution of alkaline volcanic rocks in the Rhön area of the Central European Volcanic Province.

Net accumulation rates derived from ice core stable isotope records of Pío XI glacier, Southern Patagonia Icefield

Pío XI, the largest glacier of the Southern Patagonia Icefield, reached its neoglacial maximum extent in 1994 and is one of the few glaciers in that area which is not retreating. In view of the recent warming it is important to understand glacier responses to climate changes. Due to its remoteness and the harsh conditions in Patagonia, no systematic mass balance studies have been performed. In this study we derived net accumulation rates for the period 2000–2006 from a 50 m (33.2 4 m weq) ice core collected in the accumulation area of Pío XI (2600 m a.s.l., 49°16'40"S, 73°21'14"W). Borehole temperatures indicate near temperate ice, but the average melt percent is only 16 ± 14%. Records of stable isotopes are well preserved and were used for identification of annual layers. Net accumulation rates range from 3.4–7.1 water equivalent (m weq) with an average of 5.8 m weq, comparable to precipitation amounts at the Chilean coast, but not as high as expected for the Icefield. Ice core stable isotope data correlate well with upper air temperatures and may be used as temperature proxy.

Holtite and Dumortierite from the Szklary Pegmatite, Lower Silesia, Poland

The Szklary holtite is represented by three compositional varieties: (I) Ta-bearing (up to 14.66 wt.% Ta(2)O(5)), which forms homogeneous crystals and cores within zoned crystals; (2) Ti-bearing (up to 3.82 wt.% TiO(2)), found as small domains within the core; and (3) Nb-bearing (up to 5.30 wt.% Nb(2)O(5),) forming the rims of zoned crystals. All three varieties show variable Sb+As content, reaching 19.18 wt.% Sb(2)O(3) (0.87 Sb a.p.f.u.) and 3.30 wt.% As(2)O(3) (0.22 As a.p.f.u.) in zoned Ta-bearing holtite, which constitutes the largest Sb+As content reported for the mineral. The zoning in holtite is a result of Ta-Nb fractionation in the parental pegmatite-forming melt together with contamination of the relatively thin Szklary dyke by Fe, Mg and Ti. Holtite and the As- and Sb-bearing dumortierite, which in places overgrows the youngest Nb-bearing zone, suggest the following crystallization sequence: Ta-bearing holtite -> Ti-bearing holtite -> Nb-bearing holtite -> As- and Sb-bearing, (Ta,Nb,Ti)-poor dumortierite -> As- and Sb-dominant, (Ta,Nb,Ti)-free dumortierite-like mineral (16.81 wt.% As(2)O(3) and 10.23 wt.% Sb(2)O(3)) with (As+Sb) > Si. The last phase is potentially a new mineral species, Al(6)rectangle B(Sb,As)(3)O(15). or Al(5)rectangle(2)B(Sb,As)(3)O(12)(OH)(3), belonging to the dumortierite group. The Szklary holtite shows no evidence of clustering of compositions around 'holtite I' and 'holtite II'. Instead, the substitutions of Si(4+) by Sb(3+)+As(3+) at the Si/Sb sites and of Ta(5+) by Nb(5+) or Ti(4+) at the Al(l) site suggest possible solid solutions between: (1) (Sb,As)-poor and (Sb,As)-rich holtite; (2) dumortierite and the unnamed (As+Sb)-dominant dumortierite-like mineral; and (3) Ti-bearing dumortierite and holtite, i.e. our data provide further evidence for miscibility between holtite and dumortierite, but leave open the question of defining the distinction between them. The Szklary holtite crystallized from the melt along with other primary Ta-Nb-(Ti) minerals such as columbite-(Mn), tantalite-(Mn), stibiotantalite and stibiocolumbite as the availability of Ta decreased. The origin of the parental melt can be related to anatexis in the adjacent Sowie Mountains complex, leading to widespread migmatization and metamorphic segregation in pelitic-psammitic sediments metamorphosed at similar to 390-380 Ma.

Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

Recent Increase in Black Carbon Concentrations from a Mt. Everest Ice Core Spanning 1860-2000 AD

A Mt. Everest ice core spanning 1860-2000 AD and analyzed at high resolution for black carbon (BC) using a Single Particle Soot Photometer (SP2) demonstrates strong seasonality, with peak concentrations during the winter-spring, and low concentrations during the summer monsoon season. BC concentrations from 1975-2000 relative to 1860-1975 have increased approximately threefold, indicating that BC from anthropogenic sources is being transported to high elevation regions of the Himalaya. The timing of the increase in BC is consistent with BC emission inventory data from South Asia and the Middle East, however since 1990 the ice core BC record does not indicate continually increasing BC concentrations. The Everest BC and dust records provide information about absorbing impurities that can contribute to glacier melt by reducing the albedo of snow and ice. There is no increasing trend in dust concentrations since 1860, and estimated surface radiative forcing due to BC in snow exceeds that of dust in snow. This suggests that a reduction in BC emissions may be an effective means to reduce the effect of absorbing impurities on snow albedo and melt, which affects Himalayan glaciers and the availability of water resources in major Asian rivers. Citation: Kaspari, S. D., M. Schwikowski, M. Gysel, M. G. Flanner, S. Kang, S. Hou, and P. A. Mayewski (2011), Recent increase in black carbon concentrations from a Mt. Everest ice core spanning 1860-2000 AD, Geophys. Res. Lett., 38, L04703, doi: 10.1029/2010GL046096.

Elastic Uplift in Southeast Greenland Due to Rapid Ice Mass Loss

The rapid unloading of ice from the southeastern sector of the Greenland ice sheet between 2001 and 2006 caused an elastic uplift of similar to 35 mm at a GPS site in Kulusuk. Most of the uplift results from ice dynamic-induced volume losses on two nearby outlet glaciers. Volume loss from Helheim Glacier, calculated from sequential digital elevation models, contributes about similar to 16 mm of the observed uplift, with an additional similar to 5 mm from volume loss of Kangerdlugssuaq Glacier. The remaining uplift signal is attributed to significant melt-induced ice volume loss from the ice sheet margin along the southeast coast between 62 degrees N and 66 degrees N.

Summer Temperature Trend Over the Past Two Millennia Using Air Content in Himalayan Ice

Two Himalayan ice cores display a factor-two decreasing trend of air content over the past two millennia, in contrast to the relatively stable values in Greenland and Antarctica ice cores over the same period. Because the air content can be related with the relative frequency and intensity of melt phenomena, its variations along the Himalayan ice cores provide an indication of summer temperature trend. Our reconstruction point toward an unprecedented warming trend in the 20th century but does not depict the usual trends associated with "Medieval Warm Period" (MWP), or "Little Ice Age" (LIA).

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

Evidence for a Frozen Bed, Byrd Glacier, Antarctica

Ice thickness, computed within the fjord region of Byrd Glacier on the assumptions that Byrd Glacier is in mass-balance equilibrium and that ice velocity is entirely due to basal sliding, are on average 400 m less than measured ice thicknesses along a radio-echo profile. We consider four explanations for these differences: (1) active glacier ice is separated from a zone of stagnant ice near the base of the glacier by a shear zone at depth; (2) basal melting rates are some 8 m/yr; (3) internal shear occurs with no basal sliding in much of the region above the grounding zone; or (4) internal creep and basal sliding contribute to the flow velocity in varying proportions above the grounding zone. Large gradients of surface strain rate seem to invalidate the first explanation. Computed values of basal shear stress (140 to 200 kPa) provide insufficient frictional heat to melt the ice demanded by the second explanation. Both the third and fourth explanations were examined by making simplifying assumptions that prevented a truly quantitative evaluation of their merit. Nevertheless, there is no escaping the qualitative conclusion that internal shear contributes strongly to surface velocities measured on Byrd Glacier, as is postulated in both these explanations.

On the Disintegration of Ice Shelves: The Role of Fracture

Crevasses can be ignored in studying the dynamics of most glaciers because they are only about 20 m deep, a small fraction of ice thickness. In ice shelves, however, s urface crevasses 20 m deep often reach sealevel and bottom crevasses can move upward to sea-level (Clough, 1974; Weertman, 1980). The ice shelf is fractured completely through if surface and basal crevasses meet (Barrett, 1975; Hughes, 1979). This is especially likely if surface melt water fills surface crevasses (Weertman, 1973; Pfeffer, 1982; Fastook and Schmidt, 1982). Fracture may therefore play an important role i n the disintegration of ice shelves. Two fracture criteria which can be evaluated experimentally and applied to ice shelves, are presented. Fracture is then examined for the general strain field of an ice shelf and for local strain fields caused by shear rupture alongside ice streams entering the ice shelf, fatigue rupture along ice shelf grounding lines, and buckling up-stream from ice rises. The effect of these fracture patterns on the stability of Antarctic ice shelves and the West Antarctic ice sheet is then discussed.

Seasonal Geophysical Monitoring of Biogenic Gases in a Northern Peatland: Implications for Temporal and Spatial Variability in Free Phase Gas Production Rates

A set of high resolution surface ground penetrating radar (GPR) surveys, combined with elevation rod ( to monitor surface deformation) and gas flux measurements, were used to investigate in situ biogenic gas dynamics within a northern peatland (Caribou Bog, Maine). Gas production rates were directly estimated from the time series of GPR measurements. Spatial variability in gas production was also investigated by comparing two sites with different geological and ecological attributes, showing differences and/or similarities depending on season. One site characterized by thick highly humified peat deposits (5-6 m), wooded heath vegetation and open pools showed large ebullition events during the summer season, with estimated emissions (based on an assumed range of CH(4) concentration) between 100 and 172 g CH(4) m(-2) during a single event. The other site characterized by thinner less humified peat deposits (2-3 m) and shrub vegetation showed much smaller ebullition events during the same season (between 13 and 23 g CH(4) m(-2)). A consistent period of free-phase gas (FPG) accumulation during the fall and winter, enhanced by the frozen surficial peat acting as a confining layer, was followed by a decrease in FPG after the snow/ice melt that released estimated fluxes between 100 and 200 g CH(4) m(-2) from both sites. Estimated FPG production rates during periods of biogenic gas accumulation ranged between 0.22 and 2.00 g CH(4) m(3) d(-1) and reflected strong seasonal and spatial variability associated with differences in temperature, peat soil properties, and/or depositional attributes (e. g., stratigraphy). Periods of decreased atmospheric pressure coincided with short-period increases in biogenic gas flux, including a very rapid decrease in FPG content associated with an ebullition event that released an estimated 39 and 67 g CH(4) m(-2) in less than 3.5 hours. These results provide insights into the spatial and seasonal variability in production and emission of biogenic gases from northern peatlands.

Thermal Conditions at the Bed of the Laurentide Ice Sheet in Maine During Deglaciation: Implications for Esker Formation

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average similar to 5 mm a(-1) during retreat. They decline over the next 100km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of similar to 1.2 m(3) s(-1). While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

Sliding of Ice Past an Obstacle at Engabreen, Norway

At Engabreen, Norway, an instrumented panel containing a decimetric obstacle was mounted flush With the bed surface beneath 210 m of ice. Simultaneous measurements of normal and shear stresses, ice velocity and temperature were obtained as dirty basal ice flowed past the obstacle. Our measurements were broadly consistent with ice thickness, flow conditions and bedrock topography near the site of the experiment. Ice speed 0.45 m above the bed was about 130 mm d(-1), much less than the surface velocity of 800 mm d(-1) Average normal stress on the panel was 1.0-1.6 MPa, smaller than the expected ice overburden pressure. Normal stress was larger and temperature was lower on the stoss side than on the lee side, in accord with flow dynamics and equilibrium thermodynamics. Annual differences in normal stresses were correlated with changes in sliding speed and ice-flow direction. These temporal variations may have been caused by changes in ice rheology associated with changes in sediment concentration, water content or both. Temperature and normal stress were generally correlated, except when clasts presumably collided with the panel. Temperature gradients in the obstacle indicated that regelation was negligible, consistent with the obstacle size. Melt rate was about 10% of the sliding speed. Despite high sliding speed, no significant ice/bed separation was observed in the lee of the obstacle. Frictional forces between sediment particles in the ice and the panel, estimated from Hallet's (1981) model, indicated that friction accounted for about 5% of the measured bed-parallel force. This value is uncertain, as friction theories are largely untested. Some of these findings agree with sliding theories, others do not.