Showing papers by "Geological Survey of Denmark and Greenland published in 2014"

Extensive liquid meltwater storage in firn within the Greenland ice sheet

Abstract: The accelerating loss of mass from the Greenland ice sheet is a major contribution to current sea level rise. Increased melt water runoff is responsible for half of Greenlands mass loss increase. Surface melt has been increasing in extent and intensity, setting a record for surface area melt and runoff in 2012. The mechanisms and timescales involved in allowing surface melt water to reach the ocean where it can contribute to sea level rise are poorly understood. The potential capacity to store this water in liquid or frozen form in the firn (multi-year snow layer) is significant, and could delay its sea-level contribution. Here we describe direct observation of water within a perennial firn aquifer persisting throughout the winter in the southern ice sheet,where snow accumulation and melt rates are high. This represents a previously unknown storagemode for water within the ice sheet. Ice cores, groundairborne radar and a regional climatemodel are used to estimate aquifer area (70 plue or minus 10 x 10(exp 3) square kilometers ) and water table depth (5-50 m). The perennial firn aquifer represents a new glacier facies to be considered 29 in future ice sheet mass 30 and energy budget calculations.

Policies for agricultural nitrogen management?trends, challenges and prospects for improved efficiency in Denmark

Abstract: With more than 60% of the land farmed, with vulnerable freshwater and marine environments, and with one of the most intensive, export-oriented livestock sectors in the world, the nitrogen (N) pollution pressure from Danish agriculture is severe. Consequently, a series of policy action plans have been implemented since the mid 1980s with significant effects on the surplus, efficiency and environmental loadings of N. This paper reviews the policies and actions taken and their ability to mitigate effects of reactive N (Nr) while maintaining agricultural production. In summary, the average N-surplus has been reduced from approximately 170 kg N ha?1 yr?1 to below 100 kg N ha?1 yr?1 during the past 30 yrs, while the overall N-efficiency for the agricultural sector (crop?+?livestock farming) has increased from around 20?30% to 40?45%, the N-leaching from the field root zone has been halved, and N losses to the aquatic and atmospheric environment have been significantly reduced. This has been achieved through a combination of approaches and measures (ranging from command and control legislation, over market-based regulation and governmental expenditure to information and voluntary action), with specific measures addressing the whole N cascade, in order to improve the quality of ground- and surface waters, and to reduce the deposition to terrestrial natural ecosystems. However, there is still a major challenge in complying with the EU Water Framework and Habitats Directives, calling for new approaches, measures and technologies to mitigate agricultural N losses and control N flows.

State of the Climate in 2013

Abstract: Editors note: For easy download the posted pdf of the State of the Climate for 2013 is a very low-resolution file. A high-resolution copy of the report is available by clicking here. Please be patient as it may take a few minutes for the high-resolution file to download.

Polymorphonuclear Leukocytes Restrict Growth of Pseudomonas aeruginosa in the Lungs of Cystic Fibrosis Patients

Abstract: Cystic fibrosis (CF) patients have increased susceptibility to chronic lung infections by Pseudomonas aeruginosa, but the ecophysiology within the CF lung during infections is poorly understood. The aim of this study was to elucidate the in vivo growth physiology of P. aeruginosa within lungs of chronically infected CF patients. A novel, quantitative peptide nucleic acid (PNA) fluorescence in situ hybridization (PNA-FISH)-based method was used to estimate the in vivo growth rates of P. aeruginosa directly in lung tissue samples from CF patients and the growth rates of P. aeruginosa in infected lungs in a mouse model. The growth rate of P. aeruginosa within CF lungs did not correlate with the dimensions of bacterial aggregates but showed an inverse correlation to the concentration of polymorphonuclear leukocytes (PMNs) surrounding the bacteria. A growth-limiting effect on P. aeruginosa by PMNs was also observed in vitro, where this limitation was alleviated in the presence of the alternative electron acceptor nitrate. The finding that P. aeruginosa growth patterns correlate with the number of surrounding PMNs points to a bacteriostatic effect by PMNs via their strong O2 consumption, which slows the growth of P. aeruginosa in infected CF lungs. In support of this, the growth of P. aeruginosa was significantly higher in the respiratory airways than in the conducting airways of mice. These results indicate a complex host-pathogen interaction in chronic P. aeruginosa infection of the CF lung whereby PMNs slow the growth of the bacteria and render them less susceptible to antibiotic treatment while enabling them to persist by anaerobic respiration.

