scholarly journals Data assimilation using a hybrid ice flow model

2011 ◽  
Vol 5 (2) ◽  
pp. 315-327 ◽  
Author(s):  
D. N. Goldberg ◽  
O. V. Sergienko

Abstract. Hybrid models, or depth-integrated flow models that include the effect of both longitudinal stresses and vertical shearing, are becoming more prevalent in dynamical ice modeling. Under a wide range of conditions they closely approximate the well-known First Order stress balance, yet are of computationally lower dimension, and thus require less intensive resources. Concomitant with the development and use of these models is the need to perform inversions of observed data. Here, an inverse control method is extended to use a hybrid flow model as a forward model. We derive an adjoint of a hybrid model and use it for inversion of ice-stream basal traction from observed surface velocities. A novel aspect of the adjoint derivation is a retention of non-linearities in Glen's flow law. Experiments show that in some cases, including those nonlinearities is advantageous in minimization of the cost function, yielding a more efficient inversion procedure.

2010 ◽  
Vol 4 (4) ◽  
pp. 2201-2231 ◽  
Author(s):  
D. N. Goldberg ◽  
O. V. Sergienko

Abstract. Hybrid models, or depth-integrated flow models that include the effect of both longitudinal stresses and vertical shearing, are becoming more prevalent in dynamical ice modeling. Under a wide range of conditions they closely approximate the well-known First Order stress balance, yet are of computationally lower dimension, and thus require less intensive resources. Concomitant with the development and use of these models is the need to perform inversions of observed data. Here, an inverse control method is extended to use a hybrid flow model as a forward model. We derive an adjoint of a hybrid model and use it for inversion of ice-stream basal traction from observed surface velocities. A novel aspect of the adjoint derivation is a retention of non-linearities in Glen's flow law. Experiments show that including those nonlinearities is advantageous in minimization of the cost function, yielding a more efficient inversion procedure.


2013 ◽  
Vol 6 (4) ◽  
pp. 6493-6568 ◽  
Author(s):  
R. Fischer ◽  
S. Nowicki ◽  
M. Kelley ◽  
G. A. Schmidt

Abstract. The method of elevation classes has proven to be a useful way for a low-resolution general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields, for use in one-way studies coupling GCMs and ice flow models. Past uses of elevation classes have been a cause of non-conservation of mass and energy, caused by inconsistency in regridding schemes chosen to regrid to the atmosphere vs. downscaling to the ice model. This causes problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded, and five transformations between those grids required by a typical coupled GCM–ice flow model. This paper shows how each of those transformations may be achieved in a consistent, conservative manner. These transformations are implemented in GLINT2, a library used to couple GCMs with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.


2021 ◽  
Author(s):  
Anya Schlich-Davies ◽  
Ann Rowan ◽  
Duncan Quincey ◽  
Andrew Ross ◽  
David Egholm

<p>Debris-covered glaciers in the Himalaya are losing mass more rapidly than expected. Quantifying and understanding the behaviour of these glaciers under climate change requires the use of numerical glacier models that represent the important feedbacks between debris transport, ice flow, and mass balance. However, these approaches have, so far, lacked a robust representation of the distributed mass balance forcing that is critical for making accurate simulations of ice volume change. This study forces a 3D higher-order ice flow model, with the outputs from an ensemble of distributed models of present day and future mass balance of Khumbu Glacier, Nepal. Distributed mass balance modelling, using the open access COupled Snowpack and Ice surface energy and mass balance model in PYthon (COSIPY) model (Sauter et al., 2020), was forced by three statistically downscaled climate models from the Coordinated Regional Climate Downscaling Experiment (CORDEX) project.</p><p>Climate models were selected based on their ability to reproduce observed present-day seasonality and to account for several future climate and monsoon scenarios, the latter being of particular importance for these summer-accumulation type glaciers. Two emission scenarios, RCP4.5 and RCP8.5, were also chosen to simulate glacier change to 2100. Statistical downscaling involved Quantile Mapping and Generalized Analog Regression Downscaling, and the efficacy of these approaches was informed by present day mass balance sensitivity studies. Downscaled daily climate data were trained with data from two weather stations to aid disaggregation to an hourly resolution.</p><p>The integration of the mass balance and ice flow models posed some interesting challenges. The COSIPY model was run as if Khumbu Glacier were a clean-ice glacier (with no supraglacial debris) with sub-debris ablation resolved in the ice flow model. The value of using distributed mass balance forcing is seen in the simulated present-day velocities in the Khumbu icefall, which give a better fit to remote-sensing observations than previous simulations using a simple elevation-dependent mass balance forcing. The simulated present-day glacier extent is considerably smaller than the existing glacier outline. The debris-covered tongue, known to be losing mass at an accelerating rate, is virtually absent from these results, and is indicative of a stagnant tongue that is now or very soon to be dynamically disconnected from the active upper reaches of Khumbu Glacier.</p>


