A novel convection permitting modelling framework that combines a pseudo-global warming approach with continuously forced deep soil moisture from prescribed perturbation storylines is applied in the Eastern European Alpine region and parts of the Pannonian Basin to investigate soil moisture precipitation (SMP) feedbacks on summertime precipitation and the feedbacks' role under changed climate conditions. A set of 1-year convection-permitting (3 km horizontal grid spacing) soil moisture sensitivity simulations with the regional climate model of the Consortium for Small-Scale Modelling in Climate Mode are conducted. In order to account for global warming, end-of-the-century climate change effects from four global climate models, projecting the greenhouse gas concentration scenario RCP 8.5, are imprinted. The simulations reveal that (1) the locations of precipitation events are highly sensitive to soil moisture modifications while intensities and the internal structure of precipitation events are nearly unaffected and (2) high precipitation intensities are more likely in combinations with positive temporal but distinctive (either strong positive or strong negative) spatial SMP coupling. Low precipitation intensities are in favour of combinations of negative temporal and positive spatial coupling. The analyses suggest that soil moisture at a given time acts as a guiding field for the location of the next precipitation event. Interestingly, this behaviour is independent of climate change, although the coupling strength's increase is 1.5-1.7 times larger than expected from linear climate change scaling when climate becomes 50% dryer. Finally, it is found that (1) local deviations in the climate change signal of summertime precipitation in the range of up to ±40% are caused by uncertainty in deep soil moisture in the range of ±10% and (2) these local deviations in the climate change signal are dominated by soil moisture uncertainty in future climate conditions.

Combined heat and humidity is frequently described as the main driver of human heat-related mortality, more so than dry-bulb temperature alone. While based on physiological thinking, this assumption has not been robustly supported by epidemiological evidence. By performing the first systematic comparison of eight heat stress metrics (i.e., temperature combined with humidity and other climate variables) with warm-season mortality, in 604 locations over 39 countries, we find that the optimal metric for modelling mortality varies from country to country. Temperature metrics with no or little humidity modification associates best with mortality in ~40% of the studied countries. Apparent temperature (combined temperature, humidity and wind speed) dominates in another 40% of countries. There is no obvious climate grouping in these results. We recommend, where possible, that researchers use the optimal metric for each country. However, dry-bulb temperature performs similarly to humidity-based heat stress metrics in estimating heat-related mortality in present-day climate.

A 25-year set of daily radiosonde data was used to investigate temperature and humidity inversions at Neumayer Station, coastal Dronning Maud Land, Antarctica. For the first time, inversions were studied differentiating between different synoptic conditions and different height levels. It was shown that, generally, inversions occurred on the majority (78%) of the days, with simultaneous occurrence of humidity and temperature inversions being observed on approximately two thirds of all days. Multiple inversions are common in all seasons for cyclonic and noncyclonic conditions, however, typically occur more frequently under cyclonic conditions. The seasonality of inversion occurrence and features, that is, inversion strength, depth and vertical gradients, was analysed statistically. Different formation mechanisms depending on inversion levels and prevailing weather situations are related to typical annual courses of certain inversion features. Winter maxima were found for the features that are mostly connected to the temperature close to the surface, which is mainly a result of the negative energy balance, thus influencing surface-based inversions. At the second level, both temperature and humidity inversions are often caused by advection of comparably warm and moist air masses related to the passage of cyclones and their frontal systems. Hence, maxima in several inversion features are found in spring and fall, when cyclonic activity is strongest. Monthly mean profiles of humidity and temperature inversions reveal that elevated inversions are often obscured in average profiles due to large variations in inversion height and depth.

