A new efficient 3D non-hydrostatic free-surface flow model for simulating water wave motions

A new three-dimensional, non-hydrostatic free surface flow model is presented. For simulating water wave motions over uneven bottoms, the model employs an explicit project method on a Cartesian the staggered gird system to solve the complete three-dimensional Navier–Stokes equations. A bi-conjugated gradient method with a pre-conditioning procedure is used to solve the resulting matrix system. The model is capable of resolving non-hydrostatic pressure by incorporating the integral method of the top-layer pressure treatment, and predicting wave propagation and interaction over irregular bottom by including a partial bottom-cell treatment. Four examples of surface wave propagation are used to demonstrate the capability of the model. Using a small of vertical layers (e.g. 2–3 layers), it is shown that the model could effectively and accurately resolve wave shoaling, non-linearity, dispersion, fission, refraction, and diffraction phenomena.     

Bulk particulate organic carbon in the East China Sea: Tidal influence and bottom transport

Five oceanographic surveys were carried out in the East China Sea (ECS) and Yellow Sea from 1999 to 2003. In all, seven different sections were surveyed, but one section (the PN section) was observed on every cruise. Two time-series stations were also surveyed, one located at the Changjiang River mouth, the other over the continental shelf in the PN section.We identified biogeochemical characteristics for waters close to the Changjiang Estuary and in the Kuroshio waters (KW), respectively. Resuspension is a strong feature near bottom over the ECS continental shelf, with suspended matter values 13 times higher than that for the surface. A model of particulate organic carbon (POC) dynamics based on a rectangle equation reveals that POC concentration close to the Changjiang Estuary varies with a semidiurnal period of ∼13 h, coinciding with the tidal period. The upper limit for POC residence times in the seasons we covered over the shelf are estimated to be on the order of weeks and generally increase seaward from near the Changjiang Estuary to the KW. Short POC residence times suggest that POC in the ECS is rapidly exported from euphotic waters.A nepheloid layer, observed as elevated suspended matter in near the bottom of the water column, is important in particle transfer over the shelf, especially in winter when the residual current flows mainly eastward. Cross-shelf transport of POC via the nepheloid layer is estimated to be 0.22 × 10¹² g yr⁻¹. Comparison with other work indicates that POC transport is ∼2% of the Changjiang POC input.     

Numerical simulation of Pacific water intrusions into Otsuchi Bay,northeast of Japan,with a nested-grid OGCM

A numerical simulation of Otsuchi Bay located on the northeast coast of the Honshu, the largest island of Japan, is conducted, using an ocean general circulation model (OGCM) with a nested-grid system in order to illustrate seasonal variability of the circulation in the bay. Through a year, an anticlockwise circulation is dominant in the bay, as observational studies have implied, although it is modified in the bay-mouth-half of the bay in winter. In addition, there is an intense outflow at the surface layer during spring to autumn, influenced by river water discharge. Intrusion of the Pacific water into the bay is influened by mean circulations, but it is also influenced by baroclinic tides from spring to autumn. Pacific water intrusions affected by baroclinic tides may have an impact on the environment in Otsuchi Bay.     

Underwater observations of the giant spoon worm <Emphasis Type="Italic">Ikeda taenioides</Emphasis> (Annelida: Echiura: Ikedidae) in a subtidal soft-bottom environment in northeastern Japan,which survived tsunamis of the 2011 off the Pacific Coast of Tohoku Earthquake

Tsunamis associated with the 2011 off the Pacific Coast of Tohoku Earthquake seriously disrupted the shallow marine ecosystem along a 2000 km stretch of the Pacific coast of Japan. The effects of the 2011 tsunamis on the soft-bottom benthic community have been relatively well studied in the intertidal zone, whereas tsunami effects on the subtidal benthos remain poorly understood. Here, we investigated populations of the world’s largest spoon worm Ikeda taenioides (Annelida: Echiura: Ikedidae) in subtidal zone of Funakoshi Bay, Tohoku District, northeastern Japan. Subtidal scuba-diving surveys at two sites in the bay showed extremely long proboscises frequently extending from small holes in the sandy seafloor shortly before and soon after the tsunami disturbances. Based on morphological and molecular identification, the proboscises were revealed to be parts of I. taenioides. On 30 November 2011, 265 days after the tsunami event, many large-sized individuals with >1 m long proboscises were observed; these individuals were probably not derived from post-tsunami larval recruitment but more likely survived the tsunami disturbances. This is surprising because other sympatric megabenthos (e.g. spatangoid echinoids and venerid bivalves) and seagrass beds were almost completely destroyed (although they later recovered) by the tsunamis in this bay. The burrows of I. taenioides are known to be very deep (70–90 cm), which may have sheltered them from the impacts of the tsunamis. Our observations suggest that the effects of the 2011 tsunamis on benthos in soft sediments may differ depending on their burrowing depth.     

