scholarly journals The Effect Of Rheology On Seasonal Sea-Ice Simulations

1990 ◽  
Vol 14 ◽  
pp. 340-340
Author(s):  
Bill Ip ◽  
W.D. Hibler ◽  
Greg Flato

On the seasonal time scales relevant to numerical investigations of climate, the rheology used in large-scale sea ice models significantly affects the ice thickness build-up and ice velocity fields. Plastic rheologies with a normal flow rule have been used to-date in seasonal dynamic thermodynamic simulations. These rheologies have proved useful in simulating discontinuous slip near the coast while still supplying relatively robust velocity fields in the central Arctic Basin. However, as indicated by limited numerical sensitivity studies with different types of elliptical yield curves, the amount of shear strength significantly affects the ice build-up and can possibly cause a stoppage of the ice outflow through Fram Strait. In addition to the shear strength issue, there is also the possibility that non-normal flow rule rheologies, such as the Mohr Coulomb failure criterion used in soil mechanics, may cause somewhat different types of flow patterns, especially in the Fram Strait region. However, to date no seasonal simulations with such non-normal flow rule rheologies have been carried out.In order to investigate the role of different rheologies on the large-scale flow patterns in the Arctic Basin, a more general numerical scheme than that of Hibler (1979) is developed, which allows the simulation of the dynamic thermodynamic behavior of sea ice with a wide variety of different non linear rheologies. Using this numerical scheme, comparative simulations are carried out to seasonal equilibrium with several variations of the Mohr Coulomb rheology and compared to the more standard Elliptical yield curve results. In particular, the main control Mohr Coulomb case is a capped rheology in which the shear strength is taken to be proportional to the compressive stress. In this capped case only shear flow is allowed until a maximum allowable compressive stress is reached. This cap strength is parameterized to be a function of the ice thickness and compactness. For comparison, a simulation with a very large cap strength is also carried out, and an experiment with a similar compressive cap but much lower shear strength. Overall the results are analyzed to determine the sensitivity of the ice build-up to flow rule and shear strength magnitude. In addition special attention is given to the character of the flow and stoppage (if any) through Fram Strait.

1990 ◽  
Vol 14 ◽  
pp. 340
Author(s):  
Bill Ip ◽  
W.D. Hibler ◽  
Greg Flato

On the seasonal time scales relevant to numerical investigations of climate, the rheology used in large-scale sea ice models significantly affects the ice thickness build-up and ice velocity fields. Plastic rheologies with a normal flow rule have been used to-date in seasonal dynamic thermodynamic simulations. These rheologies have proved useful in simulating discontinuous slip near the coast while still supplying relatively robust velocity fields in the central Arctic Basin. However, as indicated by limited numerical sensitivity studies with different types of elliptical yield curves, the amount of shear strength significantly affects the ice build-up and can possibly cause a stoppage of the ice outflow through Fram Strait. In addition to the shear strength issue, there is also the possibility that non-normal flow rule rheologies, such as the Mohr Coulomb failure criterion used in soil mechanics, may cause somewhat different types of flow patterns, especially in the Fram Strait region. However, to date no seasonal simulations with such non-normal flow rule rheologies have been carried out. In order to investigate the role of different rheologies on the large-scale flow patterns in the Arctic Basin, a more general numerical scheme than that of Hibler (1979) is developed, which allows the simulation of the dynamic thermodynamic behavior of sea ice with a wide variety of different non linear rheologies. Using this numerical scheme, comparative simulations are carried out to seasonal equilibrium with several variations of the Mohr Coulomb rheology and compared to the more standard Elliptical yield curve results. In particular, the main control Mohr Coulomb case is a capped rheology in which the shear strength is taken to be proportional to the compressive stress. In this capped case only shear flow is allowed until a maximum allowable compressive stress is reached. This cap strength is parameterized to be a function of the ice thickness and compactness. For comparison, a simulation with a very large cap strength is also carried out, and an experiment with a similar compressive cap but much lower shear strength. Overall the results are analyzed to determine the sensitivity of the ice build-up to flow rule and shear strength magnitude. In addition special attention is given to the character of the flow and stoppage (if any) through Fram Strait.


2015 ◽  
Vol 9 (1) ◽  
pp. 269-283 ◽  
Author(s):  
R. Lindsay ◽  
A. Schweiger

