Statistical Downscaling Methods Research Articles

Climate change impacts on magnitude and frequency of urban floods under scenario and model uncertainties

Many studies have confirmed that climate change leads to frequent urban flooding, which can lead to significant socioeconomic repercussions. However, most existing studies have not evaluated the impacts of climate change on urban flood from both event-scale and annual-scale dimensions. In addition, there are only few studies that simultaneously consider scenario and model uncertainties of climate change, and combine flood risk assessment and uncertainty analysis results to provide practical suggestions for urban drainage system management. This study uses the statistical downscaling method to calculate the design rainfall under ten rainfall return periods of four climate models and three climate change scenarios in 2040s, 2060s, and 2080s in various prefecture-level cities in China. The four climate models are HadGEM2- ES, MPI-ESM-MR, NorESM1-M and FGOALS-g2 models and the three climate change scenarios are constructed by different representative concentration pathways (RCP), namely RCP2.6, RCP4.5 and RCP8.5. On this basis, relying on the generated drainage systems using geographical information and other data, event-scale and annual-scale precipitation are combined to calculate the change ratio of annual flood volume expectation in prefecture-level cities in each future year compared with the current situation. Furthermore, the study evaluates scenario and model uncertainties of climate change, and then comprehensively integrates the flood risk and its uncertainties to provides suggestions for urban drainage system management.

Future challenges of terrestrial water storage over the arid regions of Central Asia

Since the arid regions of Central Asia (ACA) are located in the interior of Eurasia, water resources play a vital role in the stability of its ecosystem and economic development. Based on the terrestrial water storage anomaly (TWSA) of the Gravity Recovery and Climate Experiment (GRACE), we analyze the observed characteristics of the TWSA over the ACA during 2003–2014. Results indicate that the terrestrial water storage (TWS) in the region showed an overall declining trend from 2003 to 2014, and the autumn TWS in this region is the smallest compared to other seasons and exhibits a strong decreasing trend at least −4.5 cm/decade. This means water resources over the ACA are scarcer and more vulnerable in autumn. The Distance between Indices of Simulation and Observation (DISO) method is employed to evaluate the performance of the sixth phase of the Coupled Model Intercomparison Project (CMIP6) models in simulating the autumn TWSA over the ACA. Compared with observational results, the autumn TWSA values captured by CMIP6 models are larger and the declining TWS trends are weaker. Using the optimal CMIP6 models, the statistical downscaling method constrains the projection results of autumn TWSA values over the ACA using the GRACE datasets. It shows autumn TWS will continue to decrease in most parts of the ACA in the future, and water scarcity will be the most severe in Tajikistan and southwestern Kazakhstan. Under SSP126, Tajikistan’s TWSA is projected to decrease by 11.0 cm in the long term. This study reveals the current situation and possible future changes in TWS over the ACA in autumn, providing references for water resource management and sustainable development policies in this area to avoid losses caused by water scarcity.

On Snowmelt Computations in the Climate Change Context: Comparison of Single Site and Multisite Downscaling of Extreme Daily Temperature Processes

This study focuses on how multisite downscaling procedures can provide a significant improvement in the accuracy of impact assessment results as compared to single-site downscaling techniques. The study compared the accuracy of snowmelt computations given by the multisite multivariate statistical downscaling (MMSD) method and the single site downscaling (SDSM) procedure. Results of an illustrative application using daily maximum (Tmax) and minimum (Tmin) temperature series at seven stations in southern Quebec-Ontario region in Canada indicated that MMSD method was able to reproduce accurately the observed spatial and temporal variabilities of daily extreme temperature processes over different locations. MMSD approach preserved more accurately the interstation correlations of daily Tmax (0.78 and 0.7) and Tmin (0.94 and 0.84) respectively for calibration and validation period than SDSM (0.08 and 0.11 for Tmax; and 0.02 for Tmin). MMSD was found capable to produce more accurate snowmelt estimations than those given by the SDSM.