Optimization of dissolved urea measurements in coastal waters with the combination of a single reagent and high temperature

Urea is an unstable and intermediate organic nitrogenous compound present in coastal environments and is derived from the excretion of some aquatic organisms and wastewater discharges. Urea plays an important role in the nitrogen cycle, where it is utilized by algae, including diatoms. However, there are very limited relevant data on the production, consumption, and degradation of urea because of the lack of appropriate measurement techniques. The conventional method is based on the formation of a colored product when urea reacts with diacetyl monoxime in a sulfuric acid solution. We examined the optimal conditions for the formation of the colored product; specifically, we evaluated different temperatures (22–80 °C), reaction times, mixing ratios of color reagents, and sample storage times. Application of the single mixed color-developing reagent (COLDER) at 70 °C resulted in the optimal formation of the colored product within a short reaction time of 60 min. This method was then used to measure dissolved urea in different coastal environments. The concentrations detected were as follows: 0.65–0.72, 0.49–0.58, and 1.09–2.28 µM urea-N at coral reef, seagrass, and mangrove sites, respectively. Our results showed high precision (SD = 0.02, CV = 1.2%), a low detection limit (0.03 µM urea-N), and a high recovery rate (94–99%). In summary, this high-temperature procedure for urea measurements should be valuable for obtaining high-precision data that can further the understanding of urea dynamics and its role in coastal ecosystems.     

Long-term change and variation of salinity in the western North Pacific subtropical gyre revealed by 50-year long observations along 137°E

The 137°E repeat hydrographic section for 50 winters during 1967–2016 has been analyzed to examine interannual to interdecadal variations and long-term changes of salinity and temperature in the surface and intermediate layers of the western North Pacific, with a particular focus on freshening in the subtropical gyre. Rapid freshening on both isobars and isopycnals began in the mid-1990s and persisted for the last 20 years in the upper main thermocline/halocline in the western subtropical gyre. In addition, significant decadal variability of salinity existed in the subtropical mode water (STMW), as previously reported for the shallower layers. An analysis of the 144°E repeat hydrographic section during 1984–2013 supplemented by Argo profiling float data in 2014 and 2015 revealed that the freshening trend and decadal variability observed at 137°E originated in the winter mixed layer in the Kuroshio Extension (KE) region and was transmitted southwestward to 137°E 1–2 years later in association with the subduction and advection of STMW. The mechanism of these changes and variations in the source region was further investigated. In addition to the surface freshwater flux in the KE region pointed out by previous studies, the decadal KE variability in association with the Pacific Decadal Oscillation likely contributes to the decadal salinity variability through water exchange between the subtropics and the subarctic across the KE. Interdecadal change in both the surface freshwater flux and the KE state, however, failed to explain the rapid freshening for the last 20 years.     

A note on ocean surface drift with application to surface velocities measured with HF Radar

The ocean drift current consists of a (local) pure drift current generated by the interaction of wind and waves at the sea surface, to which the surface geostrophic current is added vectorially. We present (a) a similarity solution for the wave boundary layer (which has been validated through the prediction of the 10-m drag law), from which the component of pure drift current along the direction of the wind (and hence the speed factor) can be evaluated from the 10-m wind speed and the peak wave period, and (b) a similarity solution for the Ekman layers of the two fluids, which shows that under steady-state neutral conditions the pure drift current lies along the direction of the geostrophic wind, and has a magnitude 0.034 that of the geostrophic wind speed. The co-existence of these two similarity solutions indicates that the frictional properties of the coupled air-sea system are easily evaluated functions of the 10-m wind speed and the peak wave period, and also leads to a simple expression for the angle of deflection of the pure drift current to the 10 m wind. The analysis provides a dynamical model for global ocean drift on monthly and annual time scales for which the steady-state neutral model is a good approximation. In particular, the theoretical results appear to be able to successfully predict the mean surface drift measured by HF Radar, which at present is the best technique for studying the near surface velocity profile.     