Abstract. Sea ice thickness is a fundamental climate state variable that provides an integrated measure of changes in the high-latitude energy balance. However, observations of mean ice thickness have been sparse in time and space, making the construction of observation-based time series difficult. Moreover, different groups use a variety of methods and processing procedures to measure ice thickness, and each observational source likely has different and poorly characterized measurement and sampling errors. Observational sources used in this study include upward-looking sonars mounted on submarines or moorings, electromagnetic sensors on helicopters or aircraft, and lidar or radar altimeters on airplanes or satellites. Here we use a curve-fitting approach to determine the large-scale spatial and temporal variability of the ice thickness as well as the mean differences between the observation systems, using over 3000 estimates of the ice thickness. The thickness estimates are measured over spatial scales of approximately 50 km or time scales of 1 month, and the primary time period analyzed is 2000–2012 when the modern mix of observations is available. Good agreement is found between five of the systems, within 0.15 m, while systematic differences of up to 0.5 m are found for three others compared to the five. The trend in annual mean ice thickness over the Arctic Basin is −0.58 ± 0.07 m decade−1 over the period 2000–2012. Applying our method to the period 1975–2012 for the central Arctic Basin where we have sufficient data (the SCICEX box), we find that the annual mean ice thickness has decreased from 3.59 m in 1975 to 1.25 m in 2012, a 65% reduction. This is nearly double the 36% decline reported by an earlier study. These results provide additional direct observational evidence of substantial sea ice losses found in model analyses.


2018 ◽  
Author(s):  
Damien Ringeisen ◽  
Nils Hutter ◽  
Martin Losch ◽  
L. Bruno Tremblay

Abstract. Recent high resolution pan-Arctic sea ice simulations show fracture patterns (Linear Kinematic Features – LKFs) that are typical of granular materials but with intersection (fracture) angles wider than those observed from high-resolution satellite images (with a modal value of θ = 20°). In this article, We investigate the mechanism of formation and parameter dependencies of ice fracture in simple numerical bi-axial test on a 8 km x 25 km ice floe at an unprecedented resolution of 25m for two different yield curves: an elliptical (VP) and a Coulombic yield curve both with normal flow rule. In the standardVP model, the simulated angle of fracture is θ = 33.9°, compared to 20° in observations. The dependence of the angle of fracture on the ice shear strength is also contrary to that of typical granular materials with larger angle of fracture for higher shear strength – think of a wet sand castle with steeper walls than a dry sand castle. In this model, the divergence along the fracture lines (or LKFs) is entirely dictated by the ice shear strength used in the model with high shear strength resulting in convergence along LKFs and low shear strength resulting in divergence along LKFs. This is again contrary to typical granular materials where divergence (or dilation) is linked with the orientation of contacts normals that oppose the flow with divergence present for larger shear resistance and convergence for lower shear resistance. Moreover, the angle of fracture depends on the confining pressure in the uni-axial test with more convergence as the confining pressure increases, again contrary to granular material that have an angle of fracture that is independent of the confining pressure. We note that all three behaviors of the VP model are linked with the use of an associative (normal) flow rule. In the Coulombic model, the angle of fracture is smaller (θ = 23.5°), but the solution is unstable when the compressive stresses are too large because of the discontinuity between the straight limbs of the yield curve and the elliptical capping. Our results show that while the VP model gives angles of fracture that are visually correct, the bias in the magnitude of the angle of fracture and the physical dependencies of the angle of fracture on mechanical strength parameters and stress fields couple the sea ice mechanical strength parameters, the sea-ice drift, sea-ice deformation (strain-rate) field in an inconsistent way. We consider this evidence to move away from the elliptical yield curve and associative (normal) flow rule, a deformation law that is not applicable to pressure-sensitive and dilatant granular material such as sea ice.


2014 ◽  
Vol 8 (4) ◽  
pp. 4545-4580 ◽  
Author(s):  
R. Lindsay ◽  
A. Schweiger

Abstract. Sea ice thickness is a fundamental climate state variable that provides an integrated measure of changes in the high-latitude energy balance. However, observations of ice thickness have been sparse in time and space making the construction of observation-based time series difficult. Moreover, different groups use a variety of methods and processing procedures to measure ice thickness and each observational source likely has different and poorly characterized measurement and sampling biases. Observational sources include upward looking sonars mounted on submarines or moorings, electromagnetic sensors on helicopters or aircraft, and lidar or radar altimeters on airplanes or satellites. Here we use a curve-fitting approach to evaluate the systematic differences between eight different observation systems in the Arctic Basin. The approach determines the large-scale spatial and temporal variability of the ice thickness as well as the mean differences between the observation systems using over 3000 estimates of the ice thickness. The thickness estimates are measured over spatial scales of approximately 50 km or time scales of 1 month and the primary time period analyzed is 2000–2013 when the modern mix of observations is available. Good agreement is found between five of the systems, within 0.15 m, while systematic differences of up to 0.5 m are found for three others compared to the five. The trend in annual mean ice thickness over the Arctic Basin is −0.58 ± 0.07 m decade−1 over the period 2000–2013, while the annual mean ice thickness for the central Arctic Basin alone (the SCICEX Box) has decreased from 3.45 m in 1975 to 1.11 m in 2013, a 68% reduction. This is nearly double the 36% decline reported by an earlier study. These results provide additional direct observational confirmation of substantial sea ice losses found in model analyses.


1983 ◽  
Vol 29 (101) ◽  
pp. 191-195 ◽  
Author(s):  
Ronald B. Smith

