In the last years, new findings and new methods (stable isotopes of oxygen, zinc and nitrogen, 2D and 3D modelling, geometric morphometric analyses of the teeth) have enhanced our knowledge of the Neogene shark fauna and its palaeobiology. Several papers deal with the large Otodus (Megselachus) species, including the construction of a 3D model as well as insights into lifestyle and diet. In addition, skeletal remains of Carcharias gustrowensis, Carcharodon hastalis, Keasius parvus and a natural tooth set of Carcharodon hubbelli have been described in the last 13 years, and the dentition of the Neogene species Carcharoides catticus, Megachasma applegatei and Parotodus benedenii have been reconstructed. Stable isotope analyses of the teeth from the Neogene species of Araloselachus, Carcharias, Carcharodon, Galeocerdo, Hemipristris, and Mitsukurina have given insights into the trophic position of these genera during the Neogene, and shark teeth preserved near skeletal remains of prey animals (mammals) and shark bite traces on these remains provide direct evidence of trophic interactions. Tooth shape, fossil locality and palaeoenvironment have been used to better understand the taxa Carcharhinus dicelmai, Megalolamna paradoxodon, Pachyscyllium dachiardii and P. distans. Among extant species, Galeorhinus galeus can be traced back to the Eocene. The following taxa can be traced back to the Oligocene: ?Alopias superciliosus, and Rhincodon typus. Species already present in the Miocene include: Alopias vulpinus, Carcharhinus amblyrhynchoides, C. amblyrhynchos, C. albimarginatus, C. amboinensis, C. brachyurus, C. brevipinna, C. falciformis, C. glaucus, C. leucas, C. limbatus, C. longimanus, C. macloti, C. obscurus, C. perezi, C. sealei, ?Carcharodon carcharias, Centrophorus granulosus, Cetorhinus maximus, Dalatias licha, Deania calcea, Galeocerdo cuvier, , Glyphis glyphis, Heptranchias perlo, Isurus paucus, Lamna nasus, Negaprion brevirostris, Odontaspis ferox, Pseudocarcharias kamoharai, Sphyrna media, S. mokarran. First appearing in the Pliocene are: Scymnodon ringens, Somniosus rostratus, Zameus squamulosus. For some extant species (Carcharias taurus, Hexanchus griseus, Isurus oxyrinchus, Notorynchus cepedianus, Sphyrna zygaena) it is not clear if the assigned Neogene teeth represent the same species. Applying these new methods to more fossil shark taxa, a detailed search for shark fossils, as well as better knowledge of the dentition of extant species (especially those with minute-sized teeth) will further enhance knowledge of the evolution and palaeobiology of sharks.

Here is analysed the process of water re-injection in the two hydrocarbon reservoirs/fields in the Upper Miocene sandstone reservoirs, located in the western part of the Sava Depression (Croatia). Namely, those are the "A" field with "L" reservoir and the "B" field with "K" reservoir. Both currently produce hydrocarbons using a secondary recovery method, i.e. water injection (in fact re-injection of the field waters). Three regional reservoir variables had been analysed, namely porosity, permeability and injected water volumes. The number of data was small in all three cases. For porosity: reservoir “L” included 25 data, reservoir “K” 19 data; for permeability: reservoir “L” 10 data, reservoir “K” 18 data; for injected volumes of water: reservoir “L” 10 data; reservoir “K” 3 data. It defined selection of mapping algorithms mostly designed for small datasets (less than 20 points), i.e. Inverse Distance Weighting, Nearest and Natural Neighbourhood. Additionally, the Ordinary Kriging was used, but only with jack-knifed variograms, producing many “artificial points”. Results are extensively tested, using cross-validation and shape recognitions, and the Inverse Distance Weighting method is described as the most appropriate approach for mapping permeability and injected volumes in both reservoirs (“K” and “L”). The Kriging could be slightly outlined as the best approach for porosity. Obtained maps made possible application of the modified geological probability calculation as tools for prediction of successfulness of future injection (probability of 0.56). Consequently, results made possible to plan future injection more efficiently, with smaller injected volumes and same of higher hydrocarbon recovery. That could prevent useless injection, decrease number of injection wells, and save energy and funds invested in such processes.

Interpolation is procedure that depends on spatial and/or statistical properties of analysed variable(s). It is special challenging task for data that included low number of samples, like dataset with less than 20 data. This problem is especially emphasized in the subsurface geological mapping, i.e. in the cases where data are taken solely from wells. Successful solutions of such mapping problems ask for knowledge about interpolation methods designed primarily for small datasets and dataset itself. Here are compared two methods, namely Inverse Distance Weighting and Modified Shepard’s Method, applied for three variables (porosity, permeability, thickness) measured in the Neogene sandstone hydrocarbon reservoirs (Northern Croatia). The results showed that pure cross-validation is not enough condition for appropriate map selection, but also geometrical features need to be considered, for datasets with less than 20 points.

The basic principles of geological risk calculation through Probability of Success (PoS) are mostly applied for numerical estimation of additional hydrocarbon existence in proven reservoirs or potential hydrocarbon discoveries in selected geological regional subsurface volume. It can be tailored and validated for a comprehensive input dataset collected in the selected petroleum province, adapted by dividing up geological events into several probability categories and classes. The most applied categories are (existence of) reservoir rocks, traps and isolators, source rocks, migration pathways and preservation conditions for hydrocarbons The methodology results in unique probability values as multiplication of independent statistical events, which can also be applied in the assessment of a potential hydrocarbon discovery of desired minimal volume and its value in any virtual currency like risk-neutral dollars. Such methodology has been extensively developed in the last decades in the Croatian subsurface, mostly in the Croatian part of the Pannonian Basin System (CPBS). Through the adaptation of geological categories, it was also applied in hybrid, i.e., stochastical, models developed in the CPBS (Drava Depression). Stochastically estimation of porosity was already applied. As the robustness of this methodology is very high, it was also modified to estimate the influence of water-flooding in increasing oil recovery in some proven Neogene sandstone reservoirs in the CPBS (Sava Depression). This new modification is being presented to be applied to geological risk calculation, intending to assess the safety of geological environment in deep wells, where depleted radioactive fuel would be disposed, a subject of great importance. The case study encompassed the magmatic and metamorphic rocks in the pre-Neogene basement of the CPBS. For disposal purpose, these are regionally lithologies considered as to be the safest ones considering petrophysical values, water saturation, recent weathering and tectonic activity.