This work analyses three uncertainty sources affecting the observation‐based gridded data sets: station density, interpolation methodology and spatial resolution. For this purpose, we consider ...precipitation in two countries, Poland and Spain, three resolutions (0.11, 0.22 and 0.44°), three interpolation methods, both areal‐ and point‐representative implementations, and three different densities of the underlying station network (high/medium/low density). As a result, for each resolution and interpolation approach, nine different grids have been obtained for each country and inter‐compared using a variance decomposition methodology.
Results indicate larger differences among the data sets for Spain than for Poland, mainly due to the larger spatial variability and complex orography of the former region. The variance decomposition points out to station density as the most influential factor, independent of the season, the areal‐ or point‐representative implementation and the country considered, and slightly increasing with the spatial resolution. In contrast, the decomposition is stable when extreme precipitation indices are considered, in particular for the 50‐year return value.
Finally, the uncertainty due to station sub‐sampling inside a particular grid box decreases with the number of stations used in the averaging/interpolation. In the case of spatially homogeneous grid boxes, the interpolation approach obtains similar results for all the parameters, excepting the wet day frequency, independently of the number of stations. When there is a more significant internal variability in the grid box, the interpolation is more sensitive to the number of stations, pointing out to a minimum stations’ density for the target resolution (six to seven stations).
The uncertainty due to the stations density, interpolation method and resolution has been analysed for the precipitation in two countries, Poland and Spain, by means a variance decomposition analysis. The results reflect that the main factor in the development of gridded data sets is the underlying station's density for both mean and extreme precipitation. In addition, a minimum density of six to seven stations per grid box has been identified to reach an effective resolution of 0.44.
•Prediction accuracy of clay was more affected by sampling strategies.•Vis-NIR spectra had higher prediction performance compared to pXRF spectra.•Combined Vis-NIR and pXRF spectra had no improvement ...on prediction accuracy.
The recent technologies employed for rapid, cost-effective, and non-destructive prediction of soil particle size distribution (clay, sand, and silt) are becoming increasingly interesting among soil scientists. Our aims were to explore the effect of surface, profile wall, and surface + profile wall on prediction accuracy using individual and combined both soil spectra (Vis-NIR and pXRF) with machine learning algorithms for sand, silt, and clay. In total, 191 soil samples were collected from the soil surface (0–30 cm) and profile wall (1 m × 1 m) from cultivated fields in Eskisehir, Central Anatolia of Türkiye. The pXRF (0–45 keV) and Vis-NIR (350–2500 nm) spectroradiometers were used to obtain soil spectra from sieved soil samples. The prediction accuracy of each soil particle size was evaluated by 54 models to explore the predictive performance. The five machine learning algorithms (elastic net, lasso, random forest, ridge, and support vector machine-linear) were applied with calibration (70% soil samples) and validation (30% soil samples) data set for each soil particle size.Results showed the dominant clay mineral in the A and C horizons is chlorite. Moderate and high prediction accuracy for sand (R2 = 0.56–0.84) and clay (R2 = 0.61–0.80), whereas only moderate prediction accuracy for silt (R2 = 0.47–0.55) using both soil spectra in the surface, profile wall, and surface + profile wall. The highest prediction accuracy for each soil particle size was achieved in the soil profile wall using Vis-NIR spectra with elastic net, which outperformed other samplings such as individual pXRF, combined both soil spectra, and other machine learning algorithms. In addition, the prediction accuracy of clay was more affected by sampling strategies compared to sand and silt. We concluded that individual Vis-NIR spectroradiometer can be utilized to achieve the highest prediction accuracy for sand, silt, and clay ratio in semiarid ecosystems for soil surveys and land use studies.
A erosão é um grande problema ambiental e um dos principais fatores de degradação do solo, reduzindo o potencial agrícola das áreas e gerando o assoreamento de canais fluviais. Para entender o ...comportamento do solo e estabelecer práticas conservacionistas para controle da erosão é muito utilizado no modelo de estimação de perda de solo a Equação Universal da Perda de Solo (EUPS). Outro índice utilizado no estudo da erosão é o Potencial Natural de Erosão (PNE), que é baseado na EUPS, o qual, aliado a geoestatística, é capaz de produzir bons resultados. Dessa forma, buscou-se com este trabalho avaliar a distribuição espacial do potencial natural de erosão em diferentes áreas de produção de café, utilizando análises geoestatísticas e técnicas de interpolação de dados. O estudo foi realizado para quatro áreas de produção de café localizadas em Araponga - MG. O PNE foi determinado através da EUPS suplantada das variáveis referentes ao fator uso e manejo do solo (C) e fator prática conservacionista (P). Foram testadas duas metodologias, a álgebra de mapas e a interpolação do PNE, e analisou-se o erro do mapa final para os dois métodos. O fator erodibilidade (K) apresentou dependência espacial para as áreas 01, 03 e 04, com índice de dependência espacial de moderado a forte, enquanto que o PNE apresentou dependência espacial para todas as áreas de estudo, com índice de dependência espacial forte. Todas as áreas tiveram PNE classificado como muito baixo, com exceção da área 01, que apresentou maiores valores de PNE devido a maior declividade e erodibilidade encontrada na parte mais alta da área. Independentemente da área, os valores mais baixos de PNE, estão associados principalmente, a distribuição espacial das frações granulométricas, que condicionam o fator K. A forma que é calculado o mapa final pode resultar em diferentes resultados, acumulando mais ou menos erros. Dentre os dois métodos utilizados o método 2 proporcionou melhores resultados, atingindo menores erros e permitindo um processamento mais rápido.
High-spatial-resolution and long-term climate data are
highly desirable for understanding climate-related natural processes. China
covers a large area with a low density of weather stations in some ...(e.g.,
mountainous) regions. This study describes a 0.5′ (∼ 1 km)
dataset of monthly air temperatures at 2 m (minimum, maximum, and mean proxy monthly temperatures, TMPs)
and precipitation (PRE) for China in the period of 1901–2017. The dataset
was spatially downscaled from the 30′ Climatic Research Unit (CRU) time
series dataset with the climatology dataset of WorldClim using delta spatial
downscaling and evaluated using observations collected in 1951–2016 by 496
weather stations across China. Prior to downscaling, we evaluated the
performances of the WorldClim data with different spatial resolutions and
the 30′ original CRU dataset using the observations, revealing that their
qualities were overall satisfactory. Specifically, WorldClim data exhibited
better performance at higher spatial resolution, while the 30′ original CRU
dataset had low biases and high performances. Bicubic, bilinear, and
nearest-neighbor interpolation methods employed in downscaling processes
were compared, and bilinear interpolation was found to exhibit the best
performance to generate the downscaled dataset. Compared with the
evaluations of the 30′ original CRU dataset, the mean absolute error of the new dataset (i.e., of the 0.5′ dataset downscaled by bilinear interpolation) decreased by 35.4 %–48.7 % for TMPs and by 25.7 % for PRE. The root-mean-square error decreased by 32.4 %–44.9 % for TMPs and by 25.8 % for PRE. The Nash–Sutcliffe efficiency coefficients increased by
9.6 %–13.8 % for TMPs and by 31.6 % for PRE, and correlation
coefficients increased by 0.2 %–0.4 % for TMPs and by 5.0 % for PRE. The new dataset could provide detailed climatology data and annual trends of all climatic variables across China, and the results could be evaluated well using observations at the station. Although the new dataset was not evaluated before 1950 owing to data unavailability, the quality of the new
dataset in the period of 1901–2017 depended on the quality of the original
CRU and WorldClim datasets. Therefore, the new dataset was reliable, as the
downscaling procedure further improved the quality and spatial resolution of
the CRU dataset and was concluded to be useful for investigations related
to climate change across China. The dataset presented in this article has
been published in the Network Common Data Form (NetCDF) at
https://doi.org/10.5281/zenodo.3114194 for precipitation (Peng,
2019a) and https://doi.org/10.5281/zenodo.3185722 for air temperatures at 2 m
(Peng, 2019b) and includes 156 NetCDF files compressed in zip
format and one user guidance text file.
This paper describes the construction of a 0.25° × 0.25° daily soil temperature (TS) data set at six soil depths (i.e., 0, 5, 10, 15, 20, and 40 cm) based on a homogenized TS at more than 2,000 ...meteorological stations across mainland of China for 1960–2017. First, at each station, the observed TS is split into a daily climatology of 1981–2010 and a daily anomaly for 1960–2017. Then, they are respectively interpolated to 0.25° × 0.25° horizontal resolution by using the thin‐plain spline and the angular distance weight interpolation methods. Finally, a gridded multilayer TS data set is constructed from the sum of the above two gridded products. Latitude, longitude, and high‐resolution digital elevation models are explicitly incorporated into the interpolation processes. Selectivity tests indicate that the optimal interpolation distance for the TS anomaly field is 500 km, where the gridded data set may have maximum coverage of the entire Chinese mainland. The results show that the gridded TS climatology is slightly lower than the station observations. The differences between gridded and the station‐based TS varies ±3°C, with relatively large in SW and TIBET. For the gridded anomaly fields, the cross‐validation shows that the mean absolute errors between gridded and station‐based TS are generally less than 0.6°C in the YZ, SW, and SC, but they are above 1.0°C in the NW and TIBET due to the sparse station distribution and large topographic relief. From 1961 to 2016, the seasonal variability of the gridded TS is generally consistent with that of the surface air temperature. The gridded TS over the majority of land displays significantly positive trends in both summer and winter (p > 95%), indicating an overall warming land surface in China. The long‐term gridded TS data not only preserves the spatial features represented in the station observations, but also provides a spatio‐temporal continuous long‐term product of mainland China. Thus, it may be used in future studies, for example, evaluation of land surface model simulations, validation of satellite retrievals, land‐atmosphere interactions, and climate change.
The construction process of a 0.25° × 0.25° daily soil temperature data set at six soil depths (i.e., 0, 5, 10, 15, 20, and 40 cm) across mainland China for the period of 1960–2017. The data sets are based entirely on the in‐situ measurements at meteorological stations, and have been strictly homogenized. The thin‐plate spline interpolation and angular distance weight interpolation methods are respectively used to interpolate the daily climatology and its anomaly of observation from station locations to the rectangle grids.
Abstract
Most massive stars are members of a binary or a higher-order stellar system, where the presence of a binary companion can decisively alter their evolution via binary interactions. ...Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of compact-object binaries and to decipher the physical processes that lead to their formation. Here, we present
POSYDON
, a novel, publicly available, binary population synthesis code that incorporates full stellar structure and binary-evolution modeling, using the
MESA
code, throughout the whole evolution of the binaries. The use of
POSYDON
enables the self-consistent treatment of physical processes in stellar and binary evolution, including: realistic mass-transfer calculations and assessment of stability, internal angular-momentum transport and tides, stellar core sizes, mass-transfer rates, and orbital periods. This paper describes the detailed methodology and implementation of
POSYDON
, including the assumed physics of stellar and binary evolution, the extensive grids of detailed single- and binary-star models, the postprocessing, classification, and interpolation methods we developed for use with the grids, and the treatment of evolutionary phases that are not based on precalculated grids. The first version of
POSYDON
targets binaries with massive primary stars (potential progenitors of neutron stars or black holes) at solar metallicity.
A dimension reduction method called discrete empirical interpolation is proposed and shown to dramatically reduce the computational complexity of the popular proper orthogonal decomposition (POD) ...method for constructing reduced-order models for time dependent and/or parametrized nonlinear partial differential equations (PDEs). In the presence of a general nonlinearity, the standard POD-Galerkin technique reduces dimension in the sense that far fewer variables are present, but the complexity of evaluating the nonlinear term remains that of the original problem. The original empirical interpolation method (EIM) is a modification of POD that reduces the complexity of evaluating the nonlinear term of the reduced model to a cost proportional to the number of reduced variables obtained by POD. We propose a discrete empirical interpolation method (DEIM), a variant that is suitable for reducing the dimension of systems of ordinary differential equations (ODEs) of a certain type. As presented here, it is applicable to ODEs arising from finite difference discretization of time dependent PDEs and/or parametrically dependent steady state problems. However, the approach extends to arbitrary systems of nonlinear ODEs with minor modification. Our contribution is a greatly simplified description of the EIM in a finite-dimensional setting that possesses an error bound on the quality of approximation. An application of DEIM to a finite difference discretization of the one-dimensional FitzHugh-Nagumo equations is shown to reduce the dimension from 1024 to order 5 variables with negligible error over a long-time integration that fully captures nonlinear limit cycle behavior. We also demonstrate applicability in higher spatial dimensions with similar state space dimension reduction and accuracy results. PUBLICATION ABSTRACT
Celotno besedilo
Dostopno za:
CEKLJ, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Precipitation data availability plays a crucial role in many climatic, hydrological and agricultural-related applications. In this study, the use of alternative data sources (i.e. interpolation ...methods and ERA5-Land reanalysis data) was combined for improving the spatially distributed precipitation estimates at the Simeto river basin, located in Eastern Sicily (Italy). A total of 51 rain gauges were used to generate a monthly precipitation dataset for the reference period 2002–2019. Among the 6 tested interpolation methods, Natural Neighbour was the method that predicted precipitation the best at monthly level with a Distance between Indices of Simulation and Observation (DISO) of 0.51. Radial Basis Functions and Inverse Distance Weighting provided the highest precipitation accuracies, respectively, for winter and autumn (with DISO values of 0.44 and 0.50, respectively), and for spring and summer seasons (with DISO values of 0.50 and 0.63, respectively). Underestimations on the ERA5-Land precipitation estimates were observed when compared to the most accurate interpolation methods both at monthly (25%) and seasonal temporal scales (21% in winter and summer, 36% in autumn), with the exception for spring. The performance was significantly improved when the interpolation estimates were corrected with local observations (with RMSD values ranging from 35.29 mm to 26.46 mm at monthly scale, and from 23.33–55.34 mm to 23.15–37.88 mm at seasonal level). The spatial distribution of the estimation errors associated to precipitation obtained from ERA5-Land reanalysis revealed a significant positive correlation (p value <0.05) with the altitude variation in each ERA5-Land cell within the basin under study. These results confirm the good performance on the combined use of alternative precipitation data sources, while adjustments are required to reduce site-specific uncertainties due to local microclimatic conditions occurring at the basin scale.
•Monthly precipitation was estimated most accurately by deterministic interpolators.•ERA5-Land underestimated precipitation obtained by deterministic interpolators.•BIAS correction procedure improved the accuracy of the precipitation estimates.•Topography heterogeneity worsened the accuracy of ERA5-Land precipitation estimates.
AbstractIn the quest for a best approach to estimate the spatial distribution of monthly precipitation for the state of Pernambuco, Brazil, this work tests different versions (choices of parameters ...and functions) of seven interpolation methods, totaling 26 distinct interpolation schemes. The performance of each method in estimating monthly precipitation was evaluated through a set of seven measures of performance evaluation through a cross-validation procedure. It was found that the trend surface analysis yielded the best overall interpolation results, followed by natural neighbor method, inverse distance weighting, and kriging. Performance ranking of all 26 methods is provided. Moreover, in order to assess the regional performance quality of the best identified method (trend surface analysis) in estimating monthly precipitation, a spatially explicit error analysis was also performed for each of the evaluation measures. Finally, temporal evolution of the error evaluation measures of the best performing methods is presented in order to demonstrate the seasonal impact of interpolation.
Soil total phosphorus (TP) is an essential indicator to reflect the soil fertility in agricultural ecosystems. The accurate prediction of the spatial heterogeneity of TP is crucial to evaluate the ...soil productivity and quality. In this study, the interpolation methods of inverse distance weighted (IDW), radial basis functions (RBF), ordinary kriging (OK), co-kriging (COK), multiple linear regression (MLR), geographically weighted regression (GWR), regression kriging (RK), and geographically weighted regression kriging (GWRK) were used to estimate the spatial patterns of TP in four Mollisol areas with different kinds of landscapes at diverse scales. The calculation method, accuracy (mean error, mean absolute error, root mean square error and relative improvement), cost (money, labour and time), and process by which the models were established were also compared in order to ascertain the best interpolation methods to analyse and predict the spatial heterogeneity of TP in Mollisol areas. The results showed that i) generally, compared with OK, the RK and GWRK methods typically increased the simulative accuracy of the spatial distribution of TP, while IDW, COK and MLR decreased the accuracy, and RBF and GWR were not consistent at improving the accuracy in the Mollisol area; ii) the slope steepness and brightness index could be introduced to regression models as superior auxiliary variables to improve the interpolation precision (0.8%–1.9%) of TP by RK and GWRK when the study area is relatively flat; iii) the GWRK and RK methods with smaller root mean square error and higher relative improvement outperformed other methods in the region with suitable sampling evenness (>45%), while the RBF method is an optional approach when the sampling evenness is low. In summary, the GWRK and RK could be regarded as the optimal interpolation methods, despite the fact that their improvement was limited in this study. Meanwhile, when considering the cost, time and process by which the models were established, OK can also be deemed an optional interpolation method with relatively acceptable accuracy. Furthermore, the evenness and density of sample should be considered when mapping soil TP in Mollisol areas.
Display omitted
•Eight interpolation methods were used to predict TP distribution in Mollisol area.•Both of GWRK and RK were the optimal methods to predict soil TP distribution.•Brightness index as auxiliary variable can improve the interpolation precision.•Slope steepness as the best auxiliary variable can improve interpolation precision.•Interpolation precision is influenced by the evenness and density of soil sampling.