Modeling bio-geomorphological influences for offshore sandwaves

The coastal environment shows a wide range of bed patterns, for which sandwaves and sandbanks are among the most common. Less known in this context is the high benthic diversity in the coastal environment, which gives rise to the question to what extend the benthos interacts with the shape of the seabed. This paper reviews field and flume experiments on bio-geomorphological influences between benthos and sediment and tests the hypothesis that both the occurrence and the dimensions of sandwaves are dependent on the benthic diversity in the North Sea. Mathematical inclusions to account for biological activity in idealized models reveal that biota is able to influence the wavelength of sandwaves significantly, compared to the default case. More importantly, the models indicate that biota is able to induce bed patterns under conditions when the physical parameters suggest a stable flat bed and vice versa. Present model explorations indicate that future research should focus on the parameterization of subtidal biological activity on sediment dynamics and thereby on seabed patterns. Such knowledge will enable process-based modeling of the spatial and temporal variation in biological activity on seabed morphodynamics and validate the proposed modeling approach with field measurements.     

Sediment transport and resuspension due to combined motion of wave and current in the northern Adriatic Sea during a Bora event in January 2001: A numerical modelling study

The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.     

An experimental study on erodibility of intertidal sediments in the Yellow River delta

The erodibility of intertidal sediments is an important factor affecting coastal erosion.In July and October 2008,in situ measurement of erodibility of the surficial sediment were conducted using a recirculating flume at 20 tidal flat experiment sites along the seashore of the Yellow River delta.At the same time,the characteristics of sand ripples and biogenic features on the tidal flat were observed and the physical-mechanical sediment properties such as bulk density,water content,grain size distribution,plasticity,penetration resistance,shear strength and compressibility,were measured.By field measurement,it is obtained that the critical erosion shear stress of the surficial sediment on the tidal flat varies between 0.088 Pa and 0.254 Pa.The factors influencing sediment erodibility are complicated because of physical and biological reworking after the sediment deposited.There’s a positive correlation between shear strength and critical erosion shear stress.The burrowing crabs’ activities changed the sediment microtopography and made the sediment have greater roughness,and that is one possible reason for the higher erodibility.The formation of scour pits on the tidal flat correlates with the heterogeneous erodibility of the surficial sediment.     

Numerical modeling of tidal currents, sediment transport and morphological evolution in Hangzhou Bay, China

A 2D depth-averaged numerical model is set up to simulate the macro-scale hydrodynamic characteristics, sediment transport patterns and morphological evolution in Hangzhou Bay, a large macro-tidal estuary on the eastern coast of China. By incorporating the shallow water equations, the suspended sediment transport equation and the mass-balance equation for sediment; short-term hydrodynamics, sediment transport and long-term morphological evolution for Hangzhou Bay are simulated and the underlying physical mechanisms are analyzed. The model reproduces the spatial distribution patterns of suspended sediment concentration (SSC) in Hangzhou Bay, characterized by three high SSC zones and two low SSC zones. It also correctly simulates the residual flow, the residual sediment transport and the sediment accumulation patterns in Hangzhou Bay. The model results are in agreement with previous studies based on field measurements. The residual flow and the residual sediment transport are landwards directed in the northern part of the bay and seawards directed in the southern part. Sediment accumulation takes place in most areas of the bay. Harmonic analysis revealed that the tide is flood-dominant in the northern part of the bay and ebb-dominant in the southern part of the bay. The strength of the flood-dominance increases landwards along the northern Hangzhou Bay. In turn sediment transport in Hangzhou Bay is controlled by this tidal asymmetry pattern. In addition, the direction of tidal propagation in the East China Sea, the presence of the archipelago in the southeast and the funnel-shaped geometry of the bay, play important roles for the patterns of sediment transport and sediment accumulation respectively.     

Using geostatistics to elucidate temporal change in the spatial variation of aeolian sediment transport

