The removal of non-climatic biases, so-called inhomogeneities, from long climatic records needs sophistically developed statistical methods. One principle is that usually the differences between a candidate series and its neighbour series are analysed instead of directly the candidate series, in order to neutralize the possible impacts of regionally common natural climate variation on the detection of inhomogeneities. In most homogenization methods, two main kinds of time series comparisons are applied, i.e. composite reference series or pairwise comparisons. In composite reference series the inhomogeneities of neighbour series are attenuated by averaging the individual series, and the accuracy of homogenization can be improved by the iterative improvement of composite reference series. By contrast, pairwise comparisons have the advantage that coincidental inhomogeneities affecting several station series in a similar way can be identified with higher certainty than with composite reference series. In addition, homogenization with pairwise comparisons tends to facilitate the most accurate regional trend estimations. A new time series comparison method is presented here, which combines the use of pairwise comparisons and composite reference series in a way that their advantages are unified. This time series comparison method is embedded into the ACMANT homogenization method, and tested in large, commonly available monthly temperature test datasets.

Linear mixed models are frequently used for analysing heterogeneous data in a broad range of applications. The restricted maximum likelihood method is often preferred to estimate co-variance parameters in such models due to its unbiased estimation of the underlying variance parameters. The restricted log-likelihood function involves log determinants of a complicated co-variance matrix. An efficient statistical estimate of the underlying model parameters and quantifying the accuracy of the estimation requires the first derivatives and the second derivatives of the restricted log-likelihood function, i.e., the observed information. Standard approaches to compute the observed information and its expectation, the Fisher information, is computationally prohibitive for linear mixed models with thousands random and fixed effects. Customized algorithms are of highly demand to keep mixed models analysis scalable for increasing high-throughput heterogeneous data sets. In this paper, we explore how to leverage an averaged information splitting technique and dedicate matrix transform to significantly reduce computations and to accelerate computing. Together with a fill-in reducing multi-frontal sparse direct solver, the averaged information splitting approach improves the performance of the computation process.

The Cordoba Guitar Festival is one of the most important cultural events in Spain. This article analyses the musical preferences, satisfaction, attitudinal loyalty and behavioural loyalty of spectators who attended the 36th festival held in July 2016, as well as the festival’s economic impact on the city. To achieve this aim, a structural equation model (SEM) was used. The results show the goodness-of-fit of the model and indicate that the observed data fit the expected dataset.

In soil erosion estimation models, the variable with the greatest impact is rainfall erosivity (RE), which is the measurement of precipitation energy and its potential capacity to cause erosion, and erosivity density (ED), which relates RE to precipitation. The RE requires high temporal resolution records for its estimation. However, due to the limited observed information and the increasing availability of rainfall estimates based on remote sensing, recent research has shown the usefulness of using observed-corrected satellite data for RE estimation. This study evaluates the performance of a new gridded dataset of RE and ED in Peru (PISCO_reed) by merging data from the IMERG v06 product, through a new calibration approach with hourly records of automatic weather stations, during the period of 2000-2020. By using this method, a correlation of 0.7 was found between the PISCO\_reed and RE obtained by the observed data. An average annual RE for Peru of 4831 MJ·mm·ha-1·h-1 was estimated with a general increase towards the lowland Amazon regions and high values are found on the north-coast Pacific area of Peru. The spatial identification of the most risk areas of erosion, was carried out through a relationship between the ED and rainfall. Both erosivity data sets will allow us to expand our fundamental understanding and quantify soil erosion with greater precision.