Forty years (1980-2019) of reanalysis data were used to investigate climatology and trends of heat stress in the Caribbean region. Represented via the Universal Thermal Climate Index (UTCI), a multivariate thermophysiological-relevant parameter, the highest heat stress is found to be most frequent and geographically widespread during the rainy season (August, September, and October). UTCI trends indicate an increase of more than 0.2°C·decade, with southern Florida and the Lesser Antilles witnessing the greatest upward rates (0.45°C·decade). Correlations with climate variables known to induce heat stress reveal that the increase in heat stress is driven by increases in air temperature and radiation, and decreases in wind speed. Conditions of heat danger, as depicted by the heat index (HI), have intensified since 1980 (+1.2°C) and are found to occur simultaneously to conditions of heat stress suggesting a synergy between heat illnesses and physiological responses to heat. This work also includes the analysis of the record-breaking 2020 heat season during which the UTCI and HI achieved above average values, indicating that local populations most likely experienced heat stress and danger higher than the ones they are used to. These findings confirm the gradual intensification of heat stress in the Caribbean and aim to provide a guidance for heat-related policies in the region.

Seasonal precipitation forecasting is highly challenging for the northwest fringes of Europe due to complex dynamical drivers. Hybrid dynamical-statistical approaches offer potential to improve forecast skill. Here, hindcasts of mean sea level pressure (MSLP) from two dynamical systems (GloSea5 and SEAS5) are used to derive two distinct sets of indices for forecasting winter (DJF) and summer (JJA) precipitation over lead-times of 1-4 months. These indices provide predictors of seasonal precipitation via a multiple linear regression model (MLR) and an artificial neural network (ANN) applied to four Irish rainfall regions and the Island of Ireland. Forecast skill for each model, lead time, and region was evaluated using the correlation coefficient () and mean absolute error (MAE), benchmarked against (a) climatology, (b) bias corrected precipitation hindcasts from both GloSea5 and SEAS5, and (c) a zero-order forecast based on rainfall persistence. The MLR and ANN models produced skilful precipitation forecasts with leads of up to 4 months. In all tests, our hybrid method based on MSLP indices outperformed the three benchmarks (i.e., climatology, bias corrected, and persistence). With correlation coefficients ranging between 0.38 and 0.81 in winter, and between 0.24 and 0.78 in summer, the ANN model outperformed MLR in both seasons in most regions and lead-times. Forecast skill for summer was comparable to that in winter and for some regions/lead times even superior. Our results also show that climatology and persistence performed better than direct use of bias corrected dynamical outputs in most regions and lead-times in terms of MAE. We conclude that the hybrid dynamical-statistical approach developed here-by leveraging useful information about MSLP from dynamical systems-enables more skilful seasonal precipitation forecasts for Ireland, and possibly other locations in western Europe, in both winter and summer.

While weather stations generally capture near-surface ambient air temperature (Ta) at a high temporal resolution to calculate daily values (i.e., daily minimum, mean, and maximum Ta), their fixed locations can limit their spatial coverage and resolution even in densely populated urban areas. As a result, data from weather stations alone may be inadequate for Ta-related epidemiology particularly when the stations are not located in the areas of interest for human exposure assessment. To address this limitation in the Megalopolis of Central Mexico (MCM), we developed the first spatiotemporally resolved hybrid satellite-based land use regression Ta model for the region, home to nearly 30 million people and includes Mexico City and seven more metropolitan areas. Our model predicted daily minimum, mean, and maximum Ta for the years 2003-2019. We used data from 120 weather stations and Land Surface Temperature (LST) data from NASA's MODIS instruments on the Aqua and Terra satellites on a 1 × 1 km grid. We generated a satellite-hybrid mixed-effects model for each year, regressing Ta measurements against land use terms, day-specific random intercepts, and fixed and random LST slopes. We assessed model performance using 10-fold cross-validation at withheld stations. Across all years, the root-mean-square error ranged from 0.92 to 1.92 K and the ranged from .78 to .95. To demonstrate the utility of our model for health research, we evaluated the total number of days in the year 2010 when residents ≥65 years old were exposed to Ta extremes (above 30°C or below 5°C). Our model provides much needed high-quality Ta estimates for epidemiology studies in the MCM region.

