Oceans

There is increasing interest in the role that seagrasses play in storing carbon in the context of climate mitigation, but many knowledge gaps in the factors controlling this storage exist. Here, we provide a small case study that examines the role of infaunal biodiversity in influencing seagrass and the carbon stored in its sediments. A total of 25 species of invertebrate were recorded in an intertidal Zostera marina meadow, where these species were dominated by polychaete worms with no bivalves present. We find organic carbon storage (within the top 20 cm) measured by AFDW to be highly variable within a small area of seagrass meadow ranging from 2961 gC.m⁻² to 11,620 gC.m⁻² with an average (±sd) of 6460² ± 3274 gC.m⁻². Our analysis indicates that infaunal communities are significantly and negatively correlated with this sediment organic carbon. However, this effect is not as influential as hypothesised, and the relatively small sample size of the present study limits its ability to provide strong causality. Other factors, such as algal abundance, curiously had a potentially stronger influence on the carbon in the upper sediments. The increasing richness of infauna is likely reducing the build-up of organic carbon, reducing its ecosystem service role. We believe this to likely be the result of bioturbation by specific species such as Arenicola marina and Ampharete acutifrons. A change in sediment organic carbon suggests that these species could be key drivers of bioturbator-initiated redox-driven organic matter turnovers, influencing the microbial processes and remobilizing sediment compounds. Bioturbators should be considered as a limitation to C_org storage when managing seagrass C_org stocks; however, bioturbation is a natural process that can be moderated when an ecosystem is less influenced by anthropogenic change. The present study only provides small-scale correlative evidence with a range of surprising results; confirming these results within temperate seagrasses requires examining this process at large spatial scales or with targeted experiments. Full article

Perfluorooctanoic acid (PFOA) has been widely studied due to its environmental persistence and bioaccumulation potential, raising concerns about its effects on aquatic life. This research evaluates the impact of PFOA on the antioxidant defenses and stress response systems of the mussel Mytilus galloprovincialis. Mussels were exposed to three concentrations of PFOA (1, 10, and 100 µg·L⁻¹) over 28 days. Several biomarkers, including glutathione S-transferase (GST), superoxide dismutase (SOD), catalase (CAT), lipid peroxidation (LPO), total antioxidant capacity (TAC), vitellogenin (VTG), ubiquitin (UBI), and caspase-3 (CASP) were analyzed. The results suggest stress responses, particularly in animals exposed to higher concentrations, as shown by GST and SOD activities which increased according to PFOA concentrations. Additionally, oxidative stress markers such as MDA and CAT showed variable responses depending on the exposure concentration tested. This study underscores the need for further investigation into the effects of PFOA on mollusks but also the need to unveil gender-specific responses in aquatic organisms exposed to this contaminant. The concentrations of PFOA used in our research are lower than those examined in previous studies, providing crucial insights into the impacts of even minimal exposure levels. It highlights the potential of M. galloprovincialis as a bioindicator in environmental monitoring programs, providing crucial insights for environmental management and policymaking regarding regulating and monitoring PFOA in marine settings. Consequently, in a country where seafood consumption is the second largest in Europe, implementing environmental policies and regulatory measures to manage and monitor PFOA levels in marine environments is crucial. Full article

The current study looks at how the characteristics of Arabian Sea tropical cyclones (TCs) change over time. The results show that in the pre-monsoon (April–June) and the post-monsoon (October–December), the activity of TCs > 34 knots, including cyclonic storm (CS), severe cyclonic storm (SCS), very severe cyclonic storm (VSCS), extreme severe cyclonic storm (ESCS), and super cyclonic storm (Sup. CS), has significantly increased, while the tendency of TCs < 34 knots, depressions and deep depressions (Ds) over the Arabian Sea has only slightly increased. Most of the TC activity in the first two decades (1982–2001) over the Arabian Sea activated on the eastern side, while in the last two decades (2002–2021), there was an expansion toward the southwest region of the Arabian Sea, especially in the post-monsoon season. The composite analysis of environmental parameters over the Arabian Sea reveals that the negative anomalies of outgoing longwave radiation (OLR) and the positive anomalies of relative humidity at 500 hPa (RH–500 hPa) in the first decade (1982–1991) and the second decade (1992–2001) are more concentrated on the eastern side of the Arabian Sea, leading to increased activity for TCs. Decades three (2002–2011) and four (2012–2021) demonstrated a wide distribution of weak vertical wind shear (VWS) and strong convection (OLR and RH–500 hPa) over the Arabian Sea basin. This led to TCs occurring more frequently and stronger, especially in the post-monsoon season. SST over the Arabian Sea was sufficient for tropical storm activity (≥26.5 °C) for both typical seasons. Full article

