Three-Dimensional Mechanics of Yakutat Convergence in the Southern Alaskan Plate Corner

Three-dimensional numerical models are used to investigate the mechanical evolution of the southern Alaskan plate corner where the Yakutat and the Pacific plates converge on the North American plate. The evolving model plate boundary consists of Convergent, Lateral, and Subduction subboundaries with flow separation of incoming material into upward or downward trajectories forming dual, nonlinear advective thermal/mechanical anomalies that fix the position of major subaerial mountain belts. The model convergent subboundary evolves into two teleconnected orogens: Inlet and Outlet orogens form at locations that correspond with the St. Elias and the Central Alaska Range, respectively, linked to the East by the Lateral boundary. Basins form parallel to the orogens in response to the downward component of velocity associated with subduction. Strain along the Lateral subboundary varies as a function of orogen rheology and magnitude and distribution of erosion. Strain-dependent shear resistance of the plate boundary associated with the shallow subduction zone controls the position of the Inlet orogen. The linkages among these plate boundaries display maximum shear strain rates in the horizontal and vertical planes where the Lateral subboundary joins the Inlet and Outlet orogens. The location of the strain maxima shifts with time as the separation of the Inlet and Outlet orogens increases. The spatiotemporal predictions of the model are consistent with observed exhumation histories deduced from thermochronology, as well as stratigraphic studies of synorogenic deposits. In addition, the complex structural evolution of the St Elias region is broadly consistent with the predicted strain field evolution. Citation: Koons, P. O., B. P. Hooks, T. Pavlis, P. Upton, and A. D. Barker (2010), Three-dimensional mechanics of Yakutat convergence in the southern Alaskan plate corner, Tectonics, 29, TC4008, doi: 10.1029/2009TC002463.

Buckling Rate and Overhang Development at a Calving Face

Using the finite-element we have modeled the stress field near the calving face of an idealized tidewater glacier under a variety of assumptions about submarine calving-face height, subaerial calving-face height, and ice rheology These simulations all suggest that a speed maximum should be present at the calving face near the waterline. In experiments without crevassing, the decrease in horizontal velocity above this maximum culminates in a zone of longitudinal compression at the surface somewhat Up-glacier from the face. This zone of compression appears to be a consequence of the non-linear rheology of ice. It disappears when a linear rheology is assumed. Explorations of the near-surface stress field indicate that when pervasive crevassing of the surface ice is accounted for in the simulations (by rheological softening), the zone of compressive strain rates does not develop. Variations in the pattern of horizontal velocity with glacier thickness support the contention that calving rates should increase with water depth at the calving face. In addition, the height of the subaerial calving face may have an importance that is not visible ill Current field data owing to the lack of variation in height of such faces in nature. Glaciers with lower calving faces may not have sufficient tensile stress to calve actively, while tensile stresses in simulated higher faces are sufficiently high that such faces will be unlikely to build in nature.

Sliding of Ice Past an Obstacle at Engabreen, Norway

At Engabreen, Norway, an instrumented panel containing a decimetric obstacle was mounted flush With the bed surface beneath 210 m of ice. Simultaneous measurements of normal and shear stresses, ice velocity and temperature were obtained as dirty basal ice flowed past the obstacle. Our measurements were broadly consistent with ice thickness, flow conditions and bedrock topography near the site of the experiment. Ice speed 0.45 m above the bed was about 130 mm d(-1), much less than the surface velocity of 800 mm d(-1) Average normal stress on the panel was 1.0-1.6 MPa, smaller than the expected ice overburden pressure. Normal stress was larger and temperature was lower on the stoss side than on the lee side, in accord with flow dynamics and equilibrium thermodynamics. Annual differences in normal stresses were correlated with changes in sliding speed and ice-flow direction. These temporal variations may have been caused by changes in ice rheology associated with changes in sediment concentration, water content or both. Temperature and normal stress were generally correlated, except when clasts presumably collided with the panel. Temperature gradients in the obstacle indicated that regelation was negligible, consistent with the obstacle size. Melt rate was about 10% of the sliding speed. Despite high sliding speed, no significant ice/bed separation was observed in the lee of the obstacle. Frictional forces between sediment particles in the ice and the panel, estimated from Hallet's (1981) model, indicated that friction accounted for about 5% of the measured bed-parallel force. This value is uncertain, as friction theories are largely untested. Some of these findings agree with sliding theories, others do not.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Ultimate and failure strengths of normal muds and serpentinite muds of ODP Leg 125 samples

