Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Hourly values of X, Y and Z geomagnetic components recorded at Coimbra Station (IAGA code: COI) and local geomagnetic K indices for 2007 to 2014

This study examines the performance of series of two geomagnetic indices and series synthesized from a semi-empirical model of magnetospheric currents, in explaining the geomagnetic activity observed at Northern Hemipshere's mid-latitude ground-based stations. We analyse data, for the 2007 to 2014 period, from four magnetic observatories (Coimbra, Portugal; Panagyurishte, Bulgary; Novosibirsk, Russia and Boulder, USA), at geomagnetic latitudes between 40° and 50° N. The quiet daily (QD) variation is firstly removed from the time series of the geomagnetic horizontal component (H) using natural orthogonal components (NOC) tools. We compare the resulting series with series of storm-time disturbance (Dst) and ring current (RC) indices and with H series synthesized from the Tsyganenko and Sitnov (2005, doi:10.1029/2004JA010798) (TS05) semi-empirical model of storm-time geomagnetic field. In the analysis, we separate days with low and high local K-index values. Our results show that NOC models are as efficient as standard models of QD variation in preparing raw data to be compared with proxies, but with much less complexity. For the two stations in Europe, we obtain indication that NOC models could be able to separate ionospheric and magnetospheric contributions. Dst and RC series explain the four observatory H-series successfully, with values for the mean of significant correlation coefficients, from 0.5 to 0.6 during low geomagnetic activity (K less than 4) and from 0.6 to 0.7 for geomagnetic active days (K greater than or equal to 4). With regard to the performance of TS05, our results show that the four observatories separate into two groups: Coimbra and Panagyurishte, in one group, for which the magnetospheric/ionospheric ratio in QD variation is smaller, a dominantly QD ionospheric contribution can be removed and TS05 simulations are the best proxy; Boulder and Novosibirsk,in the other group, for which the ionospheric and magnetospheric contributions in QD variation can not be differentiated and correlations with TS05 series can not be made to improve. The main contributor to magnetospheric QD signal are Birkeland currents. The relatively good success of TS05 model in explaining ground-based irregular geomagnetic activity at mid-latitudes makes it an effective tool to classify storms according to their main sources. For Coimbra and Panagyurishte in particular, where ionospheric and magnetospheric daily contributions seem easier to separate, we can aspire to use the TS05 model for ensemble generation in space weather (SW) forecasting and interpretation of past SW events.

Age models and stabile isotope ratios from SO35 cores

Stable isotope measurements on the planktonic foraminifer Globigerinoides ruber (white) have been carried out on a number of selected deep-seas sediment cores from the South Lau and Norlh Fiji Basins. The d18O-curves show good correlation with the inter-ocean oraphic correlation composite d18O-record of the standard reference section (Prell et al. 1986), which, in combination with the chronostratigraphic classifications of Herterich & Sarnthein (1984, modified) and Imbrie et al. 1984), allows a detailed dating of the sedimentary sequences. The deepest layers in core no. 119 (southern Lau Basin) could be assigned to Isotope Stage 24. Measurements made on bulk carbonate in two cores show a much higher glacial-interglacial amplitude, allowing the general identification of the conventional oxygen isotope stages. The d13C-values of the benthic foraminifer Cibicidoides wuellerstorfi show progressively lighter values northwards reflecting an increasing contribution of the isotopically lighter CO2 from the remineralisation of organic matter during the general northward movement of the deep water masses. Cyclicities in the sedimentation rates were observed in core nos. 117 and 119 (both southern Lau Basin) where the interglacials exhibit higher levels than the glacials. Calculated new or export paleoproductivity show that the glacials had higher productivity in the euphotic zone. From the oxygen isotope stratigraphy, the five ash layers in core nos. 117 and 119 could be dated as about 530 ka B.P. in Stage 14, 695 ka B.P. in Stage 18, 775 ka B.P. in Stage 21, 790 ka B.P. and 825 ka B.P. in Stage 22. Carbonate dissolution occurred during stages 5, 8 and 10 to 12.

DISPRO datasets for validation of coastal hydrodynamic models

Appropriate field data are required to check the reliability of hydrodynamic models simulating the dispersion of soluble substances in the marine environment. This study deals with the collection of physical measurements and soluble tracer data intended specifically for this kind of validation. The intensity of currents as well as the complexity of topography and tides around the Cap de La Hague in the center of the English Channel makes it one of the most difficult areas to represent in terms of hydrodynamics and dispersion. Controlled releases of tritium - in the form of HTO - are carried out in this area by the AREVA-NC plant, providing an excellent soluble tracer. A total of 14 493 measurements were acquired to track dispersion in the hours and days following a release. These data, supplementing previously gathered data and physical measurements (bathymetry, water-surface levels, Eulerian and Lagrangian current studies) allow us to test dispersion models from the hour following release to periods of several years which are not accessible with dye experiments. The dispersion characteristics are described and methods are proposed for comparing models against measurements. An application is proposed for a 2 dimensions high-resolution numerical model. It shows how an extensive dataset can be used to build, calibrate and validate several aspects of the model in a highly dynamic and macrotidal area: tidal cycle timing, tidal amplitude, fixed-point current data, hodographs. This study presents results concerning the model's ability to reproduce residual Lagrangian currents, along with a comparison between simulation and high-frequency measurements of tracer dispersion. Physical and tracer data are available from the SISMER database of IFREMER (www.ifremer.fr/sismer/catal). This tool for validation of models in macro-tidal seas is intended to be an open and evolving resource, which could provide a benchmark for dispersion model validation.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Age models, organic analyses and major-element analyses of sediment cores on a transect of the tropical African rainbelt

