Snow,runoff, and mass balance modeling for the entire Mittivakkat Glacier (1998–2006), Ammassalik Island,SE Greenland

Geografisk Tidsskrift, Danish Journal of Geography 108(1):121–136, 2008

SnowModel, a physically-based snow evolution modeling system that includes four submodels—MicroMet, EnBal, SnowPack, and SnowTran-3D—was used to simulate eight full-year (1998/99 through 2005/06) evolutions of snow accumulation, blowing snow sublimation, evaporation, snow and ice surface melt, runoff, and mass changes on the entire Mittivakkat Glacier (31 km²) in southeast Greenland. Meteorological observations from two meteorological stations inside the glacier catchment were used as model input, and glaciological mass balance observations were used for model calibration (1998/99 through 2001/02) and validation (2002/03 through 2005/06) of winter snow simulations. As confirmed by observations, the spatially modeled end-of-winter snow water equivalent (SWE) accumulation increased with elevation up to 700–800 m a.s.l. in response to elevation, topography, and dominating wind direction, and maximum snow deposition occurred on the lee side of the ridge east and south of the glacier. Simulated end-of-summer cumulative runoff decreased with elevation and minimum runoff occurred on the shadowed side of the ridge east and south of the glacier. The modeled test period averaged annual mass balance was 65 mm w. eq. y^?1 or ～8% more than the observed. For the simulation period, the glacier net mass balance varies from -199 to -1,834 mm w.eq. y^?1, averaging -900 (±470) mm w.eq.y^?1. The glacier averaged annual modeled precipitation ranged from 1,299 to 1,613 mm w.eq. y^?1, evaporation and sublimation from 206 to 289 mm w.eq., and runoff from 1,531 to 2,869 mm w.eq. y^?1. The model simulated Mittivakkat Glacier net loss of900 mm w.eq. y^?1 contributes approximately 42% to the average simulated runoff of 2,140 mm w.eq. y^?1, indicating a mean specific runoff of 67.8 l s^?1 km^?2.     

3D DEM simulation of principal stress rotation in different planes of cross-anisotropic granular materials

Three-dimensional Discrete Element Method simulations have been performed to study the deformation of cross-anisotropic granular materials under principal stress rotation (PSR), for rotation planes oriented at different angles θ with respect to the bedding plane. The simulations have been conducted with a novel technique for applying specified stresses at three-dimensional boundaries. The results are qualitatively in agreement with experimental results from literature. Cumulative volume contraction is always observed under continuous PSR and increases with increasing θ. The dilatancy rate decreases with increasing number of PSR cycles, tending to zero. The noncoaxiality angle between the strain increment and the stress in the PSR plane increases with increasing number of cycles, reaching the same asymptotic value for samples of various densities and for various θ. Periodic oscillations of the dilatancy rate and noncoaxiality angle within each PSR cycle are observed with an increasing oscillation magnitude with increasing θ, due to the larger fabric anisotropy within the PSR plane. When θ = 30 or 60°, significant noncoaxial strain accumulation occurs in the plane perpendicular to the PSR plane due to the oblique angle between the PSR plane and the bedding plane, echoing the major principal fabric direction's being neither parallel nor perpendicular to the PSR plane. The macroscopic behavior of the samples is related to the microscopic parameters including coordination number and fabric anisotropy. With increasing number of cycles, the difference between normalized stress/strain/fabric increment tensors tends to become constant, with only a small lag between each pair, irrespective of θ.     

Continuous monitoring of stream δ18O and δ2H and stormflow hydrograph separation using laser spectrometry in an agricultural catchment

A portable Wavelength Scanned‐Cavity Ring‐Down Spectrometer (Picarro L2120) fitted with a diffusion sampler (DS‐CRDS) was used for the first time to continuously measure δ¹⁸O and δ²H of stream water. The experiment took place during a storm event in a wet tropical agricultural catchment in north‐eastern Australia. At a temporal resolution of one minute, the DS‐CRDS measured 2160 δ¹⁸O and δ²H values continuously over a period of 36 h with a precision of ±0.08 and 0.5‰ for δ¹⁸O and δ²H, respectively. Four main advantages in using high temporal resolution stream δ¹⁸O and δ²H data during a storm event are highlighted from this study. First, they enabled us to separate components of the hydrograph, which was not possible using high temporal resolution electrical conductivity data that represented changes in solute transfers during the storm event rather than physical hydrological processes. The results from the hydrograph separation confirm fast groundwater contribution to the stream, with the first 5 h of increases in stream discharge comprising over 70% pre‐event water. Second, the high temporal resolution stream δ¹⁸O and δ²H data allowed us to detect a short‐lived reversal in stream isotopic values (δ¹⁸O increase by 0.4‰ over 9 min), which was observed immediately after the heavy rainfall period. Third, δ¹⁸O values were used to calculate a time lag of 20 min between the physical and chemical stream responses during the storm event. Finally, the hydrograph separation highlights the role of event waters in the runoff transfers of herbicides and nutrients from this heavily cultivated catchment to the Great Barrier Reef. Copyright © 2015 John Wiley & Sons, Ltd.     

