Melting Threshold Research Articles

Enhanced two-temperature model and its application in comprehensive analysis of femtosecond laser melting of gold, copper and their alloys.

Extensive research on ultrashort laser-induced melting of noble metals like Au, Ag and Cu is available. However, studies on laser energy deposition and thermal damage of their alloys, which are currently attracting interest for energy harvesting and storage devices, are limited. This study investigates the melting damage threshold (DT) of three intermetallic alloys of Au and Cu (Au3Cu, AuCu and AuCu3) subjected to single-pulse femtosecond laser irradiation, comparing them with their constituent metals. This is accomplished by extending an earlier-developed two-temperature model (TTM)-based code with several improvements, including precise modeling of temperature-dependent optical properties and ballistic electron transport. The dynamic optical model inclusive of ballistic effects is demonstrated to reproduce the experimental DT of pure metals with minimal variation and is therefore adopted for further investigation. Our simulations reveal that the alloy films have significantly lower incipient and complete melting thresholds compared to the pure metals due to their low thermal conductivity and high electron-phonon coupling strength. Theoretical studies on varying the thickness of metal and alloy films unveil the usual trend of a rapid increase in DT up to a certain thickness, followed by a saturation region. This universal DT profile is elucidated by proposing a first-of-its-kind analytical function. Excellent agreement between the coefficients of the function with optical and electron diffusion parameters derived from the comprehensive theory proposed here reinforces the robustness of the model. The novelty of this study also lies in introducing the concept of a critical film thickness for which the entire film attains the melting temperature at its complete melting threshold.

Radiogenic heat production provides a thermal threshold for Archean cratonization process

A striking feature of the Yilgarn craton at the current erosional level is an abundance of late K-rich granites with radiogenic heat production elevated far above global crustal averages. Extrapolated back in time, the total thickness and contribution to crustal heat production and heat flow from these granites were greater, implying that the deeper crustal sources must also have had elevated radiogenic heat production. Through back-calculated and time-integrated one-dimensional thermal modeling underpinned by geological and geochemical constraints for the model crustal columns, we find that elevated radiogenic heat production provided a significant internal driver for prolonged crustal melting and eventual cratonization of the Yilgarn craton. Our results show that elevated thermal gradients driven by high heat production thermally primed the mid- and deep crust at or above the threshold for large-volume partial melting over long periods of time, as evidenced in the magmatic rock record. This would have been amplified by any additional heat that may have been provided by the mantle melting processes that punctuated the geological history. Over time, advective movement of progressively more radiogenic heat production to the shallower crust would have resulted in two complementary outcomes: progressively refractory deep crust and long-term cooling. The widespread granite “bloom” at 2650−2600 Ma records the final time at which the crust was fertile enough to melt in large volumes and the thermal gradient was hot enough to intersect the solidus. The magnitude of radiogenic heat production in the Yilgarn craton has been underestimated in previous studies, resulting in an underappreciation of the importance of its contribution to internal drivers of magmatism and ultimately cratonization.

ATP-dependent thermoring basis for the heat unfolding of the first nucleotide-binding domain isolated from human CFTR.

Traditionally, the thermostability of a protein is defined by a melting temperature, at which half of the protein is unfolded. However, this definition cannot indicate the structural origin of a heat-induced unfolding pathway. Here, the thermoring structures were studied on the ATP-dependent heat-induced unfolding of the first nucleotide-binding domain from the human cystic fibrosis transmembrane conductance regulator. The results showed that initial theoretical and experimental melting thresholds aligned well after three structural perturbations including the F508del mutation, the most common cause of cystic fibrosis. This alignment further demonstrated that the heat-induced unfolding process began with the disruption of the least-stable noncovalent interaction within the biggest thermoring along the single peptide chain. The C-terminal region, which was related to the least-stable noncovalent interaction and the ATP-dependent dimerization of two nucleotide-binding domains, emerged as a crucial determinant of the thermal stability of the isolated protein and a potential interfacial drug target to alleviate the thermal defect caused by the F508del mutation. This groundbreaking discovery significantly advances our understanding of protein activity, thermal stability, and molecular pathology.

Femtosecond Laser Ablation and Delamination of Functional Magnetic Multilayers at the Nanoscale.