Invasive marine invertebrates are increasingly recognized as a potential disturbance to coastal ecosystems. We sought to better document the taxonomic composition of subtidal communities around Long Island to obtain a baseline that can be used to monitor current and future invasions of non-indigenous species. We placed settlement blocks at 18 sites along the coast of Long Island, New York, for three months. After recovering blocks at 12 sites, we analyzed the taxonomic composition of fouling communities on the blocks. We observed 64 invertebrate and 3 algal taxa, with large variation in taxon richness among sites. Multivariate analyses revealed that although taxon composition was significantly dissimilar between north and south shores, variation in dissimilarity did not differ significantly between shores. The high variability in taxon composition observed among sites indicates that additional research is needed to expand our knowledge of invertebrate diversity in the waters surrounding Long Island. Adding more sites and replicate blocks within sites could improve future sampling designs. This research will benefit continuing efforts to monitor, manage, and prevent the establishment of marine invasive species. Full article

Understanding the ecological dynamics of zooplankton communities is crucial to precisely assessing the health of marine ecosystems and their management. Metabarcoding has contributed to a better understanding of biodiversity in marine environments. However, this methodology still requires protocol optimisation. Here, we used a complementary approach combining molecular and morphological identification methods to identify the zooplankton community inhabiting the Berlengas Archipelago, Portugal. The presence of non-indigenous species was also assessed. The results showed that the metabarcoding approach outperformed the classical morphological identification method, detecting more species with higher resolution. Nevertheless, the classical method was able to identify species that were not detected by the molecular approaches, probably due to a lack of reference data in the databases. The comparison between different molecular approaches showed that COI and bulk DNA gave better results than 18S rRNA and eDNA by detecting higher species diversity. However, complementarities were observed between them. Molecular tools also proved effective in identifying several potential non-indigenous species, identifying, for the first time, several potentially unreported NIS inhabiting the Portuguese marine ecosystems. Overall, our results confirmed the importance of combining both classical and molecular methods to obtain a more refined assessment of the zooplanktonic communities in marine environments. Full article

In Cartesian coordinates $\left(x,y,z\right),$ the gradient Richardson number $Ri$ is the ratio between the square of the buoyancy frequency N and the square of the vertical shear S, $Ri=N^2/S^2$, where $N^2=-\left(g/\rho \right)\partial \rho /\partial z$ and $S^2={\left(\partial u/\partial z\right)}^2+{\left(\partial v/\partial z\right)}^2$, with ρ potential density, $\left(u,v\right)$ the horizontal velocity components and g gravity acceleration. In isopycnic coordinates $\left(x,y,\rho \right)$, $Ri$ is expressed as the ratio between ${M^2\equiv N}^{-2}$ and the squared diapycnal shear ${S_{\rho }}^2={\left(\rho /g\right)}^2\left[{\left(\partial u/\partial \rho \right)}^2+{\left(\partial v/\partial \rho \right)}^2\right]$, $Ri=M^2/{S_{\rho }}^2$. This could suggest that a decrease (increase) in stratification brings a decrease (increase) in dynamic stability in Cartesian coordinates, but a stability increase (decrease) in isopycnic coordinates. The apparently different role of stratification arises because S and $S_{\rho }$ are related through the stratification itself, $S_{\rho }=S/N^2$. In terms of characteristic times, this is equivalent to $\tau \equiv S_{\rho }={t_o}^2/t_d$, which is interpreted as a critical dynamic time $\tau$ that equals the buoyancy period $t_o\equiv N^{-1}$ normalized by the ratio $t_d/t_o$, where $t_d=S^{-1}$ is the deformation time. Here we follow simple arguments and use field data from three different regions (island shelf break, Gulf Stream and Mediterranean outflow) to endorse the usefulness of the isopycnal approach. In particular, we define the reduced squared diapycnal shear ${{\sigma }_{\rho }}^2={S_{\rho }}^2-M^2$ and compare it with the reduced squared vertical ${\sigma }^2=S^2-N^2$, both being positive (negative) for unstable (stable) conditions. While both Ri and ${\sigma }^2$ remain highly variable for all stratification conditions, the mean ${{\sigma }_{\rho }}^2$ values approach ${S_{\rho }}^2$ with increasing stratification. Further, the field data follow the relation ${{\sigma }_{\rho }}^2=\left(1-Ri\right)/\left(N^{2}Ri\right)$, with a subcritical $Ri=0.22$ for both the island shelf break and the Mediterranean outflow. We propose ${{\sigma }_{\rho }}^2$ and ${S_{\rho }}^2$ to be good indexes for the occurrence of effective mixing under highly stratified conditions. Full article

