Search

Over the past 3 decades, inversions for ice sheet basal drag have become commonplace in glaciological modeling. Such inversions require regularization to prevent over-fitting and ensure that the structure they recover is a robust inference from the observations, confidence which is required if they are to be used to draw conclusions about processes and properties of the ice base. While L-curve analysis can be used to select the optimal regularization level, the treatment of L-curve analysis in glaciological inverse modeling has been highly variable. Building on the history of glaciological inverse modeling, we demonstrate general best practices for regularizing glaciological inverse problems, using a domain in the Filchner–Ronne catchment of Antarctica as our test bed. We show a step-by-step approach to cost function normalization and L-curve analysis. We explore the spatial and spectral characteristics of the solution as a function of regularization, and we test the sensitivity of L-curve analysis and regularization to model resolution, effective pressure, sliding nonlinearity, and the flow equation. We find that the optimal regularization level converges towards a finite non-zero limit in the continuous problem, associated with a best knowable basal drag field. Nonlinear sliding laws outperform linear sliding in our analysis, with both a lower total variance and a more sharply cornered L-curve. By contrast, geometry-based approximations for effective pressure degrade inversion performance when added to a sliding law, but an actual hydrology model may marginally improve performance in some cases. Our results with 3D inversions suggest that the additional model complexity may not be justified by the 2D nature of the surface velocity data. We conclude with recommendations for best practices in future glaciological inversions. [ABSTRACT FROM AUTHOR]

The subglacial hydrological system affects (i) the motion of ice sheets through sliding, (ii) the location of lakes at the ice margin, and (iii) the ocean circulation by freshwater discharge directly at the grounding line or (iv) via rivers flowing over land. For modeling this hydrology system, a previously developed porous-media concept called the confined–unconfined aquifer system (CUAS) is used. To allow for realistic simulations at the ice sheet scale, we developed CUAS-MPI, an MPI-parallel C/C++ implementation of CUAS (MPI: Message Passing Interface), which employs the Portable, Extensible Toolkit for Scientific Computation (PETSc) infrastructure for handling grids and equation systems. We validate the accuracy of the numerical results by comparing them with a set of analytical solutions to the model equations, which involve two types of boundary conditions. We then investigate the scaling behavior of CUAS-MPI and show that CUAS-MPI scales up to 3840 MPI processes running a realistic Greenland setup on the Lichtenberg HPC system. Our measurements also show that CUAS-MPI reaches a throughput comparable to that of ice sheet simulations, e.g., the Ice-sheet and Sea-level System Model (ISSM). Lastly, we discuss opportunities for ice sheet modeling, explore future coupling possibilities of CUAS-MPI with other simulations, and consider throughput bottlenecks and limits of further scaling. [ABSTRACT FROM AUTHOR]

Full-Stokes (FS) ice sheet models provide the most sophisticated formulation of ice sheet flow. However, their applicability is often limited due to the high computational demand and numerical challenges. To balance computational demand and accuracy, the so-called Blatter–Pattyn (BP) stress regime is frequently used. Here, we explore the dynamic consequences of using simplified approaches by solving FS and the BP stress regime applied to the Northeast Greenland Ice Stream. To ensure a consistent comparison, we use one single ice sheet model to run the simulations under identical numerical conditions. A sensitivity study to the horizontal grid resolution (from 12.8 to a resolution of 0.1 km) reveals that velocity differences between the FS and BP solution emerge below ∼ 1 km horizontal resolution and continuously increase with resolution. Over the majority of the modelling domain both models reveal similar surface velocity patterns. At the grounding line of the 79 ∘ North Glacier the simulations show considerable differences whereby the BP model overestimates ice discharge of up to 50 % compared to FS. A sensitivity study to the friction type reveals that differences are stronger for a power-law friction than a linear friction law. Model differences are attributed to topographic variability and the basal drag, in which neglected stress terms in BP become important. [ABSTRACT FROM AUTHOR]

