The paper aims to describe the preconditioning and observations of exceptionally high salinity values that were observed in summer and autumn of 2017 in the Adriatic. The observations encompassed CTD measurements carried out along the well-surveyed climatological transect in the Middle Adriatic (the Palagruža Sill, 1961–2020), Argo profiling floats and several glider missions, accompanied with satellite altimetry and operational ocean numerical model (Mediterranean Forecasting System) products. Typically, subsurface salinity maximum, with values lower than 39.0, is observed in the Southern Adriatic (usually between 200 and 400 m), related to ingressions of saltier and warmer waters originating in the eastern Mediterranean (Levantine Intermediate Water—LIW). However, seasonally strong inflow of warm and high salinity waters (S > 38.8) has been observed much closer to the surface since spring 2015. The main LIW core deepened at the same time (to 400–700 m). Such double-maxima vertical pattern was eventually disturbed by winter convection at the beginning of 2017, increasing salinities throughout the water column. A new episode of very strong inflow of high salinity waters from the Northern Ionian was observed in late winter and spring of 2017, this time restricted almost to the surface. As most of 2017 was characterized by extremely dry conditions, low riverine inputs and warmer than usual summer over the Adriatic and Northern Ionian, salinity values above the sharp and shallow (15–40 m) thermocline significantly increased. The maximum recorded salinity was 39.26, as measured by the Argo float in the Southern Adriatic. Surface salinity maximum events, but with much lower intensity, have been documented in the past. Both past events and the 2017 event were characterized by (i) concurrence with overall high salinity conditions and cyclonic or transitional phase of the Adriatic-Ionian Bimodal Oscillating System, (ii) very low river discharges preconditioning the events for a year or more, (iii) higher-than-average heat fluxes during most of the summer and early autumn periods, forming a stable warm layer above the thermocline, and (iv) higher-than-average E-P (evaporation minus precipitation) acting on this warm surface layer. Importantly, the 2017 event was also preceded by strong near-surface inflow of very saline waters from the Northern Ionian in early 2017.

Interannual variability of the primary production in the middle Adriatic Sea for the period 1961–2002 was examined and correlated to the various atmospheric and oceanographic parameters. The sequentialt-test analysis of regime shift (STARS) method and locally-weighted scatter plot smoothing (LOWESS) method were applied to the primary production, revealing the new regime with significantly different mean productivity ranging from 1980–1996. Moreover, this period with the highest primary production, consists of the two distinguished sub-periods: periods of increasing (1980–1986) and decreasing (1987–1996) primary production. Whereas in the first period the ecosystem was under the influence of warmer and nutrient richer Levantine Intermediate Water (LIW) intrusions into the Adriatic, in the second period, which started with a cold winter in 1987, the Eastern Mediterranean Transient (EMT) occurred. The EMT established a new circulation regime which prevented the LIW intrusions in the Adriatic, causing its reduced productivity. Reduced LIW inflow in the Adriatic was evidenced in the lower than normal sea temperature, salinity and oxygen concentrations below the thermocline depth. Precipitation and wind regime also arose as important local factors for the primary production variability. Our analysis connected the shifts in primary production with hemispheric and regional scale climate variations, and supports the hypothesis that atmospheric variability can trigger the ecosystem changes.

Oceanographic parameters, such as sea surface temperature, surface chlorophyll-a concentration, sea surface ice concentration, sea surface height, etc., are listed as Essential Climate Variables. Therefore, there is a crucial need for persistent and accurate measurements on a global scale. While in situ methods tend to be accurate and continuous, these qualities are difficult to scale spatially, leaving a significant portion of Earth’s oceans and seas unmonitored. To tackle this, various remote sensing techniques have been developed. One of the more prominent ways to measure the aforementioned parameters is via satellite spacecraft-mounted remote sensors. This way, spatial coverage is considerably increased while retaining significant accuracy and resolution. Unfortunately, due to the nature of electromagnetic signals, the atmosphere itself and its content (such as clouds, rain, etc.) frequently obstruct the signals, preventing the satellite-mounted sensors from measuring, resulting in gaps—missing data—in satellite recordings. One way to deal with these gaps is via various reconstruction methods developed through the past two decades. However, there seems to be a lack of review papers on reconstruction methods for satellite-derived oceanographic variables. To rectify the lack, this paper surveyed more than 130 articles dealing with the issue of data reconstruction. Articles were chosen according to two criteria: (a) the article has to feature satellite-derived oceanographic data (b) gaps in satellite data have to be reconstructed. As an additional result of the survey, a novel categorising system based on the type of input data and the usage of time series in reconstruction efforts is proposed.