Laser nanostructuring of thin films with ultrashort laser pulses is widely used for nanofabrication across various fields. A crucial parameter for optimizing and understanding the processes underlying laser processing is the absorbed laser fluence, which is essential for all damage phenomena such as melting, ablation, spallation, and delamination. While threshold fluences have been extensively studied for single compound thin films, advancements in ultrafast acoustics, magneto-acoustics, and acousto-magneto-plasmonics necessitate understanding the laser nanofabrication processes for functional multilayer films. In this work, we investigated the thickness dependence of ablation and delamination thresholds in Ni/Au bilayers by varying the thickness of the Ni layer. The results were compared with experimental data on Ni thin films. Additionally, we performed femtosecond time-resolved pump-probe measurements of transient reflectivity in Ni to determine the heat penetration depth and evaluate the melting threshold. Delamination thresholds for Ni films were found to exceed the surface melting threshold suggesting the thermal mechanism in a liquid phase. Damage thresholds for Ni/Au bilayers were found to be significantly lower than those for Ni and fingerprint the non-thermal mechanism without Ni melting, which we attribute to the much weaker mechanical adhesion at the Au/glass interface. This finding suggests the potential of femtosecond laser delamination for nondestructive, energy-efficient nanostructuring, enabling the creation of high-quality acoustic resonators and other functional nanostructures for applications in nanosciences.

Enhancing superconductivity in CoSi2 films with laser annealing

Laser annealing was employed to trigger the solid-state reaction of a thin Co film (2.5 nm) with undoped Si. A metastable disilicide layer was obtained after one laser pulse close to the melt threshold. Its diffraction pattern, relaxed lattice parameter, and residual resisitivity are consistent with the formation of the defective CsCl structure. The CoSi2 phase was found after prolonging the thermal treatment with additional pulses or rapid thermal annealing. Because CoSi is skipped in the phase sequence, CoSi2 layers are more uniform in thickness, have an increased superconductivity and a reduced formation temperature. This approach is compatible with the SALICIDE process and can be used to form smooth contacts in superconducting or regular transistors.

Atomistic simulations of oblique collisions between crystalline and amorphous SiC nanoparticles: Mechanisms of deformation and scaling coefficients of restitutions

Systematic atomistic simulations of oblique collisions between spherical crystalline, 6H and 3C, and amorphous SiC nanoparticles (NPs) are performed in the bouncing and fragmentation regimes for the particles with diameters from 8 nm to 30 nm, impact velocities from 100 ms−1 to 5000 ms−1, and in the whole range of the geometric impact parameter. The critical sticking velocity of 6H-SiC NPs reduces from ∼490 ms−1 to ∼150 ms−1 when the NP diameter increases from 10 nm to 30 nm. Below the melting threshold, the collision of crystalline NPs induces the formation of localized zones of densified material, while the remaining parts of NPs are not affected by plastic deformations. The impact of amorphous NPs results in the plastic deformation of significant parts of the NP material and much larger irreversible losses of mechanical energy. The post-collision translational and rotational velocities of crystalline NPs as well as coefficients of restitution (CORs) demonstrate weak dependence on the NP size if the NP diameter is greater than 20 nm. At constant impact velocity, the normal COR only marginally varies for head-on and oblique collisions, when the collision angle is smaller than ∼53°. The tangential center-of-mass and point-of-contact CORs, as functions of the impact parameter, demonstrate extremal behavior with the minima at collision angles from ∼27° to ∼37°. The normal COR exhibits a linear decay as a function of the impact velocity where the slope coefficient is different below and above the melting threshold. The tangential CORs attain the minimum values at the impact velocity of ∼300 ms−1. The obtained results are summarized in the form of fitting equations for CORs as functions of the impact parameter and velocity, which can be readily used to predict the post-collision translational and angular velocities of NPs in large-scale simulations of granular gas flows.

High-density nanodot structures on silicon solar cell surfaces irradiated by ultraviolet laser pulses below the melting threshold fluence

The surface morphology of silicon solar cells irradiated with KrF excimer laser pulses (λ = 248 nm, τ = 20 ns) was investigated below the experimentally observed melting threshold fluence (F th) of 0.47 J cm−2 (±20%). At laser fluences of 0.23–0.48 J cm−2 (equivalent to 0.49F th to 1.0F th), nanodot structures with a height and width of approximately 60–120 nm were periodically formed with an interdot spacing similar to the laser wavelength. The observed nanodot density (29 dots μm2) was higher than that previously obtained at longer wavelengths. Furthermore, crystallinity analysis by micro-Raman spectroscopy revealed a Raman shift of 519.56 cm−1 after irradiation (N= 1500 pulses), compared with 518.27 cm−1 prior to irradiation. A laser fluence of 0.41 J cm−2 ( = 0.87F th) was found to induce compressive stress on the silicon solar cell surface.