Little is known about the spatial and temporal scales of variation in aeolian processes. Studies that aim to investigate surface erodibility often sample aeolian sediment transport at the nodes of a regular grid of arbitrary size. Few aeolian transport investigations have the resources to obtain sufficient samples to produce reliable models for mapping the spatial variation of transport. This study reports the use of an innovative nested strategy for sampling multiple spatial scales simultaneously using 40 sediment samplers. Reliable models of the spatial variation in aeolian sediment transport were produced and used for ordinary punctual kriging and stochastic simulated annealing to produce maps for several wind erosion events over a 25 km² playa in western Queensland, Australia. The results support the existence of a highly dynamic wind erosion system that was responding to possibly cyclic variation in the availability of material and fluctuations in wind energy. The spatial scale of transport was considerably larger than the small scale expected of the factors controlling surface erodibility. Thus, it appears that transport cannot be used as a surrogate of erodibility at the scale of this investigation. Simulation maps of transport provided considerably more information than those from kriging about the variability in aeolian sediment transport and its possible controlling factors. The proposed optimal sampling strategy involves a nested approach using ca 50 samplers. Copyright © 2003 John Wiley & Sons, Ltd.     

Modeling of sequence geometry north of Gargano Peninsula by changing sediment pathways in the Adriatic Sea

Sea level changes can reorganize sediment pathways on continental shelves in ways that can alter sediment supply and the resulting sedimentary deposits. The Adriatic Sea is one place where changing sediment pathways and along-strike currents have a major impact on sequence architecture. The Adriatic Sea, the marine portion of the Apenninic and Dinaric–Hellenic foreland basins, is being filled longitudinally, similar to other active forelands, and sediment transport patterns dramatically reorganized during Quaternary sea level cycles. We investigate the dynamics of sequence formation in the central Adriatic near the Gargano Promontory and the Mid-Adriatic Deep (MAD), where four depositional sequences each recording 100-kyr glacio-eustatic cycles have been mapped. These sequences are composed primarily of progradational units separated by regional unconformities. The geometry of the units is such that the constitutive clinoforms flatten out at their seaward termination into relatively planar strata, particularly in the upper parts of each sequence. Attempts to numerically simulate the sequences using the modeling software Sequence4 were frustrated by the difficulty of flattening the clinoforms seaward of the rollover (or depositional shelf break). The clinoform flattening observed in the Adriatic sequences contrasts with clinoform and depositional shelf break development that is characteristic of both conceptual and numerical models of sequences, including the one used here. We, therefore, modified the numerical model to account for the changes in sediment pathways that occur in the Adriatic Sea. During times with high sea level, such as the present, sediment from the Po and smaller Apennine Rivers is transported southwards along the coast by marine coastal currents and storms. At times of low sea level, fluvial transport of an enlarged Po River, in which the Apennine Rivers are likely captured as tributaries, discharges directly into the MAD basin. This produces a reciprocal sediment supply pattern where the coastal dispersal at high sea level produces prograding clinoforms, but cuts off as the exposure of the northern shelf at low sea level switches supply to the fluvial system, which discharges into the 260-m deep MAD. When the model is adapted for the shift in supply, the clinoforms flatten as nearshore sediment supply decreases and is replaced by deposition in the MAD during sea level lowstands. Including these sediment supply changes as a function of sea level enabled us to obtain a good fit to the overall stratigraphic architecture, supporting conceptual depositional models based on seismic stratigraphy. Similar dramatic shifts in sediment supply and sequence architecture may also occur at other settings, such as where canyons capture fluvial systems and bypass the terrigenous sediment supply to the deep basin. The paucity of observations of the unusual geometry seen at the Adriatic margin suggests that only rarely does the shore and fluvial sediment discharge overreach the shelf edge and cut off along-strike sediment transport to continental margins.     

Chlorinated compounds and polybrominated diphenyl ethers (PBDEs) in mussels (Mytilus galloprovincialis) collected from Apulia Region coasts