A framework for testing the ability of models to project climate change and its impacts

Abstract: Models used for climate change impact projections are typically not tested for simulation beyond current climate conditions. Since we have no data truly reflecting future conditions, a key challenge in this respect is to rigorously test models using proxies of future conditions. This paper presents a validation framework and guiding principles applicable across earth science disciplines for testing the capability of models to project future climate change and its impacts. Model test schemes comprising split-sample tests, differential split-sample tests and proxy site tests are discussed in relation to their application for projections by use of single models, ensemble modelling and space-time-substitution and in relation to use of different data from historical time series, paleo data and controlled experiments. We recommend that differential-split sample tests should be performed with best available proxy data in order to build further confidence in model projections.

Greenland Ice Sheet exports labile organic carbon to the Arctic oceans

Abstract: Runoff from small glacier systems contains dissolved organic carbon (DOC) rich in protein-like, low molecular weight (LMW) compounds, designating glaciers as an important source of bioavailable carbon for downstream heterotrophic activity. Fluxes of DOC and particulate organic carbon (POC) exported from large Greenland catchments, however, remain unquantified, despite the Greenland Ice Sheet (GrIS) being the largest source of global glacial runoff (ca. 400 km3 yr−1). We report high and episodic fluxes of POC and DOC from a large (> 600 km2) GrIS catchment during contrasting melt seasons. POC dominates organic carbon (OC) export (70–89% on average), is sourced from the ice sheet bed, and contains a significant bioreactive component (9% carbohydrates). A major source of the “bioavailable” (free carbohydrate) LMW–DOC fraction is microbial activity on the ice sheet surface, with some further addition of LMW–DOC to meltwaters by biogeochemical processes at the ice sheet bed. The bioavailability of the exported DOC (26–53%) to downstream marine microorganisms is similar to that reported from other glacial watersheds. Annual fluxes of DOC and free carbohydrates during two melt seasons were similar, despite the approximately two-fold difference in runoff fluxes, suggesting production-limited DOC sources. POC fluxes were also insensitive to an increase in seasonal runoff volumes, indicating a supply limitation in suspended sediment in runoff. Scaled to the GrIS, the combined DOC (0.13–0.17 TgC yr−1 (±13 %)) and POC fluxes (mean = 0.36–1.52 TgC yr−1 (±14 %)) are of a similar order of magnitude to a large Arctic river system, and hence may represent an important OC source to the near-coastal North Atlantic, Greenland and Labrador seas.

Ice–ocean interaction and calving front morphology at two west Greenland tidewater outlet glaciers

Abstract: Warm, subtropical-originating Atlantic water (AW) has been identified as a primary driver of mass loss across the marine sectors of the Greenland Ice Sheet (GrIS), yet the specific processes by which this water mass interacts with and erodes the calving front of tidewater glaciers is frequently modelled and much speculated upon but remains largely unobserved. We present a suite of fjord salinity, temperature, turbidity versus depth casts along with glacial runoff estimation from Rink and Store glaciers, two major marine outlets draining the western sector of the GrIS during 2009 and 2010. We characterise the main water bodies present and interpret their interaction with their respective calving fronts. We identify two distinct processes of ice– ocean interaction which have distinct spatial and temporal footprints: (1) homogenous free convective melting which occurs across the calving front where AW is in direct contact with the ice mass, and (2) localised upwelling-driven melt by turbulent subglacial runoff mixing with fjord water which occurs at distinct injection points across the calving front. Throughout the study, AW at 2.8 ± 0.2°C was consistently observed in contact with both glaciers below 450 m depth, yielding homogenous, free convective submarine melting up to ∼ 200 m depth. Above this bottom layer, multiple interactions are identified, primarily controlled by the rate of subglacial freshwater discharge which results in localised and discrete upwelling plumes. In the record melt year of 2010, the Store Glacier calving face was dominated by these runoff-driven plumes which led to a highly crenulated frontal geometry characterised by large embayments at the subglacial portals separated by headlands which are dominated by calving. Rink Glacier, which is significantly deeper than Store has a larger proportion of its submerged calving face exposed to AW, which results in a uniform, relatively flat overall frontal geometry.

