Empirical Scalar Tertiary Anisotropy Regime (ESTAR)

Description

The ESTAR (Empirical Scalar Tertiary Anisotropy Regime) flow relation [Budd2013,Graham2018] is a generalized constitutive relation for polycrystalline ice in steady-state (tertiary) flow. It is a scalar power law formulation based on tertiary creep rates from laboratory experiments of ice deformation under a variety of simple shear and compression stresses. While mathematically isotropic, the ESTAR flow relation describes the deformation of ice with a flow-compatible induced anisotropy — i.e. ice that has a developed anisotropic fabric that is a function of the underlying stress regime (i.e. the relative proportion of simple shear and compression stresses). The origins of ESTAR, including the laboratory experiments than contributed to its development, its derivation, and underlying assumptions are discussed in [Budd2013] and [Graham2018].

Equations

Ice is treated as a purely viscous incompressible material [Cuffey2010], such that its material constitutive relation can be written:

where:

• is the deviatoric stress tensor (Pa)
• is the ice effective viscosity (Pa~s)
• is the strain rate tensor (s^-1)

The ESTAR flow relation viscosity can be written:

where:

• is the ice hardness or rigidity. Note that , where is the temperature-dependent flow rate parameter and is the temperature relative to the pressure dependent melting point of ice.
• is an enhancement factor that characterizes the relative proportion of simple shear and compression stresses via the shear fraction

The most notable difference between the Glen and ESTAR flow relations is realized in the form of the enhancement factor, which for the ESTAR flow relation is , given by:

Here, and are the enhancement factors above the minimum (secondary) deformation rate for isotropic ice under compression alone or simple shear alone, respectively. Laboratory evidence suggests that a suitable ratio of is [Treverrow2012]. The shear fraction characterizes the contribution of simple shear to the effective stress. The collinear nature of the ESTAR flow relation allows to be expressed equivalently in terms of stresses and strain rates. The strain rate formulation is more convenient for Stokes flow modeling, and can be written:

where (s^-1) is the magnitude of the shear strain rate on the local non-rotating shear plane. The local non-rotating shear plane contains the velocity vector and the vorticity vector associated solely with deformation, rather than local rigid body rotation. See [Graham2018] for details.

For comparison with the ESTAR viscosity, the Glen flow relation viscosity can be written:

where is a constant enhancement factor. For the standard Glen flow relation (the matice class in ISSM), ; to specify values of , the matenhancedice class can be used.

Model parameters

The parameters relevant to the ESTAR flow relation (the matestar class in ISSM) can be displayed by running:

>> md.materials

• md.materials.rheology_B: temperature-dependent flow relation parameter (NaN means no constraint)
• md.materials.rheology_Ec: compression enhancement factor
• md.materials.rheology_Es: simple shear enhancement factor
• md.materials.rheology_law: law for the temperature dependence of the rheology (None means no temperature dependence; default is Paterson)

Using the ESTAR flow relation

The ESTAR flow relation may be specified by:

>> md.materials = matestar();

In this case, values for , , and should be explicitly set.

Alternatively, the ESTAR flow relation may be specified from conversion of a Glen type relation by the following:

>> md.materials = matestar(md.materials);

The argument is the materials class of the model. This will set the same value for as for the Glen flow model default, with and .

Using the enhanced Glen flow relation

It is possible to use an alternative Glen flow relation with an explicit enhancement factor, in a similar way to the ESTAR class, as follows:

>> md.materials = matenhancedice();

in which and should be explicitly set, or as:

>> md.materials = matenhancedice(md.materials);

in which is inherited from the default Glen flow model and 1.

References