We present the results of local earthquake tomography (LET) analysis to investigate the crust and uppermost mantle structure in the northern Dinarides and southwestern Pannonian basin. Datasets of P-wave travel times are inverted to recover a three-dimensional P-wave velocity model of the survey area. Two data subsets were used in this study: (1) data recorded on 15 temporary seismic stations, which were deployed in Croatia in the framework of ALPASS-DIPS project, and (2) travel time datasets from the Croatian Seismological Survey and ORFEUS databases. The data enabled to achieve a resolution of less than a hundred kilometres in horizontal directions and a few kilometres in vertical direction in the area with good ray coverage, as is documented by the resolution tests. Velocity variations are computed on a grid using the three-dimensional nonlinear tomographic inversion method. Our study provides the first crustal three-dimensional seismic model of the studied area, and it is correlated with previous results in the survey area allowing us to infer the main crustal structures with high confidence. The velocity model reveals crustal thickening beneath the Dinarides and significant crustal thinning beneath the Pannonian basin. The Moho surface was determined on the basis of the highest velocity gradients in the vertical cross-sections. We find relatively high velocities below the northern Dinarides at shallow depths (< 10 km), and low velocities caused by deep local depressions in the Pannonian basin. A very pronounced high-velocity body is present in the transitional part between the Dinarides and the Pannonian basin at a depth range of 5–15 km. The strong velocity increase at depth of about 20 km indicates that the Dinaridic crust could be interpreted as two-layered, while the Pannonian crust is probably one-layered.

Soil gas measurements of radon (²²²Rn), CO₂, and hydrocarbon concentrations, as well as gamma-ray spectrometry, were conducted at two separate locations to estimate the measurement results for known locations of hydrocarbon accumulations in the subsurface and oil seepage on the surface. The aim of the study was to confirm the applicability of the method for identifying migration pathways (e.g., faults) and to detect possible seepages of hydrocarbons to the surface as well as to investigate possible health issue potential about the soil gas analysis results. Site A investigations were performed with a large number of sampling points to provide sufficient spatial coverage to capture the influence of subsurface lithologic variability as well as the influence of the migration pathway on the measured parameters. For the investigation of site B, sampling points were positioned to reflect the situation between the area above producing hydrocarbon fields and areas with no confirmed accumulation. The results presented show that it is possible to distinguish the near-surface lithology (gamma-ray spectrometry), characterize the migration pathway, and indicate the area of oil seepage at the surface. Areas above the known hydrocarbon accumulations generally have elevated radon concentrations and detectable heavier hydrocarbons with sporadic methane in soil gas, which contrasts with the lower radon levels and lack of detectable heavier hydrocarbons in soil gas in the area with no confirmed hydrocarbon accumulation in the subsurface.

Geophysical exploration of bauxite deposits has been carried out in the area of Posušje in Bosnia and Herzegovina, which were formed on an Upper Cretaceous carbonate substrate, whereas the hanging wall rocks can be Paleogene limestones and sedimentary clastic rocks. Karst terrains are demanding for geophysical exploration due to the relatively complex geological relationships and exceptional near-surface inhomogeneities that generate large noises and challenging conditions for taking field measurements. The fundamental question is whether geophysical research can detect exceptionally irregular karst bauxite deposits with relatively small dimensions. The basic idea is to combine several geophysical methods and a joint interpretation of several data sets to increase the efficiency of geophysical surveying in detecting complex bauxite deposits in karst terrains. Therefore, fundamental near-surface research methods, electrical tomography and seismic refraction are used. In addition, magnetometry was used to examine whether bauxite deposits yield potential magnetic anomalies that could help in detecting them. Research undertaken in the area of Posušje was carried out in the first step on already discovered and known bauxite deposits to determine whether geophysical responses correlate with the occurrence of bauxite deposits and to evaluate the effectiveness of each of the applied surface geophysical methods. Measurements were taken at several locations, and results for two micro-locations, Krstače and Mratnjača, are shown. Geophysical measurements were firstly performed on discovered bauxite deposits in order to reliably determine the possibility of identifying deposits in geophysical inverse models. Bauxite deposits were clearly recognised as characteristic geophysical responses in inversion models using both methods, electrical tomography and seismic refraction. Although the response of bauxite deposits is expressed in both models, resistivity and velocity, it is much more evident in resistivity models. The characteristic resistivity response was confirmed by the discovery of a new deposit. Therefore, the conclusion is that electrical resistivity tomography should be considered a basic method for exploring karst bauxite deposits. Seismic refraction provides a better characterisation of deposits and reduces the interpretation ambiguity. This solution can generally be applied to the problem of researching bauxite deposits in the Dinarides and similar geological models in the Mediterranean. Magnetometric measurements have shown that no magnetic anomalies could be associated with bauxite deposits, and only magnetometry was not successful in discovering bauxite deposits.