Recent years have seen a notable rise in dolphin-watching boat activities along the Algarve coast in Portugal, potentially affecting the common dolphin (Delphinus delphis) and bottlenose dolphin (Tursiops truncatus) local populations. This study examines the impact of increasing underwater noise levels from these boats on dolphin vocalizations. Field recordings were conducted from June to September 2022, analyzing dolphin whistles in various boat presence scenarios. The results indicate significant changes in whistle-frequency characteristics with boat presence, including increased start, low, and high frequencies, alongside a decrease in the number of inflection points in modulated whistles. The changes might negatively impact dolphin populations viability, underscoring the need for further research. Additionally, improved mitigation strategies may be necessary to reduce the potential negative effects of dolphin watching on cetacean communication and behavior in the Algarve region. Full article

Global climate change is profoundly impacting coral ecosystems. Rising sea surface temperatures, in particular, disrupt coral reproductive synchrony, cause bleaching, and mortality. Oculina patagonica, a temperate scleractinian coral abundant across the Mediterranean Sea, can grow at a temperature range of 10–31 °C. Studies conducted three decades ago documented this species bleaching during the summer months, the same time as its gonads mature. However, the Eastern Mediterranean Sea is experiencing some of the fastest-warming sea surface temperatures worldwide. This study repeated the year-round in situ assessment of the reproductive cycle and gonad development and correlation to summer bleaching. In addition, thermal performance of the holobiont was assessed in an ex situ thermal stress experiment. In situ monitoring revealed no temporal changes in gonad development compared to previous studies, despite sea surface warming and concurrent bleaching. Experimental thermal performance curves indicated that photosynthetic rate peaked at 23 °C, bleached coral area was significant at 29 °C, and peaked at 34 °C. With local sea surface temperature reaching 31 °C, O. patagonica is exposed beyond its bleaching threshold during the summer months in situ. Despite this, O. patagonica maintains gonad development and physiologically recovers at the end of summer demonstrating resilience to current warming trends. Full article

Wind-induced waves can lead to the partial or complete wash-over of beaches, causing erosion that impacts both the landscape and tourist infrastructure. In some regions of the world, e.g., Croatia, this process, which usually occurs during a harsh winter, has a major impact on the environment and the economy, and preventing or reducing this process is highly desirable. One of the simplest methods to reduce or prevent beach erosion is the use of innovative underwater structures designed to decrease wave energy by reducing wave height. In this study, submerged breakwaters are numerically investigated using various topologies, positions, and angles relative to the free surface. Not only is the optimal topology determined, but the most efficient arrangement of multiple breakwaters is also determined. The advantage of newly developed submerged breakwaters over traditional ones (rock-fixed piers) is that they do not require complex construction, massive foundations, or high investment costs. Instead, they comprise simple floating bodies connected to the seabed by mooring lines. This design makes them not only cheap, adaptable, and easy to install but also environmentally friendly, as they have little impact on the seabed and the environment. To evaluate wave damping effectiveness, the incompressible computational fluid dynamics (ICFD) method is used, which enables the use of a turbulence model and the possibility of accurate wave modelling. Full article

This research focuses on integrating seakeeping indices into the trajectory planning of a mother ship in order to minimize risks during UAV (unmanned aerial vehicle) takeoff and landing in challenging sea conditions. By considering vessel dynamics and environmental factors, the proposed trajectory planning algorithm computes optimal paths that prioritize the stability and safety of the ship, mitigating the impact of adverse weather on UAV operations. Specifically, the new adaptive weather routing model presented is based on a genetic algorithm. The model uses the previously evaluated response amplitude operators (RAOs) for the reference ship at different velocities and encounter angles, along with weather forecast data provided by the global wave model (GWAM). Preliminary evaluations confirm the effectiveness of the presented model in significantly improving the reliability of autonomous UAV operations from a mother ship across all encountered sea state conditions, particularly when compared with a graph-based solution. The current results clearly demonstrate that it is possible to achieve appreciable improvements in ship seakeeping performance, thereby making UAV-related operations safer. Full article