This project was carried out to assess the levels and spatial distribution of organochlorine compounds in the coastal marine environment, using mussels as bioindicators to evaluate the coastal water quality. Levels of polychlorobiphenils (PCB), chlorinated pesticides (DDT isomers, HCH isomers, Aldrin, Dieldrin, alfa-Endosulfan, Hexachlorobenzene, Pentachlorobenzene) and polybrominated diphenyl ethers (PBDEs) were determined in tissues from mussels (Mytilus galloprovincialis) collected along the Apulia Region coasts (Mediterranean Sea). Results indicate that contamination by organochlorine compounds is higher in mussels sampled in the Ionian Sea than in those from the Adriatic Sea, with PCB levels up to seven times higher in mussels from Ionian than from the Adriatic Sea. Although PCB levels were above the maximum values indicated by both European Community (EC) and National regulation in several sample sites, the PCB concentrations were particularly high in some stations, suggesting that these locations require a much specific attention. Conversely, results on the mussel contamination by PBDEs highlight their ubiquitous environmental distribution, and underline the need to establish the maximum level for these compounds in foodstuff, according to European Regulations.     

利用GRACE数据反演东海沉积物变化

GRACE卫星的成功发射为海底沉积物的监测提供了新的方法.利用2003-2014年间的GRACERL05数据,采用同期的测高数据对海面高变化进行改正,使用水文模式数据和基于均一假设的尺度因子估计方法处理泄漏误差,反演了东海地区的沉积物变化情况,并对GIA效应进行了改正.结果表明:东海入海口处沉积物的平均变化速率为5.44±0.88mm·a^-1,最大值出现在浙江沿海地区,变化速率为6~7mm·a^-1;在空间分布上,呈现河口处沉积速率大,远离河口的大洋地区沉积速率小的特征.在时空分布上均与实测数据很好的吻合.沉积物变化时间序列的周年项振为6.8cm,周年变化主要与东海泥沙扩散路径相关的海洋环流模式有关;半周年项和两周年项振幅分别为0.6cm和0.7cm,这两项变化主要与长江流域降水引起的土壤侵蚀变化有关.最后,分析讨论了本文沉积物监测方法推广到其他地区的适用性和局限性.     

Modern sediment accumulation on the Po shelf,Adriatic Sea

Sediment accumulation over the past century on the continental shelf near the Po delta varies with distance from the most active distributary channels. Near the Pila and Goro distributaries, sediment accumulation is rapid (1–4 cm yr⁻¹) and occurs in pulses. In these areas, the seabed is dominated by physical sedimentary structures that can be related to flood sedimentation. Between the two distributaries and in the southern portion of the dispersal system, sediment accumulation is slower (rates reach a minimum of 0.23 cm yr⁻¹ at ∼50 km from the Pila mouth) and steady-state, reflecting more continuous dispersal of sediment during non-flood periods. Sedimentary strata in these locations are composed of finer (clayey silt), mottled sediment. The similarity in the spatial distribution of long-term (100-yr) sediment accumulation to deposition resulting from the 2000 flood event suggests that the Po shelf is flood-dominated.About half of the sediment delivered by the Po River on a 100-yr time scale can be accounted for in the seabed deposit within ∼50 km of the Pila mouth. The remaining sediment is likely transported southward by the prevailing circulation, and this sediment coalesces with inputs from the Apennine Rivers.     

The effects of sand extraction on the macrobenthos of a relict sands area (northern Adriatic Sea): results 12 months post-extraction

Sands for the nourishment of beaches along the Emilia-Romagna coast (northern Adriatic Sea) were dredged from an offshore area characterised by relict sands formed during the last Adriatic post-glacial transgression. The short-term effects of the sand extraction on macrozoobenthic communities were investigated before, during and 1, 6 and 12 months after dredging at three impacted stations and seven control stations. Sand extraction activities did not significantly influence the granulometry and %TOC in the sediment but caused almost complete defaunation at dredging stations. Yet, just 12 months after the extraction, the recolonisation of communities at the impacted stations was at an advanced stage. Unlike other studies on the effects of extraction of marine sand, no significant settlement of opportunistic species was observed. The limited impact of the sand extraction operation on the physical characteristics of the sediment and hydrological-sedimentary characteristics in the relict sand area should aid its rapid recovery and the restoration of the original community in a short period of time (2-4 years after dredging).     