Projections of meteorological drought events in the upper Kızılırmak basin under climate change scenarios

Climate change, whose negative impacts are becoming increasingly apparent as a result of human actions, intensifies the drought problems to dangerous levels. The development of local-scale drought projections is crucial to take necessary precautions for potential risks and possible effects of drought. In this study, drought analysis was conducted in the Upper Kızılırmak Basin using the standard precipitation index (SPI) method for the near future (2020–2049), mid-century (2050–2074), and late century (2075–2099). The precipitation data required for the SPI were gathered from the data sets developed for the SSP climate change scenarios of the four chosen global climate models. Precipitation data has been made more convenient for local analysis studies with the statistical downscaling method. Forecasts have been created for the temporal variation and spatial distribution of drought events. The study findings indicate that, under the SSP 2-4.5 scenario, drought-related effects of climate change will decrease until 2100. On the other hand, the number and severity of drought events, as well as the duration of dry periods, will increase until 2100 under the SSP 5-8.5 scenario. According to the SSP 5-8.5 scenario, consisting of the most pessimistic forecasts, moderate drought will last 0–60 months, severe drought will last 0–30 months, and extreme drought will last 0–20 months in different regions of the area in the late century. The spatial distribution of droughts will differ based on the SPI index and climate change scenarios. Comparison of SPI and CZI data showed that both indices are effective in meteorological drought analyses.

Statistical downscaling of GCMs wind speed data for trend analysis of future scenarios: a case study in the Lombardy region

Near-surface wind speed is a key climatic variable, affecting many sectors, such as energy production, air pollution, and natural hazard. Lombardy region of Italy is among the European areas with lowest average wind speed, leading generally to low air quality and wind energy potential. However, it is also one of the most affected area by tornadoes in Italy. Here we investigate possible changes in wind circulation as due to prospective global warming. We analysed wind speed WS under future scenarios (SSP1-2.6 and SSP5-8.5) from six Global Climate Models (GCMs) until 2100, tuned against observed WS data. We employed a statistical downscaling method, namely Stochastic Time Random Cascade (STRC) to correct locally GCMs outputs. Three statistical tests, i.e. Linear Regression, Mann Kendall, Moving Window Average, were carried out to analyse future trends of: annual WS averages, 95th quantile (as an indicator of large WS), and the number of days of calm wind per year (NWC). The proposed STRC algorithm can successfully adjust the mean, standard deviation, and autocorrelation structure of the GCM outputs. No strong trends are found for the future. The chosen variables would all display non-stationarity, and the 95th percentile display a positive trend for most of the stations. Concerning NWC, notable discrepancies among GCMs are seen. The STRC algorithm can be used to successfully adjust GCMs outputs to reflect locally observed data and to then generate credible long-term scenarios for WSs as a tool for decision-making.

Multi-Variable Hard Physical Constraints for Climate Model Downscaling

Global Climate Models (GCMs) are the primary tool to simulate climate evolution and assess the impacts of climate change. However, they often operate at a coarse spatial resolution that limits their accuracy in reproducing local-scale phenomena. Statistical downscaling methods leveraging deep learning offer a solution to this problem by approximating local-scale climate fields from coarse variables, thus enabling regional GCM projections. Typically, climate fields of different variables of interest are downscaled independently, resulting in violations of fundamental physical properties across interconnected variables. This study investigates the scope of this problem and, through an application on temperature, lays the foundation for a framework introducing multi-variable hard constraints that guarantees physical relationships between groups of downscaled climate variables.

A new statistical downscaling approach for short‐term forecasting of summer air temperatures through a fusion of deep learning and spatial interpolation