The Magnisi peninsula is a small portion of land located near the largest Italian petrol-chemical pole of Augusta–Priolo–Melilli (40 km²), which, since the 1950s, devastated the local environment and landscape and unloaded directly into the sea an impressive quantity of pollutants. Unlike the terrestrial part of the area, where a Natural Oriented Reserve (NOR) called “Saline di Priolo” was established in the 2000s, no concrete legislative action has been implemented or proposed for the marine environment. At the same time, the fauna of the marine environment has not been studied in the same way as that of the terrestrial environment. As concerns the molluscan fauna, most of the information dates back to the 1800s. These studies were exclusively focused on the shells of some mollusks in the area. Instead, no study related to this area has ever been carried out on the group of sea slugs. This study conducted between 2022 and 2023, through snorkeling activities, allowed to provide the first faunistic list of the sea slugs of this area, together with information on the biology and ecology of these gastropods, highlighting the potential biodiversity present in this small stretch of coastline affected by high industrial pollution. Full article

The accumulation of sea-surface debris around the coastal waters of Malta poses significant ecological and environmental challenges, negatively affecting marine ecosystems and human activities. This issue is exacerbated due to the lack of an effective system tailored to predict surface-debris movement specifically for the Islands of Malta. To address this gap, a pipeline that combines a machine learning-based prediction system with a physics-based model is proposed. This pipeline uses data on historical sea-surface current velocities to forecast future conditions and visualise debris movement. Central to this system are two machine learning models trained to predict surface velocities for the next 24 h for a specific area. These predictions are then utilised in a Lagrangian model to simulate and visualise the debris movement, providing insights into future dispersion patterns. A comparative evaluation of both models using real-world data is made to determine which one performs best in this application. This method offers a tailored approach to addressing sea-surface debris around Malta by accurately predicting sea-surface current velocities and visualising debris movement, improving cleanup operations and marine conservation strategies. Full article

The expansion of the invasive Asian macroalgae Rugulopteryx okamurae along the coasts of the Azores represents a significant challenge for local marine biodiversity. A promising approach to managing the biomass produced by this alien alga is to valorize it in the context of the blue economy. This study characterizes and evaluates the potential of R. okamurae biomass for incorporation into cattle feed, with a focus on mitigating enteric methane production. The nutritional value of R. okamurae, its digestibility, and its potential as a mitigating agent for enteric methane production were analyzed in vitro. The results indicate that the inclusion of 5% R. okamurae in the diet significantly (p < 0.05) reduced accumulated methane production by 98% after 24 h of incubation. The addition of 1% algae over the same period resulted in a 38% reduction in methane production. However, a significant decrease (p < 0.05) in gas production of 57.02% and 73.5% was also observed in relation to control, with the inclusion of 1% and 5%, respectively, during 96 h. Nutritionally, R. okamurae was found to have a crude protein content of 18.68% and fiber (NDF) of 55.71% of DM. It is also worth highlighting the high content of ash (31.86%) that was identified in these brown macroalgae. In conclusion, the fresh biomass of R. okamurae could serve as a functional ingredient in cattle feed to mitigate enteric methane production, provided it is used in low percentages. However, it is important to emphasize that high concentrations in the first 12 h did not produce methane, which is also not recommended for enteric fermentation. However, before including it in animal feed, in vivo tests are needed to assess its toxicity. Full article

This paper assessed the phytoplankton structure and its relationship with the physical environment in the coastal region off Mazatlán, Mexico, in two seasons of 2022, a year in which a strong La Niña event took place: (1) the warmer (August) and (2) the transitional period to the cold phase (November), based on hydrographic data and samples collected in two systematic scientific expeditions. The results showed clear differences between both seasons. Regarding total abundance, August reached 125,200 cells L⁻¹, while November amounted to 219,900 cells L⁻¹. Regarding species composition, the diatoms Cylindrotheca closterium and Planktoniella sol were dominant in August, while Thalassionema nitzschioides and Tetramphora decussata dominated the assemblages in November. The dinoflagellate Protoperidinium punctulatum was dominant in both seasons. However, very marked differences in its abundance are reported. The differences observed in the species richness and abundance could be attributed to the physical configuration of the water column, particularly the surface temperature, which showed clear changes between both seasons. The results presented here confirmed the high phytoplankton richness (some of them with the potential to generate harmful algal blooms), abundance, and diversity values of the region, suggesting a strong relationship with the physical environment. Full article

