Understanding climate-glacier dynamics in High Mountain Asia is of critical importance to address issues including water resources, sea-level rise, mountain geodynamics, natural hazards and ecosystem sustainability. The Karakoram Himalaya is arguably the least understood region, given its extreme topography, climate-system coupling, and advancing and surge-type glaciers that exhibit complex flow patterns. Glacier fluctuations in the Karakoram Himalaya are highly variable in space and time because of numerous controlling factors, including the westerlies, the Indian summer monsoon, various teleconnections, topographic effects, glacier debris-cover characteristics, glacier dynamics, and geological conditions. The influence of the integrative coupling of forcing factors, however, has not been adequately assessed for characterizing the glaciers in the Karakoram Himalaya. Given the scarcity of in-situ data and the difficulty of conducting fieldwork on these glaciers, recent research has focused on utilizing remote sensing, geospatial technologies, and scientific modeling to obtain baseline information about the state of glaciers in the region. This review summarizes our current knowledge of glaciers, climate-glacier interaction, and topographic forcing in the Karakoram Himalaya, and demonstrates the complexities in mountain geodynamics that influence climate-glacier dynamics. Innovative analysis is also presented in support of our review and discussion.

Under global warming, Arctic sea ice has declined significantly in recent decades, with years of extremely low sea ice occurring more frequently. Recent studies suggest that teleconnections with large-scale climate patterns could induce the observed extreme sea ice loss. In this study, a probabilistic analysis of Arctic sea ice was conducted using quantile regression analysis with covariates, including time and climate indices. From temporal trends at quantile levels from 0.01 to 0.99, Arctic sea ice shows statistically significant decreases over all quantile levels, although of different magnitudes at different quantiles. At the representative extreme quantile levels of the 5th and 95th percentiles, the Arctic Oscillation (AO), the North Atlantic Oscillation (NAO), and the Pacific–North American pattern (PNA) have more significant influence on Arctic sea ice than El Niño–Southern Oscillation (ENSO), the Pacific decadal oscillation (PDO), and the Atlantic multidecadal oscillation (AMO). Positive AO as well as positive NAO contribute to low winter sea ice, and a positive PNA contributes to low summer Arctic sea ice. If, in addition to these conditions, there is concurrently positive AMO and PDO, the sea ice decrease is amplified. Teleconnections between Arctic sea ice and the climate patterns were demonstrated through a composite analysis of the climate variables. The anomalously strong anticyclonic circulation during the years of positive AO, NAO, and PNA promotes more sea ice export through Fram Strait, resulting in excessive sea ice loss. The probabilistic analyses of the teleconnections between the Arctic sea ice and climate patterns confirm the crucial role that the climate patterns and their combinations play in overall sea ice reduction, but particularly for the low and high quantiles of sea ice concentration.

Between 11 and 14 August 2017, the southern belt of the central Himalaya experienced extreme precipitation, with some stations recording more than 500 mm of accumulated rainfall, which resulted in widespread, devastating flooding. Precipitation was concentrated over the sub-Himalaya, and the established forecasting systems failed to predict the event. In this study, we evaluate the performance of six cloud microphysics schemes in the Weather Research and Forecasting (WRF) model forced with the advanced ERA5 dataset. We also examine the importance of the cumulus scheme in WRF at 3 km horizontal grid spacing in highly convective events like this. Six WRF simulations, each with one of the six different microphysics schemes with the Kain–Fritsch cumulus scheme turned off, all fail to reproduce the spatial variability of accumulated precipitation during this devastating flood-producing precipitation event. In contrast, the simulations exhibit greatly improved performance with the cumulus scheme turned on. In this study, the cumulus scheme helps to initiate convection, after which grid-scale precipitation becomes dominant. Amongst the different simulations, the WRF simulation using the Morrison microphysics scheme with the cumulus turned on displayed the best performance, with the smallest normalized root mean square error (NRMSE) of 0.25 and percentage bias (PBIAS) of −6.99%. The analysis of cloud microphysics using the two best-performing simulations reveals that the event is strongly convective, and it is essential to keep the cumulus scheme on for such convective events and capture all the precipitation characteristics showing that in regions of extreme topography, the cumulus scheme is still necessary even down to the grid spacing of at least 3 km.