Climate Time Series Research Articles

New method for statistical analysis of climate time series

After publishing four articles utilizing a new method for the statistical study of climate time series, we found it useful to provide a detailed review of the method itself, which is the primary objective of this work. Unlike the methods most commonly used by scientists analyzing such data, this new method does not seek to identify trends for explorative forecasts. Instead, it enables the detection of precise signals indicating interactions with other climate entities, thereby enhancing our understanding of the underlying phenomena. As illustrated through three example articles, the mechanisms uncovered using this method can be integrated into a mathematical model. The simulations thus obtained are more deterministic than stochastic – a significant advantage for producing high-quality forecasts in the context of global warming. Even if this was the sole application of the method, it would be sufficient to demonstrate its value. However, as a final example detailed in this work shows, reconsidering the original series using different periods (e.g., month, quarter, semester, year) can further refine our understanding of the mechanisms at play. We conclude this work by exploring the potential applicability of this method for analyzing non-climatic temporal data series.

Climate time series variability analysis of Islamabad Capital Territory using fractal dimension and Hurst exponent methods

Climate time series variability analysis of Islamabad Capital Territory using fractal dimension and Hurst exponent methods

Can residential energy systems withstand the heat? Quantifying solar photovoltaic and heat pump yields for future New Zealand climate conditions

Abstract We investigate how higher temperatures resulting from climate change impact the energy system. Specifically, we examine the cumulative effects of fluctuating solar photovoltaic (PV) generation performance, heating and cooling demand, and heat pump efficiency on such days.
To achieve this, we used the climate analogue space, which maps a given city’s future climate to an existing one. By employing climate analogues, we can predict the impact of higher temperatures by 2050, transforming Auckland, New Zealand’s climate into one akin to Sydney, Australia. This approach avoids reliance on historical weather data, which many energy system models use. We used this future climate time series as an input to a residential energy system model for Auckland, New Zealand. The residential energy system model simulates solar PV generation output via mapping of experimental data, building thermal characteristics via grey-box resistance-capacitance (RC) modelling, and hourly coefficient of performance (COP) for air source heat pumps (ASHP) via linear regression.
Our findings revealed that a future climate doubles the cooling demand and reduces the heating demand by one-third, with the heat pump demand peak load projected to be 40% higher than current demand. Although solar PV generation experiences a decrease in efficiency of 8%, there is a 40% increase in annual direct usage of ASHP. Despite the high cooling demand, the combined yearly electricity demand for heating and cooling decreased by 6.5% overall, and the system saw a 50% improvement in demand fulfilment. However, the system performance volatility at hotter-than-normal temperatures and the potential for significant energy shortfalls remain concerns. The shift from a predominantly heating to a cooling environment is a critical design condition that should be considered in energy expansion planning and future electrification.
The framework and time series developed in this work can be expanded and applied to other energy system modelling exercises.

Revitalizing temperature records: A novel framework towards continuous data reconstruction using univariate and multivariate imputation techniques

Revitalizing temperature records: A novel framework towards continuous data reconstruction using univariate and multivariate imputation techniques

A Novel Approach for Temperature Forecasting in Climate Change Using Ensemble Decomposition of Time Series

AbstractThis paper presents FORSEER (Forecasting by Selective Ensemble Estimation and Reconstruction), a novel methodology designed to address temperature forecasting under the challenges inherent to climate change. FORSEER integrates decomposition, forecasting, and ensemble methods within a modular framework. This methodology decomposes the time series into trend, seasonal, and residual components. Subsequently, multiple optimized forecast models are applied to each component. These component models are then carefully weighted and combined through an ensemble process to generate a final robust forecast. Experimental results demonstrate that FORSEER is an efficient computational forecasting methodology for complex climate time series. Furthermore, we show that FORSEER has an equivalent forecasting performance to the M4 competition champion SMYL method for temperature series. Besides, the proposed methodology has less computational complexity than SMYL, making it a more accessible and scalable option. FORSEER's modular architecture also allows flexibility when substituting techniques depending on the context of the problem, facilitating the parallel execution of independent tasks and resulting in a strategy adaptable to multiple contexts.