This study investigates the temporal trends and correlations between Saharan dust mass concentration densities (DMCD) and Sargassum concentrations (SCT) in the tropical North Atlantic. Average DMCD data for June, July, and August from 1980 to 2022, alongside SCT data for the same months from 2012 to 2022, were analyzed using Mann–Kendall tests for trends and lagged regression models to assess whether higher Saharan dust levels correlate with Sargassum outbreaks in the region. A comprehensive analysis reveals a significant upward trend in Saharan dust quantities over the study period, with the summer months of June, July, and August exhibiting consistent increases. Notably, 2018 and 2020 recorded the highest mean DMCD levels, with June showing the most significant increasing trend, peaking in 2019. These findings are consistent with previous studies indicating a continuous elevation in Saharan dust concentrations in the tropical atmosphere of the North Atlantic. Simultaneously, Sargassum concentrations also show a notable increasing trend, particularly in 2018, which experienced both peak SCT and elevated DMCD levels. Mann–Kendall tests confirm statistically significant upward trends in both Saharan dust and Sargassum concentrations. Simple linear regression and lagged regression analyses reveal positive correlations between DMCD and SCT, highlighting a temporal component with stronger associations observed in July and the overall June–July–August (JJA) period. These results underscore the potential contribution of elevated Saharan dust concentrations to the recent surge in Sargassum outbreaks in the tropical North Atlantic. Furthermore, the results from forward stepwise regression (FSR) models indicate that DMCD and chlorophyll (CHLO) are the most critical predictors of SCT for the summer months, while sea surface temperature (SST) was not a significant predictor. These findings emphasize the importance of monitoring Saharan dust and chlorophyll trends in the Eastern Caribbean, as both factors are essential for improving Sargassum modeling and prediction in the region. This study provides valuable insights into the climatic factors influencing marine ecosystems and highlights the need for integrated environmental monitoring to manage the impacts on coastal economies. Full article

The formation of dense water in the Aegean Sea is important as it affects the deep circulation and the hydrography of the Eastern Mediterranean Sea. In this study, the variability of dense water formation is investigated in relation to forcing mechanisms from 1947 to 2023 in the subbasins of the Aegean Sea, utilising in situ observations from various sources, which have been analysed in combination with satellite altimetry and reanalyses products. The analysis reveals that the Aegean Sea has been in a state of increased dense water formation since 2017 due to the combination of increased surface buoyancy loss and reduced Black Sea water inflow. Extremely high salinity has been recorded in the intermediate layers of the Aegean Sea since 2019. The anticyclonic circulation of the North Ionian gyre during 2017 and 2018 probably also contributed to the rapid transport of highly saline waters in the intermediate and, through dense water formation, the deep layers of the Aegean Sea in 2019. Until 2022, the dense waters formed during the peak of the Eastern Mediterranean Transient still occupied the bottom layers of some deep subbasins of the North and South Aegean; however, the $29.4 \mathrm{k}\mathrm{g} {\mathrm{m}}^{-3}$ isopycnal in the North Aegean and the $29.3 \mathrm{k}\mathrm{g} {\mathrm{m}}^{-3}$ isopycnal in the Southeastern Aegean have gradually deepened by 800 m, permitting the waters forming in the last ten years in the Aegean Sea to settle at ever greater depths. Temperature controls the density variability of the Cretan intermediate water up to the decadal time scale. Increased data availability since 2010 was sufficient to clarify that intrusions of dense water from the North–Central Aegean Sea contributed to the erosion of the Eastern Mediterranean transitional waters in the South Aegean Sea after 2017, as well as to raising the intermediate water masses of the South Aegean to shallower depths. The erosion of the transitional Mediterranean waters in the South Aegean Sea between 1947 and 1955 and 1973 and 1980 coincided with increased dense water formation in the North–Central Aegean Sea. During the peak of the Eastern Mediterranean Transient, the North Ionian circulation, the Black Sea water inflow, the Atlantic Multidecadal Oscillation, and the surface buoyancy fluxes favoured dense water formation in the Aegean Sea. Full article

