We observe for the first time in silico the transition to a linear regime in the primary damage production in tungsten. As the critical plasma-facing material in fusion reactors, radiation damage in tungsten has been studied extensively in experiments and simulations. Irradiation experiments routinely produce recoils in the MeV range while full atomistic modelling has been limited to a few hundred keV. Here we bridge these scales with extremely large-scale and accurate machine-learning-driven molecular dynamics simulations with recoil energies up to 2 MeV in systems up to one billion atoms. We reveal four regimes of primary damage as a function of damage energy, with a transition to a high-energy regime that deviates from all previous models. Curiously, the start of the high-energy regime coincides with the highest possible recoil energy to tungsten atoms from fusion-emitted neutrons (300 keV).

Thwaites Glacier, a large outlet glacier of the West Antarctic Ice Sheet, holds over a half meter of sea level rise equivalent. The large potential contribution to sea level is concerning given that the glacier may be vulnerable to self-sustaining processes of rapid retreat due to the retrograde bed slope that characterises much of the glacier’s bed. Such a reverse-sloping bed exists behind the relatively high ridge on which the western calving front (WCF) of the Thwaites Glacier terminus currently rests. Our study focuses on the factors that control the calving dynamics of the WCF and the ability of mélange to influence these dynamics. Employing the 3D Helsinki Discrete Element Model (HiDEM), we find that calving at this location currently occurs as rifts form and widen due to longitudinal tensile stresses associated with ice flow across the grounding line. Calving is restricted in HiDEM simulations that include a constricted mélange field that is confined within the bounds of the model domain. A thicker, constricted mélange field fully suppresses calving. These simulations show the development of robust force chains that transmit resistive forces to the Thwaites WCF. In the future, the ability for mélange to influence the calving dynamics at the WCF will depend on the degree to which it is constrained in the wide Amundsen Sea Embayment, either through binding in land-fast sea ice or jamming behind large, grounded icebergs. As such, sea-ice conditions and iceberg characteristics will need to be considered along with the presence of mélange in investigations of the future retreat of the prominently recognised Thwaites Glacier.

Radiation therapy is one of the most common cancer treatments, and dose optimization and targeting of radiation are crucial since both cancerous and healthy cells are affected. Different mathematical and computational approaches have been developed for this task. The most common mathematical approach, dating back to the late 1970's, is the linear-quadratic (LQ) model for the survival probability given the radiation dose. Most simulation models consider tissue as a continuum rather than consisting of discrete cells. While reasonable for large-scale models (e.g., human organs), continuum approaches necessarily neglect cellular-scale effects, which may play a role in growth, morphology, and metastasis of tumors. Here, we propose a method for modeling the effect of radiation on cells based on the mechanobiological CellSim3D simulation model for growth, division, and proliferation of cells. To model the effect of a radiation beam, we incorporate a Monte Carlo procedure into CellSim3D with the LQ model by introducing a survival probability at each beam delivery. Effective removal of dead cells by phagocytosis was also implemented. Systems with two types of cells were simulated: stiff slowly proliferating healthy cells and soft rapidly proliferating cancer cells. For model verification, the results were compared to prostate cancer (PC-3 cell line) data for different doses and we found good agreement. In addition, we simulated proliferating systems and analyzed the probability density of the contact forces. We determined the state of the system with respect to the jamming transition and found very good agreement with experiments.

Marine-terminating glaciers and ice shelves are notoriously complex, with a wide range of ice-dynamic and calving processes occuring in response to oceanographic, atmospheric and glaciological influences. Within this complexity, however, we can recognise order on at least two scales. First, marine ice fronts typically form vertical cliffs, reflecting competition between oversteepening (ice flow and melt-undercutting) and failure. Calving magnitude-frequency distributions have power-law form with an exponent of -1.2, characteristic of self-organising criticality (SOC). Such systems have a critical point as an attractor, such that the system converges on the failure threshold. The second scale is that of the whole ice tongue. Tidewater glaciers and ice shelves typically oscillate around stable positions for multiple years, punctuated by transitions to new quasi-stable positions. Stability is encouraged by pinning points which function as attractors at thresholds between stable and metastable states. Ice tongues may exist in metastable states for variable amounts of time, from days to decades. Factors encouraging rapid relaxation to the threshold include large stress gradients and rapid basal melt, and factors encouraging long relaxation times include low stress gradients, low melt rates, and buttressing from mélange or sea ice. Calving magnitude-frequency distributions have exponential form, reflecting the stochastic nature of calving in the metastable zone. Both scales of self-organisation emerge spontaneously from physically-based calving models such as the Helsinki Discrete Element Model (HiDEM) and the crevasse-depth (CD) calving law implemented in Elmer/Ice. Purely deterministic models, however, are not optimal for long-term simulations, especially in Antarctic contexts. We present results of preliminary simulations using a stochastic CD calving law, which opens up the possibility of a universal calving model applicable to both the Greenland and Antarctic ice sheets.

