Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Nitrogen cycle feedbacks as a control on euxinia in the mid-Proterozoic ocean

Geochemical evidence invokes anoxic deep oceans until the terminal Neoproterozoic similar to 0.55 Ma, despite oxygenation of Earth's atmosphere nearly 2 Gyr earlier. Marine sediments from the intervening period suggest predominantly ferruginous (anoxic Fe(II)-rich) waters, interspersed with euxinia (anoxic H2S-rich conditions) along productive continental margins. Today, sustained biotic H2S production requires NO3- depletion because denitrifiers outcompete sulphate reducers. Thus, euxinia is rare, only occurring concurrently with (steady state) organic carbon availability when N-2-fixers dominate the production in the photic zone. Here we use a simple box model of a generic Proterozoic coastal upwelling zone to show how these feedbacks caused the mid-Proterozoic ocean to exhibit a spatial/temporal separation between two states: photic zone NO3- with denitrification in lower anoxic waters, and N-2-fixation- driven production overlying euxinia. Interchange between these states likely explains the varying H2S concentration implied by existing data, which persisted until the Neoproterozoic oxygenation event gave rise to modern marine biogeochemistry.

Oceanic controls on the primary production of the northwest European continental shelf: model experiments under recent past conditions and a potential future scenario

In this paper we clearly demonstrate that changes in oceanic nutrients are a first order factor in determining changes in the primary production of the northwest European continental shelf on time scales of 5–10 yr. We present a series of coupled hydrodynamic ecosystem modelling simulations, using the POLCOMS-ERSEM system. These are forced by both reanalysis data and a single example of a coupled ocean-atmosphere general circulation model (OA-GCM) representative of possible conditions in 2080–2100 under an SRES A1B emissions scenario, along with the corresponding present day control. The OA-GCM forced simulations show a substantial reduction in surface nutrients in the open-ocean regions of the model domain, comparing future and present day time-slices. This arises from a large increase in oceanic stratification. Tracer transport experiments identify a substantial fraction of on-shelf water originates from the open-ocean region to the south of the domain, where this increase is largest, and indeed the on-shelf nutrient and primary production are reduced as this water is transported on-shelf. This relationship is confirmed quantitatively by comparing changes in winter nitrate with total annual nitrate uptake. The reduction in primary production by the reduced nutrient transport is mitigated by on-shelf processes relating to temperature, stratification (length of growing season) and recycling. Regions less exposed to ocean-shelf exchange in this model (Celtic Sea, Irish Sea, English Channel, and Southern North Sea) show a modest increase in primary production (of 5–10%) compared with a decrease of 0–20% in the outer shelf, Central and Northern North Sea. These findings are backed up by a boundary condition perturbation experiment and a simple mixing model.

Validation of the NEMO-ERSEM operational ecosystem model for the North West European Continental Shelf

This paper details updates to the Met Office's operational coupled hydrodynamic-ecosystem model from the 7 km Medium-Resolution Continental Shelf – POLCOMS-ERSEM (MRCS-PE) system (Siddorn et al., 2007) to the 7 km Atlantic Margin Model NEMO-ERSEM (AMM7-NE) system. We also provide a validation of the ecosystem component of the new operational system. Comparisons have been made between the model variables and available in situ, satellite and climatological data. The AMM7-NE system has also been benchmarked against the MRCS-PE system. The transition to the new AMM7-NE system was successful and it has been running operationally since March 2012 and has been providing products through MyOcean (http://www.myocean.eu.org) since that time. The results presented herein show the AMM7-NE system performs better than the MRCS-PE system with the most improvement in the model nutrient fields. The problem of nutrient accumulation in the MRCS-PE system appears to be solved in the new AMM7-NE system with nutrient fields improved throughout the domain as discussed in Sect. 4. Improvements in model chlorophyll are also seen but are more modest.

The carbonate system of the NW European shelf: sensitivity and model validation