Variability of water mass properties in the last two decades in the South Adriatic Sea with emphasis on the period 2006-2009

Spatially averaged temperature and salinity profiles from individual cruises between 1990 and 2009 were analysed to outline the temporal evolution of water mass properties in the deep convection site in the South Adriatic Pit (SAP). The long-term variability in thermohaline conditions has been explained and related to a close feedback mechanism between the Ionian and the Adriatic Sea. Prominent influences of the Eastern Mediterranean Transient are manifested in changes in the vertical temperature and salinity patterns in the South Adriatic, and the whole studied period was divided into three stages according to the main thermohaline characteristics: 1990-1995, 1995-2004 and the last period from 2005 onwards.Particular attention was given to data collected during 2006-2009, which permitted us to situate the actual thermohaline properties in the context of the decadal variability. This last period was characterised by a very low production of dense water in the northern basin during 2007, while from winter 2008 high production of North Adriatic Deep Water (NAdDW) and Adriatic Deep Water (AdDW) in the northern and southern basins, respectively, was observed. Finally, we used the Optimum Multiparameter Analysis (OMP) to identify the percentages of the different water masses contained in the SAP, and this highlighted some differences between two recent periods studied (2007 and 2008) and the production of dense waters.     

Assessment of meteorological climate models as inputs for coastal studies

Modeling studies of future changes in coastal hydrodynamics, in terms of storm surges and wave climate, need appropriate wind and atmospheric forcings, a necessary requirement for the realistic reproduction of the statistics and the resolution of small scale features. This work compares meteorological results from different climate models in the Mediterranean area, with a focus on the Adriatic Sea, in order to assess their capability to reproduce coastal meteorological features and their possibility to be used as forcings for hydrodynamic simulations. Five meteorological datasets are considered. They are obtained from two regional climate models, implemented with different spatial resolutions and setups and are downscaled from two different global climate models. Wind and atmospheric pressure fields are compared with measurements at four stations along the Italian Adriatic coast. The analysis is carried out both on simulations of the control period 1960–1990 and on the A1B Intergovernmental Panel for Climate Change scenario projections (2070–2100), highlighting the ability of each model in reproducing the statistical coastal meteorological behavior and possible changes. The importance of simulated global- and regional-scale meteorological processes, in terms of correct spatial resolution of the phenomena, is also discussed. Within the Adriatic Sea, the meteorological climate is influenced by the local orography that controls the strengthening of north-eastern katabatic winds like Bora. Results show indeed that the increase in spatial resolution provides a more realistic wind forcing for the hydrodynamic simulations. Moreover, the chosen setup and the global climate models that drive the regional downscalings appear to play an important role in reproducing correct atmospheric pressure fields. The comparison between scenario and control simulations shows a small increase in the mean atmospheric pressure values, while a decrease in mean wind speed and in extreme wind events is observed, particularly for the datasets with higher spatial resolution. Finally, results suggest that an ensemble of downscaled climate models is likely to provide the most suitable climatic forcings (wind and atmospheric pressure fields) for coastal hydrodynamic modeling.     

Elementary study of physical properties on surface sediments in Nansha Islands sea area

Elementary studies on surface sediment physical properties (physical-mechanic property, engineering features, acoustic parameters) of the Nansha Islands sea area in the South China Sea reveal that the seabed sedimentary environments are important factors of controlling sediment physical properties and the relations between them affect one another. At the same time, the relative figures of sediment physical properties and some experienced relations are given, which are beneficial to the later research and further studies. Project of the Multidisciplinary Oceanographic Team of the Chinese Academy of Sciences to the Nansha Islands in the South China Sea.     

Anthropogenic metal contamination and sapropel imprints in deep Mediterranean sediments

