Search Results

Now showing 1 - 10 of 111
  • Item
    Cryostratigraphy, sedimentology, and the late Quaternary evolution of the Zackenberg River delta, northeast Greenland
    (Katlenburg-Lindau : Copernicus, 2017-5-30) Gilbert, Graham L.; Cable, Stefanie; Thiel, Christine; Christiansen, Hanne H.; Elberling, Bo
    The Zackenberg River delta is located in northeast Greenland (74°30′ N, 20°30′ E) at the outlet of the Zackenberg fjord valley. The fjord-valley fill consists of a series of terraced deltaic deposits (ca. 2 km2) formed during relative sea-level (RSL) fall. We investigated the deposits using sedimentological and cryostratigraphic techniques together with optically stimulated luminescence (OSL) dating. We identify four facies associations in sections (4 to 22 m in height) exposed along the modern Zackenberg River and coast. Facies associations relate to (I) overriding glaciers, (II) retreating glaciers and quiescent glaciomarine conditions, (III) delta progradation in a fjord valley, and (IV) fluvial activity and niveo-aeolian processes. Pore, layered, and suspended cryofacies are identified in two 20 m deep ice-bonded sediment cores. The cryofacies distribution, together with low overall ground-ice content, indicates that permafrost is predominately epigenetic in these deposits. Fourteen OSL ages constrain the deposition of the cored deposits to between approximately 13 and 11 ka, immediately following deglaciation. The timing of permafrost aggradation was closely related to delta progradation and began following the subaerial exposure of the delta plain (ca. 11 ka). Our results reveal information concerning the interplay between deglaciation, RSL change, sedimentation, permafrost aggradation, and the timing of these events. These findings have implications for the timing and mode of permafrost aggradation in other fjord valleys in northeast Greenland.
  • Item
    Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) – Part 2: Parameter ensemble analysis
    (Katlenburg-Lindau : Copernicus, 2020) Albrecht, Torsten; Winkelmann, Ricarda; Levermann, Anders
    The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation–age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model–data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016). This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4±4.1 m (or 6.5±2.0×106km3 ), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
  • Item
    Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams
    (Katlenburg-Lindau : Copernicus, 2018) Reese, Ronja; Winkelmann, Ricarda; Gudmundsson, G. Hilmar
    Currently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.
  • Item
    Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations
    (Katlenburg-Lindau : Copernicus, 2022) Feldmann, Johannes; Reese, Ronja; Winkelmann, Ricarda; Levermann, Anders
    Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.
  • Item
    Modeling the response of Greenland outlet glaciers to global warming using a coupled flow line-plume model
    (Katlenburg-Lindau : Copernicus, 2019) Beckmann, Johanna; Perrette, Mahé; Beyer, Sebastian; Calov, Reinhard; Willeit, Matteo; Ganopolski, Andrey
    In recent decades, the Greenland Ice Sheet has experienced an accelerated mass loss, contributing to approximately 25 % of contemporary sea level rise (SLR). This mass loss is caused by increased surface melt over a large area of the ice sheet and by the thinning, retreat and acceleration of numerous Greenland outlet glaciers. The latter is likely connected to enhanced submarine melting that, in turn, can be explained by ocean warming and enhanced subglacial discharge. The mechanisms involved in submarine melting are not yet fully understood and are only simplistically incorporated in some models of the Greenland Ice Sheet. Here, we investigate the response of 12 representative Greenland outlet glaciers to atmospheric and oceanic warming using a coupled line–plume glacier–flow line model resolving one horizontal dimension. The model parameters have been tuned for individual outlet glaciers using present-day observational constraints. We then run the model from present to the year 2100, forcing the model with changes in surface mass balance and surface runoff from simulations with a regional climate model for the RCP8.5 scenario, and applying a linear ocean temperature warming with different rates of changes representing uncertainties in the CMIP5 model experiments for the same climate change scenario. We also use different initial temperature–salinity profiles obtained from direct measurements and from ocean reanalysis data. Using different combinations of submarine melting and calving parameters that reproduce the present-day state of the glaciers, we estimate uncertainties in the contribution to global SLR for individual glaciers. We also perform a sensitivity analysis of the three forcing factors (changes in surface mass balance, ocean temperature and subglacial discharge), which shows that the roles of the different forcing factors are diverse for individual glaciers. We find that changes in ocean temperature and subglacial discharge are of comparable importance for the cumulative contribution of all 12 glaciers to global SLR in the 21st century. The median range of the cumulative contribution to the global SLR for all 12 glaciers is about 18 mm (the glaciers' dynamic response to changes of all three forcing factors). Neglecting changes in ocean temperature and subglacial discharge (which control submarine melt) and investigating the response to changes in surface mass balance only leads to a cumulative contribution of 5 mm SLR. Thus, from the 18 mm we associate roughly 70 % with the glaciers' dynamic response to increased subglacial discharge and ocean temperature and the remaining 30 % (5 mm) to the response to increased surface mass loss. We also find a strong correlation (correlation coefficient 0.74) between present-day grounding line discharge and their future contribution to SLR in 2100. If the contribution of the 12 glaciers is scaled up to the total present-day discharge of Greenland, we estimate the midrange contribution of all Greenland glaciers to 21st-century SLR to be approximately 50 mm. This number adds to SLR derived from a stand-alone ice sheet model (880 mm) that does not resolve outlet glaciers and thus increases SLR by over 50 %. This result confirms earlier studies showing that the response of the outlet glaciers to global warming has to be taken into account to correctly assess the total contribution of Greenland to sea level change.
  • Item
    Atomic oxygen number densities in the mesosphere–lower thermosphere region measured by solid electrolyte sensors onWADIS-2
    (Katlenburg-Lindau : Copernicus, 2019) Eberhart, Martin; Löhle, Stefan; Strelnikov, Boris; Hedin, Jonas; Khaplanov, Mikhail; Fasoulas, Stefanos; Gumbel, Jörg; Lübken, Franz-Josef; Rapp, Markus
    Absolute profiles of atomic oxygen number densities with high vertical resolution have been determined in the mesosphere-lower thermosphere (MLT) region from in situ measurements by several rocket-borne solid electrolyte sensors. The amperometric sensors were operated in both controlled and uncontrolled modes and with various orientations on the foredeck and aft deck of the payload. Calibration was based on mass spectrometry in a molecular beam containing atomic oxygen produced in a microwave discharge. The sensor signal is proportional to the number flux onto the electrodes, and the mass flow rate in the molecular beam was additionally measured to derive this quantity from the spectrometer reading. Numerical simulations provided aerodynamic correction factors to derive the atmospheric number density of atomic oxygen from the sensor data. The flight results indicate a preferable orientation of the electrode surface perpendicular to the rocket axis. While unstable during the upleg, the density profiles measured by these sensors show an excellent agreement with the atmospheric models and photometer results during the downleg of the trajectory. The high spatial resolution of the measurements allows for the identification of small-scale variations in the atomic oxygen concentration. © Author(s) 2019. This work is distributed under the Creative Commons Attribution 4.0 License.
  • Item
    A new method of inferring the size, number density, and charge of mesospheric dust from its in situ collection by the DUSTY probe
    (Katlenburg-Lindau : Copernicus, 2019) Havnes, Ove; Antonsen, Tarjei; Baumgarten, Gerd; Hartquist, Thomas W.; Biebricher, Alexander; Fredriksen, Ã…shild; Friedrich, Martin; Hedin, Jonas
    We present a new method of analyzing measurements of mesospheric dust made with DUSTY rocket-borne Faraday cup probes. It can yield the variation in fundamental dust parameters through a mesospheric cloud with an altitude resolution down to 10 cm or less if plasma probes give the plasma density variations with similar height resolution. A DUSTY probe was the first probe that unambiguously detected charged dust and aerosol particles in the Earth's mesosphere. DUSTY excluded the ambient plasma by various biased grids, which however allowed dust particles with radii above a few nanometers to enter, and it measured the flux of charged dust particles. The flux measurements directly yielded the total ambient dust charge density. We extend the analysis of DUSTY data by using the impact currents on its main grid and the bottom plate as before, together with a dust charging model and a secondary charge production model, to allow the determination of fundamental parameters, such as dust radius, charge number, and total dust density. We demonstrate the utility of the new analysis technique by considering observations made with the DUSTY probes during the MAXIDUSTY rocket campaign in June-July 2016 and comparing the results with those of other instruments (lidar and photometer) also used in the campaign. In the present version we have used monodisperse dust size distributions.© Author(s) 2019.
  • Item
    Increasing the spatial resolution of cloud property retrievals from Meteosat SEVIRI by use of its high-resolution visible channel: Evaluation of candidate approaches with MODIS observations
    (Katlenburg-Lindau : Copernicus, 2020) Werner, Frank; Deneke, Hartwig
    This study presents and evaluates several candidate approaches for downscaling observations from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) in order to increase the horizontal resolution of subsequent cloud optical thickness (τ) and effective droplet radius (reff) retrievals from the native ≈ 3km×3km spatial resolution of the narrowband channels to ≈ 1km×1km. These methods make use of SEVIRI's coincident broadband high-resolution visible (HRV) channel. For four example cloud fields, the reliability of each downscaling algorithm is evaluated by means of collocated 1km×1km MODIS radiances, which are reprojected to the horizontal grid of the HRV channel and serve as reference for the evaluation. By using these radiances, smoothed with the modulation transfer function of the native SEVIRI channels, as retrieval input, the accuracy at the SEVIRI standard resolution can be evaluated and an objective comparison of the accuracy of the different downscaling algorithms can be made. For the example scenes considered in this study, it is shown that neglecting high-frequency variations below the SEVIRI standard resolution results in significant random absolute deviations of the retrieved τ and reff of up to ≈ 14 and ≈ 6μm, respectively, as well as biases. By error propagation, this also negatively impacts the reliability of the subsequent calculation of liquid water path (WL) and cloud droplet number concentration (ND), which exhibit deviations of up to ≈ 89gm-2 and ≈ 177cm-3, respectively. For τ , these deviations can be almost completely mitigated by the use of the HRV channel as a physical constraint and by applying most of the presented downscaling schemes. Uncertainties in retrieved reff at the native SEVIRI resolution are smaller, and the improvements from downscaling the observations are less obvious than for τ. Nonetheless, the right choice of downscaling scheme yields noticeable improvements in the retrieved reff. Furthermore, the improved reliability in retrieved cloud products results in significantly reduced uncertainties in derived WL and ND. In particular, one downscaling approach provides clear improvements for all cloud products compared to those obtained from SEVIRI's standard resolution and is recommended for future downscaling endeavors. This work advances efforts to mitigate impacts of scale mismatches among channels of multiresolution instruments on cloud retrievals. © Author(s) 2020.
  • Item
    The realization of autonomous, aircraft-based, real-time aerosol mass spectrometry in the upper troposphere and lower stratosphere
    (Katlenburg-Lindau : Copernicus, 2022) Dragoneas, Antonis; Molleker, Sergej; Appel, Oliver; Hünig, Andreas; Böttger, Thomas; Hermann, Markus; Drewnick, Frank; Schneider, Johannes; Weigel, Ralf; Borrmann, Stephan
    We report on the developments that enabled the field deployment of a fully automated aerosol mass spectrometer, especially designed for high-altitude measurements on unpressurized aircraft. The merits of the two main categories of real-time aerosol mass spectrometry, i.e. (a) single-particle laser desorption and ionization and (b) continuous thermal desorption and electron impact ionization of aerosols, have been integrated into one compact apparatus with the aim to perform in situ real-time analysis of aerosol chemical composition. The demonstrated instrument, named the ERICA (European Research Council Instrument for Chemical composition of Aerosols), operated successfully aboard the high-altitude research aircraft M-55 Geophysica at altitudes up to 20 km while being exposed to ambient conditions of very low atmospheric pressure and temperature. A primary goal of those field deployments was the in situ study of the Asian tropopause aerosol layer (ATAL). During 11 research flights, the instrument operated for more than 49 h and collected chemical composition information of more than 150 000 single particles combined with quantitative chemical composition analysis of aerosol particle ensembles. This paper presents in detail the technical characteristics of the main constituent parts of the instrument, as well as the design considerations for its integration into the aircraft and its autonomous operation in the upper troposphere and lower stratosphere (UTLS). Additionally, system performance data from the first field deployments of the instrument are presented and discussed, together with exemplary mass spectrometry data collected during those flights.
  • Item
    Photocurrent modelling and experimental confirmation for meteoric smoke particle detectors on board atmospheric sounding rockets
    (Katlenburg-Lindau : Copernicus, 2018-9-20) Giono, Gabriel; Strelnikov, Boris; Asmus, Heiner; Staszak, Tristan; Ivchenko, Nickolay; Lübken, Franz-Josef
    Characterising the photoelectron current induced by the Sun's UV radiation is crucial to ensure accurate daylight measurements from particle detectors. This article lays out the methodology used to address this problem in the case of the meteoric smoke particle detectors (MSPDs), developed by the Leibniz Institute of Atmospheric Physics in Kühlungsborn (IAP) and flown on board the PMWEs (Polar Mesosphere Winter Echoes) sounding rockets in April 2018. The methodology focuses on two complementary aspects: modelling and experimental measurements. A detailed model of the MSPD photocurrent was created based on the expected solar UV flux, the atmospheric UV absorption as a function of height by molecular oxygen and ozone, the photoelectric yield of the material coating the MSPD as a function of wavelength, the index of refraction of these materials as a function of wavelength and the angle of incidence of the illumination onto the MSPD. Due to its complex structure, composed of a central electrode shielded by two concentric grids, extensive ray-tracing calculations were conducted to obtain the incidence angles of the illumination on the central electrode, and this was done for various orientations of the MSPD in respect to the Sun. Results of the modelled photocurrent at different heights and for different materials, as well as for different orientations of the detector, are presented. As a pre-flight confirmation, the model was used to reproduce the experimental measurements conducted by Robertson et al. (2014) and agrees within an order of magnitude. An experimental setup for the calibration of the MSPD photocurrent is also presented. The photocurrent induced by the Lyman-alpha line from a deuterium lamp was recorded inside a vacuum chamber using a narrowband filter, while a UV-sensitive photodiode was used to monitor the UV flux. These measurements were compared with the model prediction, and also matched within an order of magnitude. Although precisely modelling the photocurrent is a challenging task, this article quantitatively improved the understanding of the photocurrent on the MSPD and discusses possible strategies to untangle the meteoric smoke particles (MSPs) current from the photocurrent recorded in-flight.