Two different approaches can be applied to calculate the mass transfer in a two‐phase system using computational fluid dynamics and the volume‐of‐fluid method, single‐field and two‐field approaches. The used method affects the stability of the calculation as well as the accuracy of the result. Two volume‐of‐fluid‐based approaches, namely, the generalized continuous species transfer model, a single‐field approach, and a two‐field approach are compared. The models are implemented in the OpenFOAM® framework and validated in a 1D test case. Although both approaches show great agreement with the analytical solution for large time steps, numerical instabilities of the two‐field model are revealed with a reduced time step size, which leads to false results. [ABSTRACT FROM AUTHOR]

Antarctica's ice shelves modulate the grounded ice flow, and weakening of ice shelves due to climate forcing will decrease their 'buttressing' effect, causing a response in the grounded ice. While the processes governing ice-shelf weakening are complex, uncertainties in the response of the grounded ice sheet are also difficult to assess. The Antarctic BUttressing Model Intercomparison Project (ABUMIP) compares ice-sheet model responses to decrease in buttressing by investigating the 'end-member' scenario of total and sustained loss of ice shelves. Although unrealistic, this scenario enables gauging the sensitivity of an ensemble of 15 ice-sheet models to a total loss of buttressing, hence exhibiting the full potential of marine ice-sheet instability. All models predict that this scenario leads to multi-metre (1–12 m) sea-level rise over 500 years from present day. West Antarctic ice sheet collapse alone leads to a 1.91–5.08 m sea-level rise due to the marine ice-sheet instability. Mass loss rates are a strong function of the sliding/friction law, with plastic laws cause a further destabilization of the Aurora and Wilkes Subglacial Basins, East Antarctica. Improvements to marine ice-sheet models have greatly reduced variability between modelled ice-sheet responses to extreme ice-shelf loss, e.g. compared to the SeaRISE assessments. [ABSTRACT FROM AUTHOR]

Various observational estimates indicate growing mass loss at Antarctica's margins as well as heavier precipitation across the continent. Simulated future projections reveal that heavier precipitation, falling on Antarctica, may counteract amplified iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we test how the ansatz (implementation in a mathematical framework) of the precipitation boundary condition shapes Antarctica's sea level contribution in an ensemble of ice sheet simulations. We test two precipitation conditions: we either apply the precipitation anomalies from CMIP5 models directly or scale the precipitation by the air temperature anomalies from the CMIP5 models. In the scaling approach, it is common to use a relative precipitation increment per degree warming as an invariant scaling constant. We use future climate projections from nine CMIP5 models, ranging from strong mitigation efforts to business-as-usual scenarios, to perform simulations from 1850 to 5000. We take advantage of individual climate projections by exploiting their full temporal and spatial structure. The CMIP5 projections beyond 2100 are prolonged with reiterated forcing that includes decadal variability; hence, our study may underestimate ice loss after 2100. In contrast to various former studies that apply an evolving temporal forcing that is spatially averaged across the entire Antarctic Ice Sheet, our simulations consider the spatial structure in the forcing stemming from various climate patterns. This fundamental difference reproduces regions of decreasing precipitation despite general warming. Regardless of the boundary and forcing conditions applied, our ensemble study suggests that some areas, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea, will lose ice in the future. In general, the simulated ice sheet thickness grows along the coast, where incoming storms deliver topographically controlled precipitation. In this region, the ice thickness differences are largest between the applied precipitation methods. On average, Antarctica shrinks for all future scenarios if the air temperature anomalies scale the precipitation. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies directly. The analysis reveals that the mean scaling inferred from climate models is larger than the commonly used values deduced from ice cores; moreover, it varies spatially: the highest scaling is across the East Antarctic Ice Sheet, and the lowest scaling is around the Siple Coast, east of the Ross Ice Shelf. The discrepancies in response to both precipitation ansatzes illustrate the principal uncertainty in projections of Antarctica's sea level contribution. [ABSTRACT FROM AUTHOR]

The thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheet may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs). Additionally, the highly anisotropic geometry of the 3D elements in ice sheet modelling poses a problem for stabilization approaches in advection-dominated problems. Here, we present extended enthalpy formulations within the finite-element Ice-Sheet and Sea-Level System model (ISSM) that show a better performance than earlier implementations. In a first polythermal-slab experiment, we found that the treatment of the discontinuous conductivities at the CTS with a geometric mean produces more accurate results compared to the arithmetic or harmonic mean. This improvement is particularly efficient when applied to coarse vertical resolutions. In a second ice dome experiment, we find that the numerical solution is sensitive to the choice of stabilization parameters in the well-established streamline upwind Petrov–Galerkin (SUPG) method. As standard literature values for the SUPG stabilization parameter do not account for the highly anisotropic geometry of the 3D elements in ice sheet modelling, we propose a novel anisotropic SUPG (ASUPG) formulation. This formulation circumvents the problem of high aspect ratio by treating the horizontal and vertical directions separately in the stabilization coefficients. The ASUPG method provides accurate results for the thermodynamic equation on geometries with very small aspect ratios like ice sheets. [ABSTRACT FROM AUTHOR]

Stirred tanks are widely used in chemical engineering, which is why there is a large variety of different designs. Numerical fluid dynamics is therefore often used in stirred tank studies that focus on flow behavior. In this work, a generic reaction in a stirred tank is simulated using computational fluid dynamics and the temperature control of the reactor is investigated. Design‐critical aspects are derived from the simulation. In addition, the simulation approach is very well suited for a reaction‐specific optimization of a stirred tank. [ABSTRACT FROM AUTHOR]

We present the result of the third Marine Ice Sheet Model Intercomparison Project, MISMIP+. MISMIP+ is intended to be a benchmark for ice-flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient to model buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments first tests that models are able to maintain a steady state with the grounding line located on a retrograde slope due to buttressing and then explore scenarios where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. The majority of participating models passed the first test and then produced similar responses to the loss of buttressing. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress but also the differences between numerical methods – taking a secondary role. [ABSTRACT FROM AUTHOR]

The thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheets may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs). Additionally, the highly anisotropic geometry of the 3D elements in ice sheet modelling poses a problem for stabilization approaches in advection dominated problems. Here, we present extended enthalpy formulations within the finite-element Ice Sheet System Model (ISSM) that show a better performance to earlier implementations. In a first polythermal-slab experiment, we found that the treatment of the discontinuous conductivities at the CTS with a geometric mean produce more accurate results compared to the arithmetic or harmonic mean. This improvement is particularly efficient when applied to coarse vertical resolutions. In a second ice dome experiment, we find that the numerical solution is sensitive to the choice of stabilization parameters in the well-established Streamline Upwind Petrov–Galerkin (SUPG) method. As standard literature values for the SUPG stabilization parameter are not accounting for the highly anisotropic geometry of the 3D elements in ice sheet modelling, we propose a novel Anisotropic SUPG (ASUPG) formulation. This formulation circumvents the problem of high aspect-ratio by treating the horizontal and vertical directions separately in the stabilization coefficients. The ASUPG method provides accurate results for the thermodynamic equation on geometries with very small aspect ratios like ice sheets. [ABSTRACT FROM AUTHOR]

Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity. [ABSTRACT FROM AUTHOR]

Today industrial production is performed across national boundaries in real time with great accuracy and precision using digital information communication technologies. Manufacturing in one location is possible as it receives digital signals and messages transmissions from remote activities across borders. But who are these magicians behind the curtain? They are the highly skilled and educated knowledge workers of the global, blue-collar, working class. A generation ago working in industrial settings we would have likely shared close physical proximity as well as a common language and culture. Today we share the interconnectedness of our global village across national, cultural, and linguistic boundaries. [ABSTRACT FROM AUTHOR]

Various observational estimates indicate growing mass loss at Antarctica's margins but also heavier precipitation across the continent. In the future, heavier precipitation fallen on Antarctica will counteract any stronger iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we use from nine CMIP5 models future projections, ranging from strong mitigation efforts to business-as-usual, to run an ensemble of ice-sheet simulations. We test, how the precipitation boundary condition determines Antarctica's sea-level contribution. The spatial and temporal varying climate forcings drive ice-sheet simulations. Hence, our ensemble inherits all spatial and temporal climate patterns, which is in contrast to a spatial mean forcing. Regardless of the applied boundary condition and forcing, some areas will lose ice in the future, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea. In general the simulated ice-sheet thickness grows in a broad marginal strip, where incoming storms deliver topographically controlled precipitation. This strip shows the largest ice thickness differences between the applied precipitation boundary conditions too. On average Antarctica's ice mass shrinks for all future scenarios if the precipitation is scaled by the spatial temperature anomalies coming from the CMIP5 models. In this approach, we use the relative precipitation increment per degree warming as invariant scaling constant. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies of the CMIP5 models directly. Here, the scaling factors show a distinct spatial pattern across Antarctica. Furthermore, the diagnosed mean scaling across all considered climate forcings is larger than the values deduced from ice cores. In general, the scaling is higher across the East Antarctic Ice Sheet, lower across the West Antarctic Ice Sheet, and lowest around the Siple Coast. The latter is located on the east side of the Ross Ice Shelf. [ABSTRACT FROM AUTHOR]