The Weather Research and Forecasting (WRF) model and a combination of the Regional Spectral Model (RSM) and the Japanese Meteorological Agency Non-Hydrostatic Model (NHM) were used to dynamically downscale selected CMIP5 global climate models to provide 2-km projections with hourly model output for Puerto Rico and the U.S. Virgin Islands. Two 20-year time slices were downscaled for historical (1986-2005) and future (2041-2060) periods following RCP8.5. Projected changes to mean and extreme temperature and precipitation were quantified for Holdridge life zones within Puerto Rico and for the U.S. Virgin Islands. The evaluation reveals a persistent cold bias for all islands in the U.S. Caribbean, a dry bias across Puerto Rico, and a wet bias on the windward side of mountains within the U.S. Virgin Islands. Despite these biases, model simulations show a robust drying pattern for all islands that is generally larger for Puerto Rico (25% annual rainfall reduction for some life zones) than the U.S. Virgin Islands (12% island average). The largest precipitation reductions are found during the more convectively active afternoon and evening hours. Within Puerto Rico, the model uncertainty increases for the wetter life zones, especially for precipitation. Across the life zones, both models project unprecedented maximum and minimum temperatures that may exceed 200 days annually above the historical baseline with only small changes to the frequency of extreme rainfall. By contrast, in the U.S. Virgin Islands, there is no consensus on the location of the largest drying relative to the windward and leeward side of the islands. However, the models project the largest increases in maximum temperature on the southern side of St. Croix and in higher elevations of St. Thomas and St. John.

Dynamical climate models present an initialization problem due to the poor availability of deep oceanic data, which is required for the model assimilation process. In this sense, teleconnection indices, defined from spatial and temporal patterns of climatic variables, are conceived as useful tools to complement them. In this work, the near-term climate predictability of 35 temperature and 36 precipitation time series of three cities (Barcelona, Bristol and Lisbon) was analysed using two approaches: (a) a statistical-dynamical combination of self-predictable teleconnection indices and long-term climate projections on a local scale and (b) dynamical model outputs obtained from drift-corrected decadal experiments. Fourier and wavelet analyses were used to assess the predictability of seven teleconnection indices thanks to a cross-validation process (with differentiated training and validation periods). The standardized absolute error of teleconnection-based prediction was compared with that obtained from a (9) multi-model ensemble based on the Coupled Model Intercomparison Project Phase 5. Results showed that decadal predictions at horizons between 20 and 30 years are adequate for temperature and precipitation if a teleconnection-based approach is used, while temperature is better predicted from a 5-year horizon using drift-corrected dynamical outputs.

Time series homogenization for 299 of the available precipitation records for the island of Ireland (IENet) was performed. Four modern relative homogenization methods, that is, HOMER, ACMANT, CLIMATOL and AHOPS were applied to this network of station series where contiguous intact monthly records range from 30 to 70 years within the base period 1941-2010. Break detection results are compared between homogenization methods, and coincidences with available documentary information (metadata) were analysed. The lowest (highest) number of breaks were detected with HOMER (ACMANT). Large differences of break frequency were found, namely ACMANT and AHOPS detected 8 times as many breaks than HOMER, while the break frequency with CLIMATOL was intermediate. Also, the ratio of series classified to be homogeneous varies widely between the methods. It is 85% with HOMER, 60% with CLIMATOL, 31% with AHOPS, while only 22% with ACMANT. In a further experiment, all the available time series for Ireland and Northern Ireland, (910 series) were used with ACMANT and CLIMATOL to explore the stability of break frequency for the same 299 series examined in the base experiment. While overall break frequency slightly increased (by 6-13%), the break positions notably changed for individual time series. The number of breaks changed for 59% (23%) of the series with ACMANT (CLIMATOL). For the breaks detected coincidentally by at least three methods including ACMANT and CLIMATOL in the base experiment, the second experiment confirmed the break positions in 86-87% of the breaks. The consequences of these results in relation to the reliability of statistical homogenization are discussed. Sometimes markedly different step functions provide comparable good approaches. However, the accuracy of homogenized time series cannot be related directly to the instability of break detection results.