Laser‐Annealed SiO2/Si1−xGex Scaffolds for Nanoscaled Devices, Synergy of Experiment, and Computation

Ultraviolet nanosecond laser annealing (UV‐NLA) proves to be an important technique, particularly when tightly controlled heating and melting are necessary. In the realm of semiconductor technologies, the significance of NLA grows in tandem with the escalating intricacy of integration schemes in nanoscaled devices. Silicon–germanium alloys are studied for decades for their compatibility with silicon devices. Indeed, they enable the manipulation of properties like strain, carrier mobilities, and bandgap. In this framework, they can for instance boost the performances of p‐type MOSFETs but also enable near infrared absorption and emission for applications in photodetection and photonics. Laser melting on such types of layers, however, results, up to now, in the development of extended defects and poor control over layer morphology and homogeneity. Herein, the laser melting of ≈700 nm‐thick relaxed silicon–germanium samples coated with SiO2 nanoarrays, observing the resulting material to maintain an unaltered lattice, is investigated. It is found that the geometrical parameters of the silicon oxide have an impact on the thermal budget samples, influencing melt threshold, melt depth, and germanium distribution.

Physical Regimes and Mechanisms of Picosecond Laser Fragmentation of Gold Nanoparticles in Water from X-ray Probing and Atomistic Simulations.

Laser fragmentation in liquids has emerged as a promising green chemistry technique for changing the size, shape, structure, and phase composition of colloidal nanoparticles, thus tuning their properties to the needs of practical applications. The advancement of this technique requires a solid understanding of the mechanisms of laser-nanoparticle interactions that lead to the fragmentation. While theoretical studies have made impressive practical and mechanistic predictions, their experimental validation is required. Hence, using the picosecond laser fragmentation of Au nanoparticles in water as a model system, the transient melting and fragmentation processes are investigated with a combination of time-resolved X-ray probing and atomistic simulations. The direct comparison of the diffraction profiles predicted in the simulations and measured in experiments has revealed a sequence of several nonequilibrium processes triggered by the laser irradiation. At low laser fluences, in the regime of nanoparticle melting and resolidification, the results provide evidence of a transient superheating of crystalline nanoparticles above the melting temperature. At fluences about three times the melting threshold, the fragmentation starts with evaporation of Au atoms and their condensation into small satellite nanoparticles. As fluence increases above five times the melting threshold, a transition to a rapid (explosive) phase decomposition of superheated nanoparticles into small liquid droplets and vapor phase atoms is observed. The transition to the phase explosion fragmentation regime is signified by prominent changes in the small-angle X-ray scattering profiles measured in experiments and calculated in simulations. The good match between the experimental and computational diffraction profiles gives credence to the physical picture of the cascade of thermal fragmentation regimes revealed in the simulations and demonstrates the high promise of the joint tightly integrated computational and experimental efforts.

Laser soldering of nickel plated steel sheets

The growing prominence of the electric vehicle industry, fueled by environmental concerns, has demanded innovation in various aspects of battery technologies with special emphasis on increasing the efficiency of both electric storage and its retrieval. An unexplored area of this is to identify the possibilities and limits of laser soldering. Here, we reveal the effects of surface pretreatment conditions and the amount of filler, along with the laser power and irradiation time on the characteristics of laser-soldered joints, by simultaneously evaluating the electrical and mechanical behavior of laser-soldered nickel-plated steel sheets (Hilumin®). By describing the morphological characteristics of the resolidified solder and the electrical and mechanical properties of the joints, we identify three, characteristically different morphological appearances and highlight the optimal one, where uniform and mostly void-free solder can be produced. Furthermore, we report a correlation between the threshold of upper sheet melting (either expressed as laser power or irradiation time) and joint deterioration in terms of the electrical and mechanical properties of the joint. We conclude that laser soldering can create joints with outstanding electrical conductance and adequate mechanical stability that meets the critical specifications of battery joining technologies when the surface pretreatment condition and processing parameters are properly optimized.

Atomistic-Continuum Study of an Ultrafast Melting Process Controlled by a Femtosecond Laser-Pulse Train.