AbstractReliable early forecasting of extreme summer air temperatures is essential for effectively managing and mitigating the socioeconomic damage caused by thermal disasters. Numerical weather prediction models have become valuable tools for forecasting air temperature; however, they incur high computational costs, resulting in coarse spatial resolution and systematic bias owing to imperfect parametrization. To address these problems, we developed a novel statistical downscaling and bias correction method (named DeU‐Net) for the maximum and minimum air temperature (Tmax and Tmin respectively) forecasts obtained from the Global Data Assimilation and Prediction System with a spatial resolution of 10 to 1.5 km over South Korea through the fusion of deep learning (i.e., U‐Net) and spatial interpolation. In this study, we used a methodology to decompose statistically downscaled Tmax and Tmin forecasts into temporal dynamics over South Korea and spatial fluctuations by pixels. When comparing the proposed DeU‐Net with the dynamical downscaling model (i.e., Local Data Assimilation and Prediction System) and support vector regression‐based statistical downscaling model at the seen and unseen stations for forecasting the next‐day Tmax and Tmin, DeU‐Net showed the highest spatial correlation and the lowest root‐mean‐square error in all cases. In a qualitative evaluation, DeU‐Net successfully produced a detailed spatial distribution most similar to the observations. A further comparison extending the forecast lead time to 7 days indicated that the proposed DeU‐Net is a better downscaling approach than support vector regression, regardless of the forecast lead time. These results demonstrate that bias‐corrected high spatial resolution air temperature forecasts with relatively long forecast lead times in summer can be effectively produced using the proposed model for operational forecasting.

Quantifying impacts of precipitation scenarios projected under climate change on annual probability of rainfall-induced landslides at a specific slope

Because of climate change, the intensity, duration, and frequency of future extreme rainfalls are projected to increase in many regions worldwide according to results from global climate models (GCMs). As rainfall is a major triggering factor for landslides, the frequency and risks of landslides are expected to increase. Consequently, it is crucial to quantify the impacts of climate change on landslide risks at a specific site. This study develops a statistical downscaling method based on the generalized extreme value (GEV) distribution and a probabilistic method for assessing the annual probability of rainfall-induced landslides at a specific site under projected precipitation scenarios from GCMs. The GEV-based statistical downscaling method bridges the gap between site-specific extreme rainfall (e.g., annual maximum 1-day rainfall) projection for a future scenario and GCM outputs with large spatial and temporal scales. A physics-based slope model is also used to properly depict slope failure mechanisms. The results indicate that the projected increase in extreme rainfall due to climate change leads to increases in annual frequency of heavy rainstorms and the annual slope failure probability, which are primarily controlled by the increase in the standard deviation, instead of the mean, of future extreme rainfalls.

High-resolution climate projection dataset based on CMIP6 for Peru and Ecuador: BASD-CMIP6-PE

Here, we present BASD-CMIP6-PE, a high-resolution (1d, 10 km) climate dataset for Peru and Ecuador based on the bias-adjusted and statistically downscaled CMIP6 climate projections of 10 GCMs. This dataset includes both historical simulations (1850–2014) and future projections (2015–2100) for precipitation and minimum, mean, and maximum temperature under three Shared Socioeconomic Pathways (SSP1-2.6, SSP3-7.0, and SSP5-8.5). The BASD-CMIP6-PE climate data were generated using the trend-preserving Bias Adjustment and Statistical Downscaling (BASD) method. The BASD performance was evaluated using observational data and through hydrological modeling across Peruvian and Ecuadorian river basins in the historical period. Results demonstrated that BASD significantly reduced biases between CMIP6-GCM simulations and observational data, enhancing long-term statistical representations, including mean and extreme values, and seasonal patterns. Furthermore, the hydrological evaluation highlighted the appropriateness of adjusted GCM simulations for simulating streamflow, including mean, low, and high flows. These findings underscore the reliability of BASD-CMIP6-PE in assessing regional climate change impacts on agriculture, water resources, and hydrological extremes.

Effects of systematic predictor selection for statistical downscaling of rainfall in Hawai'i