Creating accurate predictive models for drift and pack ice is crucial for a wide array of applications, from improving maritime operations to improving weather prediction and climate simulations. Traditional large-scale sea ice dynamics models rely on phenomenological ice rheology to simulate ice movements. These models are efficient on large scales but struggle to depict smaller-scale ice features. In our study, we use a new version of the HiDEM discrete element model software to examine the formation of drift and pack ice under various stress conditions. Our findings show that high-resolution size distributions of ice floes are universal and multimodal, and that compression ridges form three distinct zones. Reproducing complex characteristics of this nature in a standard rheology model is challenging, suggesting that a combination of models may be necessary for more precise predictions of sea ice dynamics. We propose a potential hybrid algorithm that integrates these approaches. Sea ice forms in cold climates and is susceptible to being easily fragmented by wind and currents, resulting in a dynamic landscape comprising solid fast ice, drift ice and pack ice. Pack ice, in particular, can pose challenges such as hindering shipping, causing damage to offshore structures, and complicating traditional fishing and hunting activities. Operational models for sea ice dynamics are currently utilized to optimize ship routes and the deployment of icebreakers. Although existing rheology-based models perform well on large scales, they encounter difficulties in capturing the finer details that are often crucial. In this study, we utilize a high-resolution Discrete Element Model computer code that is capable of simulating detailed sea ice dynamics at scales ranging from meters to kilometers. Our simulation results reveal insights that are not readily obtained from conventional large-scale models, and we explore the potential for integrating these two approaches to create a hybrid model. Ultra high-definition simulation (0.5 m elements) of sea ice fragmentation on a square kilometer scale The HiDEM model captures in fine detail the formations of leads, pressure ridge networks, and floe-size distributions The model reveals features that cannot be reproduced by rheological models, suggesting a hybrid method for prediction

The western region of the wide Thwaites Glacier terminus is characterized by a near-vertical calving front. The grounding line at this western calving front (WCF) rests on a relatively high ridge, behind which exists a reverse-sloping bed; retreat of the grounding line into this over-deepening basin could therefore expose deep calving faces that may be subject to ice-cliff failure. Here, we use the 3D Helsinki Discrete Element Model to identify the factors that control the calving dynamics in this location. We then focus on the ability of m & eacute;lange to influence these dynamics given the wide embayment in which Thwaites Glacier terminates. We find that calving along the WCF is currently influenced by ice flow across the grounding line and consequent longitudinal tensile stress and rift formation. Calving is slowed in simulations that are initiated with a highly constricted m & eacute;lange, with a thicker m & eacute;lange suppressing calving entirely. We liken the constrained simulations to a scenario in which m & eacute;lange piles behind a large grounded iceberg. In a future which may see calving become a more dominant control on the retreat of Thwaites Glacier, this type of blockage will be necessary for robust force chains to develop and transmit resistive forces to the terminus. The ability of the m & eacute;lange to hinder calving at this location will be determined by the presence and rigidity of binding land-fast sea ice and iceberg keel depths. Therefore, it is necessary to represent calving, m & eacute;lange and sea ice in a single framework to predict the fate of Thwaites Glacier. Thwaites Glacier, a major outlet glacier of the West Antarctic Ice Sheet, may be vulnerable to rapid retreat. The terminus of Thwaites Glacier is over 100 km wide; in this paper, we focus on the western terminus region, which meets the Amundsen Sea at a near-vertical ice face. We use a 3D glacier model that simulates ice fracture to identify the factors controlling iceberg calving. We also run simulations to explore the potential for m & eacute;lange, a conglomerate of variably sized ice pieces in front of a glacier, to impact calving dynamics at this location. Our simulations show that iceberg calving currently occurs through the growth of short ice-shelf extensions and the loss of resisting forces as the ice begins to float. Our simulations also demonstrate that calving is slowed when a thick immobile m & eacute;lange is present. The ability of the m & eacute;lange to apply a backforce to glaciers has been previously identified as an important brake on such calving-driven retreat. In a future that may see calving become a more dominant control on the retreat of Thwaites Glacier, this type of blockage will be necessary for the transmission of resistive forces to the terminus. Calving at Thwaites Glacier's western terminus is controlled by longitudinal tensile stresses and fracturing as basal traction is lost 3D modeling demonstrates that m & eacute;lange, if constrained, can inhibit calving at this location Future sea-ice conditions and iceberg morphology will determine if m & eacute;lange can compact and transmit substantial resistive forces