In femtosecond laser fabrication, the laser-pulse train shows great promise in improving processing efficiency, quality, and precision. This research investigates the influence of pulse number, pulse interval, and pulse energy ratio on the lateral and longitudinal ultrafast melting process using an experiment and the molecular dynamics coupling two-temperature model (MD-TTM model), which incorporates temperature-dependent thermophysical parameters. The comparison of experimental and simulation results under single and double pulses proves the reliability of the MD-TTM model and indicates that as the pulse number increases, the melting threshold at the edge region of the laser spot decreases, resulting in a larger diameter of the melting region in the 2D lateral melting results. Using the same model, the lateral melting results of five pulses are simulated. Moreover, the longitudinal melting results are also predicted, and an increasing pulse number leads to a greater early-stage melting depth in the melting process. In the case of double femtosecond laser pulses, the pulse interval and pulse energy ratio also affect the early-stage melting depth, with the best enhancement observed with a 2 ps interval and a 3:7 energy ratio. However, pulse number, pulse energy ratio, and pulse interval do not affect the final melting depth with the same total energies. The findings mean that the phenomena of melting region can be flexibly manipulated through the laser-pulse train, which is expected to be applied to improve the structural precision and boundary quality.

Surface morphological evolution and microhardness change of Cr12MoV steel by pulsed laser polishing

Laser polishing is gaining more and more attention because of its environmental friendliness and ability to process complex surfaces. It is considered as one of the technologies to replace traditional methods in the mold industry. In this work, the variation pattern of surface roughness with single pulse energy density and light-out overlap rate during nanosecond pulsed laser polishing of Cr12MoV steel was investigated. The surface roughness increased and then decreased with increasing single-pulse energy density and beam overlap rate, respectively. The lowest surface roughness was found to be 0.50 μm, which is a decrease of 79.34 % compared to the untreated sample. The effect of each parameter on the surface morphology during the laser polishing process is explained. There is minimal austenitization on the surface of the sample after laser polishing. The second phase of the remelting layer disappeared. The microhardness of the remelting layer increased by 36.17 % because the grain size was reduced by 30.65 %. This work further explains the mechanism transformation threshold for pulsed laser polishing of steel. With 95 % beam overlap, the melting and gasification energy density thresholds for Cr12MoV steel are respectively 1.269 J/cm2 and 2.937 J/cm2. The beam overlap rate of 95.5 % is the critical condition for continuous melt movement. It provides a new idea for the mechanism judgment of pulsed laser polishing.

Modeling of Short-Pulse Laser Interactions with Monolithic and Porous Silicon Targets with an Atomistic-Continuum Approach.

The acquisition of reliable knowledge about the mechanism of short laser pulse interactions with semiconductor materials is an important step for high-tech technologies towards the development of new electronic devices, the functionalization of material surfaces with predesigned optical properties, and the manufacturing of nanorobots (such as nanoparticles) for bio-medical applications. The laser-induced nanostructuring of semiconductors, however, is a complex phenomenon with several interplaying processes occurring on a wide spatial and temporal scale. In this work, we apply the atomistic-continuum approach for modeling the interaction of an fs-laser pulse with a semiconductor target, using monolithic crystalline silicon (c-Si) and porous silicon (Si). This model addresses the kinetics of non-equilibrium laser-induced phase transitions with atomic resolution via molecular dynamics, whereas the effect of the laser-generated free carriers (electron-hole pairs) is accounted for via the dynamics of their density and temperature. The combined model was applied to study the microscopic mechanism of phase transitions during the laser-induced melting and ablation of monolithic crystalline (c-Si) and porous Si targets in a vacuum. The melting thresholds for the monolithic and porous targets were found to be 0.32 J/cm2 and 0.29 J/cm2, respectively. The limited heat conduction mechanism and the absence of internal stress accumulation were found to be involved in the processes responsible for the lowering of the melting threshold in the porous target. The results of this modeling were validated by comparing the melting thresholds obtained in the simulations to the experimental values. A difference in the mechanisms of ablation of the c-Si and porous Si targets was considered. Based on the simulation results, a prediction regarding the mechanism of the laser-assisted production of Si nanoparticles with the desired properties is drawn.

Nanosecond laser processing and surface quality of SiC particle-reinforced AA2024 composite, homogeneous AA2024 aluminum alloy, and SiC: A comparative experimental study

