A review of the ocean science and technology partnership between US and Korea

In 2001, the National Oceanic and Atmospheric Administration (NOAA) and the Ministry of Maritime Affairs and Fisheries Integrated Coastal and Ocean Resources Scientific and Technical Arrangement was signed to pursue scientific and technical cooperation in integrated coastal and ocean resources management in the mutual interest of the participants. Over the last 6 years, the cooperation has been very beneficial, particularly for advancing marine science and management programs in Korea, such as Deep Sea Aquaculture, Sea Grant College Program, Marine Protected Areas, and Ocean Observations. The benefits to NOAA from the cooperation include Deep-ocean Assessment and Reporting of Tsunamis (DART) buoys, ocean and coastal observations, GEOSS program, and Global Program of Action. Also the cooperation has spillover effects on other ocean science and technological arrangement between US and Korea. This article reviews purpose, obstacles, and achievement of the cooperation and suggests recommendations for the future steps.     

44 years hindcast of sea level in the Atlantic Coast of Europe

The sea level of Northeast Atlantic Ocean is calculated for the period between 1958 and 2001 using a state-of-the-art barotropic model with a grid size of 10′ × 15′ (long × lat). The model includes astronomic effects, considering seven components of the tide, and the meteorological effects of wind and atmospheric pressure, allowing obtaining the astronomic tide, the atmospheric residuals and the non-linear addition of both components of sea level.     

A new method for the quantification of different redox-species of molybdenum (V and VI) in seawater

A new method for the direct determination of reduced and oxidized Mo species (Mo (V) and Mo (VI)) in seawater was developed and used for the first time. The method includes the complexation of Mo (V) with tartrate, solid phase extraction of the Mo (V)–tartrate complex by a XAD 7HP resin, followed by elution with acidic acetone. In this study, the eluted Mo (V) was quantified by graphite furnace atomic absorption spectrometry. The detection limit of this protocol was on the order of 0.2 nM. The analytical precision was 10% of ~ 10 nM. This method was successfully applied to the determination of Mo (V) and Mo (VI) in surface and bottom waters at the head of Peconic River Estuary. Total Mo (Mo (V) + Mo (VI)) ranged from 100–120 nM in most bottom saline waters, and 2.5–15 nM for surface fresher waters. Concentrations of Mo (V) in these environments ranged from 0 nM to ~ 15 nM, accounting for 0%–15% of the total dissolved Mo pool. The time series experiments showed that the Mo speciation changed within 1 h after the water collection, and therefore it is strongly suggested that speciation analysis be carried out within the first 15 min. However, since these are the first Mo speciation data in concentration ranges typical of normal marine and coastal waters, additional research may be required to optimize the methodology and further explore Mo cycling mechanisms.     

Tortonian-Pleistocenic oceanic features in the Southern Tyrrhenian Sea: magnetic inverse model of the Selli-Vavilov region

We show the magnetic model of the Selli-Vavilov region. The Selli Line is known as the northwestern edge of the southern Tyrrhenian Basin. The tectonic evolution of the Tyrrhenian Basin is dominated by a Tortonian-Quaternary extension through the eastward movement of the Apennine subduction system. This migration has generated a diffuse stretching of the continental crust with the emplacement of new oceanic material. This latter occurred in several localized zones where the eastward retreating of the Ionian subduction system produced a strong depletion of the crust with formation of basins and correlated spreading. Nowadays the presence of oceanic crust is confirmed through direct drilling investigation but a complete mapping of the oceanic crustal distribution is still lacking. The Selli-Vavilov region shows a differentiated crustal setting where seamount structures, the oceanic basement portions and continental crust blocks are superimposed. To this aim, a 2D inversion of the magnetic data of this region was conducted to define buried structures. The magnetic susceptibility pattern was computed by solving the least squares problem of the misfit between the predicted and real data for separated wavebands. This method produced two 2D models of the high and low frequency fields of the Selli-Vavilov region. The two apparent susceptibility maps provide different information for distinct ranges of depth. The results of the inversions were also combined with seismic data of the Selli region highlighting the position of the highly magnetized buried bodies. The results confirm a role for the Selli Line as a deep crustal boundary dividing the Sardinian passive domain from the easternmost active region where different oceanic structures are located. The Selli Line has worked as a detachment fault system which has moved eastward. Finally, the Selli-Vavilov region may be interpreted as a tectonic result due to a passive asymmetrical rift occurred between the Tortonian and Pliocene.     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

Mesozooplankton biomass and grazing responses to Cyclone Opal, a subtropical mesoscale eddy

