The problem begins with scale. To study how climate change has reshaped Alpine catchments since the mid-twentieth century, researchers often use the Weather Research and Forecasting model, or WRF, to dynamically downscale coarse global reanalysis data to a much finer grid. In this study, the team took the Twentieth Century Reanalysis Version 3, a global dataset with a spatial resolution of roughly 100 kilometers, and downscaled it to 2 kilometers with hourly output for the period 1969 to 2015. That is a remarkable achievement, but even 2 kilometers is too coarse to capture the fine-grained reality of a glacierized catchment, where snowfields, rock debris, and ice interact with the atmosphere in ways that vary from one slope to the next.

When the researchers evaluated the WRF output against observations from four meteorological stations spanning nearly 900 meters of elevation in the upper Kauner Valley, they found a pronounced cold bias in the diagnosed 2-meter air temperature at elevations above roughly 2600 meters. The culprit was the model’s coarse representation of snow- and ice-covered surfaces. In the WRF land surface scheme, grid cells above about 2900 meters were classified as permanent Snow or Ice, and the cooling influence of those bright, frozen surfaces bled into the simulated near-surface air temperature. The result was air that was too cold in exactly the places where snowmelt and glacier melt are most sensitive to temperature, a distortion that would ripple through every downstream calculation.

Conventional bias correction, such as empirical quantile mapping, adjusts model output by matching its statistical distribution to that of observations. But that approach requires long, high-quality observational records, and such records are notoriously scarce at high elevations in the Alps. Stations above 2700 meters are few, gappy, and prone to sensor failures, and their statistical properties may not be stationary over decades of rapid warming. Rather than lean on sparse data, Hofmeister and his colleagues turned to the model’s own atmosphere for the answer.

Their method exploits a fundamental property of the atmosphere: temperature generally decreases with height according to a lapse rate. In the WRF simulations, the researchers found that at the lowest atmospheric model levels, temperatures were strongly influenced by the underlying land cover, with snow- and ice-covered grid cells showing markedly colder values and considerable scatter. But higher up, between the first and twentieth vertical levels, the temperature-elevation relationship became almost perfectly linear, with very little scatter. At these heights, the model atmosphere was essentially free of surface contamination, governed by a uniform domain-wide lapse-rate structure. The lapse rate at the first model level averaged minus 0.78 degrees Celsius per 100 meters in a test period, close to the free-atmospheric value of minus 0.65, and it decreased steadily with height until it approached zero at the highest levels.

The correction workflow is a three-step procedure applied to hourly time series. First, the team computed spatial averages of hourly temperature and geopotential height across the 49 grid cells covering roughly 200 square kilometers of the Kauner Valley domain, at the first and twentieth vertical levels. Averaging in space reduces the dependence of the lapse rate on land cover type. Second, they calculated an hourly temperature lapse rate from the difference between the two levels, capturing the sub-daily variability of the atmospheric temperature gradient rather than assuming a fixed value. Third, they used that hourly lapse rate to extrapolate the temperature at the first atmospheric model level, roughly 17 meters above the ground on average, down to any station or catchment elevation. Because the reconstruction starts from the free atmosphere rather than the surface, it sidesteps the surface-related biases entirely. The mean lapse rate derived between the two levels over 1973 to 2015 was minus 0.57 degrees Celsius per 100 meters, slightly steeper in summer than in winter.

Temperature alone was not enough. Because snowmelt and glacier melt models rely on energy balance calculations in which latent heat fluxes play a substantial role, the team also applied a sequential, thermodynamically consistent correction of relative humidity, adjusting moisture to match the corrected temperatures so that the physical relationship between the two variables was preserved. Notably, they left shortwave radiation and wind speed untouched, since physics-based hydrological models such as WaSiM, the Water Flow and Balance Simulation Model used in this study, already account for topographic effects on radiation and can represent wind-driven snow redistribution on their own.

The validation was rigorous and multi-layered. Against station observations, the corrected temperatures showed substantially better agreement than the uncorrected output, with the improvement most dramatic at the highest station, Vernagt Pegel at 2640 meters, where the maximum distance between observed and modeled temperature distributions dropped from 0.13 to 0.03. The consequences for hydrological modeling were striking. When WaSiM was forced with uncorrected WRF temperatures, it could not reproduce the observed dynamics of snow accumulation and melt, glacier mass balance, snow cover, or streamflow: the cold bias suppressed melt so severely that simulated annual streamflow was nearly 50 percent too low, and the simulated glacier mass balance stayed unrealistically positive during a period when Alpine glaciers were actually shrinking. With the corrected temperatures, the model matched the observed mean annual streamflow almost exactly, and the Kling-Gupta efficiency of daily streamflow improved from 0.42 to 0.71. Adding the humidity correction pushed performance further, raising the efficiency to 0.78 and improving the representation of late-summer flows during peak glacier melt, when the Fagge carries its maximum meltwater load.