Exotic pollen as an indicator of variable atmospheric circulation over the Labrador Sea region during the mid to late Holocene

Variability in the abundance of exotic (non‐native) pollen in sediment cores has long been considered as a potential proxy for changing atmospheric circulation, but the difficulty of gaining sufficient total exotic pollen and the incomplete understanding of atmospheric pollen transport patterns has hindered its application. In light of recent advances in the study of pollen transport, we present an exotic pollen record from two fjord sediment cores taken from the west (Placentia Bay, Newfoundland) and east (Narsaq Sund, Greenland) Labrador Sea as a basis for studying variations in regional atmospheric circulation. The two cores cover the last ca. 5500 years and indicate a shift in dominant spring/summer air masses at ca. 2000 (southern Greenland) and 3000 cal a BP (Newfoundland) transporting reduced concentrations of pollen from southerly and south‐westerly vegetation zones. This may suggest a shift away from more dominantly zonal atmospheric circulation (a feature of positive North Atlantic Oscillation years) to more frequent meridional circulation. These results support sea ice/sea‐surface temperature proxy reconstructions from Newfoundland, investigated as part of the same project, which also suggest increased winter atmospheric circulation during the early part of the time period studied. In this region, more positive North Atlantic Oscillation years, and therefore more zonal atmospheric circulation, are associated with increased atmospheric circulation in both the winter and the summer seasons. Copyright © 2011 John Wiley & Sons, Ltd.     

Temperature in lowland Danish streams: contemporary patterns,empirical models and future scenarios

Continuous temperature measurements at 11 stream sites in small lowland streams of North Zealand, Denmark over a year showed much higher summer temperatures and lower winter temperatures along the course of the stream with artificial lakes than in the stream without lakes. The influence of lakes was even more prominent in the comparisons of colder lake inlets and warmer outlets and led to the decline of cold‐water and oxygen‐demanding brown trout. Seasonal and daily temperature variations were, as anticipated, dampened by forest cover, groundwater input, input from sewage plants and high downstream discharges. Seasonal variations in daily water temperature could be predicted with high accuracy at all sites by a linear air‐water regression model (r²: 0·903–0·947). The predictions improved in all instances (r²: 0·927–0·964) by a non‐linear logistic regression according to which water temperatures do not fall below freezing and they increase less steeply than air temperatures at high temperatures because of enhanced heat loss from the stream by evaporation and back radiation. The predictions improved slightly (r²: 0·933–0·969) by a multiple regression model which, in addition to air temperature as the main predictor, included solar radiation at un‐shaded sites, relative humidity, precipitation and discharge. Application of the non‐linear logistic model for a warming scenario of 4–5 °C higher air temperatures in Denmark in 2070‐2100 yielded predictions of temperatures rising 1·6–3·0 °C during winter and summer and 4·4–6·0 °C during spring in un‐shaded streams with low groundwater input. Groundwater‐fed springs are expected to follow the increase of mean air temperatures for the region. Great caution should be exercised in these temperature projections because global and regional climate scenarios remain open to discussion. Copyright © 2006 John Wiley & Sons, Ltd.     

The recession of the Inland Ice margin during the Holocene climatic optimum in the Jakobshavn Isfjord area of West Greenland

Recent subsurface mapping of parts of the Greenland Inland Ice margin in the region of Jakobshavn Isbræ indicates that the fjord system in the period of at least 2700–4700 calendar yr B.P. was more ice free than at present, and that the front of the glacier was at least 15 km behind the present position. The ¹⁴C-dating of subfossils brought to the present ice margin fit with the climatic records from ice cores and confirm the favourable conditions for Greenland's first settlers, the Sarqaq people, who arrived in the region about 400 yr ago to find hunting grounds 10–20% larger than the present.     

Gamma-burst studies using long-duration balloon flights in the Arctic

A new detector for accurate localization of the cosmic gamma-burst sources is presented. The detector is intended for satellite applications, but balloon flights will be used in the development phase. The advantages and limitations of Arctic balloon flights are discussed based on a trial flight in August 1979.Paper presented at the Symposium on Cosmic Gamma-Ray Bursts, held at Toulouse, France, 26–29 November, 1979.     

Oxygen isotope and palaeotemperature records from six Greenland ice‐core stations: Camp Century,Dye‐3, GRIP,GISP2, Renland and NorthGRIP