Sediment cores from the deep Balearic basin and the Cretan Sea provide evidence for the accumulation of Cd, Pd and Zn in the top few centimeters of the abyssal Mediterranean sea-bottom. In both cores, 206Pb/207Pb profiles confirm this anthropogenic impact with less radiogenic imprints toward surface sediments. The similarity between excess 210Pb accumulated in the top core and the 210Pb flux suggests that top core metal inventories reasonably reflect long-term atmospheric deposition to the open Mediterranean. Pb inventory in the western core for the past 100 years represents 20-30% of sediment coastal inventories, suggesting that long-term atmospheric deposition determined from coastal areas has to be used cautiously for mass balance calculations in the open Mediterranean. In the deeper section of both cores, Al normalized trace metal profiles suggest diagenetic remobilization of Fe, Mn, Cu and, to a lesser extent, Pb that likely corresponds to sapropel event S1.     

Comparing efficacy of different taxonomic resolutions and surrogates in detecting changes in soft bottom assemblages due to coastal defence structures

Sandy shores on the West coast of the North Adriatic Sea are extensively protected by different types of defence structures to prevent coastal erosion. Coastal defence schemes modify the hydrodynamic regime, the sediment structure and composition thus affecting the benthic assemblages. This study examines the effectiveness in detecting changes in soft bottom assemblages caused by coastal defence structures by using different levels of taxonomic resolution, polychaetes and/or bivalves as surrogates and different data transformations. A synoptic analyses of three datasets of subtidal benthic macrofauna used in studies aimed at assessing the impact of breakwaters along the North Adriatic coast has been done. Analyses of similarities and correlations between distance matrices were done using matrices with different levels of taxonomic resolution, and with polychaetes or bivalves data alone. Lentidium mediterraneum was the most abundant species in all datasets. Its abundance was not consistently related to the presence of defence structures. Moreover, distribution patterns of L. mediterraneum were masking the structure of the whole macrofaunal assemblages. Removal of L. mediterraneum from the datasets allowed the detection of changes in benthic assemblages due to coastal defences. Analyses on different levels of taxonomic resolution showed that the level of family maintained sufficient information to detect the impacts of coastal defence structures on benthic assemblages. Moreover, the outcomes depended on the transformation used. Patterns of distribution of bivalves, used as surrogates, showed low correlations with the patterns of the total macrofaunal species assemblages. Patterns of polychaetes, if identified to the species or genus level showed higher correlations with the whole dataset. However, the identification of polychaetes to species and genus level is as costly as the identification of all macrobenthic taxa at family level.This study provided additional evidences that taxonomic sufficiency is a useful tool in environmental monitoring, also in investigations on the impacts of coastal defence structures on subtidal macrofauna. The use of coarser taxonomic level, being time-efficient, would allow improving sampling designs of monitoring programs by increasing replication in space and time and by allowing long term monitoring studies.     

Seabed erosion and deposition related to the typhoon activity of the past millennium on the southeast coast of China

Catastrophic events often interrupt long-term Earth surface processes. In coastal areas, although millennial-scale trends of climatic and sea-level changes determine the trajectory of sedimentary landform evolution, storm and/or tsunami activity can cause abrupt changes in depositional conditions that may alter the long-term sedimentary processes. Here, we report a sedimentary hiatus that is widely observed from the late Holocene sedimentary sequence at the seabed along the southeast China coast. This hiatus was discovered by close temporal sedimentary and radiocarbon dating analyses of a seabed sedimentary sequence in the mouth area of the Pearl River estuary. The results suggest that a couple of metres in the middle to late Holocene sediment at the seabed were eroded by a catastrophic event happening c. 1000–800 cal. Years BP. In theory, a mega-tsunami generated from the Manila Trench could have caused such erosion at the seabed, but there is a lack of direct geological and historical evidence to support such a hypothesis. Much more likely, a super-typhoon struck the coast and caused the erosion. This hypothesis is strongly supported by the region's historical and geological records, which suggest a period characteristic of intense typhoons ranging from the Medieval Warm Period to the climate transition phase (c. 1000–600 cal. Years BP). During the subsequent Little Ice Age, deposition of sandy sediment continued, suggesting frequent but weaker typhoon activity. Over the past two centuries the deposition of sandy sediment and gravels began, implying the beginning of a phase of intensifying typhoon conditions, possibly a result of the recent warming climate. © 2020 John Wiley & Sons, Ltd.     