Crucially for climate research, the correction preserved the long-term warming signal. Comparing the corrected product with homogenized observations at Obergurgl, a station in a neighboring valley with records stretching back to 1851, the team found that both the observations and the corrected WRF temperatures showed a significant warming trend of about 0.05 to 0.06 degrees Celsius per year over 1973 to 2015, and both detected the well-documented abrupt Alpine warming shift in May 1987. A 42-year simulation of the catchment reproduced the observed transition from near-equilibrium glacier conditions in the 1980s to continuous mass loss thereafter, along with a significant shortening of snow cover duration after 1985. Because the method relies on lapse rates derived from the climate model itself rather than on local observations, it is transferable to other glacierized mountain regions where ground data are lacking. As the authors note, future work should test the approach in catchments with different climatic and topographic settings, but the message is already clear: physically consistent corrections of model output can rescue hydrological simulations in the world’s most water-critical landscapes, without erasing the climate signal that makes those simulations worth running in the first place.

Subject of Research: Physically based bias correction of regional climate model temperatures for hydrological simulation in a glacierized Alpine catchment

Article Title: Technical note: A physically based correction of WRF temperature using atmospheric lapse rates for high-resolution hydrological simulations in a glacierized catchment

Article References: Hofmeister, F., Fan, X., Pfeiffer, M., Marzeion, B., Schaefli, B., & Chiogna, G. (2026). Technical note: A physically based correction of WRF temperature using atmospheric lapse rates for high-resolution hydrological simulations in a glacierized catchment. Hydrology and Earth System Sciences, 30(19), 6265-6288. https://doi.org/10.5194/hess-30-6265-2026

Image Credits: AI Generated

DOI: 10.5194/hess-30-6265-2026

Keywords: WRF model, temperature lapse rate, bias correction, glacierized catchment, hydrological modeling, Alps, snowmelt, glacier mass balance, streamflow, regional climate model, relative humidity, climate change

The effort behind the data was anything but routine. High Mountain Asia contains the largest glacierized area outside the polar regions, yet field research there is notoriously difficult. Low barometric pressure, extreme weather and steep terrain challenge both personnel and instruments, and the region’s stations are frequently disabled by riming, heavy snowfall, battery depletion and communication failures. The Everest stations were installed as part of the National Geographic and Rolex Perpetual Planet Everest Expeditions, with Sherpa climbers bolting custom aluminium tripods designed by Campbell Scientific directly to rock and guying them to additional anchor points. The Balcony station at 8,430 metres was toppled by extreme winds and anchor failure in January 2020, prompting a return expedition in May 2022 that installed the new Bishop Rock station at 8,810 metres and upgraded the South Col site.

Each station carries a carefully engineered suite of sensors. Air temperature and relative humidity are measured with Vaisala HMP155A-L5-PT probes housed in naturally ventilated 14-plate solar radiation shields, while barometric pressure is recorded with Vaisala PTB110 and PTB210 sensors inside the datalogger enclosures. Wind speed and direction initially came from R. M. Young 05108-45 anemometers, but all four units installed at the two highest sites failed in the extreme conditions, and the team replaced them with a polycarbonate version of the sensor plus a custom Richards C5C anemometer and pitot tube made by the Mount Washington Observatory. Radiation components are captured with Hukseflux pyranometers and pyrgeometers, and precipitation is measured at the two lowest stations with OTT Pluvio2 weighing gauges fitted with double-Alter wind shields to reduce the wind-induced undercatch that plagues snowfall measurement in exposed mountain terrain.

The engineering details reveal how altitude dictated design choices. At the lower stations, sensors sit 2 metres above the ground and batteries live inside the datalogger box; at the highest sites, sensors were mounted lower, at about 1.5 metres, to reduce both wind loading and the leverage, or torque, that a taller pole would exert on the tripod base during severe gusts. Batteries at the upper stations were placed in separate insulated boxes to protect them from extreme cold, and both boxes were bolted to rock to cut wind drag. The lower stations weigh roughly 60 kilograms each, the upper ones about 52 kilograms including the experimental pitot tube. Data are logged on Campbell Scientific CR1000X dataloggers at intervals from 10 minutes to daily, and the published archive provides quality-controlled hourly values in Coordinated Universal Time.