Univariate and multivariate imputation methods evaluation for reconstructing climate time series data: A case study of Mosul station-Iraq

Comprehensive climate time series data is indispensable for monitoring the impacts of climate change. However, observational datasets often suffer from data gaps within their time series, necessitating imputation to ensure dataset integrity for further analysis. This study evaluated six univariate and multivariate imputation methods to infill missing values. These methods were applied to complete the subsets of time series data for precipitation, temperature, and relative humidity from Mosul station spanning 1980–2013. Artificial gaps of 5%, 10%, 20%, and 30% missing observations were introduced under scenarios of missing completely at random (MCAR) missing at random (MAR), and missing not at random (MNAR). Evaluation metrics including RMSE and Kling-Gupta Efficiency were utilized for performance evaluation. Results revealed that seasonal decomposition was the most effective univariate imputation method across all variables. For the multivariate imputation, kNN demonstrated superior performance in infilling the precipitation missing data under MCAR, while norm.predict exhibited optimal performance in the temperature missing data under all missing scenarios. Moreover, missForest was identified as the most suitable method for infilling missing relative humidity data. This study's methodology offers insights into selecting appropriate imputation methods for other climate stations, thereby enhancing the accuracy of the climate change effects analysis.

ANOVA (Benova) correction in relative homogenization: Why it is indispensable

AbstractThis paper reviews the role of ANOVA correction model in the homogenization of climatic time series. In the present context ANOVA has only weak connection to its original meaning (analysis of variance), so we propose the new name “Benova” to replace the confusing old name. In the linear model of Benova corrections (hereafter Benova) the information of statistically detected inhomogeneities and metadata are jointly considered for all time series of a given climatic region. Benova has indisputable advantages on the accuracy of homogenization results, and this has both theoretical and practical evidence. The study presents two principal versions of Benova: in simple Benova the climate signal is presumed to be spatially invariant, while in weighted Benova the spatial variation of climate is considered. In Benova models usually only breaks (i.e., sudden shifts of section mean) are considered, but this restriction has practical reasons, rather than theoretical limits, and the study shows an extended version of the method with which trend‐like inhomogeneity biases can also be removed. Benova can be used for a group of time series covering varied time periods, the operations can be performed in any time resolution, and statistical characteristics others than climatic means can also be homogenized by the method. Benova can be used together with any break detection method. The study discusses the likely reasons of the relatively slow spread in practical application. PRODIGE was the first homogenization method which used Benova corrections. Until now, all homogenization methods including Benova corrections include also the same kind break detection method, i.e., penalized maximum likelihood method with step function fitting. The brief descriptions of two modern methods of this method family, that is, ACMANT and Bart methods, are also provided.

Relative Homogenization of Climatic Time Series

Homogenization of the time series of observed climatic data aims to remove non-climatic biases caused by technical changes during the history of the climate observations. The spatial redundancy of climate information helps to recognize station-specific inhomogeneities with statistical methods, but the correct detection and removal of inhomogeneity biases is generally not easy for the combined effects of individual inhomogeneities. In a homogenization procedure, several time series of a given climatic variable observed in one climatic region are usually homogenized together via a large number of spatial comparisons between them. Such procedures are called relative homogenization. A relative homogenization procedure may include one or more homogenization cycles where a cycle includes the steps of time series comparison, inhomogeneity detection and corrections for inhomogeneities, and they may include other steps like the filtering of outlier values or spatial interpolations for infilling data gaps. Relative homogenization methods differ according to the number and content of the individual homogenization cycles, the procedure for the time series comparisons, the statistical inhomogeneity detection method, the way of the inhomogeneity bias removal, among other specifics. Efficient homogenization needs the use of tested statistical methods to be included in partly or fully automated homogenization procedures. Due to the large number and high variety of homogenization experiments fulfilled in the Spanish MULTITEST project (2015–2017), its method comparison test results are still the most informative about the efficiencies of homogenization methods in use. This study presents a brief review of the advances in relative homogenization, recalls some key results of the MULTITEST project, and analyzes some theoretical aspects of successful homogenization.

