Observational reconstructions indicate a contemporary increase in coastal ocean CO₂ uptake. However, the mechanisms and their relative importance in driving this globally intensifying absorption remain unclear. Here we integrate coastal carbon dynamics in a global model via regional grid refinement and enhanced process representation. We find that the increasing coastal CO₂ sink is primarily driven by biological responses to climate-induced changes in circulation (36%) and increasing riverine nutrient loads (23%), together exceeding the ocean CO₂ solubility pump (41%). The riverine impact is mediated by enhanced export of organic carbon across the shelf break, thereby adding to the carbon enrichment of the open ocean. The contribution of biological carbon fixation increases as the seawater capacity to hold CO₂ decreases under continuous climate change and ocean acidification. Our seamless coastal ocean integration advances carbon cycle model realism, which is relevant for addressing impacts of climate change mitigation efforts.

Bottom trawling represents the most widespread anthropogenic physical disturbance to seafloor sediments on continental shelves. While trawling-induced changes to benthic ecology have been widely recognized, the impacts on long-term organic carbon storage in marine sediments remains uncertain. Here we combined datasets of sediment and bottom trawling for a heavily trawled region, the North Sea, to explore their potential mutual dependency. A pattern emerges when comparing the surface sediment organic carbon-to-mud ratio with the trawling intensity represented by the multi-year averaged swept area ratio. The organic carbon-to-mud ratio exhibits a systematic response to trawling where the swept area ratio is larger than 1 yr⁻¹. Three-dimensional physical–biogeochemical simulation results suggest that the observed pattern is attributed to the correlated dynamics of mud and organic carbon during transport and redeposition in response to trawling. Both gain and loss of sedimentary organic carbon may occur in weakly trawled areas, whereas a net reduction of sedimentary organic carbon is found in intensely trawled grounds. Cessation of trawling allows restoration of sedimentary carbon stock and benthic biomass, but their recovery occurs at different timescales. Our results point out a need for management of intensely trawled grounds to enhance the CO₂ sequestration capacity in shelf seas.

Coupled physical–biological models usually resolve only parts of the trophic food chain; hence, they run the risk of neglecting relevant ecosystem processes. Additionally, this imposes a closure term problem at the respective “ends” of the trophic levels considered. In this study, we aim to understand how the implementation of higher trophic levels in a nutrient–phytoplankton–zooplankton–detritus (NPZD) model affects the simulated response of the ecosystem using a consistent NPZD–fish modelling approach (ECOSMO E2E) in the combined North Sea–Baltic Sea system. Utilising this approach, we addressed the above-mentioned closure term problem in lower trophic ecosystem modelling at a very low computational cost; thus, we provide an efficient method that requires very little data to obtain spatially and temporally dynamic zooplankton mortality. On the basis of the ECOSMO II coupled ecosystem model we implemented one functional group that represented fish and one group that represented macrobenthos in the 3-D model formulation. Both groups were linked to the lower trophic levels and to each other via predator–prey relationships, which allowed for the investigation of both bottom-up processes and top-down mechanisms in the trophic chain of the North Sea–Baltic Sea ecosystem. Model results for a 10-year-long simulation period (1980–1989) were analysed and discussed with respect to the observed patterns. To understand the impact of the newly implemented functional groups for the simulated ecosystem response, we compared the performance of the ECOSMO E2E to that of a respective truncated NPZD model (ECOSMO II) applied to the same time period. Additionally, we performed scenario tests to analyse the new role of the zooplankton mortality closure term in the truncated NPZD and the fish mortality term in the end-to-end model, which summarises the pressure imposed on the system by fisheries and mortality imposed by apex predators. We found that the model-simulated macrobenthos and fish spatial and seasonal patterns agree well with current system understanding. Considering a dynamic fish component in the ecosystem model resulted in slightly improved model performance with respect to the representation of spatial and temporal variations in nutrients, changes in modelled plankton seasonality, and nutrient profiles. Model sensitivity scenarios showed that changes in the zooplankton mortality parameter are transferred up and down the trophic chain with little attenuation of the signal, whereas major changes in fish mortality and fish biomass cascade down the food chain.