2019 ◽  
Vol 65 (254) ◽  
pp. 1000-1010 ◽  
Author(s):  
Michael A. Imhof ◽  
Denis Cohen ◽  
Julien Seguinot ◽  
Andy Aschwanden ◽  
Martin Funk ◽  
...  

AbstractModelling paleo-glacier networks in mountain ranges on the millennial timescales requires ice flow approximations. Hybrid models calculating ice flow by combining vertical shearing (shallow ice approximation) and longitudinal stretching (shallow shelf approximation) have been applied to model paleo-glacier networks on steep terrain, yet their validity has not yet been assessed quantitatively. Moreover, hybrid models consistently yield higher ice thicknesses than Last Glacial Maximum geomorphological reconstructions in the European Alps. Here, we compare results based on the hybrid Parallel Ice Sheet Model (PISM) and the Stokes model Elmer/Ice on the Rhine Glacier, a catchment of the former European Alpine Icefield. For PISM, we also test two magnitudes of flux limitation in a scheme that reduces shearing velocities. We find that the flux limitation typically used in PISM yields significantly reduced shearing speeds and increases ice thicknesses by up to 500 m, partly explaining previous overestimations. However, reducing the ice flux limitation allows the hybrid model to minimize this mismatch and captures sliding speeds, ice thicknesses, ice extent and basal temperatures in close agreement with those obtained with the Stokes model.


2011 ◽  
Vol 57 (206) ◽  
pp. 1003-1016 ◽  
Author(s):  
Surendra Adhikari ◽  
Shawn J. Marshall

AbstractIn a two-dimensional (plane strain) glacier domain, gravity-driven ice flow is balanced by basal drag and the resistance associated with longitudinal stress gradients. The plane strain Stokes model accommodates both these resistances, whereas several simpler models only account for basal drag. Solving the Stokes equations is numerically challenging and computationally expensive, but simpler models may lead to unrealistic dynamical behaviour. Here, we propose a factor which can be introduced in shear-deformational flow models to yield results comparable to those from the plane strain Stokes model. As this factor adapts simpler models to capture the effects of missing dynamics, i.e. longitudinal stress gradients, we refer to it as the longitudinal stress (L-)factor. We assess the usefulness of this factor for idealized domains with complex basal topography and evolving geometry. We apply the model to Haig Glacier, Canadian Rockies, in order to present an illustration of how simulations of glacier response to climate forcing can be improved through the introduction of the L-factor in a shear-deformational flow model.


2016 ◽  
Vol 62 (231) ◽  
pp. 147-158 ◽  
Author(s):  
BRENT MINCHEW ◽  
MARK SIMONS ◽  
HELGI BJÖRNSSON ◽  
FINNUR PÁLSSON ◽  
MATHIEU MORLIGHEM ◽  
...  