The ocean plays an important role in regulating the climate, acting as a sink for carbon dioxide, perturbing the carbonate system and resulting in a slow decrease of seawater pH. Understanding the dynamics of the carbonate system in shelf sea regions is necessary to evaluate the impact of Ocean Acidification (OA) in these societally important ecosystems. Complex hydrodynamic and ecosystem coupled models provide a method of capturing the significant heterogeneity of these areas. However rigorous validation is essential to properly assess the reliability of such models. The coupled model POLCOMS–ERSEM has been implemented in the North Western European shelf with a new parameterization for alkalinity explicitly accounting for riverine inputs and the influence of biological processes. The model has been validated in a like with like comparison with North Sea data from the CANOBA dataset. The model shows good to reasonable agreement for the principal variables, physical (temperature and salinity), biogeochemical (nutrients) and carbonate system (dissolved inorganic carbon and total alkalinity), but simulation of the derived variables, pH and pCO2, are not yet fully satisfactory. This high uncertainty is attributed mostly to riverine forcing and primary production. This study suggests that the model is a useful tool to provide information on Ocean Acidification scenarios, but uncertainty on pH and pCO2 needs to be reduced, particularly when impacts of OA on ecosystem functions are included in the model systems.

Potential Climate Change Effects on the Habitat of Antarctic Krill in the Weddell Quadrant of the Southern Ocean

Antarctic krill is a cold water species, an increasingly important fishery resource and a major prey item for many fish, birds and mammals in the Southern Ocean. The fishery and the summer foraging sites of many of these predators are concentrated between 0 degrees and 90 degrees W. Parts of this quadrant have experienced recent localised sea surface warming of up to 0.2 degrees C per decade, and projections suggest that further widespread warming of 0.27 degrees to 1.08 degrees C will occur by the late 21st century. We assessed the potential influence of this projected warming on Antarctic krill habitat with a statistical model that links growth to temperature and chlorophyll concentration. The results divide the quadrant into two zones: a band around the Antarctic Circumpolar Current in which habitat quality is particularly vulnerable to warming, and a southern area which is relatively insensitive. Our analysis suggests that the direct effects of warming could reduce the area of growth habitat by up to 20%. The reduction in growth habitat within the range of predators, such as Antarctic fur seals, that forage from breeding sites on South Georgia could be up to 55%, and the habitat's ability to support Antarctic krill biomass production within this range could be reduced by up to 68%. Sensitivity analysis suggests that the effects of a 50% change in summer chlorophyll concentration could be more significant than the direct effects of warming. A reduction in primary production could lead to further habitat degradation but, even if chlorophyll increased by 50%, projected warming would still cause some degradation of the habitat accessible to predators. While there is considerable uncertainty in these projections, they suggest that future climate change could have a significant negative effect on Antarctic krill growth habitat and, consequently, on Southern Ocean biodiversity and ecosystem services.

Real-time sea surface air temperature retrieval over Southern Ocean using AMSR-E data

The AMSR-E satellite data and in-situ data were applied to retrieve sea surface air temperature (Ta) over the Southern Ocean. The in-situ data were obtained from the 24~(th) -26~(th) Chinese Antarctic Expeditions during 2008-2010. First, Ta was used to analyze the relativity with the bright temperature (Tb) from the twelve channels of AMSR-E, and no high relativity was found between Ta and Tb from any of the channels. The highest relativity was 0.38 (with 23.8 GHz). The dataset for the modeling was obtained by using in-situ data to match up with Tb, and two methods were applied to build the retrieval model. In multi-parameters regression method, the Tbs from 12 channels were used to the model and the region was divided into two parts according to the latitude of 50°S. The retrieval results were compared with the in-situ data. The Root Mean Square Error (RMS) and relativity of high latitude zone were 0.96℃and 0.93, respectively. And those of low latitude zone were 1.29 ℃ and 0.96, respectively. Artificial neural network (ANN) method was applied to retrieve Ta.The RMS and relativity were 1.26 ℃ and 0.98, respectively.

The Ocean Colour Climate Change Initiative: III. A round-robin comparison on in-water bio-optical algorithms