Because raw readings from such a hostile environment are riddled with artefacts, the team applied a multistage quality control procedure. Relative humidity values exceeding 100 percent were capped at that physical limit, and artificially low readings caused by the logger computing humidity with respect to water rather than ice at sub-zero temperatures were corrected using Buck’s equations for saturation vapour pressure. Periods of zero wind speed combined with zero directional variability were flagged as sensor freezing and removed rather than mistaken for calm conditions. Night-time incoming shortwave radiation below 7 watts per square metre was set to zero, and an albedo-based correction recalculated compromised radiation values whenever fresh snow or rime on the upward-facing sensor pushed the apparent surface albedo above the realistic threshold of 0.95.

The resulting climatology, though based on a short record, paints a vivid seasonal picture. Mean annual temperatures were 4.1 degrees Celsius at Phortse, minus 3.1 degrees at Base Camp and minus 10.2 degrees at Camp II. At South Col, July was the warmest month with a mean of minus 12.2 degrees and a striking diurnal range of 10.2 degrees, while February averaged minus 29.7 degrees. Precipitation clearly delineates the seasons: winters are predominantly dry, amounts build through the pre-monsoon, and the June-to-September monsoon delivers 72 percent of annual precipitation at Phortse and 77 percent at Base Camp. The mean annual precipitation gradient between the two sites, separated by roughly 1,500 metres, was minus 107 millimetres per kilometre, but it weakens by almost half during the monsoon, indicating that relative precipitation drops with elevation are far steeper outside the monsoon season.

Temperature gradients proved equally revealing. Between Phortse and Base Camp the mean temperature gradient was minus 4.8 degrees Celsius per kilometre, least negative in winter and most negative in the pre-monsoon. When the four stations with minimal data gaps were combined, the gradients became more strongly negative, at minus 6.5 degrees per kilometre in the pre-monsoon, minus 5.7 in the monsoon and minus 6.0 annually. The authors caution that the temperature-altitude relationship in the Khumbu region is fundamentally non-linear, which limits direct comparison with lapse rates calculated across different elevational ranges, and they therefore provide a non-linear equation for the gradient between Phortse and South Col. Radiation measurements add further nuance: maximum daily incoming longwave radiation during the monsoon exceeds 350 watts per square metre at Phortse but stays below about 250 at South Col, likely reflecting the colder, drier, less cloudy atmosphere aloft.

A central contribution of the study is a rigorous comparison with ERA5, the fifth-generation global reanalysis produced by the European Centre for Medium-Range Weather Forecasts that many researchers use as a stand-in for sparse mountain observations. Extracted from the nearest grid point at the 350 hectopascal pressure level, ERA5 temperatures tracked observed variability at South Col reasonably well, with coefficients of determination above 0.6 in both 2019 and 2022, and even higher agreement of 0.85 with the short Balcony record in 2019. But the reanalysis systematically underestimated air temperature, showed far less diurnal variability because the pressure level represents the free atmosphere rather than the surface, and consistently overestimated mean wind speeds, with mean absolute errors of 4.8 and 5.3 metres per second in 2019 and 2022. Relative humidity comparisons carried mean absolute errors above 20 percent, though 10-day running means from both datasets clearly captured the monsoon onset at South Col on 1 July 2019 and 14 June 2022.

The practical implications reach well beyond atmospheric science. Observations from South Col reveal mean winds that can exceed 30 metres per second and gusts above 60 metres per second, conditions capable of blowing mountaineers off their feet and inducing cold injuries. Pressure data from the network have already shown that oxygen availability at the summit varies on synoptic timescales, meaning the apparent elevation of Everest, how high the mountain would feel without supplemental oxygen, can shift by almost 750 metres, and a winter ascent without bottled oxygen may at times be impossible. The team anticipates that ERA5, once its time-varying biases are corrected with empirical-statistical or machine-learning approaches, could be used to gap-fill and extend the intermittent records from Everest’s upper slopes, enabling hyper-local forecasts that help expeditions identify optimal climbing windows.

The archive also underpins research on the region’s fragile cryosphere and its role as a freshwater source for downstream communities. It has already been used to estimate surface energy balances at the summit and at South Col Glacier, revealing a high-altitude ice system acutely sensitive to changes in effective precipitation because of extremely high insolation and its responsiveness to albedo variations. Combined with the longer-running EvK2CNR and GLACIOCLIM networks at lower elevations, the new data allow quantification of elevational gradients in key meteorological variables across roughly five vertical kilometres, covering the entire glacierized range of the Khumbu region, information essential for distributed glacier and hydrological modelling. And because meteorological measurements from the highest reaches only began in 2022, the authors note that scientists are still at the very beginnings of exploring the weather of this extreme environment, with questions the data might answer likely to grow rapidly across disciplines far beyond the climate sciences.