Serpentinite seamounts in the Mariana forearc have been explained as diapirs rising from the Benioff zone. This hypothesis predicts that the serpentinites should have low strengths as well as low densities relative to the surrounding rocks. Drilling during Leg 125 showed that the materials forming Conical Seamount in the Mariana forearc and Torishima Forearc Seamount in the Izu-Bonin forearc are water-charged serpentinite muds of density <2 g/cm**3. Wykeham-Farrance torsion-vane tests showed that they are plastic solids with a rheology that bears many similarities to the idealized Cam clay soil model and is well described by critical-state soil mechanics. The serpentinite muds have ultimate strengths of 1.3 to 273.7 kPa and yield strengths of approximately 1.0 to 50 kPa. These muds thus are orders of magnitude weaker than salt and are, in fact, comparable in density and strength to common deep-sea clay muds. Such weak and low-density materials easily become diapiric. Serpentinite muds from the summit of Conical Seamount are weaker and more ductile than those on its flanks or on Torishima Forearc Seamount. Moreover, the summit muds do not contain the pronounced pure- and simple-shear fabrics that characterize those on the seamount flanks. The seamounts are morphologically similar to shield volcanoes, and anastomosing serpentinite debris flows descending from their summits are similar in map view to pahoehoe flows. These morphologic features, together with the physical properties of the muds and their similarities to other oceanic muds and the geochemistry of the entrained waters, suggest that many forearc serpentinite seamounts are gigantic (10-20 km wide, 1.5-2.0 km high) mud volcanoes that formed by the eruption of highly liquid serpentinite muds. Torishima Forearc Seamount, which is blanketed by more ìnormalî pelagic/volcaniclastic sediment, has probably been inactive since the Miocene. Conical Seamount, which seems to consist entirely of serpentinite mud and is venting fresh water of unusual chemistry from its summit, is presently active. Muds from the flanks of Conical Seamount are stronger and more brittle than those from the summit site, and muds from Torishima Forearc Seamount are stronger yet; this suggests that the serpentinite debris flows are compacted and dewatered as they mature. The shear fabrics probably result from downslope creep and flow, but may also be inherited.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Two-way traveltimes of internal layers and ice thickness between Dome C - Vostok - Dome A

Chinese scientists will start to drill a deep ice core at Kunlun station near Dome A in the near future. Recent work has predicted that Dome A is a location where ice older than 1 million years can be found. We model flow, temperature and the age of the ice by applying a three-dimensional, thermomechanically coupled full-Stokes model to a 70 × 70 km**2 domain around Kunlun station, using isotropic non-linear rheology and different prescribed anisotropic ice fabrics that vary the evolution from isotropic to single maximum at 1/3 or 2/3 depths. The variation in fabric is about as important as the uncertainties in geothermal heat flux in determining the vertical advection which in consequence controls both the basal temperature and the age profile. We find strongly variable basal ages across the domain since the ice varies greatly in thickness, and any basal melting effectively removes very old ice in the deepest parts of the subglacial valleys. Comparison with dated radar isochrones in the upper one third of the ice sheet cannot sufficiently constrain the age of the deeper ice, with uncertainties as large as 500 000 years in the basal age. We also assess basal age and thermal state sensitivities to geothermal heat flux and surface conditions. Despite expectations of modest changes in surface height over a glacial cycle at Dome A, even small variations in the evolution of surface conditions cause large variation in basal conditions, which is consistent with basal accretion features seen in radar surveys.

An investigation of the performance of a new Mechanical thrombectomy device using Bond Graph modelling: application to the extraction of blood clots in the middle cerebral artery

A number of thrombectomy devices using a variety of methods have now been developed to facilitate clot removal. We present research involving one such experimental device recently developed in the UK, called a ‘GP’ Thrombus Aspiration Device (GPTAD). This device has the potential to bring about the extraction of a thrombus. Although the device is at a relatively early stage of development, the results look encouraging. In this work, we present an analysis and modeling of the GPTAD by means of the bond graph technique; it seems to be a highly effective method of simulating the device under a variety of conditions. Such modeling is useful in optimizing the GPTAD and predicting the result of clot extraction. The aim of this simulation model is to obtain the minimum pressure necessary to extract the clot and to verify that both the pressure and the time required to complete the clot extraction are realistic for use in clinical situations, and are consistent with any experimentally obtained data. We therefore consider aspects of rheology and mechanics in our modeling.