The largest uncertainty in future sea-level rise is loss of ice from the Greenland and Antarctic Ice Sheets. Ice shelves, freely floating platforms of ice that fringe the ice sheets, play a crucial role in restraining discharge of grounded ice into the ocean through buttressing. However, since the 1990s, several ice shelves have thinned, retreated, and collapsed. If this pattern continues, it could expose thick cliffs that become structurally unstable and collapse in a process called marine ice cliff instability (MICI). However, the feedbacks between calving, retreat, and other forcings are not well understood. Here we review observed modes of calving from ice shelves and marine-terminating glaciers, and their relation to environmental forces. We show that the primary driver of calving is long-term internal glaciological stress, but as ice shelves thin they may become more vulnerable to environmental forcing. This vulnerability—and the potential for MICI—comes from a combination of the distribution of preexisting flaws within the ice and regions where the stress is large enough to initiate fracture. Although significant progress has been made modeling these processes, theories must now be tested against a wide range of environmental and glaciological conditions in both modern and paleo conditions. ▪Ice shelves, floating platforms of ice fed by ice sheets, shed mass in a near-instantaneous fashion through iceberg calving.▪Most ice shelves exhibit a stable cycle of calving front advance and retreat that is insensitive to small changes in environmental conditions.▪Some ice shelves have retreated or collapsed completely, and in the future this could expose thick cliffs that could become structurally unstable called ice cliff instability.▪The potential for ice shelf and ice cliff instability is controlled by the presence and evolution of flaws or fractures within the ice.

Forecasts of sea-ice motion and fragmentation are of vital importance for all human interactions with sea ice, ranging from those involving indigenous hunters to shipping in polar regions. Sea-ice models are also important for simulating long-term changes in a warming climate. Here, we apply the Helsinki Discrete Element Model (HiDEM), originally developed for glacier calving, to sea-ice breakup and dynamics. The code is highly optimized to utilize high-end supercomputers to achieve an extreme time and space resolution. Simulated fracture patterns and ice motion are compared with satellite images of the Kvarken region of the Baltic Sea from March 2018. A second application of HiDEM involves ice ridge formation in the Gulf of Riga. With a few tens of graphics processing units (GPUs), the code is capable of reproducing observed ice patterns that in nature may take a few days to form; this is done over an area of ∼100km×100km, with an 8 m resolution, in computations lasting ∼10 h. The simulations largely reproduce observed fracture patterns, ice motion, fast-ice regions, floe size distributions, and ridge patterns. The similarities and differences between observed and computed ice dynamics and their relation to initial conditions, boundary conditions, and applied driving forces are discussed in detail. The results reported here indicate that the HiDEM has the potential to be developed into a detailed high-resolution model for sea-ice dynamics at short timescales, which, when combined with large-scale and long-term continuum models, may form an efficient framework for forecasts of sea-ice dynamics.

Tissue growth kinetics and interface dynamics depend on the properties of the tissue environment and cell-cell interactions. In cellular environments, substrate heterogeneity and geometry arise from a variety factors, such as the structure of the extracellular matrix and nutrient concentration. We used the CellSim3D model, a kinetic cell division simulator, to investigate the growth kinetics and interface roughness dynamics of epithelial tissue growth on heterogeneous substrates with varying topologies. The results show that the presence of quenched disorder has a clear effect on the colony morphology and the roughness scaling of the interface in the moving interface regime. In a medium with quenched disorder, the tissue interface has a smaller interface roughness exponent, α, and a larger growth exponent, β. The scaling exponents also depend on the topology of the substrate and cannot be categorized by well-known universality classes.

The universality of interfacial roughness in growing epithelial tissue has remained a controversial issue. Kardar-Parisi-Zhang (KPZ) and molecular beam epitaxy (MBE) universality classes have been reported among other behaviors including a total lack of universality. Here, we simulate tissues using the cellsim3d kinetic division model for deformable cells to investigate cell-colony scaling. With seemingly minor model changes, it can reproduce both KPZ- and MBE-like scaling in configurations that mimic the respective experiments. Tissue growth with strong cell-cell adhesion in a linear geometry is KPZ like, while weakly adhesive tissues in a radial geometry are MBE like. This result neutralizes the apparent scaling controversy.

