Wilkes Land and Totten Glacier (TG) in East Antarctica (EA) have been losing ice mass significantly since 1989. There is a lack of knowledge of long-term mass balance in the region which hinders the estimation of its contribution to global sea level rise. Here we show that this acceleration trend in TG has occurred since the 1960s. We reconstruct ice flow velocity fields of 1963–1989 in TG from the first-generation satellite images of ARGON and Landsat-1&4, and build a five decade-long record of ice dynamics. We find a persistent long-term ice discharge rate of 68 ± 1 Gt/y and an acceleration of 0.17 ± 0.02 Gt/y² from 1963 to 2018, making TG the greatest contributor to global sea level rise in EA. We attribute the long-term acceleration near grounding line from 1963 to 2018 to basal melting likely induced by warm modified Circumpolar Deep Water. The speed up in shelf front during 1973–1989 was caused by a large calving front retreat. As the current trend continues, intensified monitoring in the TG region is recommended in the next decades.

Antarctic ice velocity maps describe the ice flow dynamics of the ice sheet and are one of the primary components used to estimate the Antarctic mass balance and contribution to global sea level changes. In comparison to velocity maps derived from recent satellite images of monthly to weekly time spans, historical maps, from before the 1990s, generally cover longer time spans, e.g., over 10 years, due to the scarce spatial and temporal coverage of earlier satellite image data. We found velocity overestimations (OEs) in such long-span maps that can be mainly attributed to velocity gradients and time span of the images used. In general, they are less significant in slow-flowing grounded regions with low spatial accelerations. Instead, they take effect in places of high ice dynamics, for example, near grounding lines and often in ice shelf fronts. Velocities in these areas are important for estimating ice sheet mass balance and analyzing ice shelf instability. We propose an innovative Lagrangian velocity-based method for OE correction without the use of field observations or additional image data. The method is validated by using a set of ground truth velocity maps for the Totten Glacier and Pine Island Glacier which are produced from high-quality Landsat 8 images from 2013 to 2020. Subsequently, the validated method is applied to a historical velocity map of the David Glacier region from images from 1972–1989 acquired during Landsat 1, 4, and 5 satellite missions. It is demonstrated that velocity overestimations of up to 39 m a⁻¹ for David Glacier and 195 m a⁻¹ for Pine Island Glacier can be effectively corrected. Furthermore, temporal acceleration information, e.g., on basal melting and calving activities, is preserved in the corrected velocity maps and can be used for long-term ice flow dynamics analysis. Our experiment results in the Pine Island Glacier (PIG) show that OEs of a 15-year span can reach up to 1300 m a⁻¹ along the grounding line and cause an overestimated grounding line (GL) flux of 11.5 Gt a⁻¹ if not corrected. The magnitudes of the OEs contained in both velocity and mass balance estimates are significant. When used alongside recent velocity maps of 1990s–2010s, they may lead to underestimated long-term changes for assessment and forecast modeling of the global climate change impact on the Antarctic ice sheet. Therefore, the OEs in the long-span historical maps must be seriously examined and corrected. We recommend that overestimations of more than the velocity mapping uncertainty (1σ) be corrected. This velocity overestimation correction method can be applied to the production of regional and ice-sheet-wide historical velocity maps from long-term satellite images.

The seasonal deposition and sublimation of CO₂ constitute a major element in Martian volatile cycles. We reprocess the Mars Orbiter Laser Altimeter (MOLA) data and apply co‐registration procedures to obtain spatio‐temporal variations in levels of the Seasonal North Polar Cap (SNPC). The maximum level over the Residual North Polar Cap (RNPC) is 1.3 m, approximately half of that at the south pole (2.5 m). However, the maximum level in the dune fields at Olympia Undae can be up to 3.8 m. Furthermore, off‐season decreases up to 3 m during the northern winter at Olympia Undae are observed. These are likely due to metamorphism effects accentuated by the reduced snowfall at this period. Meanwhile, off‐season increases of up to 2 m during the northern spring are noted, the cause of which remains to be explored. The volume of the SNPC peaks at the end of northern winter and is estimated to be approximately 9.6 × 10¹² m³, which is 2% more than that of the Seasonal South Polar Cap. The bulk density of the SNPC can go through phased decreases in accordance with phased accumulation at northern high‐latitudes. These findings can put important constraints on the Martian volatile cycling models.

Fractures are an important indicator of ice shelf instability and predecessors of calving events. It is of great importance to identify fractures and understand the formation and propagation processes of fractures. However, automated methods of fracture identification also erroneously extract non‐fracture features, which are mixed with the automatically extracted fractures, negatively impacting the accuracy of the results. In this study, a framework for fracture extraction with glaciological property‐based constraints was proposed to extract fractures and filter out non‐fracture features. The extraction module consists of image preprocessing and linear feature extraction. The constraints were based on glaciological property analysis, including constraints based on physical properties of fractures and structural properties of ice shelves. The framework was applied to Mosaic of Antarctica and Landsat 8 Operational Land Imager images to extract fractures on the Filchner–Ronne Ice Shelf. The results suggest that most of the non‐fracture features were successfully filtered out, accounting for about half of all linear features. Upon comparison, it was found that the proposed method detected ∼97% of the rifts and ∼42% of the crevasses detected by manual extraction. The spatial pattern and formation mechanisms of the fractures were analyzed based on strain maps. The formation of crevasses is closely related to shear margin and longitudinal tension, while the large rifts near the ice shelf front are likely to be related to the strikes of lateral ice flow tributaries. The constraint strategy can be further improved by utilizing more glaciological properties to obtain more‐accurate fracture extraction.

As a critical parameter, thermal conductivity directly determines the heat transfer and temperature variation within rocks, which can lead to mechanical damage and chemical corrosion. Consequently, understanding the thermal conductivity of stone heritage is vital for assessing their deterioration mechanisms and developing effective conservation strategies. This study obtained sandstone samples from the Yungang Grottoes and subjected them to freeze–thaw cycle experiments to generate weathered sandstone samples. Subsequently, the thermal conductivity of these samples was measured under both dry and water-saturated state using the transient plane source method. To analyze the relationship between air–water saturation, porosity, and thermal conductivity, a saturation influence coefficient was introduced. Thereafter, the effectiveness and applicability of 13 commonly used thermal conductivity mixing law prediction models were evaluated based on experimental data. The results suggested that the influence of water saturation on the thermal conductivity of rocks varies with porosity, and water saturation significantly enhances the thermal conductivity of weathered sandstone. Among the 13 common models, the Geometric mean model was found to be more accurate than other models, with superior performance in both dry (MAE, RMSE, MAPE are 0.148, 0.214, 5.59% respectively) and water-saturated (MAE, RMSE, MAPE are 0.244, 0.170, 8.4% respectively) state. The Albert model demonstrates a good fit in the dry state, whereas the Walsh model (with maximum effect), Ribaud model, and Huang model also exhibit good fitting efficacy in the water-saturated state. This study provides a solid foundation for better predicting the thermal conductivity of weathered stone heritage and developing effective preventive conservation strategies.