Historical precipitation records are fundamental for the management of water resources, yet rainfall observations typically span 100-150 years at most, with considerable uncertainties surrounding earlier records. Here, we analyse some of the longest available precipitation records globally, for England and Wales, Scotland and Ireland. To assess the credibility of these records and extend them further back in time, we statistically reconstruct (using independent predictors) monthly precipitation series representing these regions for the period 1748-2000. By applying the Standardized Precipitation Index at 12-month accumulations (SPI-12) to the observed and our reconstructed series we re-evaluate historical meteorological droughts. We find strong agreement between observed and reconstructed drought chronologies in post-1870 records, but divergence in earlier series due to biases in early precipitation observations. Hence, the 1800s decade was less drought prone in our reconstructions relative to observations. Overall, the drought of 1834-1836 was the most intense SPI-12 event in our reconstruction for England and Wales. Newspaper accounts and documentary sources confirm the extent of impacts across England in particular. We also identify a major, "forgotten" drought in 1765-1768 that affected the British-Irish Isles. This was the most intense event in our reconstructions for Ireland and Scotland, and ranks first for accumulated deficits across all three regional series. Moreover, the 1765-1768 event was also the most extreme multi-year drought across all regional series when considering 36-month accumulations (SPI-36). Newspaper and other sources confirm the occurrence and major socio-economic impact of this drought, such as major rivers like the Shannon being fordable by foot. Our results provide new insights into historical droughts across the British Irish Isles. Given the importance of historical droughts for stress-testing the resilience of water resources, drought plans and supply systems, the forgotten drought of 1765-1768 offers perhaps the most extreme benchmark scenario in more than 250-years.

Glaciers are of key importance to freshwater supplies in the Himalayan region. Their growth or decline is among other factors determined by an interaction of 2-m air temperature (TAS) and precipitation rate (PR) and thereof derived positive degree days (PDD) and snow and ice accumulation (SAC). To investigate determining factors in climate projections, we use a model ensemble consisting of 36 CMIP5 general circulation models (GCMs) and 13 regional climate models (RCMs) of two Asian CORDEX domains for two different representative concentration pathways (RCP4.5 and RCP8.5). First, we downsize the ensemble in respect to the models' ability to correctly reproduce dominant circulation patterns (i.e., the Indian summer monsoon [ISM] and western disturbances [WDs]) as well as elevation-dependent trend signals in winter. Within this evaluation, a newly produced data set for the Indus, Ganges and Brahmaputra catchments is used as observational data. The reanalyses WFDEI, ERA-Interim, NCEP/NCAR and JRA-55 are used to further account for observational uncertainty. In a next step, remaining TAS and PR data are bias corrected applying a new bias adjustment method, scale distribution mapping, and subsequently PDD and SAC computed. Finally, we identify and quantify projected climate change effects. Until the end of the century, the ensemble indicates a rise of PDD, especially during summer and for lower altitudes. Also TAS is rising, though the highest increases are shown for higher altitudes and between December and April (DJFMA). PRs connected to the ISM are projected to robustly increase, while signals for PR changes during DJFMA show a higher level of uncertainty and spatial heterogeneity. However, a robust decline in solid precipitation is projected over our research domain, with the exception of a small area in the high mountain Indus catchment where no clear signal emerges.

Melting snow and glacier ice in the Himalaya forms an important source of water for people downstream. Incoming longwave radiation (LW) is an important energy source for melt, but there are only few measurements of LW at high elevation. For the modelling of snow and glacier melt, the LW is therefore often represented by parameterizations that were originally developed for lower elevation environments. With LW measurements at eight stations in three catchments in the Himalaya, with elevations between 3,980 and 6,352 m.a.s.l., we test existing LW parameterizations. We find that these parameterizations generally underestimate the LW, especially in wet (monsoon) conditions, where clouds are abundant and locally formed. We present a new parameterization based only on near-surface temperature and relative humidity, both of which are easy and inexpensive to measure accurately. The new parameterization performs better than the parameterizations available in literature, in some cases halving the root-mean-squared error. The new parameterization is especially improving existing parameterizations in cloudy conditions. We also show that the choice of longwave parameterization strongly affects melt calculations of snow and ice.