AbstractWhile there has been rapid advancement in the development and application of statistical downscaling methods for climate projection, determining the best predictive large‐scale climate information for the targeted local climate variable remains a challenge. The choice of predictor variables is one of the most influential steps of model development and has the potential to lead to varying results, contributing to the total uncertainty of future projections. Despite being a well‐known problem, predictor selection often does not receive adequate attention and the development of a straightforward and feasible prescreening process is needed to provide guidance to simple (e.g., linear regression) and complex (e.g., machine learning) modeling tools alike. In this project, Akaike information criterion (AIC) and leave one out cross validation are used to evaluate sets of predictor variables common in statistical downscaling models of precipitation (e.g., temperature, geopotential height, moisture transport). Because thousands of predictor sets were found to be competitive in their statistical skill in projecting rainfall, results suggest that an ensemble of predictor sets should be used to account for the resulting variance associated with predictor selection. This ensemble approach is applied to make improvements to the most recent statistical downscaling model developed for rainfall projection in Hawai'i. An ensemble validation study is performed by projecting future wet season rainfall in Hawai'i for scenario RCP4.5 using 17 CMIP5 GCMs with ensembles of highly ranked and poorly ranked predictor sets. Results show significantly less variance and improved agreement in the projected sign of precipitation change for the highly ranked predictor set ensemble. In conclusion, it is recommended that statistical downscaling models are run in an ensemble‐mode with multiple combinations of predictors instead of using a single model, allowing for additional quantification of model uncertainty and a “safety net” for when the true physically best‐suited predictor set is unknown.

Precipitation Projection in Cambodia Using Statistically Downscaled CMIP6 Models

The consequences of climate change are arising in the form of many types of natural disasters, such as flooding, drought, and tropical cyclones. Responding to climate change is a long horizontal run action that requires adaptation and mitigation strategies. Hence, future climate information is essential for developing effective strategies. This study explored the applicability of a statistical downscaling method, Bias-Corrected Spatial Disaggregation (BCSD), in downscaling climate models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and then applied the downscaled data to project the future condition of precipitation pattern and extreme events in Cambodia. We calculated four climate change indicators, namely mean precipitation changes, consecutive dry days (CDD), consecutive wet days (CWD), and maximum one-day precipitation (rx1day) under two shared socioeconomic pathways (SSPs) scenarios, which are SSP245 and SSP585. The results indicated the satisfactory performance of the BCSD method in capturing the spatial feature of orographic precipitation in Cambodia. The analysis of downscaled CMIP6 models shows that the mean precipitation in Cambodia increases during the wet season and slightly decreases in the dry season, and thus, there is a slight increase in annual rainfall. The projection of extreme climate indices shows that the CDD would likely increase under both climate change scenarios, indicating the potential threat of dry spells or drought events in Cambodia. In addition, CWD would likely increase under the SSP245 scenario and strongly decrease in the eastern part of the country under the SSP585 scenario, which inferred that the wet spell would have happened under the moderate scenario of climate change, but it would be the opposite under the SSP585 scenario. Moreover, rx1day would likely increase over most parts of Cambodia, especially under the SSP585 scenario at the end of the century. This can be inferred as a potential threat to extreme rainfall triggering flood events in the country due to climate change.

Integrating GRACE/GRACE Follow-On and Wells Data to Detect Groundwater Storage Recovery at a Small-Scale in Beijing Using Deep Learning