Satellite-derived remote-sensing reflectance (Rrs) can be used for mapping biogeochemically relevant variables, such as the chlorophyll concentration and the Inherent Optical Properties (IOPs) of the water, at global scale for use in climate-change studies. Prior to generating such products, suitable algorithms have to be selected that are appropriate for the purpose. Algorithm selection needs to account for both qualitative and quantitative requirements. In this paper we develop an objective methodology designed to rank the quantitative performance of a suite of bio-optical models. The objective classification is applied using the NASA bio-Optical Marine Algorithm Dataset (NOMAD). Using in situRrs as input to the models, the performance of eleven semi-analytical models, as well as five empirical chlorophyll algorithms and an empirical diffuse attenuation coefficient algorithm, is ranked for spectrally-resolved IOPs, chlorophyll concentration and the diffuse attenuation coefficient at 489 nm. The sensitivity of the objective classification and the uncertainty in the ranking are tested using a Monte-Carlo approach (bootstrapping). Results indicate that the performance of the semi-analytical models varies depending on the product and wavelength of interest. For chlorophyll retrieval, empirical algorithms perform better than semi-analytical models, in general. The performance of these empirical models reflects either their immunity to scale errors or instrument noise in Rrs data, or simply that the data used for model parameterisation were not independent of NOMAD. Nonetheless, uncertainty in the classification suggests that the performance of some semi-analytical algorithms at retrieving chlorophyll is comparable with the empirical algorithms. For phytoplankton absorption at 443 nm, some semi-analytical models also perform with similar accuracy to an empirical model. We discuss the potential biases, limitations and uncertainty in the approach, as well as additional qualitative considerations for algorithm selection for climate-change studies. Our classification has the potential to be routinely implemented, such that the performance of emerging algorithms can be compared with existing algorithms as they become available. In the long-term, such an approach will further aid algorithm development for ocean-colour studies.

Large bio-geographical shifts in the north-eastern Atlantic Ocean: From the subpolar gyre, via plankton, to blue whiting and pilot whales

Pronounced changes in fauna, extending from the English Channel in the south to the Barents Sea in the north-east and off Greenland in the north-west, have occurred in the late 1920s, the late 1960s and again in the late 1990s. We attribute these events to exchanges of subarctic and subtropical water masses in the north-eastern North Atlantic Ocean, associated with changes in the strength and extent of the subpolar gyre. These exchanges lead to variations in the influence exerted by the subarctic or Lusitanian biomes on the intermediate faunistic zone in the north-eastern Atlantic. This strong and persistent bottom-up bio-physical link is demonstrated using a numerical ocean general circulation model and data on four trophically connected levels in the food chain – phytoplankton, zooplankton, blue whiting, and pilot whales. The plankton data give a unique basin-scale depiction of these changes, and a long pilot whale record from the Faroe Islands offers an exceptional temporal perspective over three centuries. Recent advances in simulating the dynamics of the subpolar gyre suggests a potential for predicting the distribution of the main faunistic zones in the north-eastern Atlantic a few years into the future, which might facilitate a more rational management of the commercially important fisheries in this region.

Macrophysiology of Calanus finmarchicus in the North Atlantic Ocean

Copepods represent the major part of the dry weight of the mesozooplankton in pelagic ecosystems and therefore have a central role in the secondary production of the North Atlantic Ocean. The calanoid copepod species Calanus finmarchicus is the main large copepod in subarctic waters of the North Atlantic, dominating the dry weight of the mesozooplankton in regions such as the northern North Sea and the Norwegian Sea. The objective of this work was to investigate the relationships between both the fundamental and realised niches of C. finmarchicus in order to better understand the future influence of global climate change on the abundance, the spatial distribution and the phenology of this key-structural species. Based on standardised Principal Component Analyses (PCAs), a macroecological approach was applied to determine factors affecting the spatial distribution of C. finmarchicus and to characterise its realised niche. Second, an ecophysiological model was used to calculate the Potential Egg Production Rate (PEPR) of C. finmarchicus and the centre of its fundamental niche. Relationships between the two niches were then investigated by correlation analysis. We found a close relationship between the fundamental and realised niches of C. finmarchicus at spatial, monthly and decadal scales. While the species is at the centre of its niche in the subarctic gyre, our joint macroecological and macrophysiological analyses show that it is at the edge of its niche in the North Sea, making the species in this region more vulnerable to temperature changes.