We investigate factors influencing European winter (DJFM) air temperatures for the period 1979-2015 with the focus on changes during the recent period of rapid Arctic warming (1998-2015). We employ meteorological reanalyses analysed with a combination of correlation analysis, two pattern clustering techniques, and back-trajectory airmass identification. In all five selected European regions, severe cold winter events lasting at least 4 days are significantly correlated with warm Arctic episodes. Relationships during opposite conditions of warm Europe/cold Arctic are also significant. Correlations have become consistently stronger since 1998. Large-scale pattern analysis reveals that cold spells are associated with the negative phase of the North Atlantic Oscillation (NAO-) and the positive phase of the Scandinavian (SCA+) pattern, which in turn are correlated with the divergence of dry-static energy transport. Warm European extremes are associated with opposite phases of these patterns and the convergence of latent heat transport. Airmass trajectory analysis is consistent with these findings, as airmasses associated with extreme cold events typically originate over continents, while warm events tend to occur with prevailing maritime airmasses. Despite Arctic-wide warming, significant cooling has occurred in northeastern Europe owing to a decrease in adiabatic subsidence heating in airmasses arriving from the southeast, along with increased occurrence of circulation patterns favouring low temperature advection. These dynamic effects dominated over the increased mean temperature of most circulation patterns. Lagged correlation analysis reveals that SCA- and NAO+ are typically preceded by cold Arctic anomalies during the previous 2-3 months, which may aid seasonal forecasting.

Globally, few precipitation records extend to the 18th century. The England Wales Precipitation (EWP) series is a notable exception with continuous monthly records from 1766. EWP has found widespread use across diverse fields of research including trend detection, evaluation of climate model simulations, as a proxy for mid-latitude atmospheric circulation, a predictor in long-term European gridded precipitation data sets, the assessment of drought and extremes, tree-ring reconstructions and as a benchmark for other regional series. A key finding from EWP has been the multi-centennial trends towards wetter winters and drier summers. We statistically reconstruct seasonal EWP using independent, quality-assured temperature, pressure and circulation indices. Using a sleet and snow series for the UK derived by Profs. Gordon Manley and Elizabeth Shaw to examine winter reconstructions, we show that precipitation totals for pre-1870 winters are likely biased low due to gauge under-catch of snowfall and a higher incidence of snowfall during this period. When these factors are accounted for in our reconstructions, the observed trend to wetter winters in EWP is no longer evident. For summer, we find that pre-1820 precipitation totals are too high, likely due to decreasing network density and less certain data at key stations. A significant trend to drier summers is not robustly present in our reconstructions of the EWP series. While our findings are more certain for winter than summer, we highlight (a) that extreme caution should be exercised when using EWP to make inferences about multi-centennial trends, and; (b) that assessments of 18th and 19th Century winter precipitation should be aware of potential snow biases in early records. Our findings underline the importance of continual re-appraisal of established long-term climate data sets as new evidence becomes available. It is also likely that the identified biases in winter EWP have distorted many other long-term European precipitation series.

Despite the importance of snow in alpine regions, little attention has been given to the homogenization of snow depth time series. Snow depth time series are generally characterized by high spatial heterogeneity and low correlation among the time series, and the homogenization thereof is therefore challenging. In this work, we present a comparison between two homogenization methods for mean seasonal snow depth time series available for Austria: the standard normal homogeneity test (SNHT) and HOMOP. The results of the two methods are generally in good agreement for high elevation sites. For low elevation sites, HOMOP often identifies suspicious breakpoints (that cannot be confirmed by metadata and only occur in relation to seasons with particularly low mean snow depth), while the SNHT classifies the time series as homogeneous. We therefore suggest applying both methods to verify the reliability of the detected breakpoints. The number of computed anomalies is more sensitive to inhomogeneities than trend analysis performed with the Mann-Kendall test. Nevertheless, the homogenized dataset shows an increased number of stations with negative snow depth trends and characterized by consecutive negative anomalies starting from the late 1980s and early 1990s, which was in agreement with the observations available for several stations in the Alps. In summary, homogenization of snow depth data is possible, relevant and should be carried out prior to performing climatological analysis.