Benthic community response to habitat variation: A case of study from a natural protected area,the Celestun coastal lagoon

Little information currently exists on spatial and temporal benthic community variations in tropical coastal lagoons. Here, the benthic community response to habitat variation in the Celestun coastal lagoon, northwest Yucatan peninsula, was seasonally examined during the 1994–1995 climatic cycle into a grid of 12 sampling sites distributed along the salinity gradient of the lagoon. Habitat variation was assessed through physical factors associated both to the water column (e.g. salinity) and the bottom sediment (e.g. sand, silt and clay fractions). The benthic community response was assessed through species diversity measures and abundance. Under the influence of climatic seasonality, variations in habitat conditions followed by changes in the benthic community characteristics were expected. Results from two-way ANOVAs showed that for the period of study, Celestun lagoon was more heterogeneous along the spatial axis of variability than along the temporal one. Multiple regression analysis showed that salinity was spatially the main factor influencing the benthic community characteristics. Temporally, the sediment characteristics were observed to exert significant effects on the species diversity characteristics but not on abundance. Other variables assessed (dissolved oxygen, pH, temperature and water column transparency) exhibited no significant covariance with species diversity and abundance. Since generated from historical data, these results have the potential to be useful as a benchmark to the establishment of monitoring programs in the light of the increasing anthropogenic pressure on the natural resources of the lagoon and surrounding coastal area.     

Measurement of Enzymatic Activity in Deep Marine Sediments Using the Semiquantitative APIZym Test System Messung enzymatischer Aktivität in tiefen Meeressedimenten mit Hilfe des APIZym-Testsystems

APIZym strips were used for measuring the activity of 19 enzymatic reactions in marine sediment samples from the continental slope of the Celtic Sea in 135…1630 m depth. The activities of phosphatases were dominant followed by esterases, peptidases, and glycosidases. The activities were influenced by the sediment depth within the cores and the depth of the station. The results indicate that the APIZym system is able to resolve vertical stratification and spatial variability of enzymes possessing activity rates exceeding 10 nmol substrate h⁻¹ cm⁻³ in sediment but not of those with lower rates which for example would be found in deep-sea sediments.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.     

Satellite-derived variability of the Aegean Sea ecohydrodynamics

Eight years of AVHRR-derived sea surface temperature (SST) and SeaWiFS-derived surface chlorophyll (Chl) data (1998–2005) are used to investigate key processes affecting the spatial and temporal variability of the two parameters in the Aegean Sea. Seasonal mean SST and Chl maps are constructed using daily data to study seasonal dynamics whereas empirical orthogonal function (EOF) and correlational analysis is applied to the 8-day composite SST and Chl anomaly time-series in order to study the variability and co-variability of the two parameters from subseasonal to interannual time-scales. The seasonal mean fields show that Black Sea cold and chlorophyll-rich waters enter through the Dardanelles Strait and they are accumulated in the north-eastern part of the Aegean Sea, steered by the Samothraki anticyclone. Large chlorophyll concentrations are encountered in the hydrological front off the Dardanelles Strait as well as in coastal areas affected by large riverine/anthropogenic nutrient loads. The SST seasonal mean patterns reveal strong cooling that is associated with upwelling along the eastern boundary of the basin during summer due to strong northerly winds, a process which is not present in the surface chlorophyll climatology. The Chl dataset presents much stronger sub-seasonal variability than SST, with large variations in the phase and strength of the phytoplankton seasonal cycles. EOF analysis of the anomaly time-series shows that SST non-seasonal variability is controlled by synoptic weather variations and anomalies in the north–south wind-stress component regulating the summer coastal upwelling regime. Mean SST and Chl patterns, and their associated variations, are not closely linked implying that Black Sea and riverine inputs mainly control the intra-annual and interannual variability of the surface chlorophyll in the Aegean Sea rather than mixing and/or upwelling processes.