The forty-year-old Long-Term Ecological Research MareChiara (LTER-MC) program started on 26 January 1984, with fortnightly oceanographic sampling until 1991 and then, from 1995, with weekly sampling up to the present time. LTER-MC produced >150 publications that have been cited by thousands of other studies. In this scoping review, we analyzed this corpus using a semantic approach based on topic modeling, a machine-driven procedure to identify and map topics and their interactions. Understanding the causes behind the evolution of scientific topics, their emergence, splitting, hybridization, or merging within a scientific community is an important step in science policy in managing collaborative research and bringing it into the future. Across different topics, mainly represented by studies on Natural History, Biodiversity, Phenology, Life Cycles, and Community Ecology, the LTER-MC work expanded the knowledge on planktonic organisms, describing in detail their lifestyles and delineating their relationships with environmental conditions. In presenting these results, the potential strengths, weaknesses, opportunities, and threats connected to the overall scientific dimension of LTER-MC are discussed. Finally, the upcoming effort is envisioned in reinforcing internal collaboration to integrate basic and applied research around scientific investigations suitable for establishing a stronger interaction between science and policy, as indicated by the United Nations Decade of Ocean Science for Sustainable Development. Full article

The impacts of climate change-related stressors are becoming more noticeable in the ocean, particularly in coastal marine ecosystems. Yet limited information still exists on the physiological state and ecological resilience of marine fish species, especially during their early life stages (i.e., larvae and juveniles). The present study investigated the effects of chronic exposure to seawater warming (OW; ΔT = +4 °C) and acidification (OA; ΔpH = −0.3 pH units, equivalent to pCO₂~1000 µatm), acting alone or combined (OWA), on juvenile zebra seabream (Diplodus cervinus) physiological resilience, considering distinct levels of biological organization (i.e., biochemical, cell, organ and individual levels). After 60 days of exposure, both stressors, in isolation or combination, significantly decreased specific growth rate (−11% in OW, −42% in OA and −49% in OWA) and leukocyte counts (from −29% in OA and OWA up to −37% in OW) in relation to the control treatment. In addition, a decreased Fulton’s condition index (K) was observed under warming and acidification in combination (−35% in OWA). At the cell level, OW, OA and OWA triggered different biomarker responses in D. cervinus (i.e., up-regulation, down-regulation, or absence of significant effect). In general, the results are suggestive of an antagonistic effect when warming and acidification are combined. OWA yielded the highest integrated biomarker response (IBR) index value in the whole organism, muscle, brain and gills of D. cervinus juveniles, therefore suggesting that the effects of these stressors are more severe when they act together. The distinct patterns observed in each stress scenario highlight the importance of carrying out further studies adjusted to the specificities of different regions, i.e., accounting not only for the type and degree of severity of environmental stressors already felt and/or projected for that specific area, but also the physiological plasticity of species that inhabit a particular ecosystem. The gathered knowledge will allow one to determine the vulnerability of particular marine species and geographic areas and, most importantly, to draw up effective and tailor-made conservation strategies to overcome climate change impacts. Full article

Shipping is a vital component of global trade. However, shipping activities have negative environmental impacts, including emissions to air and water. The Baltic Sea is severely affected by eutrophication, mainly due to nutrient inputs, especially nitrogen and phosphorus. Understanding the contributions of various nutrient sources is essential for informed environmental regulation. The aim of this research is to quantify the nutrient content of black and grey water discharged from cargo ships in the Baltic Sea in order to assess their contribution to the total nutrient load. Specifically, this research examines the nitrogen and phosphorus loadings from ships calling the port of HaminaKotka and addresses key questions regarding nutrient generation, discharge proportions, emission shares by ship type, and their importance compared to other sources. Using a methodology based on ship data and nutrient production estimates, this study found that 2545 cargo ships generated a total of 781 kg of nitrogen and 134 kg of phosphorus in their effluents during their voyages to the port in 2021. However, only a small fraction (0.5%) of the cargo ships discharged their wastewater at port reception facilities. This study concludes that nutrient discharges from cargo ships’ wastewaters are relatively small, contributing less than 0.06% phosphorus and 0.01% nitrogen to the total load in the area. This result indicates that the impact of cargo ships’ wastewater on the eutrophication of the Baltic Sea is smaller than previously thought. The methodology used in this study can be applied in other ports and regions to assess nutrient discharges from cargo ships’ wastewaters on a wider scale. Further research is recommended to assess the effectiveness of on-board treatment systems, the impact of other contaminants in wastewater, and the development of port reception facilities to facilitate proper wastewater management in the maritime trade. Full article