Abstract The derivation of a constitutive law to describe large-scale sea-ice deformation commonly uses the so-called “flow rule”. This method introduces a mathematically imposed relationship between shear strength and dilatation which is not based on physical postulates. In this note a more general procedure is described which uses the Reiner–Rivlin equation as a starting point. The method is illustrated by deriving a particular constitutive law and applying it to a simple problem of sea ice blown against a shoreline.


2019 ◽  
Vol 13 (4) ◽  
pp. 1167-1186 ◽  
Author(s):  
Damien Ringeisen ◽  
Martin Losch ◽  
L. Bruno Tremblay ◽  
Nils Hutter

Abstract. Recent high-resolution pan-Arctic sea ice simulations show fracture patterns (linear kinematic features or LKFs) that are typical of granular materials but with wider fracture angles than those observed in high-resolution satellite images. Motivated by this, ice fracture is investigated in a simple uni-axial loading test using two different viscous–plastic (VP) rheologies: one with an elliptical yield curve and a normal flow rule and one with a Coulombic yield curve and a normal flow rule that applies only to the elliptical cap. With the standard VP rheology, it is not possible to simulate fracture angles smaller than 30∘. Further, the standard VP model is not consistent with the behavior of granular material such as sea ice because (1) the fracture angle increases with ice shear strength; (2) the divergence along the fracture lines (or LKFs) is uniquely defined by the shear strength of the material with divergence for high shear strength and convergent with low shear strength; (3) the angle of fracture depends on the confining pressure with more convergence as the confining pressure increases. This behavior of the VP model is connected to the convexity of the yield curve together with use of a normal flow rule. In the Coulombic model, the angle of fracture is smaller (θ=23∘) and grossly consistent with observations. The solution, however, is unstable when the compressive stress is too large because of non-differentiable corners between the straight limbs of the Coulombic yield curve and the elliptical cap. The results suggest that, although at first sight the large-scale patterns of LKFs simulated with a VP sea ice model appear to be realistic, the elliptical yield curve with a normal flow rule is not consistent with the notion of sea ice as a pressure-sensitive and dilatant granular material.


1983 ◽  
Vol 29 (101) ◽  
pp. 191-195
Author(s):  
Ronald B. Smith

AbstractThe derivation of a constitutive law to describe large-scale sea-ice deformation commonly uses the so-called “flow rule”. This method introduces a mathematically imposed relationship between shear strength and dilatation which is not based on physical postulates. In this note a more general procedure is described which uses the Reiner–Rivlin equation as a starting point. The method is illustrated by deriving a particular constitutive law and applying it to a simple problem of sea ice blown against a shoreline.


1991 ◽  
Vol 15 ◽  
pp. 17-25 ◽  
Author(s):  
Chi F. Ip ◽  
William D. Hibler ◽  
Gregory M. Flato

A generalized numerical model which allows for a variety of non-linear rheologies is developed for the seasonal simulation of sea-ice circulation and thickness. The model is used to investigate the effects (such as the role of shear stress and the existence of a flow rule) of different rheologies on the ice-drift pattern and build-up in the Arctic Basin. Differences in local drift seem to be closely related to the amount of allowable shear stress. Similarities are found between the elliptical and square cases and between the Mohr-Coulomb and cavitating fluid cases. Comparisons between observed and simulated buoy drift are made for several buoy tracks in the Arctic Basin. Correlation coefficients to the observed buoy drift range from 0.83 for the cavitating fluid to 0.86 for the square rheology. The average ratio of buoy-drift distance to average model-drift distance for several buoys is 1.15 (square), 1.18 (elliptical), 1.30 (Mohr-Coulomb) and 1.40 (cavitating fluid).


2019 ◽  
Author(s):  
Maciej Miernecki ◽  
Lars Kaleschke ◽  
Nina Maaß ◽  
Stefan Hendricks ◽  
Sten Schmidl Søbjrg

Abstract. Sea ice thickness measurements with L-band radiometry is a technique which allows daily, weather-independent monitoring of the polar sea ice cover. The sea-ice thickness retrieval algorithms relay on the sensitivity of the L-band brightness temperature to sea-ice thickness. In this work, we investigate the decimetre-scale surface roughness as a factor influencing the L-band emissions from sea ice. We used an airborne laser scanner to construct a digital elevation model of the sea ice surface. We found that the probability density function of surface slopes is exponential for a range of degrees of roughness. Then we applied the geometrical optics, bounded with the MIcrowave L-band LAyered Sea ice emission model in the Monte Carlo simulation to simulate the effects of surface roughness. According to this simulations, the most affected by surface roughness is the vertical polarization around Brewster's angle, where the decrease in brightness temperature can reach 8 K. The vertical polarization for the same configuration exhibits a 4 K increase. The near-nadir angles are little affected, up to 2.6 K decrease for the most deformed ice. Overall the effects of large-scale surface roughness can be expressed as a superposition of two factors: the change in intensity and the polarization mixing. The first factor depends on surface permittivity, second shows little dependence on it. Comparison of the brightness temperature simulations with the radiometer data does not yield definite results.


Sign in / Sign up

Export Citation Format

Share Document