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.

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ∼4 lenses per square degree down to i(AB) = 24.8. This is equivalent to ∼60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

Abstract. The standard viscous–plastic (VP) sea ice model with an elliptical yield curve and a normal flow rule has at least two issues. First, it does not simulate fracture angles below 30∘ in uni-axial compression, in contrast with observations of linear kinematic features (LKFs) in the Arctic Ocean. Second, there is a tight, but unphysical, coupling between the fracture angle, post-fracture deformation, and the shape of the yield curve. This tight coupling was identified as the reason for the overestimation of fracture angles. In this paper, these issues are addressed by removing the normality constraint on the flow rule in the standard VP model. The new rheology is tested in numerical uni-axial loading tests. To this end, an elliptical plastic potential – which defines the post-fracture deformations, or flow rule – is introduced independently of the elliptical yield curve. As a consequence, the post-fracture deformation is decoupled from the mechanical strength properties of the ice. We adapt Roscoe's angle theory, which is based on observations of granular materials, to the context of sea ice modeling. In this framework, the fracture angles depend on both yield curve and plastic potential parameters. This new formulation predicts accurately the results of the numerical experiments with a root-mean-square error below 1.3∘. The new rheology allows for angles of fracture smaller than 30∘ in uni-axial compression. For instance, a plastic potential with an ellipse aspect ratio smaller than 2 (i.e., the default value in the standard viscous–plastic model) can lead to fracture angles as low as 22∘. Implementing an elliptical plastic potential in the standard VP sea ice model requires only small modifications to the standard VP rheology. The momentum equations with the modified rheology, however, are more difficult to solve numerically. The independent plastic potential solves the two issues with VP rheology addressed in this paper: in uni-axial loading experiments, it allows for smaller fracture angles, which fall within the range of satellite observations, and it decouples the angle of fracture and the post-fracture deformation from the shape of the yield curve. The orientation of the post-fracture deformation along the fracture lines (convergence and divergence), however, is still controlled by the shape of the plastic potential and the location of the stress state on the yield curve. A non-elliptical plastic potential would be required to change the orientation of deformation and to match deformation statistics derived from satellite measurements.

Abstract. The standard viscous-plastic (VP) sea ice model with an elliptical yield curve and normal flow rule does not simulate fracture angles below 30° in uni-axial compression, in stark contrast with observations of Linear Kinematic Features (LKFs) in the Arctic Ocean. In this paper, we remove the normality constraint in the standard VP model and study its impact on the fracture angle in a simple uni-axial compressive loading test. To this end, we introduce a plastic potential independent of the yield curve that defines the post-fracture deformations or flow rule. The numerical experiments show that the fracture angle strongly depends on the flow rule details. For instance, a plastic potential with an ellipse aspect ratio smaller than that of the standard ellipse gives fracture angles that are as low as 22°. A newly adapted theory – based on one developed from observations of granular material – predicts numerical simulations of the fracture angles for plastic materials with a normal or non-normal flow rule with a root-mean-square error below 1.3°. Implementing an elliptical plastic potential in the standard VP sea ice model requires only minor modifications. The modified rheology, however, takes longer to solve numerically for a fixed level of numerical convergence. In conclusion, the use of a plastic potential addresses several issues with the standard VP rheology: the fracture angle can be reduced to values within the range of satellite observations and it can be decoupled from the exact shape of the yield curve. Furthermore, a different plastic potential function will be required to change the post-fracture deformation along the fracture lines (convergence or divergence) and to make the fracture angle independent on the confining pressure (as in observations).

No abstract