The distribution of rainfall in tropical Africa is controlled by the African rainbelt**1, which oscillates on a seasonal basis. The rainbelt has varied on centennial to millennial timescales along with changes in Northern Hemisphere high-latitude climate**2, 3, 4, 5, the Atlantic meridional overturning circulation**6 and low-latitude insolation**7 over the past glacial-interglacial cycle. However, the overall dynamics of the African rainbelt remain poorly constrained and are not always consistent with a latitudinal migration**2, 4, 5, 6, as has been proposed for other regions**8, 9. Here we use terrestrially derived organic and sedimentary markers from marine sediment cores to reconstruct the distribution of vegetation, and hence rainfall, in tropical Africa during extreme climate states over the past 23,000 years. Our data indicate that rather than migrating latitudinally, the rainbelt contracted and expanded symmetrically in both hemispheres in response to changes in climate. During the Last Glacial Maximum and Heinrich Stadial 1, the rainbelt contracted relative to the late Holocene, which we attribute to a latitudinal compression of atmospheric circulation associated with lower global mean temperatures**10. Conversely, during the mid-Holocene climatic optimum, the rainbelt expanded across tropical Africa. In light of our findings, it is not clear whether the tropical rainbelt has migrated latitudinally on a global scale, as has been suggested**8,9.

Age models of northern North Atlantic sediment cores

Benthic foraminiferal stable isotope records from four high-resolution sediment cores, forming a depth transect between 1237 m and 2303 m on the South Iceland Rise, have been used to reconstruct intermediate and deep water paleoceanographic changes in the northern North Atlantic during the last 21 ka (spanning Termination I and the Holocene). Typically, a sampling resolution of ~100 years is attained. Deglacial core chronologies are accurately tied to North Greenland Ice Core Project (NGRIP) ice core records through the correlation of tephra layers and changes in the percent abundance of Neogloboquadrina pachyderma (sinistral) with transitions in NGRIP. The evolution from the glacial mode of circulation to the present regime is punctuated by two periods with low benthic d13C and d18O values, which do not lie on glacial or Holocene water mass mixing lines. These periods correlate with the late Younger Dryas/Early Holocene (11.5-12.2 ka) and Heinrich Stadial 1 (14.7-16.8 ka) during which time freshwater input and sea-ice formation led to brine rejection both locally and as an overflow exported from the Nordic seas into the northern North Atlantic, as earlier reported by Meland et al. (2008). The export of brine with low ?13C values from the Nordic seas complicates traditional interpretations of low d13C values during the deglaciation as incursions of southern sourced water, although the spatial extent of this brine is uncertain. The records also reveal that the onset of the Younger Dryas was accompanied by an abrupt and transient (~200-300 year duration) decrease in the ventilation of the northern North Atlantic. During the Holocene, Iceland-Scotland Overflow Water only reached its modern flow strength and/or depth over the South Iceland Rise by 7-8 ka, in parallel with surface ocean reorganizations and a cessation in deglacial meltwater input to the North Atlantic.

Age models, sedimentology and geochemistry from sediment cores in Drake Passage throughflow

The Drake Passage (DP) is the major geographic constriction for the Antarctic Circumpolar Current (ACC) and exerts a strong control on the exchange of physical, chemical, and biological properties between the Atlantic, Pacific, and Indian Ocean basins. Resolving changes in the flow of circumpolar water masses through this gateway is, therefore, crucial for advancing our understanding of the Southern Ocean's role in global ocean and climate variability. Here, we reconstruct changes in DP throughflow dynamics over the past 65,000 y based on grain size and geochemical properties of sediment records from the southernmost continental margin of South America. Combined with published sediment records from the Scotia Sea, we argue for a considerable total reduction of DP transport and reveal an up to ~40% decrease in flow speed along the northernmost ACC pathway entering the DP during glacial times. Superimposed on this long-term decrease are high-amplitude, millennial-scale variations, which parallel Southern Ocean and Antarctic temperature patterns. The glacial intervals of strong weakening of the ACC entering the DP imply an enhanced export of northern ACC surface and intermediate waters into the South Pacific Gyre and reduced Pacific-Atlantic exchange through the DP ("cold water route"). We conclude that changes in DP throughflow play a critical role for the global meridional overturning circulation and interbasin exchange in the Southern Ocean, most likely regulated by variations in the westerly wind field and changes in Antarctic sea ice extent.