The international endeavour to retrieve a continuous ice core, which spans the middle Pleistocene climate transition ca. 1.2–0.9 Myr ago, encompasses a multitude of field and model-based pre-site surveys. We expand on the current efforts to locate a suitable drilling site for the oldest Antarctic ice core by means of 3-D continental ice-sheet modelling. To this end, we present an ensemble of ice-sheet simulations spanning the last 2 Myr, employing transient boundary conditions derived from climate modelling and climate proxy records. We discuss the imprint of changing climate conditions, sea level and geothermal heat flux on the ice thickness, and basal conditions around previously identified sites with continuous records of old ice. Our modelling results show a range of configurational ice-sheet changes across the middle Pleistocene transition, suggesting a potential shift of the West Antarctic Ice Sheet to a marine-based configuration. Despite the middle Pleistocene climate reorganisation and associated ice-dynamic changes, we identify several regions conducive to conditions maintaining 1.5 Myr (million years) old ice, particularly around Dome Fuji, Dome C and Ridge B, which is in agreement with previous studies. This finding strengthens the notion that continuous records with such old ice do exist in previously identified regions, while we are also providing a dynamic continental ice-sheet context. [ABSTRACT FROM AUTHOR]

The East Antarctic Ice Sheet contains a wealth of information that can be extracted from its internal architecture such as distribution of age, past flow features, and surface and basal properties. Airborne radar surveys can sample this stratigraphic archive across broad areas. Here, we identify and trace key horizons across several radar surveys to obtain the stratigraphic information. We transfer the age–depth scales from ice cores to intersecting radar data. We then propagate these age scales across the ice sheet using the high fidelity continuity of the radar horizons. In Dronning Maud Land, including Dome Fuji, we mapped isochrones with ages of 38 and 74 ka. In the central region of East Antarctica around Dome Concordia, Vostok and Dome Argus, we use isochrone ages of 38, 48, 90 and 161 ka. Taking together both regions, we provide isochrone depths traced along a combined profile length of more than 40 000 km and discuss uncertainties of the obtained stratigraphy, as well as factors important to consider for further expansion. This data set is the most extensive distribution of internal horizons in East Antarctica to date. The isochrone depths presented in this study are available on PANGAEA (10.1594/PANGAEA.895528;). [ABSTRACT FROM AUTHOR]

Subglacial hydrology plays an important role in ice sheet dynamics as it determines the sliding velocity. It also drives freshwater into the ocean, leading to undercutting of calving fronts by plumes. Modeling subglacial water has been a challenge for decades. Only recently have new approaches been developed such as representing subglacial channels and thin water sheets by separate layers of variable hydraulic conductivity. We extend this concept by modeling a confined–unconfined aquifer system (CUAS) in a single layer of an equivalent porous medium (EPM). The advantage of this formulation is that it prevents unphysical values of pressure at reasonable computational cost. We performed sensitivity tests to investigate the effect of different model parameters. The strongest influence of model parameters was detected in terms of governing the opening and closure of the system. Furthermore, we applied the model to the Northeast Greenland Ice Stream, where an efficient system independent of seasonal input was identified about 500 km downstream from the ice divide. Using the effective pressure from the hydrology model, the Ice Sheet System Model (ISSM) showed considerable improvements in modeled velocities in the coastal region. [ABSTRACT FROM AUTHOR]

