Continuum sea ice rheology determined from subcontinuum mechanics

[1] A method is presented to calculate the continuum-scale sea ice stress as an imposed, continuum-scale strain-rate is varied. The continuum-scale stress is calculated as the area-average of the stresses within the floes and leads in a region (the continuum element). The continuum-scale stress depends upon: the imposed strain rate; the subcontinuum scale, material rheology of sea ice; the chosen configuration of sea ice floes and leads; and a prescribed rule for determining the motion of the floes in response to the continuum-scale strain-rate. We calculated plastic yield curves and flow rules associated with subcontinuum scale, material sea ice rheologies with elliptic, linear and modified Coulombic elliptic plastic yield curves, and with square, diamond and irregular, convex polygon-shaped floes. For the case of a tiling of square floes, only for particular orientations of the leads have the principal axes of strain rate and calculated continuum-scale sea ice stress aligned, and these have been investigated analytically. The ensemble average of calculated sea ice stress for square floes with uniform orientation with respect to the principal axes of strain rate yielded alignment of average stress and strain-rate principal axes and an isotropic, continuum-scale sea ice rheology. We present a lemon-shaped yield curve with normal flow rule, derived from ensemble averages of sea ice stress, suitable for direct inclusion into the current generation of sea ice models. This continuum-scale sea ice rheology directly relates the size (strength) of the continuum-scale yield curve to the material compressive strength

Fracture of andesite in the brittle and brittle-ductile transition regimes

see Abstract volum

Continuum sea ice rheology determined from subcontinuum mechanics

[1] A method is presented to calculate the continuum-scale sea ice stress as an imposed, continuum-scale strain-rate is varied. The continuum-scale stress is calculated as the area-average of the stresses within the floes and leads in a region (the continuum element). The continuum-scale stress depends upon: the imposed strain rate; the subcontinuum scale, material rheology of sea ice; the chosen configuration of sea ice floes and leads; and a prescribed rule for determining the motion of the floes in response to the continuum-scale strain-rate. We calculated plastic yield curves and flow rules associated with subcontinuum scale, material sea ice rheologies with elliptic, linear and modified Coulombic elliptic plastic yield curves, and with square, diamond and irregular, convex polygon-shaped floes. For the case of a tiling of square floes, only for particular orientations of the leads have the principal axes of strain rate and calculated continuum-scale sea ice stress aligned, and these have been investigated analytically. The ensemble average of calculated sea ice stress for square floes with uniform orientation with respect to the principal axes of strain rate yielded alignment of average stress and strain-rate principal axes and an isotropic, continuum-scale sea ice rheology. We present a lemon-shaped yield curve with normal flow rule, derived from ensemble averages of sea ice stress, suitable for direct inclusion into the current generation of sea ice models. This continuum-scale sea ice rheology directly relates the size (strength) of the continuum-scale yield curve to the material compressive strength

Surface structure and stability of the Larsen C ice shelf, Antarctic Peninsula

A structural glaciological description and analysis of surface morphological features of the Larsen C ice shelf, Antarctic Peninsula, is derived from satellite images spanning the period 1963-2007. The data are evaluated in two time ranges: a comparison of a 1963 satellite image photomosaic with a modern digital mosaic compiled using 2003/04 austral summer data; and an image series since 2003 showing recent evolution of the shelf. We map the ice-shelf edge, rift swarms, crevasses and crevasse traces, and linear longitudinal structures (called 'flow stripes' or 'streak lines'). The latter are observed to be continuous over distances of up to 200km from the grounding line to the ice-shelf edge, with little evidence of changes in pattern over that distance. Integrated velocity measurements along a flowline indicate that the shelf has been stable for similar to 560years in the mid-shelf area. Linear longitudinal features may be grouped into 12 units, each related to one or a small group of outlet feeder glaciers to the shelf. We observe that the boundaries between these flow units often mark rift terminations. The boundary zones originate upstream at capes, islands or other suture areas between outlet glaciers. In agreement with previous work, our findings imply that rift terminations within such suture zones indicate that they contain anomalously soft ice. We thus suggest that suture zones within the Larsen C ice shelf, and perhaps within ice shelves more generally, may act to stabilize them by reducing regional stress intensities and thus rates of rift lengthening

Developing a Dynamic Web-GIS Based Landslide Early Warning System for the Chittagong Metropolitan Area, Bangladesh