A decade (2002–2012) of supraglacial lake volume estimates across Russell Glacier, West Greenland

Abstract: Supraglacial lakes represent an ephemeral storage buffer for meltwater runoff and lead to significant, yet short-lived, episodes of ice-flow acceleration by decanting large meltwater and energy fluxes into the ice sheet's hydrological system. Here, a methodology for calculating lake volume is used to quantify storage and drainage across Russell Glacier, West Greenland, between 2002 and 2012. Using 502 MODIS scenes, water volume at ~200 seasonally occurring lakes was derived using a depth–reflectance relationship, which was independently calibrated and field validated against lake bathymetry. The inland expansion of lakes is strongly correlated with air temperature: during the record melt years of 2010 and 2012, lakes formed and drained earlier, attaining their maximum volume 38 and 20 days earlier than the 11 yr mean, as well as occupying a greater area and forming at higher elevations (> 1800 m) than previously. Despite occupying under 2% of the study area, lakes delay the transmission of up to 7–13% of the bulk meltwater discharged. Although the results are subject to an observational bias caused by periods of cloud cover, we estimate that across Russell Glacier, 28% of supraglacial lakes drain rapidly (< 4 days). Clustering of such events in space and time suggests a synoptic trigger mechanism. Further, we find no evidence to support a unifying critical size or depth-dependent drainage threshold.

Persistent flow acceleration within the interior of the Greenland ice sheet

Abstract: We present surface velocity measurements from a high-elevation site located 140 km from the western margin of the Greenland ice sheet, and ~ 50 km into its accumulation area. Annual velocity increased each year from 51.78 ± 0.01 m yr−1 in 2009 to 52.92 ± 0.01 m yr−1 in 2012—a net increase of 2.2%. These data also reveal a strong seasonal velocity cycle of up to 8.1% above the winter mean, driven by seasonal melt and supraglacial lake drainage. Sole et al. (2013) recently argued that ice motion in the ablation area is mediated by reduced winter flow following the development of efficient subglacial drainage during warmer, faster, summers. Our data extend this analysis and reveal a year-on-year increase in annual velocity above the equilibrium line altitude, where despite surface melt increasing, it is still sufficiently low to hinder the development of efficient drainage under thick ice.

Assimilation of SMOS-derived soil moisture in a fully integrated hydrological and soil-vegetation-atmosphere transfer model in Western Denmark

Abstract: Real surface soil moisture retrieved from the Soil Moisture and Ocean Salinity (SMOS) satellite is downscaled and assimilated in a fully integrated hydrological and soil-vegetation-atmosphere transfer (MIKE SHE SW-ET) model using a bias aware ensemble transform Kalman filter (Bias-ETKF). Satellite-derived soil moisture assimilation in a catchment scale model is typically restricted by two challenges: (1) passive microwave is too coarse for direct assimilation and (2) the data tend to be biased. The solution proposed in this study is to disaggregate the SMOS bias using a higher resolution land cover classification map that was derived from Landsat thermal images. Using known correlations between SMOS bias and vegetation type, the assimilation filter is adapted to calculate biases online, using an initial bias estimate. Real SMOS-derived soil moisture is assimilated in a precalibrated catchment model in Denmark. The objective is to determine if any additional gains can be achieved by SMOS surface soil moisture assimilation beyond the optimized model. A series of assimilation experiments were designed to (1) determine how effectively soil moisture corrections propagate downward in the soil column, (2) compare the efficacy of in situ versus SMOS assimilation, and (3) determine how soil moisture assimilation affects fluxes and discharge in the catchment. We find that assimilation of SMOS improved R2 soil moisture correlations in the upper 5 cm compared to a network of 30 in situ sensors for most land cover classes. Assimilation also brought modest gains in R2 at 25 cm depth but slightly degraded the correlation at 50 cm depth. Assimilation overcorrected discharge peaks.