Groundwater depletion is adversely affecting Beijing’s ecology and environment. However, the effective execution of the South-to-North Water Diversion Project’s middle route (SNDWP-MR) is anticipated to mitigate Beijing’s groundwater depletion. Here, we propose a robust hybrid statistical downscaling method aimed at enhancing the capability of the Gravity Recovery and Climate Experiment (GRACE) to detect the small-scale groundwater storage anomaly (GWSA) in Beijing. We used three deep learning (DL) methods to reconstruct the 0.5° × 0.5° terrestrial water storage anomaly (TWSA) between 2004 and 2021. Moreover, multiple processing strategies were used to downscale the GWSA to 0.25° from 2004 to 2021 by integrating wells and GRACE/GRACE follow-on data from the optimal DL model. Additionally, we analyzed the spatiotemporal evolution trends of GW in Beijing before and after the implementation of the SNDWP-MR. The results show that the long short-term memory model delivers optimal performance in the TWSA reconstruction of Beijing, with the correlation coefficient (CC), Nash–Sutcliffe coefficient (NSE), and root mean square error (RMSE) being 0.98, 0.96, and 10.19 mm, respectively. The GWSA before and after downscaling is basically consistent with wells data, but the CC and RMSE of downscaling the GWSA from 2004 to 2021 are improving by 34% and 31%, respectively. Before the SNDWP-MR (2004–2014), the trend of GWSA in Beijing was −17.68 ± 4.46 mm/y, with a human contribution of 69.30%. After SNDWP-MR (2015–2021), GWSA gradually increased by 10.00 mm per year, with the SNDWP-MR accounting for 18.30%. This study delivers a technical innovation reference for dynamically monitoring a small-scale GWSA from GRACE/GRACE-FO data.

A Deep Learning Method for Statistical Downscaling of CLDAS Relative Humidity with Different Sources of Data: Sensitivity Analysis

A Deep Learning Method for Statistical Downscaling of CLDAS Relative Humidity with Different Sources of Data: Sensitivity Analysis

Investigating transformer‐based models for spatial downscaling and correcting biases of near‐surface temperature and wind‐speed forecasts

AbstractHigh‐resolution and accurate prediction of near‐surface weather parameters based on numerical weather prediction (NWP) models is essential for many downstream and real‐world applications. Traditional dynamical or statistical downscaling methods are insufficient to derive high‐resolution data from operational NWP forecasts, making it essential to devise new approaches. In recent years, an increasing number of researchers have explored the implementations of deep learning (DL) based models for spatial downscaling, motivated by the similarity between the super‐resolution (SR) problem in computer vision (CV) and downscaling. Furthermore, while transformer‐based models have become state‐of‐the‐art models for many SR tasks, they are rarely applied for downscaling of weather forecasts or climate projections. This study adapted transformer‐based models such as SwinIR and Uformer to downscale the temperature at 2 m () and wind speed at 10 m () over Eastern Inner Mongolia, encompassing the area from 39.6–46°N latitude and 111.6–118°E longitude. We used high‐resolution forecast (HRES) data from the European Centre for Medium‐range Weather Forecast (ECMWF) with a spatial resolution of 0.1° as the input and gridded observation data from the China Meteorological Administration (CMA) Land Data Assimilation System (CLDAS) at a spatial resolution of 0.01° as the target. Given that the models use observation data rather than a coarse‐grained version of forecast data as the target, they accomplish both bias correction and spatial downscaling. The results demonstrate that the performance of SwinIR and Uformer is superior to that of two convolutional neural network (CNN) based models (UNet and RCAN). Additionally, we introduced a novel module to extract features of varying resolution from the high‐resolution topography data and applied a multiscale feature fusion module to merge features of different scales, contributing to further enhancement of Uformer's performance.

STATISTICAL DOWNSCALING USING REGRESSION NONPARAMETRIC OF FOURIER SERIES-POLYNOMIAL LOCAL OF CLIMATE CHANGE

Indonesia is a tropical country that is vulnerable to the impacts of climate change. Climate change causes an effect on the level of comfort (heat stress) that can affect the level of human immunity, one of the indices to calculate the level of human comfort (heat stress) is the Thermal Humidity Index (THI). Climate change scenarios modeled in Earth System Models (ESMs). ESM has a coarse resolution and is subject to considerable bias. This research is using secondary data. The data source used in this study comes from the Coupled Model Intercomparison Project (CMIP5). This research will focus on projected heat stress which is calculated based on THI with the temperature and humidity variables. Therefore, in this research to reduce the bias correction method used Statistical Downscaling (SD) and nonparametric regression. The results of the bias correction using the Statistical Downscaling (SD) method and Nonparametric Regression Fourier-Polynomial Local Series in this study the R-square value for Relative Humidity yields 95% and for Temperature yields 94%. The projection of climate change based on the value of the Temperature Humidity Index (THI) in Indonesia in the category of 50% of the population of Indonesians feeling comfortable conditions occurred in 2006-2059. Then the population of citizens in Indonesia felt uncomfortable conditions occurred in 2060 to 2100 with a THI value of 27.0730°C - 27.7800°C.