This study performed comparative experiments on ablation behavior and surface quality obtained by nanosecond laser processing of silicon carbide (SiC) particle-reinforced 2024 aluminum alloy matrix composite (SiCp/AA2024), homogeneous SiC, and 2024 aluminum alloy. A nanosecond laser with a pulse fluence range of 0.68 ∼ 6.83 J/cm2 was used at irradiate samples of these materials, which were then ablated by various numbers of scanning passes (1 to 10) at a fluence of 6.83 J/cm2. The material melting could be achieved only for homogeneous AA2024. The ablation depth for SiC exceeded that of the composite, and the former could be ablated with a pulse fluence exceeding 3.42 J/cm2. The composite ablation threshold was higher than 4.10 J/cm2. Due to the thermal transition between the reinforcement and matrix, the melting threshold of the composite’s AA2024 matrix was lower than that of homogeneous AA2024 alloy. Moreover, the ablation threshold of SiC particles exceeded that of homogeneous SiC. The ablated surface roughness for the homogeneous material showed a sudden increase because of melting or oxidization, while that of the composite exhibited a gradual increase. For homogeneous SiC, oxidization accompanied by ablation and the accumulation of oxidization products on the surface deteriorated the laser pulse processing quality, especially in multipass scanning. In contrast, the dispersed distribution of SiC particles in the composite improved the composite’s machining feasibility.

Numerical study of ultrashort laser-induced microjet formation on the metal film based on the Navier–Stokes equation

A two-temperature model (TTM) for the electron-phonon thermal equilibrium is used to determine the heat distribution and laser fluence threshold for melting a thin metal film coated on a glass substrate and irradiated by an ultrashort laser pulse. This study proposes a novel model based on the Navier–Stokes equation to explain the formation of jet-shaped structures in the film's molten region. By solving this equation and obtaining the temporal evolution of the velocity distribution and displacement in the molten region, the Marangoni convection effect can be numerically demonstrated, and the circular motion of the fluid can describe the formation of a jet-shaped structure in the central region of the radiation. The results are compared to those obtained by numerically solving the thermo-elastoplastic equations, and also, to the previously reported experimental results to ensure the accuracy of the microjet height calculated by the Navier–Stokes equation. Good agreement is observed, particularly when the temperature of the irradiated area is significantly over the film's melting temperature. In addition, several calculations are performed for various pulse fluences. In both models, increasing the pulse fluences leads to an increase in the height of microjets.

Glacier projections sensitivity to temperature-index model choices and calibration strategies

AbstractThe uncertainty of glacier change projections is largely influenced by glacier models. In this study, we focus on temperature-index mass-balance (MB) models and their calibration. Using the Open Global Glacier Model (OGGM), we examine the influence of different surface-type dependent degree-day factors, temporal climate resolutions (daily, monthly) and downscaling options (temperature lapse rates, temperature and precipitation corrections) for 88 glaciers with in-situ observations. Our findings indicate that higher spatial and temporal resolution observations enhance MB gradient representation due to an improved calibration. The addition of surface-type distinction in the model also improves MB gradients, but the lack of independent observations limits our ability to demonstrate the added value of increased model complexity. Some model choices have systematic effects, for example weaker temperature lapse rates result in smaller projected glaciers. However, we often find counter balancing effects, such as the sensitivity to different degree-day factors for snow, firn and ice, which depends on how the glacier accumulation area ratio changes in the future. Similarly, using daily versus monthly climate data can affect glaciers differently depending on the shifting balance between melt and solid precipitation thresholds. Our study highlights the importance of considering minor model design differences to predict future glacier volumes and runoff accurately.

Orientation-dependent phase transition pathways of single-crystal nickel over large shock range