The Hawai‘i Climate Data Portal (HCDP)

Abstract The Hawai‘i Climate Data Portal (HCDP) is designed to facilitate streamlined access to a wide variety of climate data and information for the State of Hawai‘i. Prior to the development of the HCDP, gridded climate products and point datasets were fragmented, outdated, not easily accessible, and not available in near–real time. To address these limitations, HCDP researchers developed the cyberinfrastructure necessary to 1) operationalize data acquisition and product production in a near-real-time environment and 2) make data and products easily accessible to a wide range of users. The HCDP hosts several high-resolution (250 m) gridded products including monthly rainfall and daily temperature (maximum, minimum, and mean), station data, and gridded future projections of rainfall and temperature. HCDP users can visualize both gridded and point data, create and download custom maps, and query station and gridded data for export with relative ease. The “virtual station” feature allows users to create a climate time series at any grid point. The primary objective of the HCDP is to promote sharing and access to data and information to streamline research activities, improve awareness, and promote the development of tools and resources that can help to build adaptive capacities. The HCDP products have the potential to serve a wide range of users including researchers, resource managers, city planners, engineers, teachers, students, civil society organizations, and the broader community.

Machine learning based on functional principal component analysis to quantify the effects of the main drivers of wheat yields

Machine learning based on functional principal component analysis to quantify the effects of the main drivers of wheat yields

Modeling deficit irrigation water demand of maize and potato in Eastern Germany using ERA5-Land reanalysis climate time series

Modeling deficit irrigation water demand of maize and potato in Eastern Germany using ERA5-Land reanalysis climate time series

Effects of major storms over a spit-tidal flat interaction system (Punta Rasa-Samborombón, Argentina)

Sandy coastal areas are very dynamic systems in which morphological changes occur over different time scales from hours to decades. However, it has been widely reported that major storms are the main responsible of the most significant changes in short to medium time scales. Major storms have been defined using a variety of environmental variables, but they are normally associated with high values of wave height, duration, return period and direction. Here, we aim to characterize types of major storms and to categorize associated morphological impacts over a complex coastal system. The study site, known as Punta Rasa, is located in the Samborombón bay in the outer part of the Río de La Plata estuary (Argentina) and corresponds to a zone of interaction between a large sandy spit and a backwash tidal flat system. Methods combine statistics of wave climate time-series, analysis of wave energy using nearshore numerical modelling (SWAN) and comparison of pre- and post-storm morphological changes by means of shoreline change detection and satellite images derived indexes (CoastSat Toolkit and NDWI index respectively). Results allowed to characterize four types of major storms impacting the study area: High-Energy Storms (HES), defined by an average storm wave below the 1% exceedance (>2.6m), Long-Lived Storms (LLS) represented by an exceedance of the 1% of Du (>60h), Storm Groups (SG) in which storm return period is <6days and Northeastern moderate storms (NMS) defined by their eastern, onshore oriented direction. Under HES and NMS storms erosional areas are dominant over depositional, causing shoreline retreat, a growth of the end-spit and the increase on sand deposition on the back-barrier areas. Under LLS and SG storms, the morphological impact varies alongshore, with multiple erosional hotspots found along the shoreline accompanied by a general flattens of the end-spit system.

Developing spatially explicit and stochastic measures of ecological departure

Background Ecological departure is a metric applied to mapped ecological systems measuring dissimilarity between the distributions of observed and expected proportions of non-stochastic reference vegetation classes within an area. Aims We created spatially explicit measures of ecological departure incorporating stochasticity for each ecological system and all ecological systems from a central Nevada, USA, landscape. Methods Spatially explicit ecological departures were estimated from a radius from each pixel governed by a distance-decay function within a moving window. Variability was introduced by simulating replicate climate time series for each spatial reference condition and calculating departure per replicate. Key results Single-system spatial ecological departure was high and extensive, except for one area of low-elevation groundwater-dependent systems. Variance of spatial ecological departure was extensively low, except in areas of lower ecological departure, despite vegetation differences among replicates. The multiple-system ecological departure exhibited lower values. Conclusions Spatial ecological departure is warranted for efficient land management as results were concordant between non-spatial and spatial metrics; however, rapid coding languages will be required. Implications Spatially explicit ecological departure of both single and multiple systems facilitate localised vegetation and wildlife habitat management and land protection decisions.