In subsurface geological mapping, it is more than advisable to compare different solutions obtained with neural and other algorithms. Here, for such comparison, we used the previously published and well-prepared dataset of subsurface data collected from the Bjelovar Subdepression, a 2900 km² large regional macrounit in the Croatian part of the Pannonian Basin System. Data on depth were obtained for the youngest (the shallowest) Lonja Formation (Pliocene, Quaternary) and mapped using neural network (NN), inverse distance weighting (IDW), and ordinary kriging (OK) algorithms. The obtained maps were compared based on square error (using k-fold cross-validation) and the visual interpretation of isopaches. Two other algorithms were also tested, namely, random forest (RF) and extreme gradient boosting (XGB) algorithms, but they were rejected as inappropriate for this purpose solely based on the visuals of the obtained maps, which did not follow any interpretable geological structures. The results showed that NN is a highly adjustable method for interpolation, with adjustment for numerous hyperparameters. IDW showed its strength as one of the classical interpolators, and its results are always located close to the top if several methods are compared. OK is the relative winner, showing the flexibility of variogram analysis regarding the number of data points and possible clustering. The presented variogram model, even with a relatively high sill and occasional nugget effect, can be well fitted into OK, giving better results than other methods when applied to the presented area and datasets. This was not surprising because kriging is a well-established method used exclusively for interpolation. In contrast, NN and machine learning algorithms are used in many fields, and these algorithms, particularly the fitting of hyperparameters in NN, simply cannot be the best solution for all.

Based on the statistical analysis of the time series of data, the influence of the change of Sava River temperature and geothermal anomalies on the changes of the groundwater temperature of the Zagreb and Samobor-Zaprešić aquifers is...

Gamma-ray parameters are widely used to screen organic-rich intervals, yet the influence of lithofacies heterogeneity on the reliability of gamma-ray-based TOC predictions remains incompletely evaluated. This study aims to (1) assess the relationships between gamma-ray parameters and TOC, (2) evaluate how lithofacies heterogeneity influences predictive reliability, and (3) establish practical screening criteria for identifying organic-rich intervals in two Middle Miocene marine successions from the Croatian part of the Pannonian Basin System. Outcrop gamma-ray spectrometry measurements (K, U, Th) and dose rate (DR) were paired with co-located laboratory TOC analyses (Voćin n = 45; Podsused n = 96). Pearson correlation, multiple and simple linear regression, logistic regression for TOC ≥ 1 wt.%, and threshold analysis were applied to evaluate relationships between gamma-ray parameters and organic matter enrichment. U and DR show moderate positive relationships with TOC in both sections, whereas K and Th are not statistically significant predictors. Regression and threshold analyses identify U as the primary predictor of TOC and show more stable and systematic predictive behaviour in the homogeneous Podsused succession than in the heterogeneous Voćin section. The results indicate that lithofacies heterogeneity primarily influences predictive stability rather than the existence of TOC–gamma-ray relationships. The study demonstrates that lithofacies heterogeneity is an important factor influencing the reliability of gamma-ray-based TOC screening and provides a framework for assessing the suitability of uranium and dose-rate proxies in different sedimentary settings.