The Baltoscandian margin detrital zircon signatures of the central Scandes

Abstract: In central parts of the Scandinavian Caledonides, detrital zircon signatures provide evidence of the change in character of the Baltoscandian crystalline basement, from the characteristic Late Palaeoproterozoic granites of the Transscandinavian Igneous Belt (TIB, c. 1650–1850 Ma) in the foreland Autochthon to the typical, mainly Mesoproterozoic-age profile (c. 950–1700 Ma) of the Sveconorwegian Orogen of southwestern Scandinavia in the hinterland. Late Ediacaran to Early Cambrian shallow-marine Vemdal quartzites of the Jamtlandian Nappes (Lower Allochthon) provide strong bimodal signatures with TIB (1700–1800 Ma) and Sveconorwegian, sensu stricto (900–1150 Ma) ages dominant. Mid-Ordovician turbidites (Norraker Formation) of the Lower Allochthon in Sweden, sourced from the west, have unimodal signatures dominated by Sveconorwegian ages with peaks at 1000–1100 Ma, but with subordinate components of older Mesoproterozoic zircons (1200–1650 Ma). Latest Ordovician shallow-marine quartzites also yield bimodal signatures, but are more dispersed than in the Vemdal quartzites. In the greenschist facies lower parts of the Middle Allochthon, the Fuda (Offerdal Nappe) and Sarv Nappe signatures are either unimodal or bimodal (950–1100 and/or 1700–1850 Ma), with variable dominance of the younger or older group, and subordinate other Mesoproterozoic components. In the overlying, amphibolite to eclogite facies lower part of the Seve Nappe Complex, where the metasediments are dominated by feldspathic quartzites, calcsilicate-rich psammites and marbles, most units have bimodal signatures similar to the Sarv Nappes, but more dispersed; one has a unimodal signature very similar to the Ordovician turbidites of the Jamtlandian Nappes. In the overlying Upper Allochthon, Lower Koli (Baltica-proximal, Virisen Terrane), Late Ordovician quartzites provide unimodal signatures dominated by Sveconorwegian ages (sensu stricto). Further north in the Scandes, previously published zircon signatures in quartzites of the Lower Allochthon are similar to the Vemdal quartzites in Jamtland. Data from the Kalak Nappes at 70°N are in no way exotic to the Sveconorwegian Baltoscandian margin. They do show a Timanian influence (ages of c. 560–610 Ma), as would be expected from the palinspastic reconstructions of the nappes. Thus the detrital zircon signatures reported here and published elsewhere provide supporting evidence for a continuation northwards of the Sveconorwegian Orogen in the Neoproterozoic, from type areas in the south, along the Baltoscandian margin of Baltica into the high Arctic.

Increasing meltwater discharge from the Nuuk region of the Greenland ice sheet and implications for mass balance (1960-2012)