Age models, ostracods and trace elements of five gravity cores from the NW Black Sea

New Mg/Ca, Sr/Ca, and published stable oxygen isotope and 87Sr/86Sr data obtained on ostracods from gravity cores located on the northwestern Black Sea slope were used to infer changes in the Black Sea hydrology and water chemistry for the period between 30 to 8 ka B.P. (calibrated radiocarbon years). The period prior to 16.5 ka B.P. was characterized by stable conditions in all records until a distinct drop in d18O values combined with a sharp increase in 87Sr/86Sr occurred between 16.5 and 14.8 ka B.P. This event is attributed to an increased runoff from the northern drainage area of the Black Sea between Heinrich Event 1 and the onset of the Bølling warm period. While the Mg/Ca and Sr/Ca records remained rather unaffected by this inflow; they show an abrupt rise with the onset of the Bølling/Allerød warm period. This rise was caused by calcite precipitation in the surface water, which led to a sudden increase of the Sr/Ca and Mg/Ca ratios of the Black Sea water. The stable oxygen isotopes also start to increase around 15 ka B.P., although in a more gradual manner, due to isotopically enriched meteoric precipitation. While Sr/Ca remains constant during the following interval of the Younger Dryas cold period, a decrease in the Mg/Ca ratio implies that the intermediate water masses of the Black Sea temporarily cooled by 1-2°C during the Younger Dryas. The 87Sr/86Sr values drop after the cessation of the water inflow at 15 ka B.P. to a lower level until the Younger Dryas, where they reach values similar to those observed during the Last Glacial Maximum. This might point to a potential outflow to the Mediterranean Sea via the Sea of Marmara during this period. The inflow of Mediterranean water started around 9.3 ka B.P., which is clearly detectable in the abruptly increasing Mg/Ca, Sr/Ca, and 87Sr/86Sr values. The accompanying increase in the d18O record is less pronounced and would fit to an inflow lasting ~100 a.

Age models and summer sea surface temperature and winter sea ice concentration for the EPILOG-LGM time slice in the Pacific Southern Ocean

Sea surface temperatures and sea-ice extent are the most critical variables to evaluate the Southern Ocean paleoceanographic evolution in relation to the development of the global carbon cycle, atmospheric CO2 variability and ocean-atmosphere circulation. In contrast to the Atlantic and the Indian sectors, the Pacific sector of the Southern Ocean has been insufficiently investigated so far. To cover this gap of information we present diatom-based estimates of summer sea surface temperature (SSST) and winter sea-ice concentration (WSI) from 17 sites in the polar South Pacific to study the Last Glacial Maximum (LGM) at the EPILOG time slice (19,000-23,000 cal. years BP). Applied statistical methods are the Imbrie and Kipp Method (IKM) and the Modern Analog Technique (MAT) to estimate temperature and sea-ice concentration, respectively. Our data display a distinct LGM east-west differentiation in SSST and WSI with steeper latitudinal temperature gradients and a winter sea-ice edge located consistently north of the Pacific-Antarctic Ridge in the Ross sea sector. In the eastern sector of our study area, which is governed by the Amundsen Abyssal Plain, the estimates yield weaker latitudinal SSST gradients together with a variable extended winter sea-ice field. In this sector, sea-ice extent may have reached sporadically the area of the present Subantarctic Front at its maximum LGM expansion. This pattern points to topographic forcing as major controller of the frontal system location and sea-ice extent in the western Pacific sector whereas atmospheric conditions like the Southern Annular Mode and the ENSO affected the oceanographic conditions in the eastern Pacific sector. Although it is difficult to depict the location and the physical nature of frontal systems separating the glacial Southern Ocean water masses into different zones, we found a distinct temperature gradient in latitudes straddled by the modern Southern Subtropical Front. Considering that the glacial temperatures north of this zone are similar to the modern, we suggest that this represents the Glacial Southern Subtropical Front (GSSTF), which delimits the zone of strongest glacial SSST cooling (>4K) to its North. The southern boundary of the zone of maximum cooling is close to the glacial 4°C isotherm. This isotherm, which is in the range of SSST at the modern Antarctic Polar Front (APF), represents a circum-Antarctic feature and marks the northern edge of the glacial Antarctic Circumpolar Current (ACC). We also assume that a glacial front was established at the northern average winter sea ice edge, comparable with the modern Southern Antarctic Circumpolar Current Front (SACCF). During the glacial, this front would be located in the area of the modern APF. The northward deflection of colder than modern surface waters along the South American continent leads to a significant cooling of the glacial Humboldt Current surface waters (4-8K), which affects the temperature regimes as far north as into tropical latitudes. The glacial reduction of ACC temperatures may also result in the significant cooling in the Atlantic and Indian Southern Ocean, thus may enhance thermal differentiation of the Southern Ocean and Antarctic continental cooling. Comparison with temperature and sea ice simulations for the last glacial based on numerical simulations show that the majority of modern models overestimate summer and winter sea ice cover and that there exists few models that reproduce our temperature data rather well.