ClimateDT: A Global Scale-Free Dynamic Downscaling Portal for Historic and Future Climate Data

Statistical downscaling of climate data has been widely described in the literature, with the aim of improving the reliability of local climatic parameters from coarse-resolution (often &gt;20 km) global datasets. In this article, we present ClimateDT, a dynamic downscaling web tool for monthly historical and future time series at a global scale. The core of ClimateDT is the 1 km 1981–2010 climatology from CHELSA Climate (version 2.1), where the CRU-TS layers for the period 1901-current are overlayed to generate a historic time series. ClimateDT also provides future scenarios from CMIP5 using UKCP18 projections (rcp2.6 and rcp8.5) and CMIP6 using 5 GCMs, also available on the CHELSA website. The system can downscale the grids using a dynamic approach (scale-free) by computing a local environmental lapse rate for each location as an adjustment for spatial interpolation. Local predictions of temperature and precipitation obtained by ClimateDT were compared with climate time series assembled from 12,000 meteorological stations, and the Mean Absolute Error (MAE) and the explained variance (R2) were used as indicators of performance. The average MAEs for monthly values on the whole temporal scale (1901–2022) were around 1.26 °C for the maximum monthly temperature, 0.80 °C for the average monthly temperature, and 1.32 °C for the minimum monthly temperature. Regarding monthly total precipitation, the average MAE was 19 mm. As for the proportion of variance explained, average R2 values were always greater than 0.95 for temperatures and around 0.70 for precipitation due to the different degrees of temporal autocorrelation of precipitation data across time and space, which makes the estimation more complex. The elevation adjustment resulted in very accurate estimates in mountainous regions and areas with complex topography and substantially improved the local climatic parameter estimations in the downscaling process. Since its first release in November 2022, more than 1300 submissions have been processed. It takes less than 2 min to calculate 45 locations and around 8 min for the full dataset (512 records).

A statistical method for the attribution of change‐points in segmented Integrated Water Vapor difference time series

AbstractMany segmentation or change‐point detection methods for homogenizing climate time series compare candidate station data with reference data to eliminate common climate signals and more efficiently detect spurious, non‐climatic changes. One drawback is that it is difficult to decide whether the detected change‐point is due to the candidate series or to the reference. A so‐called attribution procedure is typically applied in a post‐processing step for each detected change‐point. This article describes a new statistical method for the attribution of change‐points detected in Global Navigation Satellite System (GNSS) minus reanalysis series of integrated water vapour. It requires at least one nearby station with similar GNSS and reanalysis data. Six series of differences are formed from the four base series (BS) and are tested for a significant jump at the time of the change‐point detected in the candidate station. The six test results are analysed with a statistical predictive rule to attribute the change‐point to one, or several, of the four BS. Original aspects of our method are: (1) the significance test, which is based on a generalized linear regression approach, taking both heteroscedasticity and autocorrelation into account; (2) the predictive rule, which uses a machine learning method and is constructed from the test results obtained with the real data by using a resampling strategy. Four popular machine learning methods have been compared using cross‐validation and the best one was applied to a real data set (49 main stations with 114 change‐points). The results depend on the choice of the test significance level and the aggregation method combining the prediction results when several nearby stations are available. We find that 62% of the change‐points are attributed to GNSS, 19% to the reanalysis, and 10% are due to coincident detections.

Improved estimation of semiparametric dynamic copula models with filtered nonstationarity

Improved estimation of semiparametric dynamic copula models with filtered nonstationarity

Hazard Analysis and Vulnerability Assessment of Cultural Landscapes Exposed to Climate Change-Related Extreme Events: A Case Study of Wachau (Austria)