Snowfall in the six basins of the Catskill/Delaware Watershed in south-central New York State historically contributes roughly 20-30% of the water resources derived from the watershed for use in the New York City water supply. The watershed regularly experiences snowfall from three distinctive weather patterns: coastal mid-latitude cyclones, overrunning systems, and lake-effect or Great Lakes enhanced storms. Using synoptic weather classification techniques, these distinct regional atmospheric patterns impacting the watershed are isolated and analysed in conjunction with daily snowfall observations from 1960 to 2009 to allow the influence of each synoptic weather pattern on snowfall to be evaluated independently. Results indicate that snowfall-producing events occur on average approximately 63 days/year, or once every 4 days during the October-May season, leading to an average of 213 cm/year of snowfall within the watershed. Snowfall from Great Lakes enhanced storms and overrunning systems contribute nearly equally to seasonal totals, representing 38 and 39%, respectively. Coastal mid-latitude cyclones, while producing the highest amount of snowfall per event on average, contribute only 16% to the watershed average total snowfall. Predicted climate change is expected to impact snowfall differently depending on the specific atmospheric pattern producing the snow. As such, quantifying the contribution of snowfall to the watershed by synoptic pattern can inform future water management and reservoir operation practices for the New York City Water Supply Management System.

Winter windstorms are known to be among the most dangerous and loss intensive natural hazards in Europe. In order to gain a better understanding of their variability and driving mechanisms, this study analyses the temporal variability which is often referred to as serial or seasonal clustering. This is realized by developing a statistical model relating the winter storm counts to known teleconnection patterns affecting European weather and climate conditions (e.g., North Atlantic Oscillation [NAO], Scandinavian pattern [SCA], etc.). The statistical model is developed via a stepwise Poisson regression approach that is applied to windstorm counts and large-scale indices retrieved from the ERA-20C reanalysis. Significant large-scale drivers accountable for the inter-annual variability of storms for several European regions are identified and compared. In addition to the SCA and the NAO which are found to be the essential drivers for most areas within the European domain, other teleconnections (e.g., East Atlantic pattern) are found to be more significant for the inter-annual variability in certain regions. Furthermore, the statistical model allows an estimation of the expected number of storms per winter season and also whether a season has the characteristic of being what we define an active or inactive season. The statistical model reveals high skill particularly over British Isles and central Europe; however, even for regions with less frequent storm events (e.g., southern and eastern Europe) the model shows adequate positive skill. This feature could be of specific interest for the actuarial sector.

Precipitation patterns over Europe are largely controlled by atmospheric cyclones embedded in the general circulation of the mid-latitudes. This study evaluates the climatologic features of precipitation for selected regions in central Europe with respect to cyclone track types for 1959-2015, focusing on large-scale heavy precipitation. The analysis suggests that each of the cyclone track types is connected to a specific pattern of the upper level atmospheric flow, usually characterized by a major trough located over Europe. A dominant upper level cut-off low (COL) is found over Europe for strong continental (CON) and van Bebber's type (Vb) cyclones which move from the east and southeast into central Europe. Strong Vb cyclones revealed the longest residence times, mainly due to circular propagation paths. The central European cyclone precipitation climate can largely be explained by seasonal track-type frequency and cyclone intensity; however, additional factors are needed to explain a secondary precipitation maximum in early autumn. The occurrence of large precipitation totals for track events is strongly related to the track type and the region, with the highest value of 45% of all Vb cyclones connected to heavy precipitation in summer over the Czech Republic and eastern Austria. In western Germany, Atlantic winter cyclones are most relevant for heavy precipitation. The analysis of the top 50 precipitation events revealed an outstanding heavy precipitation period from 2006 to 2011 in the Czech Republic, but no gradual long-term change. The findings help better understand spatio-temporal variability of heavy precipitation in the context of floods and may be used for evaluating climate models.