Here we present results from a long-term model simulation of the 3-D coupled ecosystem model ECOSMO II for a North Sea and Baltic Sea set-up. The model allows both multi-decadal hindcast simulation of the marine system and specific process studies under controlled environmental conditions. Model results have been analysed with respect to long-term multi-decadal variability in both physical and biological parameters with the help of empirical orthogonal function (EOF) analysis. The analysis of a 61-year (1948–2008) hindcast reveals a quasi-decadal variation in salinity, temperature and current fields in the North Sea in addition to singular events of major changes during restricted time frames. These changes in hydrodynamic variables were found to be associated with changes in ecosystem productivity that are temporally aligned with the timing of reported regime shifts in the areas. Our results clearly indicate that for analysing ecosystem productivity, spatially explicit methods are indispensable. Especially in the North Sea, a correlation analysis between atmospheric forcing and primary production (PP) reveals significant correlations between PP and the North Atlantic Oscillation (NAO) and wind forcing for the central part of the region, while the Atlantic Multi-decadal Oscillation (AMO) and air temperature are correlated to long-term changes in PP in the southern North Sea frontal areas. Since correlations cannot serve to identify causal relationship, we performed scenario model runs perturbing the temporal variability in forcing condition to emphasize specifically the role of solar radiation, wind and eutrophication. The results revealed that, although all parameters are relevant for the magnitude of PP in the North Sea and Baltic Sea, the dominant impact on long-term variability and major shifts in ecosystem productivity was introduced by modulations of the wind fields.

The north Atlantic subpolar gyre (SPG) has been widely implicated as the source of large-scale changes in the subpolar marine environment. However, inconsistencies between indices of SPG-strength have raised questions about the active role SPG-strength and size play in determining water properties in the eastern subpolar North Atlantic (ENA). Here, by analyzing various SPG indices derived from observations and a global coupled model, we show that the choice of the SPG index dictates the interpretation of SPG strength-salinity relationship in the ENA. Variability in geostrophic currents derived from observed hydrography and model based Lagrangian trajectories reveal zonal shifts of advective pathways in the ENA and meridional shifts in the western intergyre region. Such shifts in advective pathways are manifestations of variability in the size and strength of the SPG, and they impact salinity by modulating the proportion of subpolar and subtropical waters reaching the ENA. SPG indices based on subsurface density and principal component analysis of sea surface height variability capture these shifts in advective pathways, and are therefore best suited to describe SPG-salinity relationship in the ENA. Our results establish the dynamical constraints on the choice of the SPG index and emphasize that SPG indices should be cautiously interpreted.

Simulations of a Regional Climate Model (RCM) driven by identical lateral boundary conditions but initialized at different times exhibit the phenomenon of so-called internal model variability (or in short, Internal Variability—IV), which is defined as the inter-member spread between members in an ensemble of simulations. Our study investigates the effects of air-sea coupling on IV of the regional atmospheric model COSMO-CLM (CCLM) of the new regional coupled system model GCOAST-AHOI (Geesthacht Coupled cOAstal model SysTem: Atmosphere, Hydrology, Ocean and Sea Ice). We specifically address physical processes parameterized in CCLM, which may cause a large IV during an extreme event, and where this IV is affected by the air-sea coupling. Two six-member ensemble simulations were conducted with GCOAST-AHOI and the stand-alone CCLM (CCLM_ctr) for a period of 1 September–31 December 2013 over Europe. IV is expressed by spreads within the two sets of ensembles. Analyses focus on specific events during this period, especially on the storm Christian occurring from 27 to 29 October 2013 in northern Europe. Results show that simulations of CCLM_ctr vary largely amongst ensemble members during the storm. By analyzing two members of CCLM_ctr with opposite behaviors, we found that the large uncertainty in CCLM_ctr is caused by a combination of two factors (1) uncertainty in parameterization of cloud-radiation interaction in the atmospheric model. and (2) lack of an active two-way air-sea interaction. When CCLM is two-way coupled with the ocean model, the ensemble means of GCOAST-AHOI and CCLM_ctr are relatively similar, but the spread is reduced remarkably in GCOAST-AHOI, not only over the ocean where the coupling is done but also over land due to the land-sea interactions.