A major challenge within material science is the proper modeling of force transmission through fragmenting materials under compression. A particularly demanding material is sea ice, which on small scales is an anisotropic material with quasibrittle characteristics under failure. Here we use the particle-based model HiDEM and laboratory-scale experiments on saline ice to develop a material model for fragmenting ice. The material behavior of the HiDEM model-ice, and the experiments are compatible on force transmission and fragmentation if: (i) the typical HiDEM glacier-scale particle size of meters is brought down to millimeters corresponding to the grain size of the laboratory ice, (ii) the often used HiDEM lattice structure is replaced by a planar random structure with an anisotropy in the direction normal to the randomized plane, and (iii) the instant tensile and bending failure criterion, used in HiDEM on glacier scale, is replaced by a cohesive softening failure potential for energy dissipation. The main outcomes of this exercise is that many of the, more or less, traditional ice modeling schemes are proven to be incomplete. In particular, local crushing of ice is not valid as a generic failure mode for fragmented ice under compression. Rather, shear failure, as described by Mohr-Coulomb theory is demonstrated to be the dominant failure mode.

<p>Below ice shelves, complex interactions between the ice and the ocean are at stake that have large implications for future sea level rise. Basal melting from the ocean is recognised to have large impacts on the stability. Many studies focus on theses interactions in coupled models at different spatio-temporal scales. However, most of them consider the basal topography of the shelf as smooth ignoring its irregular state due to basal crevassing or channel-like features. We propose to investigate the impact of these features on basal melt and ice shelf stability by using a discrete particule model and an ocean model applied at the ice shelf of Thwaites glacier.</p>

Universality of interfacial roughness in growing epithelial tissue has remained a controversial issue. Kardar-Parisi-Zhang (KPZ) and Molecular Beam Epitaxy (MBE) universality classes have been reported among other behaviors including total lack of universality. Here, we utilize a kinetic division model for deformable cells to investigate cell-colony scaling. With seemingly minor model changes, it can reproduce both KPZ- and MBE-like scaling in configurations that mimic the respective experiments. This result neutralizes the apparent scaling controversy. It can be speculated that this diversity in growth behavior is beneficial for efficient evolution and versatile growth dynamics.

Ice shelves play a key role in the dynamics of marine ice sheets by buttressing grounded ice and limiting rates of ice flux to the oceans. In response to recent climatic and oceanic change, ice shelves fringing the West Antarctic Ice Sheet (WAIS) have begun to fragment and retreat, with major implications for ice-sheet stability. Here, we focus on the Thwaites Eastern Ice Shelf (TEIS), the remaining pinned floating extension of Thwaites Glacier. We show that TEIS has undergone a process of fragmentation in the last 5 years, including brittle failure along a major shear zone, formation of tensile cracks on the main body of the shelf, and a release of tabular bergs on both the eastern and western flanks. Simulations with the Helsinki Discrete Element Model (HiDEM) show that this pattern of failure is associated with high backstress from a submarine pinning point at the distal edge of the shelf. We show that a significant zone of shear, upstream of the main pinning point, developed in response to the rapid acceleration of the shelf between 2002 and 2006, seeding damage on the shelf. Subsequently, basal melting and positive feedback between damage and strain rates weakened TEIS, allowing damage to accumulate. Thus, although backstress on TEIS has likely diminished over time as the pinning point shrunk, accumulation of damage has ensured that the ice in the shear zone remained the weakest link in the system. Experiments with the BISICLES ice-sheet model indicate that additional damage to or unpinning of TEIS is unlikely to trigger significantly increased ice loss from WAIS, but the calving response to the loss of TEIS remains highly uncertain. It is widely recognised that ice-shelf fragmentation and collapse can be triggered by hydrofracturing and/or unpinning from ice-shelf margins or grounding points. Our results indicate a third mechanism, backstress triggered failure, that can occur if and when an ice shelf is no longer able to withstand stress imposed by pinning points. In most circumstances, pinning points are essential for ice-shelf stability, but as ice shelves thin and weaken, the concentration of backstress in damaged ice upstream of a pinning point may provide the seeds of their demise.

Abstract Iceberg calving strongly controls glacier mass loss, but the fracture processes leading to iceberg formation are poorly understood due to the stochastic nature of calving. The size distributions of icebergs produced during the calving process can yield information on the processes driving calving and also affect the timing, magnitude, and spatial distribution of ocean fresh water fluxes near glaciers and ice sheets. In this study, we apply fragmentation theory to describe key calving behaviours, based on observational and modelling data from Greenland and Antarctica. In both regions, iceberg calving is dominated by elastic-brittle fracture processes, where distributions contain both exponential and power law components describing large-scale uncorrelated fracture and correlated branching fracture, respectively. Other size distributions can also be observed. For Antarctic icebergs, distributions change from elastic-brittle type during ‘stable’ calving to one dominated by grinding or crushing during ice shelf disintegration events. In Greenland, we find that iceberg fragment size distributions evolve from an initial elastic-brittle type distribution near the calving front, into a steeper grinding/crushing-type power law along-fjord. These results provide an entirely new framework for understanding controls on iceberg calving and how calving may react to climate forcing.