Subject of Research: High-altitude weather station observations across the Mount Everest region of Nepal

Article Title: Weather station data from the Mount Everest region, Nepal: 3810–8810 m above sea level

Article References: Khadka, A., Perry, L. B., Matthews, T., Sherpa, T. G., Shrestha, C. B., Shrestha, D., Aryal, D., Tuladhar, S., Pradhananga, N., Kayastha, D., Raichle, B., Athans, P., Sherpa, D. Y., Garrett, K., Wheeler, G., Young, T., & Elmore, A. (2026). Weather station data from the Mount Everest region, Nepal: 3810–8810 m above sea level. Earth System Science Data, 18(10), 7253-7267. https://doi.org/10.5194/essd-18-7253-2026

Image Credits: AI Generated

DOI: 10.5194/essd-18-7253-2026

Keywords: Mount Everest, automatic weather stations, Himalaya, climate data, ERA5 reanalysis, temperature lapse rate, monsoon, glaciology, high-altitude meteorology, Khumbu, quality control, open data

Zhejiang presents an ideal and demanding test case. This coastal province in southeastern China is home to roughly 66.7 million people and an economy exceeding 9 trillion CNY in 2024, yet nearly 75 percent of its land is covered by hills and mountains, with only about 20 percent plains and a sliver of rivers and lakes. Rapid urbanization compounds the challenge, since weather stations are sparse, unevenly distributed, and subject to relocations, instrumentation changes, and gaps in the historical record. Ground observations offer accuracy but limited spatial coverage; satellite land surface temperature products offer detail but generally begin only in the early 2000s and are vulnerable to cloud cover, terrain shading, and atmospheric interference; reanalysis products such as ERA5 provide continuity but at coarse spatial resolutions, typically around 0.25 degrees or coarser, far too blunt to resolve the fine thermal texture of mountainous landscapes.

The team’s solution is an elegant two-part construction they call a spatial-background residual framework. First, they built a high-resolution climatological baseline from an existing hourly, 1-kilometer gridded temperature dataset covering 2008 to 2018, which had itself been developed using the INCA data-fusion framework with reanalysis fields and dense automatic weather station observations. This baseline serves purely as a spatial background, encoding how temperature varies across the terrain on each calendar day of the year. Second, daily temperature residuals—the departures of each station observation from that climatological expectation—were calculated for every meteorological station across the full 1961 to 2020 period. Because these residuals are computed directly from observed temperatures, they inherently preserve the long-term warming trend and interannual variability, while the baseline contributes the terrain-driven spatial detail. Summing the interpolated residual field with the baseline yields the finished reconstruction.

A critical methodological question was how best to interpolate those daily residuals across space. The researchers evaluated three schemes representing different levels of topographic correction and complexity: plain inverse distance weighting, or IDW, which relies only on spatial proximity; a lapse-rate-adjusted version of IDW, in which station temperatures are first corrected to grid-cell elevation using a fixed adiabatic lapse rate of 6.0 degrees Celsius per kilometer; and a multiple linear regression incorporating longitude, latitude, and elevation as predictors. Using leave-one-out cross-validation, in which each station is successively withheld and predicted from the others, the team assessed performance with mean absolute error, root-mean-square error, and the coefficient of determination across decades, seasons, and elevation zones.

The verdict was striking: the simplest method won. IDW without any lapse-rate correction consistently delivered the lowest errors and highest skill, achieving the best performance at roughly 68 percent of stations for daily minimum temperature, 70.7 percent for mean temperature, and 74.7 percent for maximum temperature. In a representative example from 1971 to 1980, IDW reconstructed minimum temperatures with a mean absolute error of just 0.73 degrees Celsius and an R-squared of 0.987, comfortably beating both rivals. The reason lies in the behavior of the lapse rate itself. Analysis of 60 years of observations revealed that near-surface temperature lapse rates in Zhejiang are strongly non-stationary: they peak in summer, with minimum-temperature lapse rates exceeding 7.0 degrees Celsius per kilometer in mountainous areas during July and August, yet collapse toward zero or even turn negative in winter lowlands, where temperature inversions prevail. Applying a fixed correction therefore risks systematic, elevation-related biases—a caution with implications well beyond Zhejiang.