This article aims to develop a Web-GIS based landslide early warning system (EWS) for the Chittagong Metropolitan Area (CMA), Bangladesh, where, in recent years, rainfall-induced landslides have caused great losses of lives and property. A method for combining static landslide susceptibility maps and rainfall thresholds is proposed by introducing a purposely-build hazard matrix. To begin with, eleven factor maps: soil permeability; surface geology; landcover; altitude; slope; aspect; distance to stream; fault line; hill cut; road cut; and drainage network along with a detailed landslide inventory map were produced. These maps were used, and four methods were applied: artificial neural network (ANN); multiple regressions; principal component analysis; and support vector machine to produce landslide susceptibility maps. After model validation, the ANN map was found best fitting and was classified into never warning, low, medium, and high susceptibility zones. Rainfall threshold analysis (1960–2017) revealed consecutive 5-day periods of rainfall of 71–282 mm could initiate landslides in CMA. Later, the threshold was classified into three rainfall rates: low rainfall (70–160 mm), medium rainfall (161–250 mm), and high rainfall (>250 mm). Each landslide was associated with a hazard class (no warning vs. warning state) based on the assumption that the higher the susceptibility, the lower the rainfall. Finally, the EWS was developed using various libraries and frameworks that is connected with a reliable online-based weather application programming interface. The system is publicly available, dynamic, and replicable to similar contexts and is able to disseminate alerts five days in advance via email notifications. The proposed EWS is novel and the first of its kind in Bangladesh, and can be applied to mitigate landslide disaster risks

The 2016 Kumamoto Earthquakes: Cascading Geological Hazards and Compounding Risks

A sequence of two strike-slip earthquakes occurred on April 14 and 16, 2016 in the intraplate region of Kyushu Island, Japan, apart from subduction zones, and caused significant damage and disruption to the Kumamoto region. The analyses of regional seismic catalog and available strong motion recordings reveal striking characteristics of the events, such as migrating seismicity, earthquake surface rupture, and major foreshock-mainshock earthquake sequences. To gain valuable lessons from the events, a UK Earthquake Engineering Field Investigation Team (EEFIT) was dispatched to Kumamoto, and earthquake damage surveys were conducted to relate observed earthquake characteristics to building and infrastructure damage caused by the earthquakes. The lessons learnt from the reconnaissance mission have important implications on current seismic design practice regarding the required seismic resistance of structures under multiple shocks and the seismic design of infrastructure subject to large ground deformation. The observations also highlight the consequences of cascading geological hazards on community resilience. To share the gathered damage data widely, geo-tagged photos are organized using Google Earth and the kmz file is made publicly available.The work is funded by the EPSRC grant (EP/I01778X/1) for the Earthquake Engineering Field Investigation Team (EEFIT). The financial supports for industrial members (GC, LH, LK, and RM) are provided by Arup, Mott MacDonald, and Willis

Computations of Viking Lander Capsule Hypersonic Aerodynamics with Comparisons to Ground and Flight Data

Comparisons are made between the LAURA Navier-Stokes code and Viking Lander Capsule hypersonic aerodynamics data from ground and flight measurements. Wind tunnel data are available for a 3.48 percent scale model at Mach 6 and a 2.75 percent scale model at Mach 10.35, both under perfect gas air conditions. Viking Lander 1 aerodynamics flight data also exist from on-board instrumentation for velocities between 2900 and 4400 m/sec (Mach 14 to 23.3). LAURA flowfield solutions are obtained for the geometry as tested or flown, including sting effects at tunnel conditions and finite-rate chemistry effects in flight. Using the flight vehicle center-of-gravity location (trim angle approx. equals -11.1 deg), the computed trim angle at tunnel conditions is within 0.31 degrees of the angle derived from Mach 6 data and 0.13 degrees from the Mach 10.35 trim angle. LAURA Mach 6 trim lift and drag force coefficients are within 2 percent of measured data, and computed trim lift-to-drag ratio is within 4 percent of the data. Computed trim lift and drag force coefficients at Mach 10.35 are within 5 percent and 3 percent, respectively, of wind tunnel data. Computed trim lift-to-drag ratio is within 2 percent of the Mach 10.35 data. Using the nominal density profile and center-of-gravity location, LAURA trim angle at flight conditions is within 0.5 degrees of the total angle measured from on-board instrumentation. LAURA trim lift and drag force coefficients at flight conditions are within 7 and 5 percent, respectively, of the flight data. Computed trim lift-to-drag ratio is within 4 percent of the data. Computed aerodynamics sensitivities to center-of-gravity location, atmospheric density, and grid refinement are generally small. The results will enable a better estimate of aerodynamics uncertainties for future Mars entry vehicles where non-zero angle-of-attack is required

Fracture and damage localization in volcanic edifice rocks from El Hierro, Stromboli and Tenerife

© 2018 The Author(s). We present elastic wave velocity and strength data from a suite of three volcanic rocks taken from the volcanic edifices of El Hierro and Tenerife (Canary Islands, Spain), and Stromboli (Aeolian Islands, Italy). These rocks span a range of porosity and are taken from volcanoes that suffer from edifice instability. We measure elastic wave velocities at known incident angles to the generated through-going fault as a function of imposed strain, and examine the effect of the damage zone on P-wave velocity. Such data are important as field measurements of elastic wave tomography are key tools for understanding volcanic regions, yet hidden fractures are likely to have a significant effect on elastic wave velocity. We then use elastic wave velocity evolution to calculate concomitant crack density evolution which ranges from 0 to 0.17: highest values were correlated to the damage zone in rocks with the highest initial porosity