As part of E-Flux III cruise studies in March 2005, plankton net collections were made to assess the effects of a cyclonic cold-core eddy (Cyclone Opal) on the biomass and grazing of mesozooplankton. Mesozooplankton biomass in the central region of Cyclone Opal, an area of uplifted nutricline and a subsurface diatom bloom, averaged 0.80±0.24 and 1.51±0.59 g DW m⁻², for day and night tows, respectively. These biomass estimates were about 80% higher than control (OUT) stations, with increases more or less proportionately distributed among size classes from 0.2 to >5 mm. Though elevated relative to surrounding waters south of the Hawaiian Islands (Hawai’i lee), total biomass and size distribution in Cyclone Opal were almost exactly the same as contemporary measurements made at Stn. ALOHA, 100 km north of the islands, by the HOT (Hawaii Ocean Time-series) Program. Mesozooplankton biomass and community composition at the OUT stations were also similar to ALOHA values from 1994 to 1996, preceding a recent decadal increase. These comparisons may therefore provide insight into production characteristics or biomass gradients associated with decadal changes at Stn. ALOHA. Gut fluorescence estimates were higher in Opal than in ambient waters, translating to grazing impacts of 0.11±0.02 d⁻¹ (IN) versus 0.03±0.01 d⁻¹ (OUT). Over the depth-integrated euphotic zone, mesozooplankton accounted for 30% of the combined grazing losses of phytoplankton to micro- and meso-herbivores in Opal, as compared to 13% at control stations. Estimates of active export flux by migrating zooplankton averaged 0.81 mmol C m⁻² d⁻¹ in Cyclone Opal and 0.37 mmol C m⁻² d⁻¹ at OUT stations, 53% and 24%, respectively, of the carbon export measured by passive sediment traps. Migrants also exported 0.18 mmol N m⁻² d⁻¹ (117% of trap N flux) in Cyclone Opal compared to 0.08 mmol N m⁻² d⁻¹ (51% of trap flux) at control stations. Overall, the food-web importance of mesozooplankton increased in Cyclone Opal both in absolute and relative terms. Diel migrants provided evidence for enhanced export flux in the eddy that was missed by sediment trap and ²³⁴Th techniques, and migrant-mediated flux was the major export term in the observed bloom-perturbation response and N mass balance of the eddy.     

Population structure of the sicklefin lemon shark Negaprion acutidens within the Curieuse Marine National Park,Seychelles

This study produced the first data on population parameters for sicklefin lemon shark Negaprion acutidens at the granitic Seychelles Islands, adding to a limited catalogue of species information. Juveniles at a site in the Curieuse Marine National Park were tagged over two periods, in 2014–2015 and 2015–2016, with measurements of length (precaudal, fork and total) and weight, sex, and state of the umbilical opening recorded. Recapture data were used to calculate growth rates and population size. Mean annual growth rate was 5.40 cm y^–1. A mean size at birth of 62.5 cm TL is in line with previous assessments of the species. Population-size estimates were similar for the two sampling periods (comprising 311 and 255 individuals) and higher than anticipated. Movements to outside the lagoon were recorded for early ages, and this, combined with weight loss and a reduction in condition factor as displayed in several individuals, suggests that the study area alone could not support the large number of sharks pupped yearly.     

The WOCE-era 3-D Pacific Ocean circulation and heat budget

To address questions concerning the intensity and spatial structure of the three-dimensional circulation within the Pacific Ocean and the associated advective and diffusive property flux divergences, data from approximately 3000 high-quality hydrographic stations collected on 40 zonal and meridional cruises have been merged into a physically consistent model. The majority of the stations were occupied as part of the World Ocean Circulation Experiment (WOCE), which took place in the 1990s. These data are supplemented by a few pre-WOCE surveys of similar quality, and time-averaged direct-velocity and historical hydrographic measurements about the equator.An inverse box model formalism is employed to estimate the absolute along-isopycnal velocity field, the magnitude and spatial distribution of the associated diapycnal flow and the corresponding diapycnal advective and diffusive property flux divergences. The resulting large-scale WOCE Pacific circulation can be described as two shallow overturning cells at mid- to low latitudes, one in each hemisphere, and a single deep cell which brings abyssal waters from the Southern Ocean into the Pacific where they upwell across isopycnals and are returned south as deep waters. Upwelling is seen to occur throughout most of the basin with generally larger dianeutral transport and greater mixing occurring at depth. The derived pattern of ocean heat transport divergence is compared to published results based on air–sea flux estimates. The synthesis suggests a strongly east/west oriented pattern of air–sea heat flux with heat loss to the atmosphere throughout most of the western basins, and a gain of heat throughout the tropics extending poleward through the eastern basins. The calculated meridional heat transport agrees well with previous hydrographic estimates. Consistent with many of the climatologies at a variety of latitudes as well, our meridional heat transport estimates tend toward lower values in both hemispheres.     

Infaunal burrow ventilation and pore-water transport in muddy sediments

The ventilation of burrows by tube-dwelling benthos is understood to be important in determining rates of exchange of solutes between the sediment and overlying water. However, few models have attempted to link the burrow ventilation behavior of tube-dwelling organisms with their geochemical consequences. The classic cylinder model of bioirrigation in muddy sediments (Aller, R.C., 1980. Quantifying solute distributions in the bioturbated zone of marine sediments by defining an average microenvironment. Geochimica et Cosmochimica Acta 44, 1955–1965) links pore-water processes and burrow sizes and distributions in the sediment by assuming that burrows are fully flushed. The equivalence between the cylinder model and the more commonly used one-dimensional non-local exchange model depends upon this assumption. However, this assumption has seldom been tested in the field. We have extended the cylinder model of bioirrigation to include burrow ventilation activities of organisms. Burrow ventilation is modeled as a simple non-local exchange of burrow water with overlying water. Model simulations indicate that burrow ventilation has a large effect on vertical profiles and fluxes of solute tracers. We collected data on burrow geometry in the field by CT-scanning freshly collected sediment cores. At the same study sites, we measured activity profiles of ²²²Rn, a naturally occurring radionuclide tracer of pore-water transport. With model geometry independently constrained, we tested the model by comparing our model-predicted profiles with measured profiles. Our results demonstrate that burrows in the field are not fully flushed. Our estimated burrow ventilation rates compare favorably with previous laboratory measurements. The inclusion of realistic burrow ventilation in this pore-water transport model strongly affects modeled solute profiles and fluxes. We demonstrate how model parameters can be determined from field samples and present a model that more realistically simulates pore-water transport processes in muddy sediments.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.