Local ice caps in Finderup Land,North Greenland,survived the Holocene Thermal Maximum

Local glaciers and ice caps (GICs) comprise only ~5.4% of the total ice volume, but account for ~14–20% of the current ice loss in Greenland. The glacial history of GICs is not well constrained, however, and little is known about how they reacted to Holocene climate changes. Specifically, in North Greenland, there is limited knowledge about past GIC fluctuations and whether they survived the Holocene Thermal Maximum (HTM, ~8 to 5 ka). In this study, we use proglacial lake records to constrain the ice‐marginal fluctuations of three local ice caps in North Greenland including Flade Isblink, the largest ice cap in Greenland. Additionally, we have radiocarbon dated reworked marine molluscs in Little Ice Age (LIA) moraines adjacent to the Flade Isblink, which reveal when the ice cap was smaller than present. We found that outlet glaciers from Flade Isblink retreated inland of their present extent from ~9.4 to 0.2 cal. ka BP. The proglacial lake records, however, demonstrate that the lakes continued to receive glacial meltwater throughout the entire Holocene. This implies that GICs in Finderup Land survived the HTM. Our results are consistent with other observations from North Greenland but differ from locations in southern Greenland where all records show that the local ice caps at low and intermediate elevations disappeared completely during the HTM. We explain the north–south gradient in glacier response as a result of sensitivity to increased temperature and precipitation. While the increased temperatures during the HTM led to a complete melting of GICs in southern Greenland, GICs remained in North Greenland probably because the melting was counterbalanced by increased precipitation due to a reduction in Arctic sea‐ice extent and/or increased poleward moisture transport.     

Studying Hydrodynamic Instability Using Shock-Tube Experiments

The hydrodynamic instability, which develops on the contact surface between two fluids, has great importance in astrophysical phenomena such as the inhomogeneous density distribution following a supernova event. In this event acceleration waves pass across a material interface and initiate and enhance unstable conditions in which small perturbations grow dramatically. In the present study, an experimental technique aimed at investigating the above-mentioned hydrodynamic instability is presented. The experimental investigation is based on a shock-tube apparatus by which a shock wave is generated and initiates the instability that develops on the contact surface between two gases. The flexibility of the system enables one to vary the initial shape of the contact surface, the shock-wave Mach number, and the density ratio across the contact surface. Three selected sets of shock-tube experiments are presented in order to demonstrate the system capabilities: (1) large-initial amplitudes with low-Mach-number incident shock waves; (2) small-initial amplitudes with moderate-Mach-number incident shock waves; and (3) shock bubble interaction. In the large-amplitude experiments a reduction of the initial velocity with respect to the linear growth prediction was measured. The results were compared to those predicted by a vorticity-deposition model and to previous experiments with moderate- and high-Mach number incident shock waves that were conducted by others. In this case, a reduction of the initial velocity was noted. However, at late times the growth rate had a 1/t behavior as in the small-amplitude low-Mach number case. In the small-amplitude moderate-Mach number shock experiments a reduction from the impulsive theory was noted at the late stages. The passage of a shock wave through a spherical bubble results in the formation of a vortex ring. Simple dimensional analysis shows that the circulation depends linearly on the speed of sound of the surrounding material and on the initial bubble radius.     

Observations of Jupiter at 26.3 MHz using a large array

A 640 element phase-steerable dipole array has been used to make highly sensitive observations of the planet Jupiter during the 1973 apparition. The satellite Io is found to have very little influence at the low flux levels, whereas the definition of sources A and B appears to be relatively flux independent. A two-dimensional analysis of the data in the Jupiter-Io plane has revealed considerable source B activity at low intensities which is not influenced by Io.     

Fracture and bedding plane control on groundwater flow in a chalk aquitard

Shallow groundwater in the northern Negev desert of Israel flows preferentially through a complex system of discontinuities. These discontinuities intersect what would otherwise be a massive, low-conductivity, high-porosity Eocene chalk. Vertical fractures and horizontal bedding planes were observed and mapped along approximately 1,200 m of scanline, 600 m of core and 30 two-dimensional trace planes. A bimodal distribution of size exists for the vertical fractures which occur as both single-layer fractures and multi-layer fractures. A bimodal distribution of log transmissivity was observed from slug tests conducted in packed-off, vertical intervals within the saturated zone. The different flow characteristics between the horizontal bedding planes and vertical-type fractures appear to be the cause of the bimodality. Two distinct conceptual models (discrete fracture network) were developed based on the fracture orientation, size, intensity and transmissivity statistics derived from field data. A correlation between fracture size and hydraulic aperture was established as the basis for calibrating the simulated model transmissivity to the field observations. This method of defining transmissivity statistically based on prior information is shown to be a reasonable and workable alternative to the usual conjecture approach towards defining transmissivity in a fractured-rock environment.     