Marine ice-cliff instability could accelerate ice loss from Antarctica, and according to some model predictions could potentially contribute >1 m of global mean sea level rise by 2100 at current emission rates. Regions with over-deepening basins >1 km in depth (e.g., the West Antarctic Ice Sheet) are particularly susceptible to this instability, as retreat could expose increasingly tall cliffs that could exceed ice stability thresholds. Here, we use a suite of high-fidelity glacier models to improve understanding of the modes through which ice cliffs can structurally fail and derive a conservative ice-cliff failure retreat rate parameterization for ice-sheet models. Our results highlight the respective roles of viscous deformation, shear-band formation, and brittle-tensile failure within marine ice-cliff instability. Calving rates increase non-linearly with cliff height, but runaway ice-cliff retreat can be inhibited by viscous flow and back force from iceberg mélange.

Projections of future ice sheet mass loss and thus sea level rise rely on the parametrization of iceberg calving in ice sheet models. The interconnection between submarine melt-induced undercutting and calving is still poorly understood, which makes predicted contributions of tidewater glaciers to sea level rise uncertain. Here, we compare detailed 3-D simulations of fracture initiation obtained with the Helsinki Discrete Element Model (HiDEM) to observations, prior to a major calving event at Bowdoin Glacier, Northwest Greenland. Observations of a plume surfacing at the calving location suggest that local melt-undercutting influenced the size of the major calving event. Therefore, several experiments are conducted with various local and distributed (front-wide) undercut geometries. Although the number of undercut experiments is limited by computational requirements, one of the conjectured undercut geometries reproduces the crevasse leading to the observed major calving event in great detail. Our simulations show that undercutting leads to initiation of wider fractures more than 100 m upstream of the terminus, well-beyond the directly undercut region. When combining a moderate distributed undercut with local amplified undercuts at the two observed plumes, fracture initiation also increases in between the local undercuts. Thus, our results agree with previous studies suggesting the existence of a “calving amplifier” effect by submarine melt, both upglacier and across-glacier. Consequently, the simulations show the potentially large impact of submarine melt-induced undercutting on iceberg size.

Rapid grounding line retreat at marine-terminating glaciers could expose ice cliffs with heights greater than those on observational record. However, the finite strength of ice places a limit on the height of subareal cliffs. It is proposed that marine ice-cliff instability (MICI) will begin once a stable height threshold is exceeded. If a glacier is situated over a retrograde slope, as is the case for Thwaites Glacier and much of the West Antarctic Ice Sheet, MICI can be expected to accelerate as retreat progresses and increasingly tall and unstable ice cliffs are formed. This is consequential for global sea level rise, yet large uncertainties remain in the prediction of MICI retreat rates. We investigate MICI by pairing the full Stokes continuum model Elmer/Ice and the Helsinki Discrete Element Model (HiDEM). Viscous flow, simulated in Elmer/Ice, is found to be a necessary pre-condition for MICI collapse. Forward advance and bulging lead to ice-front instability and pervasive crevassing in HiDEM. This culminates in full-thickness calving events. We do not observe calving at ice faces prior to viscous deformation. HiDEM simulations that implement viscous flow (HiDEM-ve) also show forward advance and waterline bulging, similar to the Elmer/Ice simulations. However, the importance of granular shear is highlighted by pronounced shear bands and patterns of surface lowering in HiDEM-ve output. These results emphasize the importance and complexity of viscous and brittle process interaction during MICI. A simulation matrix of grounded termini shows that calving frequency and magnitude increase with the thickness of the calving front. The time required for viscous flow to recreate unstable conditions is influenced by thickness as well as ice temperature and basal friction. Simulations of buoyant termini are seen to calve through basal-crevassing and block-rotation, as opposed to incising surface-crevasses. Lastly, we observe that buttressing mélange can suppress retreat rate if a sufficient resistive force is delivered to the calving front. A physically-based law for MICI retreat rate is derived from our simulation matrix; this calving rate law can be incorporated into large-scale ice sheet models to constrain projections of Antarctic retreat and associated global sea level rise. Our results will also be used to investigate the future retreat of Thwaites Glacier, which is vulnerable to MICI due to a retreating grounding line, fragile floating ice shelf, and precarious positioning above an overdeepening basin.