Weather-pattern, or weather-type, classifications are a valuable tool in many applications as they characterize the broad-scale atmospheric circulation over a given region. This study analyses the aspects of regional UK precipitation and meteorological drought climatology with respect to a new set of objectively defined weather patterns. These new patterns are currently being used by the Met Office in several probabilistic forecasting applications driven by ensemble forecasting systems. Weather pattern definitions and daily occurrences are mapped to Lamb weather types (LWTs), and parallels between the two classifications are drawn. Daily precipitation distributions are associated with each weather pattern and LWT. Standardized precipitation index (SPI) and drought severity index (DSI) series are calculated for a range of aggregation periods and seasons. Monthly weather-pattern frequency anomalies are calculated for SPI wet and dry periods and for the 5% most intense DSI-based drought months. The new weather-pattern definitions and daily occurrences largely agree with their respective LWTs, allowing comparison between the two classifications. There is also broad agreement between weather pattern and LWT changes in frequencies. The new data set is shown to be adequate for precipitation-based analyses in the UK, although a smaller set of clustered weather patterns is not. Furthermore, intra-pattern precipitation variability is lower in the new classification compared to the LWTs, which is an advantage in this context. Six of the new weather patterns are associated with drought over the entire UK, with several other patterns linked to regional drought. It is demonstrated that the new data set of weather patterns offers a new opportunity for classification-based analyses in the UK.

Flexible spatio-temporal models are widely used to create reliable and accurate estimates for precipitation climatologies. Most models are based on square root transformed monthly or annual means, where a normal distribution seems to be appropriate. This assumption becomes invalid on a daily time scale as the observations involve large fractions of zero observations and are limited to non-negative values. We develop a novel spatio-temporal model to estimate the full climatological distribution of precipitation on a daily time scale over complex terrain using a left-censored normal distribution. The results demonstrate that the new method is able to account for the non-normal distribution and the large fraction of zero observations. The new climatology provides the full climatological distribution on a very high spatial and temporal resolution, and is competitive with, or even outperforms existing methods, even for arbitrary locations.

Sub-daily rainfall extremes may be associated with flash flooding, particularly in urban areas but, compared with extremes on daily timescales, have been relatively little studied in many regions. This paper describes a new, hourly rainfall dataset for the UK based on ∼1600 rain gauges from three different data sources. This includes tipping bucket rain gauge data from the UK Environment Agency (EA), which has been collected for operational purposes, principally flood forecasting. Significant problems in the use of such data for the analysis of extreme events include the recording of accumulated totals, high frequency bucket tips, rain gauge recording errors and the non-operation of gauges. Given the prospect of an intensification of short-duration rainfall in a warming climate, the identification of such errors is essential if sub-daily datasets are to be used to better understand extreme events. We therefore first describe a series of procedures developed to quality control this new dataset. We then analyse ∼380 gauges with near-complete hourly records for 1992-2011 and map the seasonal climatology of intense rainfall based on UK hourly extremes using annual maxima, n-largest events and fixed threshold approaches. We find that the highest frequencies and intensities of hourly extreme rainfall occur during summer when the usual orographically defined pattern of extreme rainfall is replaced by a weaker, north-south pattern. A strong diurnal cycle in hourly extremes, peaking in late afternoon to early evening, is also identified in summer and, for some areas, in spring. This likely reflects the different mechanisms that generate sub-daily rainfall, with convection dominating during summer. The resulting quality-controlled hourly rainfall dataset will provide considerable value in several contexts, including the development of standard, globally applicable quality-control procedures for sub-daily data, the validation of the new generation of very high-resolution climate models and improved understanding of the drivers of extreme rainfall.