The present paper aims to study the Wachau Valley in Austria as a representative Cultural Landscape under threat from extreme hydrometeorological hazards linked to climate change. The primary objective is to investigate the impacts and assess the vulnerability associated with the events of heavy rain and flooding. The methodology employed consists of an investigation of recorded past events impacting the Wachau; a vulnerability ranking system; a climate time series analysis based on earth observation products; and future hazard maps at territorial level, developed with outputs from regional and global climate models. The investigation we carried out provides a vulnerability assessment of two terraced areas with a surface of about 10,000 m2 in total, characterized by the presence of dry stone walls, with different state of conservation in the Municipality of Krems (Wachau). In addition, climate projections at territorial level for the extreme climate indices R20mm, R95pTOT, and R×5day—selected for investigating the likelihood of increases/decreases in events of heavy rain and large basin flooding—are provided, with a spatial resolution of ~12 km for the near and far future (2021–2050; 2071–2100) under stabilizing (RCP 4.5) and pessimistic (RCP 8.5) scenarios. The results indicate a general increase for the three indices in the studied areas during the far future under the pessimistic scenario, suggesting a heightened risk of heavy rain and flooding. These findings aim to inform policymakers and decision-makers in their development of strategies for safeguarding cultural heritage. Furthermore, they serve to assist local stakeholders in enhancing their understanding of prioritizing interventions related to preparedness, emergency response, and recovery.

Analysis of Pros and Cons in Using the Water–Energy–Food Nexus Approach to Assess Resource Security: A Review

The water–energy–food (WEF) nexus is drawing much attention in scholarly literature as a novel alternative to address complex resources and achieve resource security. The aim of this study is to analyze and review existing nexus studies to investigate the current status of nexus research worldwide. This study used a narrative review approach to provide a comprehensive overview on the WEF nexus using a variety of databases. It is indicated that the majority of studies in Asia and Africa focused on the water–energy–food (WEF) nexus. China and Brazil had the largest nexus research. Based on the existing literature, most of attention has been paid to food production. However, food consumption patterns and dietary change are rarely evaluated, and there is a lack of study on impacts of dietary change on the WEF nexus. Moreover, there is a lack of frameworks for the evaluation of the WEF nexus under dietary change scenarios. The major challenge of the nexus approach is data availability in crop production, which can be solved by using remote sensing data. There is a lack of standard and conceptual frameworks for nexus assessment and, then, an essential need to provide a new holistic and standard approach that be applicable worldwide to increase connections between researchers and decision makers, as well as the applicability of nexus approaches. Future research must couple the development of a holistic standard approach with experimental tests in different areas, involving interdisciplinary research groups able to carry out all the experimental activities, the numerical simulations, and the statistical analyses of climatic time series (in a climate change perspective) indispensable to demonstrate the real benefits of using a WEF-derived nexus approach.

Ground deformation monitoring via PS-InSAR time series: An industrial zone in Sacco River Valley, central Italy

Ground deformation monitoring via PS-InSAR time series: An industrial zone in Sacco River Valley, central Italy

Agricultural Yield Responses to Climate Variabilities in West Africa: A Food Supply and Demand Analysis

Agricultural productivity is expected to decrease under changing climate conditions that correspond to the stability of West African food systems. Although numerous studies have evaluated impacts of climate variability on crop yields, many uncertainties are still associated with climate extremes as well as the rapid population growth and corresponding dietary lifestyle. Here, we present a food supply and demand analysis based on the relationship between climate change, crop production, and population growth in three sites from southwestern Burkina Faso to southwestern Ghana. Climate and agricultural time series were analyzed by using boxplots mixed with a Mann–Kendall trend test and Sen’s slope. Food balance sheets were calculated by estimating the demand using a population growth model linked to food supply with local consumption patterns. We found almost insignificant rainfall and temperature trends for both sites in the Sudano-Guinean savannah. Conversely, the climate regime of southwestern Ghana revealed a strong significant increasing temperature over time. Crop yield trends demonstrated that maize and sorghum were significantly enhanced in both study areas of the Sudano-Guinean savannah. Southwestern Ghana depicted a different crop pattern where cassava and plantain showed a strong upward yield trend. The grouped food balance sheets across the regions illustrated a surplus for the Sudano-Guinean savannah while southwestern Ghana exhibited a deficit. Despite the growing yield of various crops, food demand is outpacing regional production.