Elastic properties and stability of coexisting 3T and 2M 1 phengite polytypes

The elastic properties of coexisting natural 3T and 2M ₁ phengite samples (Cima Pal, Sesia Zone; Val Savenca; Western Alps, Italy) with similar chemical compositions have been studied by room temperature–high pressure powder diffraction, using synchrotron radiation on the ID9A beam-line at ESRF (Grenoble, France). The P–V curves have been modelled by the Birch–Murnaghan model; a third-order expansion fitted to the experimental data yields for 3T and 2M ₁ K ₀=60.4(±0.7) GPa, K′=5.79(±0.11) at V ₀=703.8851 ų, and K ₀=57.3(±1.0) GPa, K′=6.97(±0.24) at V ₀=938.8815 ų, respectively. The relative stability of 3T vs. 2M ₁ has been explored as a function of pressure and temperature in terms of configuration and deformation contributions to the Gibbs energy, using the elastic properties determined here and other thermodynamic parameters from earlier investigations. The results presented agree with the hypothesis of stability of the 3T polytype in the high pressure regime.     

Flow and transport predictions during multi-borehole tests in fractured chalk using discrete fracture network models

Multi-borehole pumping and tracer tests on the 10 to 100-m scale were conducted in a fractured chalk aquitard in the Negev Desert, Israel. Outcrop and core fracture surveys, as well as slug tests in packed-off intervals, were carried out at this site to obtain the parameters needed for construction of a stochastic discrete fracture network (DFN). Calibration of stochastic DFNs directly to the multiple borehole test data was inadequate. Instead, two equivalent deterministic DFN flow models were used: the vertical-fractures (VF) model, consisting of only vertical fractures, and the fractures’ intersections (INT) model, consisting of vertical and horizontal fractures with enhanced transmissivity at their intersections. Both models were calibrated against the multi-borehole response of one pumping test and their predictions were tested against three other independent pumping tests. The average accuracies of all transient drawdown predictions of the VF and INT models were 65 and 66%, respectively. In contrast to this equality in average drawdown predictions of both models, the INT model predicted better important breakthrough curve features (e.g., first and peak arrival times), than the VF model. This result is in line with previously assumed channeled flow, derived from analytical analysis of these pumping and tracer tests. Ronit Nativ, deceased, may her memory be blessed.     

Structure and compressibility of synthetic ZnAl2O4 (gahnite) under high-pressure conditions, from synchrotron X-ray powder diffraction

 The structural behavior of synthetic gahnite (ZnAl₂O₄) has been investigated by X-ray powder diffraction at high pressure (0–43 GPa) and room temperature, on the ID9 beamline at ESRF. The equation of state of gahnite has been derived using the models of Birch–Murnaghan, Vinet and Poirier–Tarantola, and the results have been mutually compared (the elastic bulk modulus and its derivatives versus P determined by the third-order Birch–Murnaghan equation of state are K ₀=201.7(±0.9) GPa, K ^′ ₀=7.62(±0.09) and K ^″ ₀=−0.1022 GPa⁻¹ (implied value). The compressibilities of the tetrahedral and octahedral bond lengths [0.00188(8) and 0.00142(5) GPa⁻¹ at P=0, respectively], and the␣polyhedral volume compressibilities of the four-␣and␣sixfold coordination sites [0.0057(2) and 0.0041(2) GPa⁻¹ at P=0, respectively] are discussed. Received: 15 January 2001 / Accepted: 23 April 2001     

Synchrotron X-ray powder diffraction study of natural P2 /n-omphacites at high-pressure conditions

 Synchrotron X-ray powder diffraction experiments at high pressure conditions (0.0001–13 GPa) were performed at ESRF (Grenoble-F), on the beamline ID9, to investigate the bulk elastic properties of natural P2/n-omphacites, with quasi-ideal composition. The monoclinic cell parameters a, b, c and β were determined as a function of pressure, and their compressibility coefficients are 0.00277(7), 0.00313(8), 0.00292(5) and 0.00116(4) GPa⁻¹, respectively. The third-order Birch-Murnaghan equation of state was used to interpolate the experimental P−V data, obtaining K ₀=116.6(±2.5) GPa and K′₀=6.03(±0.60). K ₀ was also determined by means of the axial and angular compressibilities [122.5(±1.7) GPa], and of the finite Lagrangian strain theory [121.5(±1.0) GPa]. The discrepancies on K ₀ are discussed in the light of a comparison between techniques to determine the bulk modulus of crystalline materials from static compression diffraction data. Received: 22 February 2000 / Accepted: 10 July 2000