Recovery Glacier reaches far into the East Antarctic Ice Sheet. Recent projections point out that its dynamic behavior has a considerable impact on future Antarctic ice loss (Golledge et al., 2017, https://doi.org/10.1002/2016GL072422). Subglacial lakes are thought to play a major role in the initiation of the rapid ice flow (Bell et al., 2007, https://doi.org/10.1038/nature05554). Satellite altimetry observations have even suggested several actively filling and draining subglacial lakes beneath the main trunk (B. E. Smith et al., 2009, https://doi.org/10.3189/002214309789470879). We present new data of the geometry of this glacier and investigate its basal properties employing radio‐echo sounding. Using ice sheet modeling, we were able to constrain estimates of radar absorption in the ice, but uncertainties remain large. The magnitude of the basal reflection coefficient is thus still poorly known. However, its spatial variability, in conjunction with additional indicators, can be used to infer the presence of subglacial water. We find no clear evidence of water at most of the previously proposed lake sites. Especially, locations, where altimetry detected active lakes, do not exhibit lake characteristics in radio‐echo sounding. We argue that lakes far upstream the main trunk are not triggering enhanced ice flow, which is also supported by modeled subglacial hydrology. Key Points: We find no clear evidence of water at most of the previously proposed lake sitesLocations, where altimetry detected active lakes, do not exhibit lake characteristics in RESLakes far upstream the main trunk are not triggering enhanced ice flow [ABSTRACT FROM AUTHOR]

We introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. The aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961-1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961-1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation-surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm until the year 2300, including our simulated additional sea level rise due to elevation-surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers. [ABSTRACT FROM AUTHOR]

The East Antarctic Ice Sheet contains a wealth of information that can be extracted from its internal architecture such as distribution of age, past flow features and surface and basal properties. Airborne radar surveys can sample this stratigraphic archive across broad areas. Here, we identify and trace key horizons across several radar surveys to obtain the stratigraphic information. We transfer the age-depth scales from ice cores to intersecting radar data. We then propagate these age scales across the ice sheet using the high fidelity continuity of the radar horizons. In Dronning Maud Land, including Dome Fuji, we mapped isochrones with ages of 38 ka and 74 ka. In the central region of East Antarctica around Dome Concordia, Vostok, and Dome Argus, we use isochrone ages of 38 ka, 48 ka, 90 ka, and 161 ka. Taking together both regions, we provide isochrone depths traced along a combined profile length of more than 40,000 km and discuss uncertainties of the obtained stratigraphy, as well as factors important to consider for further expansion. This dataset the most extensive distribution of internal horizons in East Antarctica to date. The isochrone depths are available on PANGAEA: https://doi.pangaea.de/10.1594/PANGAEA.895528. [ABSTRACT FROM AUTHOR]

Ice-stream dynamics are strongly controlled by processes taking place at the ice/bed interface where subglacial water both lubricates the base and saturates any existing, underlying sediment. Large parts of the former Eurasian ice sheet were underlain by thick sequences of soft, marine sediments and many areas are imprinted with geomorphological features indicative of fast flow and wet basal conditions. Here, we study the effect of subglacial water on past Eurasian ice-sheet dynamics by incorporating a thin-film model of basal water flow into the ice-sheet model SICOPOLIS and use it to better represent flow in temperate areas. The adjunction of subglacial hydrology results in a smaller icesheet building up over time and generally faster ice velocities, which consequently reduces the total area fraction of temperate basal ice and ice streaming areas. Minima in the hydraulic pressure potential, governing water flow, are used as indicators for potential locations of past subglacial lakes and a probability distribution of lake existence is presented based on estimated lake depth and longevity. [ABSTRACT FROM AUTHOR]

As ice flows over a subglacial lake, the drop in bed resistance leads to an increase in ice velocities and a draw down of isochrones and cold ice. The ice surface flattens as it adjusts to the lack of resisting forces at the base. The rapid transition in velocity induces changes in ice viscosity and releases deformation energy that can raise the temperature locally. Recent studies of Antarctic subglacial lakes indicate that many lakes experience very fast and possibly episodic drainage, during which the lake size is rapidly reduced as water flows out. Questions that arise are what effect this would have on internal layers within the ice and whether such past drainage events could be inferred from isochrone structures downstream. Here, we study the effect of a subglacial lake on ice dynamics as well as the influence that such short timescale drainage would have on the internal layers of the ice. To this end, we use a full Stokes, polythermal ice flow model. An enthalpy-gradient method is used to account for the evolution of temperature and water content within the ice. We find that a rapid transition between slow-moving ice outside the lake, and full sliding over the lake, can release considerable amounts of deformational energy, with the potential to form a temperate layer at depth in the transition zone. In addition, we provide an explanation for a characteristic surface feature commonly seen at the edges of subglacial lakes, a hummocky surface depression in the transition zone between little to full sliding. We also conclude that rapid changes in the horizontal extent of subglacial lakes and slippery patches, compared to the average ice column velocity, can create a traveling wave at depth within the isochrone structure that transfers downstream with the advection of ice, thus indicating the possibility of detecting past drainage events with ice penetrating radar. [ABSTRACT FROM AUTHOR]

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates.We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales. [ABSTRACT FROM AUTHOR]

The article focuses on a study conducted by the authors related to models of plan-view marine ice-sheets. It states that they conducted perturbation experiments which specified spatial variation in basal sliding and helped in evolution of grounding lines. It discusses intercomparison of the ice2sea MISMIP3d model. It highlights that examination of normal stresses perpendicular to the grounding lines reveals extension across the grounding lines.

The article presents a study that compares the Brunt Ice Shelf/Stancomb-Wills Ice Tongue system in Antarctica. The present flow regime of the ice shelf is simulated with two different models and equate the results with feature tracking. Two observed features were also integrated into the models with two different approaches, and show the impacts of variations in numerical values for viscosity and shear strength. A major result is that implementation of both kinds of the rift head to the same effects on the whole velocity field.

Lifecycle costs are frequently discussed as an aspect of operating costs. To do this, however, can be misleading for what use are low energy costs if the central components of the pump cannot stand up to arduous and demanding daily operation? Thomas G. Kleiner, business unit manager for Crane Process Flow Technologies GmbH in Düsseldorf, Germany, believes that when it comes to innovation in pumping, peristaltic pumps are undervalued. [Copyright &y& Elsevier]

• CFD simulations of pure substance condensation on annular low finned tubes. • Excellent agreement of simulation results with measurements. • Parameter study on fin slope using CFD. • Determination of condensate film along fin flank. • Increased smearing of flooding angle for a lower fin slope. • Optimum heat transfer for a medium fin slope. The use of annular low-finned tubes in tube bundle condensers greatly increases the efficiency. This enables enhanced heat coupling within chemical plants, reducing the overall CO 2 emission and power consumption. Due to the complex geometry of these tubes, no generalized condensation model is present so far. In this study, we use highly resolved computational fluid dynamics simulations to investigate pure substance condensation on said tubes with a condensation model which is independent of empirical parameters. Within these simulations, the condensate film is fully resolved and heat transfer coefficients are calculated providing the complete information about the condensation process on annular low-finned tubes for the first time. Additionally, a parameter study for the incline of the annular fin is provided. Therefore, computational fluid dynamics is used to predictively evaluate the influence of a single fin parameter, the incline of the fin, of annular low-finned tubes, for the first time. The simulations provide information about the film thickness along the fin flank and the flooding behavior of the tubes. An accurate fluid dynamic behavior of the two-phase flow is gained. Furthermore, the resulting heat transfer coefficients stand in excellent agreement to experimental data. [ABSTRACT FROM AUTHOR]

Ice shelf basal melt driven by oceanic heat is governing grounding line migration and ice dynamics of the Antarctic Ice Sheet (AIS). The recent ice sheet model intercomparison initMIP-Antarctica as part of ISMIP6 revealed a large spread in the simulated AIS response to ocean warming between the contributing models, partly based on models choices of parameterisation of ice-shelf basal melt. For the next phase of ISMIP6, the community agreed on using a prescribed ocean forcing parameterisation in combination with pre-processed forcing datasets including the observed present-day state of the ocean. Here we compare the simulations performed with the Parallel Ice Sheet Model (PISM) based on the different types of thermal forcing from the ocean. In order to assess the response to short time perturbations of oceanic heat content, we apply idealised perturbation as anomalies to the present-day ocean state on a decadal to centennial timescale. In the analysis, we focus on the regions of significant response. We evaluate surface elevation trends and ice sheet spreading rates along representative flow lines to investigate the ice sheet response in the context of dynamic thinning or rather thickening. This is serving as a baseline to understand the response to ocean forcing from GCMs with considerable decadal variability. [ABSTRACT FROM AUTHOR]

Recovery Glacier drains about 8% of the Antarctic ice sheet and feeds the Filchner Ice Shelf, that is suggested to be prone to changes in the oceanic forcing in future. Recent studies have used radar data to infer basal properties in particular with respect to the existence of subglacial water and its role in the genesis of the ice stream. Here, we take a different approach and use inverse modelling to infer the basal friction coefficient and basal drag. For this purpose, we apply the ice sheet model ISSM in higher order Blatter-Pattyn approximation using unstructured grids in high resolution in the main trunk and moderate resolution in the surrounding region ranging as far as adjacent ice streams. The temperature field is adopted from a paleo-spin up developed using PISM. Optimisation of the cost function is based on linear and logarithmic misfit of surface velocity and a regularisation of the friction parameter. We conducted a sensitivity study to various parameters, including the enhancement factor in the flow law to assess the role of deformation versus sliding, as well as initial parameters testing the robustness of the minimum. The general pattern of the distribution does not depend on minor variations in the parameters giving a high reliability of the solution. The main trunk is as expected found to have a smooth bed beside the ice falls, whereas we find a clear transition into a low friction regime just upstream the formerly suggested large lakes. Finally, we compare this to inferences of basal properties based on radar data. [ABSTRACT FROM AUTHOR]

Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18-17 ka and 6-5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland's climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately., Competing Interests: NO authors have competing interests.

Microstructures from deep ice cores reflect the dynamic conditions of the drill location as well as the thermodynamic history of the drill site and catchment area in great detail. Ice core parameters (crystal lattice-preferred orientation (LPO), grain size, grain shape), mesostructures (visual stratigraphy) as well as borehole deformation were measured in a deep ice core drilled at Kohnen Station, Dronning Maud Land (DML), Antarctica. These observations are used to characterize the local dynamic setting and its rheological as well as microstructural effects at the EDML ice core drilling site (European Project for Ice Coring in Antarctica in DML). The results suggest a division of the core into five distinct sections, interpreted as the effects of changing deformation boundary conditions from triaxial deformation with horizontal extension to bedrock-parallel shear. Region 1 (uppermost approx. 450 m depth) with still small macroscopic strain is dominated by compression of bubbles and strong strain and recrystallization localization. Region 2 (approx. 450-1700 m depth) shows a girdle-type LPO with the girdle plane being perpendicular to grain elongations, which indicates triaxial deformation with dominating horizontal extension. In this region (approx. 1000 m depth), the first subtle traces of shear deformation are observed in the shape-preferred orientation (SPO) by inclination of the grain elongation. Region 3 (approx. 1700-2030 m depth) represents a transitional regime between triaxial deformation and dominance of shear, which becomes apparent in the progression of the girdle to a single maximum LPO and increasing obliqueness of grain elongations. The fully developed single maximum LPO in region 4 (approx. 2030-2385 m depth) is an indicator of shear dominance. Region 5 (below approx. 2385 m depth) is marked by signs of strong shear, such as strong SPO values of grain elongation and strong kink folding of visual layers. The details of structural observations are compared with results from a numerical ice sheet model (PISM, isotropic) for comparison of strain rate trends predicted from the large-scale geometry of the ice sheet and borehole logging data. This comparison confirms the segmentation into these depth regions and in turn provides a wider view of the ice sheet.This article is part of the themed issue 'Microdynamics of ice'., (© 2016 The Authors.)

With two-photon absorption induced polymerization arbitrary three dimensional nano- and microstructures can be patterned directly into photoresists. We report on the fabrication of a low threshold organic semiconductor distributed feedback laser using the technique of two-photon absorption induced polymerization. A surface grating with 400 nm periodicity and 40 nm height modulation was fabricated by two-photon absorption induced polymerization in the organic-inorganic hybrid material ORMOCER. With structuring several stacked layers acting as a planar basis for the nanostructure microscopic substrate tilt can be compensated simply. This enabled us to uniformly nano-structure the surface grating over an area of 200 x 200 microm(2).