The errors that do remain follow clear and intelligible patterns. Reconstruction accuracy improved steadily from the 1960s onward as station density grew, and summer months outperformed winter months because spatial temperature gradients are weaker in warm weather. Low-elevation areas below 400 meters consistently yielded smaller errors than high-elevation zones, where complex terrain and sparse instrumentation conspire against interpolation. Spatially, larger uncertainties cluster in the mountainous southwest, including parts of Lishui and western Wenzhou, while the plains around Hangzhou, Shaoxing, and Jinhua show excellent agreement, with most stations achieving R-squared values above 0.90 and many above 0.95. Across the entire 60-year span, the annual mean error for all three temperature variables stayed within plus or minus 0.1 degrees Celsius, with no systematic drift across decades—a testament to the temporal stability of the method.

Perhaps the most compelling validation came from real disasters. The team tested ZJ-DAT against two extreme events from 2007, using more than a thousand automatic weather stations as independent ground truth while deliberately excluding any stations that had contributed to the reconstruction. During the cold wave of 4 to 9 March 2007, ZJ-DAT tracked the south-to-north advance of the cold air, accurately reproducing the observed cold centers around Lishui, with R-squared values of 0.65 to 0.79 and root-mean-square errors of 1.01 to 1.67 degrees Celsius. By comparison, the CDAT national dataset managed only moderate agreement, while ERA5 performed poorly, with near-zero or negative correlations and errors approaching 3 degrees Celsius. The heat wave of 30 June to 10 July 2007 told the same story: ZJ-DAT best resolved the core hot zones above 37 degrees Celsius over Jinhua, Shaoxing, and Ningbo and the inland-coastal thermal contrast, while CDAT smoothed away local extremes and ERA5 drifted with warm biases and excessive homogenization. Case studies of cold and heat events in January and July 2020 at four environmentally distinct stations—an island, a mountain site, and two plain stations—further confirmed the reconstruction’s fidelity, with discrepancies generally under 2 degrees Celsius.

Beyond validation, the dataset enabled a first-of-its-kind hazard assessment for the province. Using Gumbel distribution analysis of return periods, the researchers mapped the intensity of extreme cold and heat expected at 5-, 20-, and 50-year recurrence intervals. The results expose stark geographic contrasts in climate risk. Extreme low-temperature hazards concentrate in the northwestern inland regions, where 50-year minimum temperatures plunge below minus 15 degrees Celsius, while the southeastern coast stays comparatively mild. Extreme heat hazards show the opposite pattern, dominated by low-altitude basins in central and northern Zhejiang, where 50-year maximum temperatures climb above 43 to 44 degrees Celsius—figures that carry sobering weight given projections of accelerating heatwave duration under global warming. These maps, grounded in kilometer-scale temperature data rather than coarse reanalysis, offer planners a far sharper picture of where adaptation investments are most needed.

The authors are candid about limitations. Anchoring the reconstruction to a climatology drawn from 2008 to 2018 means the reference field does not reflect the climate state of earlier decades, though because it functions only as a spatial scaffold while temporal signals come from station residuals, warming trends and variability remain intact. The team suggests that future refinements could employ temporally adaptive reference fields. The broader significance, however, is clear: ZJ-DAT demonstrates that a simple, computationally efficient interpolation of station anomalies, layered onto a modern high-resolution climatology, can outperform more elaborate schemes in complex terrain—provided the scheme is chosen with local lapse-rate physics in mind. The framework, and the publicly available dataset released through Zenodo, is designed to be transferable to other topographically complex, observation-limited regions, offering a practical foundation for climate monitoring, extreme-event risk assessment, and adaptation planning as the planet continues to warm.

Subject of Research: High-resolution daily near-surface air temperature reconstruction for Zhejiang Province, China, from 1961 to 2020 using statistical residual interpolation

Article Title: A high-resolution daily temperature reconstruction for Zhejiang Province during 1961–2020 using statistical residual interpolation

Article References: Li, Y., Guo, H., Dong, M., Wu, J., Deng, F., Chen, Y., & Chen, F. (2026). A high-resolution daily temperature reconstruction for Zhejiang Province during 1961–2020 using statistical residual interpolation. Theoretical and Applied Climatology, 157(10), Article 632. https://doi.org/10.1007/s00704-026-06521-3

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06521-3

Keywords: Zhejiang Province, air temperature reconstruction, inverse distance weighting, spatial interpolation, temperature lapse rate, extreme temperature events, heat wave, cold wave, climate monitoring, hazard assessment, ERA5, high-resolution gridded dataset