Projection of temperature-related mortality among the elderly under advanced aging and climate change scenario

Advanced aging is expected to become a major social problem in China during the second half of the 21st century. Current projections of temperature-related mortality in the context of advanced aging are inadequate, and may underestimate the risks posed by global warming on people aged 90+ years. The present study addresses this issue in Nantong City, which was the city in China with the highest aging and advanced aging rates in the 2000, 2010, and 2020 population censuses. Based on 27 global climate models from the Coupled Model Intercomparison Project Phase 6 and statistical downscaling methods, the impacts of climate change and advanced aging on future temperature-related elderly mortality were explored. Our results indicate that global warming will continue to increase the proportion of people who die from non-optimal temperatures even without considering the impact of advanced aging. Moreover, a higher warming range led to a more significant increase in net-temperature-related mortality and advanced aging is likely to increase heat-related mortality and offset the decline in cold-related mortality. Our study demonstrates a “1 + 1 > 2” effect between advanced aging and climate change, under the four shared socioeconomic pathway climate change scenarios considered here. These findings contribute to a better understanding of the impact of climate change on elderly health and facilitate the development of more effective adaptive strategies for advanced aging societies.

Projected heat wave increasing trends over China based on combined dynamical and multiple statistical downscaling methods

Extensive investigations on the projection of heat waves (HWs) were conducted on the basis of coarse-resolution global climate models (GCMs). However, these investigations still fail to characterise the future changes in HWs regionally over China. PRECIS dynamical downscaling with a horizontal resolution of 25 km × 25 km was employed on the basis of GCM-HadCM3 to provide reliable projections on HWs over the Chinese mainland, and six statistical downscaling methods were used for bias correction under RCP4.5 and RCP8.5 scenarios. The multi-method ensemble (MME) of the top three dynamical downscaling methods with good performance was used to project future changes. Results showed that PRECIS primarily replicated the detailed spatiotemporal pattern of HWs. However, PRECIS overestimated the HWs in the Northwest and Southeast and expanded the areas of HWs in the Northeast and Southwest. Three statistical downscaling methods (quantile mapping, CDF-t and quantile delta mapping) demonstrated good performance in improving PRECIS simulation for reproducing HWs. By contrast, parametric-based trend-preserving approaches such as scaled distribution mapping and ISI-MIP are outperformed by the three aforementioned methods in downscaling HWs, particularly in the high latitudes of China. Based on MME projections, at the end of the 21st century, the national average of the number of HW days each year, the length of the longest HW event in the year and the extreme maximum temperature in HW will increase by 3 times, 1 time and 1.3 °C, respectively, under the RCP4.5 scenario, whilst that under the RCP8.5 scenario will increase by 8 times, 3 times and 3.7 °C, respectively, relative to 1986–2005. The Northwest is regionally projected to suffer long and hot HWs, whilst the South and Southeast will experience frequent consecutive HWs. Thus, HWs projected by the combined dynamical and statistical downscaling method are highly reliable in projecting HWs over China.

Comparison between statistical and dynamical downscaling of rainfall over the Gwadar‐Ormara basin, Pakistan