AbstractThe mechanical properties of glacier beds play a fundamental role in regulating the sensitivity of glaciers to environmental forcing across a wide range of timescales. Glaciers are commonly underlain by deformable till whose mechanical properties and influence on ice flow are not well understood but are critical for reliable projections of future glacier states. Using synoptic-scale observations of glacier motion in different seasons to constrain numerical ice flow models, we study the mechanics of the bed beneath Hofsjökull, a land-terminating ice cap in central Iceland. Our results indicate that the bed deforms plastically and weakens following incipient summertime surface melt. Combining the inferred basal shear traction fields with a Coulomb-plastic bed model, we estimate the spatially distributed effective basal water pressure and show that changes in basal water pressure and glacier accelerations are non-local and non-linear. These results motivate an idealized physical model relating mean basal water pressure and basal slip rate wherein the sensitivity of glacier flow to changes in basal water pressure is inversely related to the ice surface slope.


Author(s):  
Nataliya Stoyanets ◽  
◽  
Mathias Onuh Aboyi ◽  

The article defines that for the successful implementation of an innovative project and the introduction of a new product into production it is necessary to use advanced technologies and modern software, which is an integral part of successful innovation by taking into account the life cycle of innovations. It is proposed to consider the general potential of the enterprise through its main components, namely: production and technological, scientific and technical, financial and economic, personnel and actual innovation potential. Base for the introduction of technological innovations LLC "ALLIANCE- PARTNER", which provides a wide range of support and consulting services, services in the employment market, tourism, insurance, translation and more. To form a model of innovative development of the enterprise, it is advisable to establish the following key aspects: the system of value creation through the model of cooperation with partners and suppliers; creating a value chain; technological platform; infrastructure, determine the cost of supply, the cost of activities for customers and for the enterprise as a whole. The system of factors of influence on formation of model of strategic innovative development of the enterprise is offered. The expediency of the cost of the complex of technological equipment, which is 6800.0 thousand UAH, is economically calculated. Given the fact that the company plans to receive funds under the program of socio-economic development of Sumy region, the evaluation of the effectiveness of the innovation project, the purchase of technological equipment, it is determined that the payback period of the project is 3 years 10 months. In terms of net present value (NPV), the project under study is profitable. The project profitability index (PI) meets the requirements for a positive decision on project implementation> 1.0. The internal rate of return of the project (IRR) also has a positive value of 22% because it exceeds the discount rate.


2012 ◽  
Vol 9 (1) ◽  
pp. 175-180
Author(s):  
Yu.D. Chashechkin

According to the results of visualization of streams, the existence of structures in a wide range of scales is noted: from galactic to micron. The use of a fundamental system of equations is substantiated based on the results of comparing symmetries of various flow models with the usage of theoretical group methods. Complete solutions of the system are found by the methods of the singular perturbations theory with a condition of compatibility, which determines the characteristic equation. A comparison of complete solutions with experimental data shows that regular solutions characterize large-scale components of the flow, a rich family of singular solutions describes formation of the thin media structure. Examples of calculations and observations of stratified, rotating and multiphase media are given. The requirements for the technique of an adequate experiment are discussed.


1986 ◽  
Vol 51 (11) ◽  
pp. 2489-2501
Author(s):  
Benitto Mayrhofer ◽  
Jana Mayrhoferová ◽  
Lubomír Neužil ◽  
Jaroslav Nývlt

A model is derived for a multi-stage crystallization with cross-current flows of the solution and the crystals being purified. The purity of the product is compared with that achieved in the countercurrent arrangement. A suitable function has been set up which allows the cross-current and countercurrent flow models to be compared and reduces substantially the labour of computation for the countercurrent arrangement. Using the recrystallization of KAl(SO4)2.12 H2O as an example, it is shown that, when the cross-current and countercurrent processes are operated at the same output, the countercurrent arrangement is more advantageous because its solvent consumption is lower.


Sign in / Sign up

Export Citation Format

Share Document