Karst aquifers are important sources of thermal and groundwater in many parts of the world, such as the Alpine–Dinaric–Carpathian region in Europe. The Upper Triassic dolomites are regionally recognized thermal and groundwater aquifers but also hydrocarbon reservoirs. They are characterized by predominantly fractured porosity, but the actual share of depositional and diagenetic porosity is rarely investigated. In this research, we presented the geometric characterization of the measured microporosity of the Upper Triassic dolomites of the Žumberak Mts (Croatia), through thin-section image processing and particle analysis techniques. Pore parameters were analyzed on microphotographs of impregnated thin sections in scale. A total of 2267 pores were isolated and analyzed. The following parameters were analyzed: pore area, pore perimeter, circularity, aspect ratio (AR), roundness, solidity, Feret AR, compactness, and fractal dimension. Furthermore, porosity was calculated based on the pore portion in each image. The effective porosity on rock samples was determined using saturation and buoyancy techniques as an accompanying research method. We analyzed distributions of each parameter, their correlation, and most of the parameters are characterized by an asymmetric or asymmetric normal distribution. Parameters that quantify pore irregularities have similar distributions, and their values indicate the high complexity of the pore geometry, which can significantly impact permeability.

Fractured aquifers, especially dolomites, are important hydrocarbon reservoirs and sources of thermal and groundwater in many parts of the world, especially in the Alpine-Dinaric-Carpathian region. The most dominant porosity type is fracture porosity, which acts as the preferential fluid pathway in the subsurface, thus strongly controlling fluid flow. Outcrops provide valuable information for the characterization of fracture networks. Dolomite rock properties and structural and diagenetic processes result in fractured systems that can be considered fractals. The fracture network was analyzed on 14 vertical outcrops in 35 digitized photographs. The values of the fractal dimensions varied slightly by the software and method used, but the trends were consistent, which confirms that all methods are valid. Small values of fractal dimension indicate the dominance of a few small or large fractures, and high values of fractal dimension result from a combination of large numbers of small fractures accompanied by a few large fractures. The mean value of the fractal dimension for analyzed fracture networks was 1.648. The results indicate that the fracture network of the Upper Triassic dolomites can be approximated by fractal distribution and can be considered a natural fractal, and values can be extrapolated to higher and lower scales (1D and 3D).

The Organic Rankine Cycle (ORC) is commonly accepted as a viable technology to convert from low to medium temperature geothermal energy into electrical energy. In practice, the reference technology for converting geothermal energy to electricity is the subcritical simple ORC system. Over time, geothermal ORC plants with more complex configurations (architectures) have been developed. In the open literature, a large number of advanced architectures or configurations have been introduced. An analysis of the scientific literature indicates that there is some confusion regarding the terminology of certain advanced ORC system architectures. A new categorization of advanced configurations has been proposed, with a special emphasis on the application of geothermal energy. The basic division of advanced plant configurations is into dual-pressure and dual-stage ORC systems. In this study, the real potential of advanced ORC architectures or configurations to improve performance as compared with the simple ORC configuration was explored. The research was conducted for a wide range of geothermal heat source temperatures (from 120 °C to 180 °C) and working fluids. Net power output improvements as compared with the basic subcritical simple ORC (SORC) configuration were examined. The ability to produce net power with different ORC configurations depends on the magnitude of the geothermal fluid temperature and the type of working fluid. At a lower value of geothermal fluid temperature (120 °C), the most net power of 18.71 (kW/(kg/s)) was realized by the dual-pressure ORC (DP ORC configuration) with working fluid R1234yf, while the double stage serial-parallel ORC configuration with a low-temperature preheater in a high-temperature stage ORC (DS parHTS LTPH ORC) generated 18.51 (kW/(kg/s)) with the working fluid combination R1234yf/R1234yf. At 140 °C, three ORC configurations achieved similar net power values, namely the simple ORC configuration (SORC), the DP ORC configuration, and the DS parHTS LTPH ORC configuration, which generated 31.03 (kW/(kg/s)) with R1234yf, 31.07 (kW/(kg/s)) with R1234ze(E), and 30.96 (kW/(kg/s)) with R1234ze(E)/R1234yf, respectively. At higher values of geothermal fluid temperatures (160 °C and 180 °C) both the SORC and DP ORC configurations produced the highest net power values, namely 48.58 (kW/(kg/s)) with R1234ze(E), 67.23 (kW/(kg/s)) with isobutene for the SORC configuration, and 50.0 (kW/(kg/s)) with isobutane and 69.67 (kW/(kg/s)) with n-butane for the the DP ORC configuration.

The Dinaric karst in the north differs from the rest of the karst in Croatia in terms of karstification depth. The infiltrating precipitation drains in cascades from deeply karstified mountainous areas to the shallow or fluviokarst, forming the tributaries of the Kupa River. Time series analyses were conducted on a 5-year dataset to elucidate the hydrogeological conceptual model of the area and clarify disparate findings from tracer tests under varying hydrological conditions. The flow duration curve, autocorrelation functions, and recession curves were used to evaluate the spring discharge variability, the karstification degree, and the karst aquifer’s size. The crosscorrelation function and temperature dynamics were employed to assess the spring’s response to recharge and the hydrogeological system behavior. Comparative analysis with previous studies was conducted to contextualize the obtained results. The research outcomes delineated several key findings: (i) the deep karst zone is less developed than the shallow karst zone; (ii) groundwater exchange is significantly faster in shallow karst; (iii) groundwater divides in the Kapela Mountain are zonal; (iv) the homogenization of groundwater occurs during periods of high water levels; (v) fast water exchange transpires without concurrent groundwater temperature homogenization; and (vi) a definition of the boundary between deep and fluviokarst in Croatia.

This study presents the findings derived from a three-year monitoring programme focusing on the hydrochemical composition and stable isotope signatures (δ18O, δ2H) within the vulnerable karst groundwater system in Croatia's northern Dinaric karst region. Covering an area of 1980 km2 across Kapela Mountain and its foothills, this groundwater system falls within the Kupa River basin and encompasses the catchment areas of the main springs across two spring levels of the Dobra, Mrežnica, and Slunjčica Rivers (namely the Vitunj, Zagorska Mrežnica, Dretulja, Veliko Vrelo, Gojak, Tounjčica, Mrežnica, and Slunjčica springs). Given the complex hydrogeological framework, prior studies have not extensively addressed the hydrochemical characteristics of this region, thereby necessitating a comprehensive investigation to elucidate system dynamics. This paper outlines the first investigation employing stable isotopic analyses within the karst groundwater system of Kapela Mt. The main results and conclusions of the research include: (i) the aquifers across Kapela Mt. drain mainly limestones, (ii) water–rock interaction and carbonate weathering are significant contributors to water geochemistry, (iii) there is a slight human impact on the Gojak and Tounjčica spring waters, (iv) the mean residence time of water in the observed aquifer is up to 1.5 years, (v) the new LMWL was calculated for Kapela Mt. based on a three-year dataset (June 2018 – May 2021), (vi) in the north Dinaric karst, the predominant origin of precipitation is from the Mediterranean air mass, (vii) Velebit Mt. has a strong influence on the precipitation isotopic composition of the study area.

This paper presents and derives mathematical formulas that describe radial flow to a pumping well. Firstly, the general partial differential equation describing groundwater flow is derived using the principle of mass conservation. Then a fully penetrating pumping well is introduced into consideration, and the groundwater flow equation is adjusted accordingly. The cases of a confined and an unconfined aquifer are considered separately, as well as the cases of transient and steady-state flow. For each of these cases, the solution to the corresponding differential equation is derived and explained.