Oxygen isotope variations spanning the last glacial cycle and the Holocene derived from ice‐core records for six sites in Greenland (Camp Century, Dye‐3, GRIP, GISP2, Renland and NorthGRIP) show strong similarities. This suggests that the dominant influence on oxygen isotope variations reflected in the ice‐sheet records was regional climatic change. Differences in detail between the records probably reflect the effects of basal deformation in the ice as well as geographical gradients in atmospheric isotope ratios. Palaeotemperature estimates have been obtained from the records using three approaches: (i) inferences based on the measured relationship between mean annual δ¹⁸O of snow and of mean annual surface temperature over Greenland; (ii) modelled inversion of the borehole temperature profile constrained either by the dated isotopic profile, or (iii) by using Monte Carlo simulation techniques. The third of these approaches was adopted to reconstruct Holocene temperature variations for the Dye 3 and GRIP temperature profiles, which yields remarkably compatible results. A new record of Holocene isotope variations obtained from the NorthGRIP ice‐core matches the GRIP short‐term isotope record, and also shows similar long‐term trends to the Dye‐3 and GRIP inverted temperature data. The NorthGRIP isotope record reflects: (i) a generally stronger isotopic signal than is found in the GRIP record; (ii) several short‐lived temperature fluctuations during the first 1500 yr of the Holocene; (iii) a marked cold event at ca. 8.2 ka (the ‘8.2 ka event’); (iv) optimum temperatures for the Holocene between ca. 8.6 and 4.3 ka, a signal that is 0.6‰ stronger than for the GRIP profile; (v) a clear signal for the Little Ice Age; and (vi) a clear signal of climate warming during the last century. These data suggest that the NorthGRIP stable isotope record responded in a sensitive manner to temperature fluctuations during the Holocene. Copyright © 2001 John Wiley & Sons, Ltd.     

Recent volcano-ice interaction and outburst flooding in a Mars polar cap re-entrant

Formation of chasms in the polar ice caps of Mars has been attributed to meltwater outburst floods, but the cause of melting has remained uncertain. In a cap re-entrant enveloping Abalos Colles, west of Casma Boreale in the north polar cap, we have found possible evidence of recent volcano-ice interaction and outburst flooding. In this paper we demonstrate that these two mechanisms can have acted together to form or expand the Abalos re-entrant. Flat-topped ridges and circular rims protruding above the ice cap surface in the re-entrant apex may be lava ridges and volcano craters, and can have caused melting of 3.3 to 7.7×¹⁰³ km³ of ice. The surrounding cap surface appears to have subsided and the likely volume of missing ice matches the melt estimate. Outburst flooding from this area may have reached peak discharges of 0.3 to according to scour patterns in one of the re-entrant channels. This required ponding of melt water during lava eruption and catastrophic release through a sub- or englacial melt water tunnel, the collapse of which has left a chasm in the ice cap margin. The flood features are geologically recent, and volcano-ice interaction may have occurred within the last 20,000 years.     

Characterization of the reactivity of riverine heterogeneous sediments using a facies-based approach; the Rhine–Meuse delta (The Netherlands)

Large groundwater resources are found in densely populated lowland areas, which consist often of young unconsolidated and reduced sediments. When anthropogenic activities lead to oxygenation of the aquifer, breakdown of the main reduced fractions, i.e. sedimentary organic matter (SOM) and pyrite, could lead to severe groundwater deterioration such as acidification, heavy metal mobilization, and increased hardness. The characterization of the reactive properties of these sediments is important in predicting groundwater deterioration, but is often complicated by the high degree of heterogeneity of these sediments. In this study, the potential reduction capacity (PRC, based on SOM and pyrite content), the potential buffer capacity (PBC, based on carbonate content), potential acidification capacity (PAC, based on the potential acid production by sulfide oxidation), and the measured reduction capacity (MRC) of five facies, which are typical of the riverine sediments in the Rhine–Meuse delta (The Netherlands) were determined. A universal facies-classification model was used to classify the deposits into more homogeneous sub-units based on lithologic and geogenic properties, with a further sub-division into oxic or anoxic redox environment based upon groundwater data and field observations. The bulk chemical data show strong variation across facies for the median values of PRC (186–9093 mmol O₂ kg⁻¹), PBC (17–132 mmol O₂ kg⁻¹), and PAC (36–1530 mmol H⁺ kg⁻¹). The MRC was measured as reactivity to molecular O₂ exposure and was 0.5–567.3 mmol O₂ kg⁻¹. Steady-state oxidation rates were in the wide range of 0.001–10.355 mmol O₂ kg⁻¹ day⁻¹ but were typically about 3–8 times faster in fine facies than in coarse facies. Both the PRC and MRC depend strongly on grain size, but also on the syn/post-depositional environment and redox conditions. The main part of the PRC consists of SOM, but pyrite reactivity is higher than SOM reactivity as shown by the relative depletion of pyrite in oxic subfacies and the preferential oxidation during the oxidation experiments. Some facies are very prone to acidification because the PAC is higher than the PBC, but the oxidation experiments also show that acidification could already start before the PRC is fully exhausted. This study, is one of the few that combines bulk chemical data, groundwater data, and reactivity measurements and shows that a facies-based approach is a practical tool in characterizing the reactivity of heterogeneous deposits.