Exponential leap-forward gradient scheme for determining the isothermal layer depth from profile data

Two distinct layers usually exist in the upper ocean. The first has a near-zero vertical gradient in temperature (or density) from the surface and is called the isothermal layer (or mixed layer). Beneath that is a layer with a strong vertical gradient in temperature (or density), called the thermocline (or pycnocline). The isothermal layer depth (ILD) or mixed layer depth (MLD) for the same profile varies depending on the method used to determine it. Also, whether they are subjective or objective, existing methods of determining the ILD do not estimate the thermocline (pycnocline) gradient. Here, we propose a new exponential leap-forward gradient (ELG) method of determining the ILD that retains the strengths of subjective (simplicity) and objective (gradient change) methods and avoids their weaknesses (subjective methods are threshold-sensitive and objective methods are computationally intensive). This new method involves two steps: (1) the estimation of the thermocline gradient G _th for an individual temperature profile, and (2) the computation of the vertical gradient by averaging over gradients using exponential leap-forward steps. Such averaging can filter out noise in the profile data. Five existing methods of determining the ILD (difference, gradient, maximum curvature, maximum angle, and optimal linear fitting methods) as well as the proposed ELG method were verified using global expendable bathythermograph (XBT) temperature and conductivity–temperature–depth (CTD) datasets. Among all the methods considered, the ELG method yielded the highest skill score and the lowest Shannon information entropy (i.e., the lowest uncertainty).     

Morphology of Shatsky Rise oceanic plateau from high resolution bathymetry

Newly collected, high resolution multi-beam sonar data are combined with previous bathymetry data to produce an improved bathymetric map of Shatsky Rise oceanic plateau. Bathymetry data show that two massifs within Shatsky Rise are immense central volcanoes with gentle flank slopes declining from a central summit. Tamu Massif is a slightly elongated, dome-like volcanic edifice; Ori Massif is square shaped and smaller in area. Several down-to-basin normal faults are observed on the western flank of the massifs but they do not parallel the magnetic lineations, indicating that these faults are probably not related to spreading ridge faulting. Moreover, the faults are observed only on one side of the massifs, which is contrary to expectations from a mechanism of differential subsidence around the massif center. Multi-beam data show many small secondary cones with different shapes and sizes that are widely-distributed on Shatsky Rise massifs, which imply small late-stage magma sources scattered across the surface of the volcanoes in the form of lava flows or explosive volcanism. Erosional channels occur on the flanks of Shatsky Rise volcanoes due to mass wasting and display evidence of down-slope sediment movement. These channels are likely formed by sediments spalling off the edges of summit sediment cap.     

The potential hydrothermal systems unexplored in the Southwest Indian Ocean

Deep-sea hydrothermal vents possess complex ecosystems and abundant metallic mineral deposits valuable to human being. On-axial vents along tectonic plate boundaries have achieved prominent results and obtained huge resources, while nearly 90% of the global mid-ocean ridge and the majority of the off-axial vents buried by thick oceanic sediments within plates remain as relatively undiscovered domains. Based on previous detailed investigations, hydrothermal vents have been mapped along five sections along the Southwest Indian Ridge (SWIR) with different bathymetry, spreading rates, and gravity features, two at the western end (10°–16°E Section B and 16°–25°E Section C) and three at the eastern end (49°–52°E Section D, 52°–61°E Section E and 61°–70°E Section F). Hydrothermal vents along the Sections B, C, E and F with thin oceanic crust are hosted by ultramafic rocks under tectonic-controlled magmatic-starved settings, and hydrothermal vents along the Section D are associated with exceed magmatism. Limited coverage of investigations is provided along the 35°–47°E SWIR (between Marion and Indomed fracture zones) and a lot of research has been done around the Bouvet Island, while no hydrothermal vents has been reported. Analyzing bathymetry, gravity and geochemical data, magmatism settings are favourable for the occurrence of hydrothermal systems along these two sections. An off-axial hydrothermal system in the southern flank of the SWIR that exhibits ultra-thin oceanic crust associated with an oceanic continental transition is postulated to exist along the 100-Ma slow-spreading isochron in the Enderby Basin. A discrete, denser enriched or less depleted mantle beneath the Antarctic Plate is an alternative explanation for the large scale thin oceanic crust concentrated on the southern flank of the SWIR.