Abstract: We assess the runoff and surface mass balance (SMB) of the Greenland ice sheet in the Nuuk region (southwest) using output of two regional climate models (RCMs) evaluated by observations. The region encompasses six glaciers that drain into Godthabsfjord. RCM data (1960-2012) are resampled to a high spatial resolution to include the narrow (relative to the native grid spacing) glacier trunks in the ice mask. Comparing RCM gridded results with automatic weather station (AWS) point measurements reveals that locally models can underestimate ablation and overestimate accumulation by up to tens of per cent. However, comparison with lake discharge indicates that modelled regional runoff totals are more accurate. Model results show that melt and runoff in the Nuuk region have doubled over the past two decades. Regional SMB attained negative values in recent high-melt years. Taking into account frontal ablation of the marine-terminating glaciers, the region lost 10-20 km 3 w.e. a -1 in 2010-12. If 2010 melting prevails during the remainder of this century, a low-end estimate of sea-level rise of 5 mm is expected by 2100 from this relatively small section (2.6%) of the ice sheet alone.

From volcanic plains to glaciated peaks : BURIAL, uplift and exhumation history of southern East Greenland after opening of the NE Atlantic

Abstract: In southern East Greenland (68-70°N), voluminous flood basalts erupted onto a largely horizontal lava plain near sea level at the Paleocene-Eocene transition when sea-floor spreading started in the ...

Eocene to Miocene igneous activity in NE Greenland: northward younging of magmatism along the East Greenland margin

Abstract: Radiometric dating by the 40 Ar– 39 Ar incremental heating method was carried out on lavas, sills, dykes, and a central intrusion from NE Greenland. Eighteen samples gave acceptable crystallization ages. Lavas of both Lower and Upper Plateau Lava Series gave ages in the range 55.5–53.5 Ma and cannot be constrained to better than 56–53 Ma. Sills and dykes from Traill O to Shannon, with compositions fairly similar to those of the lavas, gave ages of 55.1–51.3 Ma, contemporaneous with and slightly younger than the lavas. Alkaline lavas on inland nunataks have ages of 53–50 Ma, and the Kap Broer Ruys intrusive centre has an age of 48.7 ± 0.5 Ma. An alkaline sill on Hvalroso is much younger at 20.3 ± 0.1 Ma. There are no pre-breakup lavas onshore NE Greenland. We surmise that the hot mantle of the Iceland plume arrived and melted extensively beneath the northern basins only at the time of breakup around 55 Ma. Post-breakup intrusive events in NE Greenland coincided with plate-tectonic events such as reorganization, uplift and opening in the north. The Hvalroso sill represents a local small melting event that may be related to coeval opening of the Lena Trough. Supplementary materials: Details of the source data, results, and the compositions and locations of the dated samples, are available at www.geolsoc.org.uk/SUP18738.

Theory of sampling (TOS) versus measurement uncertainty (MU): A call for integration

Abstract: We assess current approaches to measurement uncertainty (MU) with respect to the complete ensemble of sources affecting the measurement process, in particular the extent to which sampling errors as set out in the Theory of Sampling (TOS) are appropriately considered in the GUM and EURACHEM/CITAC guides. All pre-analysis sampling steps play an important, often dominant role in the total uncertainty budget, thereby critically affecting the validity of MU estimates, but most of these contributions are not included in the current MU framework. The TOS constitutes the only complete theoretical platform for dealing appropriately with the entire pathway from field sample to test portion. We here propose a way to reconcile the often strongly felt differences between MU and TOS. There is no need to debate terminology, as both TOS and MU can be left with their current usages.

Large-scale 3-D modeling by integration of resistivity models and borehole data through inversion

Abstract: We present an automatic method for parameterization of a 3-D model of the subsurface, integrating lithological information from boreholes with resistivity models through an inverse optimization, with the objective of further detailing of geological models, or as direct input into groundwater models. The parameter of interest is the clay fraction, expressed as the relative length of clay units in a depth interval. The clay fraction is obtained from lithological logs and the clay fraction from the resistivity is obtained by establishing a simple petrophysical relationship, a translator function, between resistivity and the clay fraction. Through inversion we use the lithological data and the resistivity data to determine the optimum spatially distributed translator function. Applying the translator function we get a 3-D clay fraction model, which holds information from the resistivity data set and the borehole data set in one variable. Finally, we use k -means clustering to generate a 3-D model of the subsurface structures. We apply the procedure to the Norsminde survey in Denmark, integrating approximately 700 boreholes and more than 100 000 resistivity models from an airborne survey in the parameterization of the 3-D model covering 156 km 2 . The final five-cluster 3-D model differentiates between clay materials and different high-resistivity materials from information held in the resistivity model and borehole observations, respectively.