The response of single crystals to increasing shock strength involves complex phase transition processes from solid-solid to solid-liquid phases. However, it is not yet fully understood how these phase transitions occur along different crystallographic orientations, and what the dominant micro-mechanisms are at different shock stages. To fill this knowledge gap, non-equilibrium molecular dynamics simulations are performed to explore the dynamic phase-transition characteristics of face-centered cubic Ni single crystal over a large shock strength range of 0–350 GPa. The orientation-dependent phase transition in single-crystal Ni is examined through three typical shock directions: [001], [011], and [111], and a phase-transition locus map is established. The results show that the shock compression of single-crystal Ni undergoes four response stages: (I) elasticity, (II) plasticity, (III) partial melting, and (IV) complete melting. The shock orientation plays a crucial role in the shock wave structure, Hugoniot, microstructure evolution, and shock melting behavior. The [001] shock direction maintains a single-wave structure and exhibits a shear stress superelasticity. However, the [011]/[111] orientations undergo a single-double-single wave transformation, and an interesting U-shape transition wave is identified along the [111] orientation due to its unique progressive plasticity. Superheating and cold melting phenomena are observed along the [001] and [011]/[111] orientations, respectively. The differences in the critical melting temperature and threshold shock pressure along different orientations are up to ∼3000 K and ∼80 GP. Atomic-structure analysis reveals that due to the Bain path, FCC→BCC (BCT) is the dominant solid-solid phase transition along the [001] orientation, and the dislocation slips are limited in a narrow shock window. Melting nuclei emerge, grow, and coalesce homogeneously on the BCC (BCT) substrate, leading to the superheating phenomenon. In contrast, a dynamic process involving shock melting and recrystallization is observed along the [011]/[111] orientations, which finally results in the coexisting of the solid and liquid phases. Heterogeneous melting, marked by mixed phases, is the most distinctive feature in the [011]/[111] orientations, which is responsible for the cold melting phenomenon. Our study shed light on the shock-induced phase transitions of FCC single crystal, which may offer insight into understanding the complex polycrystalline.

Melting Thresholds of Materials Irradiated with a Wide Class of Pulsed Electron Beams

Based on the proposed criterion of the type of heating, a classification of the sources of pulsed electron beams was carried out, both to obtain a better understanding of the nature of the thermal processes occurring under irradiation and to predict their suitability for certain applications. The melting thresholds of materials were calculated over a wide ranges of accelerating voltages and pulse durations. On the basis of calculations, a refractoriness series was proposed for metals for surface–volume pulsed heating.

Dependence of impact regime boundaries on the initial temperatures of projectiles and targets

Towards higher impact velocities, ballistic events are increasingly determined by the material temperatures. Related effects might range from moderate thermal softening to full phase transition. In particular, it is of great interest to quantify the conditions for incipient or full melting of metals during impact interactions, which result in a transition from still strength-affected to hydrodynamic material behavior. In this work, we investigate to which extent the respective melting thresholds are also dependent on the initial, and generally elevated, temperatures of projectiles and targets before impact. This is studied through the application of a model developed recently by the authors to characterize the transition regime between high-velocity and hypervelocity impact, for which the melting thresholds of materials were used as the defining quantities. The obtained results are expected to be of general interest for ballistic application cases where projectiles or targets are preheated. Such conditions might result, for example, from aerodynamic forces acting onto a projectile during atmospheric flight, explosive shaped-charge-jet formation or armor exposure to environmental conditions. The performed analyses also broaden the scientific understanding of the relevance of temperature in penetration events, generally known since the 1960s, but often not considered thoroughly in impact studies.

Characteristics of Surface “Melt Potential” over Antarctic Ice Shelves based on Regional Atmospheric Model Simulations of Summer Air Temperature Extremes from 1979/80 to 2018/19

Abstract We calculate a regional surface “melt potential” index (MPI) over Antarctic ice shelves that describes the frequency (MPI-freq; %) and intensity (MPI-int; K) of daily maximum summer temperatures exceeding a melt threshold of 273.15 K. This is used to determine which ice shelves are vulnerable to melt-induced hydrofracture and is calculated using near-surface temperature output for each summer from 1979/80 to 2018/19 from two high-resolution regional atmospheric model hindcasts (using the MetUM and HIRHAM5). MPI is highest for Antarctic Peninsula ice shelves (MPI-freq 23%–35%, MPI-int 1.2–2.1 K), lowest (2%–3%, <0 K) for the Ronne–Filchner and Ross ice shelves, and around 10%–24% and 0.6–1.7 K for the other West and East Antarctic ice shelves. Hotspots of MPI are apparent over many ice shelves, and they also show a decreasing trend in MPI-freq. The regional circulation patterns associated with high MPI values over West and East Antarctic ice shelves are remarkably consistent for their respective region but tied to different large-scale climate forcings. The West Antarctic circulation resembles the central Pacific El Niño pattern with a stationary Rossby wave and a strong anticyclone over the high-latitude South Pacific. By contrast, the East Antarctic circulation comprises a zonally symmetric negative Southern Annular Mode pattern with a strong regional anticyclone on the plateau and enhanced coastal easterlies/weakened Southern Ocean westerlies. Values of MPI are 3–4 times larger for a lower temperature/melt threshold of 271.15 K used in a sensitivity test, as melting can occur at temperatures lower than 273.15 K depending on snowpack properties.