AbstractThis paper evaluated and compared the performance of a statistical downscaling method and a dynamical downscaling method to simulate the spatial–temporal rainfall distribution. Outputs from RegCM4 Regional Climate Model (RCM) and the CanESM2 Atmosphere–Ocean General Circulation Model (AOGCM) were selected for the data scarce Gwadar‐Ormara basin, Pakistan. The evaluation was based on the climatological average and standard deviation for historic (1971–2000) and future (2041–2070) time periods under Representative Concentration Pathways (RCP) 4.5 and 8.5 scenarios. The performance evaluation showed that statistical downscaling is preferred to simulate and project rainfall patterns in the study area. Additionally, the Statistical DownScaling Model (SDSM) showed low R2 values in calibration and validation of the simulations with respect to observed data for the historic period. Overall, SDSM generated satisfactory results in simulating the monthly rainfall cycle of the entire basin. In this study, RegCM4 showed large rainfall errors and missed one rainfall season in the historic period. This study also explored whether the grid‐based rainfall time series of the Asian Precipitation—Highly Resolved Observational Daily Integration Towards Evaluation (APHRODITE) dataset could be used to enlarge and complement the sample of in situ observed rainfall time series. A spatial correlogram was used for observed and APHRODITE rainfall data to assess the consistency between the two data sources, which resulted in rejecting APHRODITE data. For the future time period (2041–2070) under RCPs 4.5 and 8.5 scenarios, rainfall projections did not show significant difference for both downscaling approaches. This may relate to the driving model (CanESM2 AOGCM) and not necessarily suggests poor performance of downscaling; either statistical or dynamical. Hence, the study recommends evaluating a multi‐model ensemble including other GCMs and RCMs for the same area of study.

Development of statistical downscaling methods for the assessment of rainfall characteristics under climate change scenarios

Abstract The objective of this research was to develop a statistical downscaling approach in the Phetchaburi River Basin, Thailand, consisting of two main processes: predictor selection and relationship construction between predictors and local rainfall. Super predictor (SP) and stepwise regression (SR) were employed with principal component analysis (PCA) and the statistical downscaling model (SDSM) was applied later. Four statistical models were finally used: M1 (SP–PCA), M2 (SP–SDSM), M3 (SR–PCA), and M4 (SR–SDSM) with 26 large circulation indices generated by CanESM2 (CMIP5). Finally, the characteristics of extreme rainfall events were observed under three climate change scenarios during three different periods (2020–2040, 2041–2070, and 2071–2100). The results revealed that rainfall and geostrophic airflow velocity were the best predictors for rainfall downscaling, followed by divergence, meridional velocity, and relative humidity. Three objective functions (R2, NSE, and RMSE) were applied to evaluate model performance. The M4 model presented the highest performance while M1 showed the lowest skill. The average annual rainfall specifically increased compared with the historical rainfall for RCP2.6, RCP4.5, and RCP8.5 scenarios in future periods. Very to extremely wet years determined by the standardized anomaly index (SAI) occurred more often in the far-future while severe to extremely dry years frequently occurred in the mid- and far-future.

Impact assessment of climate change on energy performance and thermal load of residential buildings in New Zealand

While it is evident that climate change will have an impact on the energy demand for heating and cooling in buildings, the exact extent of this impact is not yet fully understood. Quantification of future cooling and heating need in buildings provides a basis for taking appropriate measures for building climate change adaptation. The focus of this study is to examine how future climate change scenarios will impact the heating and cooling of residential buildings across different climatic regions in New Zealand. The future weather data under changing climate were generated for six climatic zones of New Zealand employing the statistical downscaling method. The study used various climate change scenarios, which represent concentration pathways (RCPs), to generate weather data. Specifically, the RCP8.5 and RCP4.5 scenarios were employed in the building performance simulations for different prototypes of residential buildings. The results showed there would be a significant change in the thermal performance of residential buildings, with a noticeable increase in cooling load and a decrease in heating load. These changes include a maximum thermal load change of 3 kWh/m2 in Auckland by 2090, 2.7 kWh/m2 in Hamilton, 8.3 kWh/m2 in Wellington, 4.2 kWh/m2 in Rotorua, 11 kWh/m2 in Christchurch, and 11.6 kWh/m2 in Queenstown. The warmer climatic zones are expected to change from a heating dominated to a cooling-dominated zone. The results indicated the importance of considering present and future climatic conditions in design and establishing a foundation for actions for the resilience of buildings to climate change.