Glacier dynamics at Helheim and Kangerdlugssuaq glaciers, southeast Greenland, since the Little Ice Age

Abstract: . Observations over the past decade show significant ice loss associated with the speed-up of glaciers in southeast Greenland from 2003, followed by a deceleration from 2006. These short-term, episodic, dynamic perturbations have a major impact on the mass balance on the decadal scale. To improve the projection of future sea level rise, a long-term data record that reveals the mass balance beyond such episodic events is required. Here, we extend the observational record of marginal thinning of Helheim and Kangerdlugssuaq glaciers from 10 to more than 80 years. We show that, although the frontal portion of Helheim Glacier thinned by more than 100 m between 2003 and 2006, it thickened by more than 50 m during the previous two decades. In contrast, Kangerdlugssuaq Glacier underwent minor thinning of 40–50 m from 1981 to 1998 and major thinning of more than 100 m after 2003. Extending the record back to the end of the Little Ice Age (prior to 1930) shows no thinning of Helheim Glacier from its maximum extent during the Little Ice Age to 1981, while Kangerdlugssuaq Glacier underwent substantial thinning of 230 to 265 m. Comparison of sub-surface water temperature anomalies and variations in air temperature to records of thickness and velocity change suggest that both glaciers are highly sensitive to short-term atmospheric and ocean forcing, and respond very quickly to small fluctuations. On century timescales, however, multiple external parameters (e.g. outlet glacier shape) may dominate the mass change. These findings suggest that special care must be taken in the projection of future dynamic ice loss.

Challenges in conditioning a stochastic geological model of a heterogeneous glacial aquifer to a comprehensive soft data set

Abstract: In traditional hydrogeological investigations, one geological model is often used based on subjective inter- pretations and sparse data availability. This deterministic approach usually does not account for any uncertainties. Stochastic simulation methods address this problem and can capture the geological structure uncertainty. In this study the geostatistical software TProGS is utilized to simulate an en- semble of realizations for a binary (sand/clay) hydrofacies model in the Norsminde catchment, Denmark. TProGS can incorporate soft data, which represent the associated level of uncertainty. High-density (20 m◊ 20 m◊ 2 m) airborne geo- physical data (SkyTEM) and categorized borehole data are utilized to define the model of spatial variability in horizontal and vertical direction, respectively, and both are used for soft conditioning of the TProGS simulations. The category proba- bilities for the SkyTEM data set are derived from a histogram probability matching method, where resistivity is paired with the corresponding lithology from the categorized borehole data. This study integrates two distinct data sources into the stochastic modeling process that represent two extremes of the conditioning density spectrum: sparse borehole data and abundant SkyTEM data. In the latter the data have a strong spatial correlation caused by its high data density, which trig- gers the problem of overconditioning. This problem is ad- dressed by a work-around utilizing a sampling/decimation of the data set, with the aim to reduce the spatial correlation of the conditioning data set. In the case of abundant condi- tioning data, it is shown that TProGS is capable of repro- ducing non-stationary trends. The stochastic realizations are validated by five performance criteria: (1) sand proportion, (2) mean length, (3) geobody connectivity, (4) facies proba- bility distribution and (5) facies probability-resistivity bias. In conclusion, a stochastically generated set of realizations soft-conditioned to 200 m moving sampling of geophysical data performs most satisfactorily when balancing the five performance criteria. The ensemble can be used in subse- quent hydrogeological flow modeling to address the pre- dictive uncertainty originating from the geological structure uncertainty.