Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations

Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching , or even higher.In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES.This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3×10³Re3×10⁴, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re=10⁵ and 10⁶, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67% when Re is increased from Re=10³ to Re=10⁶, within the computation period, with the most significant increase taking place from Re=3×10³ to Re=10⁵.     

Contributions to the dynamics of the Antarctic Circumpolar Current (A. C. C.) II. Integral relationship

Non-dimensional equations of motion are derived for the A.C.C. of the barotropic mode, including the bottom friction and the horizontal eddy viscosity. Integration of the vorticity equation along a streamline leads to the zeroth order stream function which is dependent only on depth divided by Coriolis parameter. Integration of the momentum equation along a streamline yields the relation between the momentum input by wind stress and its dissipation by the bottom friction and by the horizontal eddy viscosity. This relation determines the magnitude of the stream function. It explains differences in the total transport of the A.C.C. obtained byBryan andCox (1972), though it gives only one third of the total transport obtained byKamenkovich (1972) with his vertical eddy viscosity of 10²cm² s^?1. With 1 cm² s^?1 of this viscosity,Bryan andCox obtained the transport of about 650 or less than 32×10⁶m³s^?1 for constant or variable depth models, respectively. The higher transport is mainly due to broadening of the width of the A.C.C., whereas the lower value is due to its narrowing and meandering which in turn make the horizontal eddy viscosity more effective (by exercising friction on both sides of the A.C.C.) and the wind stress input smaller than the almost zonal streamlines for constant depth. In the Appendix dynamics of the bottom boundary layer is treated to give rational estimates of the bottom stress in terms of the geostrophic flow and is compared to the recent observations of the benthic boundary current in the Straits of Florida and off San Diego.     

Offshore Absolute Calibration of Space-Borne Radar Altimeters

Absolute calibration of sea level measurements collected from space-borne radar altimeters is usually performed with respect to collocated sea level in situ records from tide gauges or GPS buoys (Ménard et al. 1994 Ménard, Y., Jeansou, E. and Vincent, P. 1994. Calibration of the TOPEX-Poseidon altimeters at Lampedusa: Additional results at Harves. J. Geophys Res., 99(C12): 24487–24504. http://dx.doi.org/10.1029%2F94JC01300 [Google Scholar]; Haines et al. 1996 Haines, B. J., Christensen, E. J., Norman, R. A., Parke, M. E., Born, G. H. and Gill, S. K. 1996. Altimeter calibration and geophysical monitoring from collocated measurements at the Harvest oil platform. EOS Trans. Suppl., 77(22): W16 [Google Scholar]; Bonnefond et al. 2003; Haines et al. 2003 Haines, B. J., Dong, D., Born, G. H. and Gill, S. K. 2003. The Harvest experiment: Monitoring Jason-1 and TOPEX/Poseidon from a California offshore platform. Mar. Geod., 26: 239–259. [Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]; Schum et al. 2003 Schum, C. K., Yi, Y., Cheng, K., Kuo, C., Braun, A., Calmant, S. and Chambers, D. 2003. Calibration of Jason-1 Altimeter over Lake Erie. Mar. Geod., 26: 335–354.  [Google Scholar]; Watson et al. 2003 Watson, C., Coleman, R., White, N., Church, J. and Govind, R. 2003. Absolute calibration of TOPEX/ Poseidon and Jason-1 using GPS buoys in Bass Strait, Australia. Mar. Geod., 26: 285–304. [Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]; Watson et al. 2004 Watson, C., White, N., Coleman, R., Church, J., Morgan, P. and Govind, R. 2004. TOPEX/Poseidon and Jason-1: Absolute calibration in Bass Strait, Australia. Mar. Geod., 27: 107–131. http://dx.doi.org/10.1080%2F01490410490465373[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). Such a method allows regular and long-term control of altimetric systems with independent records. However, this approach is based on a single, geographically dependent point. In order to obtain more significant and accurate bias and drift estimates, there is a strong interest in multiplying the number of calibration opportunities. This article describes a method, called the “offshore method” that was developed to extend the single-point approach to a wider regional scale. The principle is to compare altimeter and tide gauge sea level data not only at the point of closest approach of an overflying pass, but also at distant points along adjacent satellite passes. However, connecting sea level satellite measurements with more distant in situ data requires a more accurate determination of the geoid and mean ocean dynamic topography slopes, and also of the ocean dynamical changes. In this demonstration experiment, 10 years of averaged TOPEX/Poseidon mean sea level profiles are used to precisely determine the geoid and the mean ocean circulation slope. The Mog2d barotropic ocean model (Carerre et Lyard 2003 Carrère, L. and Lyard, F. 2003. Modelling the barotropic response of the global ocean to atmospheric wind and pressure forcing-comparisons with observations. GRL, 30(6): 1275 [Google Scholar]) is used to improve our estimate of the ocean dynamics term. The method is first validated with Jason-1 data, off Corsica, where the dedicated calibration site of Senetosa provides independent reference data. The method is then applied to TOPEX/Poseidon on its new orbit and to Geosat Follow On. The results demonstrate that it is feasible to make altimeter calibrations a few tens to hundreds of kilometers away from a dedicated site, as long as accurate mean sea level altimeter profiles can be used to ensure the connection with reference tide gauges.     

Artificial neural network modeling of cross-shore profile on sand beaches: The coast of the province of Valencia (Spain)

The paper describes the training, validation, testing, and application of models of artificial neural networks (ANN) for computing the cross-shore beach profile of the sand beaches of the province of Valencia (Spain). Sixty ANN models were generated by modifying both the input variables as the number of neurons in the hidden layer. The input variables consist of wave data and sedimentological data. To select and evaluate the performance of the optimal model, the following parameters were used: R², absolute error, mean absolute percentage error, and percentage relative error. Finally, the results are compared with the numerical model proposed by Aragonés et al. (2016b Aragonés, L., Y. Villacampa, F. J. Navarro-González, and I. López. 2016b. Numerical modelling of the equilibrium profile in Valencia (Spain). Ocean Engineering, 123:164–73. doi:10.1016/j.oceaneng.2016.07.036[Crossref], [Web of Science ®] , [Google Scholar]) for the equilibrium profile in the study area. The results show a mean absolute error of 0.21?m compared to 0.33?m Aragones’ model, significantly improving the results of the numerical model in the bar area around de Valencia Port. In addition, when comparing the results with other methods currently used (Dean’s or Vellinga formulation), the errors of these compared to ANN are of the order of 167 and 1538% higher, respectively.     

Validation of Practice Oriented Models and Influence of Soil Stiffness on Lateral Pile Response Due to Kinematic Loading

Earthquakes usually cause both inertial and kinematic loading of pile foundations. Inertial loading and lateral pile response can be predicted to a certain extent by lateral pile load tests for important engineering projects. However, there is no accepted in situ assessment method for pile behavior due to kinematic loading. Therefore, it is important for geotechnical engineering practice to validate the practice oriented models and examine the influence of soil strength parameters on the performance of kinematically loaded piles due to lateral spreading. For this purpose, a well-documented case of full-scale lateral spreading test conducted at Port of Tokachi in Japan by Ashford et al. (2006 Ashford , S. A. , T. Juirnarongrit , T. Sugano , and M. Hamada . 2006 . Soil-pile response to blast-induced lateral spreading. I: Field test . Journal of Geotechnical and Geoenvironmental Engineering 132 ( 2 ): 152 – 162 . [Google Scholar]) is selected as a primary reference in terms of test data and site conditions. Four p-y analyses are performed using a special purpose software based on finite difference technique. The first two investigate the employment of different p-y curves suggested in the literature. The third one investigates the influence of variations in internal friction angle and undrained shear strength. In the fourth analysis, liquefied shear strength and residual shear strength approaches are utilized for the liquefiable layers and their effects are investigated. In all analyses, mobilized lateral pressures in non-liquefied layers are also scrutinized. The results of the analyses are discussed and compared with the measured values in the field. Finally, conclusions are given based on the analyses and discussions.     

Deltaic soil behavior classification using CPTu records—Proposed approach and applied to fifty-four case histories

For the past four decades, the CPT has played a key role in onshore and offshore soil investigations. One of the main applications of cone penetration test (CPT) is the soil behavioral classification. Most of the developed methods for soil identification using CPT and CPTu (piezocone) data are well categorized for common soils, such as clays, silts, and sands. Soils with low resistance or more compressibility generally involve problems in geotechnical engineering practice and construction projects. Consequently, these unusual deposits require further evaluation and more detailed data. Five major groups of problematic soils including: liquefiable, sensitive, peaty, collapsible, and expansive soils have been considered in this study. One hundred and forty CPT and CPTu test records were collected from fifteen countries. Sixty-one of the records are related to difficult soils. A brief comparison is performed for currently used soil behavioral classification charts, such as by Campanella et al. (1985 Campanella, R. G., P. K. Robertson, D. Gillespie, and J. Greig. 1985. Recent developments in in-situ testing of soils. Proceedings of 11th International Conference on Soil Mechanics and Foundation Engineering, ICSMFE, San Francisco, Vol. 2, 849–54. [Google Scholar]), Robertson (1990 Robertson, P. K. 1990. Soil classification using the cone penetration test. Canadian Geotechnical Journal 27 (1):151–58. doi:10.1139/t90–014[Crossref], [Web of Science ®] , [Google Scholar]), Jefferies and Davies (1991 Jefferies, M. G., and M. P. Davies. 1991. Soil classification using the cone penetration test: Discussion. Canadian Geotechnical Journal 28 (1):173–76. doi:10.1139/t91–023[Crossref], [Web of Science ®] , [Google Scholar]) and Eslami and Fellenius (1997 Eslami, A., and B. H. Fellenius. 1997. Pile capacity by direct CPT and CPTu methods applied to 102 case histories. Canadian Geotechnical Journal 34 (6):886–904. doi:10.1139/cgj-34–6-886[Crossref], [Web of Science ®] , [Google Scholar]). Analysis based on CPT data indicates that a few commonly used charts recognize relatively well problematic deposits. However, further studies are needed to increase the accuracy and capability of methods. Existing charts have some problems due to the limitations of the nature of rectangular charts based on two axes. A new format of classification chart, i.e., triangular form containing cone tip resistance (q_c), sleeve friction (f_s), and pore pressure (u₂) is proposed for soil identification which can be realized in practice. The proposed chart with more accuracy and less scattering of data than the previous charts is able to identify soil types particularly for deltaic soils.     

Ventilation revealed by the observation of dissolved oxygen concentration south of the Kuroshio Extension during 2012–2013

From the moored buoy observation at \(33.9^{\circ }\)N, \(144.9^{\circ }\)E south of the Kuroshio Extension (KE), we obtained Eulerian time series of dissolved oxygen concentration (DO) at 200, 400, and 600 m depths from June 2012 to March 2013. We observed ventilation by meso- and submesoscale processes that transport water southward across the KE jet. First, the cyclonic mesoscale eddy in June 2012 substantially depressed DO at depths of 400 and 600 m but maintained DO at 200 m, suggesting near-surface lateral transport of high-DO water derived from the north of the KE. Second, subduction of high-DO (>230 \({\upmu }\)mol kg\(^{-1}\)) water to a depth of 600 m was observed from early February to March 2013, associated with a mesoscale/submesoscale meandering of the KE jet. In mid-March 2013, shipboard hydrographic data were collected where the water mass at the mooring site would be advected by the eastward current on the southern flank of the KE. Based on these data, the subduction event was identified as an intrusion of an anomalously thick water mass from approximately 400–900 dbar. Ventilation of the subtropical mode water at a depth of 200 m around a subsurface DO maximum layer was detected as a rapid rise in DO in January 2013. This occurred after a significant seasonal decrease in DO at a rate of \(-0.130 \pm 0.007\) \({\upmu }\)mol kg\(^{-1}\) day\(^{-1}\) from July to December 2012.     

Ibiza Absolute Calibration Experiment: Survey and Preliminary Results

Within the framework of a project comprising part of the Spanish Space Program related to the JASON-1 CNES (Centre National d'Etudes Spatiales)/NASA (National Aeronautics and Space Administration) mission, a campaign was conducted from June 9–17, 2003, on the Absolute Calibration Site of the island of Ibiza. The objective was to determine the local marine geoid slope under the ascending (187) and descending (248) Jason-1 ground tracks, in order to allow a better extrapolation of the open-ocean altimetric data with on-shore tide gauge locations, and thereby improve the overall precision of the calibration process. For this we have used a catamaran with two GPS antennas onboard, following the Corsica/Senetosa design (Bonnefond et al. 2003a Bonnefond, P., Exertier, P., Laurain, O., Menard, Y., Orsoni, A., Jeansou, E., Haines, B., Kubitschek, D. and Born, G. 2003a. Leveling Sea Surface using a GPS catamaran. Marine Geodesy, 26(3–4): 319–334. [Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). Five GPS reference stations were deployed in order to reduce the distance between the areas covered by the catamaran and the fixed GPS receiver used in the kinematic process. The geodetic activities (e.g., GPS, leveling) have enabled the building of a very accurate (few mm) network in a reference frame compatible with the satellite altimetry missions (ITRF 2000). The GPS kinematic data were processed using two different software programmes, allowing checking of the consistency of the solutions. If the standard deviation of the differences (3.3 cm) is close to the kinematic process precision, they exhibit some large values (up to 14 cm). These large discrepancies have been reduced using a weighting based on the crossover differences. Inasmuch as the distances between the tide gauges and the areas covered by the GPS catamaran were becoming large, we have used the MOG2D ocean model (Carrère and Lyard 2003 Carrère, L. and Lyard, F. 2003. Modelling the barotropic response of the global ocean to atmospheric wind and pressure forcing—comparisons with observations. Geophys. Res. Letters, 30(6) [Google Scholar]) to correct the sea surface from tides. In the farthest areas, the crossover differences show an improvement by a factor of two. Finally, we also present preliminary results on Jason-1 altimeter calibration using the derived marine geoid. From this analysis, the altimeter bias is estimated to be 120 ± 5 mm. The quality of this first result validates the whole GPS-based marine geoid processing, for which the accuracy is estimated to be better than 3 cm rms at crossovers.     

Closure to “anisotropic hydraulic conductivity and critical hydraulic gradient of a crushed sandstone–mudstone particle mixture”

A closure is presented to the discussion of Anisotropic hydraulic conductivity and critical hydraulic gradient of a crushed sandstone–mudstone particle mixture by Salahou, Jiao, and Yousif (2017 Fener, M., and N. Yesiller. 2013. Vertical pore structure profile of a compacted clayey soil. Engineering Geology 166:204–15. doi:10.1016/j.enggeo.2013.08.001.[Crossref], [Web of Science ®] , [Google Scholar]). The relationship between the critical hydraulic gradient and the uniformity coefficient can be limited to two curves. The internal erosion mechanism of soils with different particle grade curves is different. The maximum diameter of fine particles of different particle grade curves can be calculated by a lot of methods. The density of the original paper is the average density. The anisotropic hydraulic conductivity is caused by uneven density distribution. The hydraulic conductivities are tested on the saturated soils. The saturated hydraulic conductivity can be regard as a constant value for a saturated soil.     

Gravitational Sedimentation Behavior of Sensitive Marine Ariake Clays

It has been well documented that natural marine Ariake clays are sensitive clays. In this study, extensive data of marine Ariake clays are obtained to investigate the gravitational compression behavior for sensitive clays. Analysis results indicate that the compression behavior of remolded Ariake clays is not different from that of other remolded/reconstituted soils. But natural Ariake clays do not follow the gravitational compression pattern reported by Skempton (1970) Skempton, A. W. 1970. “The consolidation of clays by gravitational compaction”. In Q. J. Geol. Soc 373–411.  [Google Scholar] for natural sedimentary soils. At a given value of effective overburden pressure, the void ratios of natural Ariake clays are almost independent of liquid limits. Most natural Ariake clays lie above the sedimentation compression line proposed by Burland (1990) Burland, J. B. 1990. On the compressibility and shear strength of natural clays. Gèotechnique, 40: 329–378. [Crossref], [Web of Science ®] , [Google Scholar]. When the liquid limit is larger than 90% and the ratio of natural water content over liquid limit ranges 0.8–1.1, the natural Ariake clays lie around the sedimentation compression line. In addition, the natural Ariake clay with higher value of the ratio of natural water content over liquid limit lies above the natural Ariake clay with lower value of the ratio of natural water content over liquid limit. Salt removal is the most probable cause for such a phenomenon.     

Backscattering from a sandy seabed measured by a calibrated multibeam echosounder in the 190–400 kHz frequency range

A SeaBat T50 calibration that combines measurements in a test tank with data from numerical models is presented. The calibration is assessed with data obtained from a series of tests conducted over a sandy seabed outside the harbor of Santa Barbara, California (April 2016). The tests include different tone-burst durations, sound pressure levels, and receive gains in order to verify that the estimated seabed backscattering strength \((S_b)\) is invariant to sonar settings. Finally, \(S_b\)-estimates obtained in the frequency range from 190 kHz in steps of 10 kHz up to 400 kHz, and for grazing angles from \(20^\circ\) up to 90\(^\circ\) in bins of width \(5^\circ ,\) are presented. The results are compared with results found in the literature.     

Simulation of Seawater Intrusion Range in Coastal Aquifer Using the FEMWATER Model for Disaster Information

The FEMWATER model is a three-dimensional finite element model that can simulate both groundwater flow and mass transport at the same time, and so by using this FEMWATER model seawater intrusion was simulated and the results were compared to observation data. By using the observation data, a conceptual model was made of 3 layers based on the type of soil, and the hydrogeological parameter was selected. There are two types of dispersivities: longitudinal dispersivity and transverse dispersivity. Longitudinal dispersivity was calculated by Xu and Eckstein's method (1995 Xu , M. and Y. Eckstein . 1995 . Use of weighted least-squares method in evaluation of the relationship between dispersivity and field scale . Ground Water 33 ( 6 ): 905 – 908 .[Crossref], [Web of Science ®] , [Google Scholar]) and the transverse dispersivity was put at 3/10 of the longitudinal dispersivity for the simulation. The electrical conductivity suggested by Freeze and Cherry (1979 Freeze , R. A. and J. A. Cherry , J. A. 1979 . Groundwater . Eaglewood Cliffs , NJ : Prentice-Hall . [Google Scholar]) and the correlation with total dissolved solids, TDS, was used to classify the state of the water. Upon analyzing the depths at 15 m, 30 m, 60 m, and 90 m, occurrence of seawater intrusion near the ground surface was limited to the vicinity of the coastline, but as the depths became deeper, seawater intrusion expanded inland. From these results we can come to the conclusion that in order to prevent seawater intrusion there needs to be sufficient preliminary examination.     

The self-thinning exponent in overcrowded stands of the mangrove, Kandelia obovata, on Okinawa Island, Japan

Weller??s allometric model assumes that the allometric relationships of mean area occupied by a tree $ \bar{s} $ , i.e., the reciprocal of population density $ \rho $ , $ \bar{s}\left( { = {1 \mathord{\left/ {\vphantom {1 {\rho = g_{\varphi } \cdot \bar{w}^{\varphi } }}} \right. \kern-0em} {\rho = g_{\varphi } \cdot \bar{w}^{\varphi } }}} \right) $ , mean tree height $ \bar{H}\left( { = g_{\theta } \cdot \bar{w}^{\theta } } \right) $ , and mean aboveground mass density $ \bar{d}\left( { = g_{\delta } \cdot \bar{w}^{\delta } } \right) $ to mean aboveground mass $ \bar{w} $ hold. Using the model, the self-thinning line $ \left( {\bar{w} = K \cdot \rho^{ - \alpha } } \right) $ of overcrowded Kandelia obovata stands in Okinawa, Japan, was studied over 8?years. Mean tree height increased with increasing $ \bar{w} $ . The values of the allometric constant $ \theta $ and the multiplying factor $ g_{\theta } $ are 0.3857 and 2.157?m?kg^???, respectively. The allometric constant $ \delta $ and the multiplying factor $ g_{\delta } $ are ?0.01673 and 2.685?m^?3?kg^1???, respectively. The $ \delta $ value was not significantly different from zero, showing that $ \bar{d} $ remains constant regardless of any increase in $ \bar{w} $ . The average of $ \bar{d} $ , i.e., biomass density $ \left( {{{\bar{w} \cdot \rho } \mathord{\left/ {\vphantom {{\bar{w} \cdot \rho } {\bar{H}}}} \right. \kern-0em} {\bar{H}}}} \right) $ , was 2.641?±?0.022?kg?m^?3, which was considerably higher than 1.3?C1.5?kg?m^?3 of most terrestrial forests. The self-thinning exponent $ \alpha \left( { = {1 \mathord{\left/ {\vphantom {1 {\varphi = }}} \right. \kern-0em} {\varphi = }}{1 \mathord{\left/ {\vphantom {1 {\left\{ {1 - \left( {\theta + \delta } \right)} \right\}}}} \right. \kern-0em} {\left\{ {1 - \left( {\theta + \delta } \right)} \right\}}}} \right) $ and the multiplying factor $ K\left( { = \left( {g_{\theta } \cdot g_{\delta } } \right)^{\alpha } } \right) $ were estimated to be 1.585 and 16.18?kg?m^?2??, respectively. The estimators $ \theta $ and $ \delta $ are dependent on each other. Therefore, the observed value of $ \theta + \delta $ cannot be used for the test of the hypothesis that the expectation of the estimator $ \theta + \delta $ equals 1/3, i.e., $ \alpha = {3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2} $ , or 1/4, i.e., $ \alpha = {4 \mathord{\left/ {\vphantom {4 3}} \right. \kern-0em} 3} $ . The $ \varphi $ value was 0.6310, which is the same as the reciprocal of the self-thinning exponent of 1.585, and was not significantly different from 2/3 (t?=?1.860, df?=?191, p?=?0.06429), i.e., $ \alpha = {3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2} $ . Thus the self-thinning exponent is not significantly different from 3/2 based on the simple geometric model. On the other hand, the self-thinning exponent was significantly different from 3/4 (t?=?6.213, df?=?191, p?=?3.182?×?10^?9), i.e., $ \alpha = {4 \mathord{\left/ {\vphantom {4 3}} \right. \kern-0em} 3} $ . Therefore, the self-thinning exponent is significantly different from 4/3 based on the metabolic model.     

Method for comparing fish growth curves

The method of comparing growth curves that is described in this paper consists of describing growth for the ith group by the equation Li(t) = ai + bi r^t, which is linear in the two parameters ai and bi, and subsequently using the analysis of variance to compare the curves. The advantages of the linear growth equation over the von Bertalanffy equation, l i (t) ⁼ ai + bi r^t i, are easier analysis and a much simpler interpretation of the differences in the parameters in terms of growth rate and maximum size. An example comparing the growth of male and female red gurnard Chelidonichthys kumu from four areas within the Hauraki Gulf is given.     

Characterizing air–sea CO2 exchange dynamics during winter in the coastal water off the Hugli-Matla estuarine system in the northern Bay of Bengal, India

The distribution of the fugacity of CO₂ ( $ f_{{{\text{CO}}_{ 2} }} $ ) and air–sea CO₂ exchange were comprehensively investigated in the outer estuary to offshore shallow water region (lying adjacent to the Sundarban mangrove forest) covering an area of ~2,000 km² in the northern Bay of Bengal during the winter. A total of ten sampling surveys were conducted between 1 December, 2011 and 21 February, 2012. Physico-chemical variables like sea surface temperature (SST), salinity, pH, total alkalinity (TAlk), dissolved inorganic carbon (DIC) and in vivo chlorophyll-a along with atmospheric variables were measured in order to study their role in controlling the CO₂ flux. Surface water $ f_{{{\text{CO}}_{ 2} }} $ ranged between 111 and 459 μatm which correlated significantly with the SST (r = 0.71, p < 0.001, n = 62). Neither DIC nor TAlk showed any linear relationship with varying salinity in the estuarine mixing zone, demonstrating the significant presence of non-carbonate alkalinity. An overall net biological control on the surface $ f_{{{\text{CO}}_{ 2} }} $ distribution was established during the study, although no significant correlation was found between chlorophyll-a and $ f_{{{\text{CO}}_{ 2} }} $ (water). The shallow water region studied was mostly under-saturated with CO₂ and acted as a sink for atmospheric CO₂. The difference between surface water and atmospheric $ f_{{{\text{CO}}_{ 2} }} $ ( $ \Updelta f_{{{\text{CO}}_{ 2} }} $ ) ranged from ?274 to 69 μatm, with an average seaward flux of ?10.5 ± 12.6 μmol m^?2 h^?1. The $ \Updelta f_{{{\text{CO}}_{ 2} }} $ and hence the air–sea CO₂ exchange was primarily regulated by the variation in sea surface $ f_{{{\text{CO}}_{ 2} }} $ , since atmospheric $ f_{{{\text{CO}}_{ 2} }} $ varied over a comparatively narrow range of 361.23–399.05 μatm.     

Effects of stratification and mesoscale eddies on Kuroshio path variation south of Japan

Bimodality of the Kuroshio current path south of Japan is investigated, focusing on the effects of stratification and mesoscale eddies. For this purpose, wind-driven numerical experiments are executed in barotropic and two-layered ocean models. Stratification has two effects on the path selection of the Kuroshio south of Japan. First, it makes an alongshore path stable at intermediate wind stress strength τ₀ by arresting an eddy southeast of Kyushu. This enables an alongshore path to appear in the entire experimental range of τ₀. Second, the upper limit of τ₀ which allows a meandering path decreases from ( in the Sverdrup transport at the Tokara Strait) to () as Δρ/ρ₀ increases from 2.0×10^-3 to 4.0×10^-3. While an anticyclonic eddy imposed upstream (southeast of Kyushu) can cause the transition from an alongshore to a meandering path, it occurs most easily when (). The transition from a meandering to an alongshore path requires an eddy imposed downstream (east of the meandering segment) which suppresses redevelopment of the meandering segment and breaks the balance between the advective and beta effects. Applicability of the results to previously observed path variations is discussed.     

Energy Exchange between Jets of the Antarctic Circumpolar Current and Synoptic Eddies in the Drake Passage and Scotia Sea

We study the energy exchange between jets of the Antarctic Circumpolar Current (ACC) and synoptic eddies generated by them in the surface layer of the ocean in the Drake Passage and Scotia Sea based on 22-year-long satellite altimetry time series from the French CLS Agency (DT Global–MADT–Upd product, http://www.aviso.altimetry.fr) under the assumption, based on observations, that each of the jets is confined between two fixed contour lines of the absolute dynamic topography of the ocean. We calculate and analyze the 22-year evolution of the kinetic energy of each ACC jet and cyclonic and anticyclonic eddies generated by it. We demonstrate the fundamental dependence of fluctuations in jet energy on the phase of their meander and eddy formation, as well as their back absorption by jets. We calculate the mean and extreme energetic characteristics of jets and eddies and compare the jets in terms of the intensity of the generated eddies.     

Nonlinear Creep and Swelling Behavior of Bentonite Mixed with Different Sand Contents Under Oedometric Condition

Due to the viscous nature, especially the swelling property, bentonite is selected as buffer and backfill material for marine radioactive waste disposal. This article presents experimental results from one-dimensional oedometric compression tests with multi-stage loading and unloading/reloading on specimens of bentonite mixed with different sand contents. The sand contents are 50%, 60%, 70%, 80%, and 90%, respectively, in terms of the dry mass of sand over the dry mass of the mixture. Influences of sand contents on time-dependent stress-strain behavior of the bentonite-sand mixture under 1-D straining are examined with special attention to the nonlinear creep and swelling behavior. It is found that the bentonite-sand mixtures exhibit distinct nonlinear creep and swelling behavior that cannot be neglected. Moreover, a nonlinear creep function proposed by Yin (1999 Yin , J.-H. 1999 . Non-linear creep of soils in oedometer tests . Geotechnique 49 ( 5 ): 699 – 707 .[Crossref], [Web of Science ®] , [Google Scholar]) for Hong Kong Marine Deposits has been modified to describe the nonlinear creep behavior of the bentonite-sand mixture.     

Nutricline shoaling in the eastern Pacific warm pool during the last two glacial maxima

Isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs) and alkenones were analyzed in sediment samples retrieved from Ocean Drilling Program Site 1241 covering the last 150000 years to understand the hydrological evolution of the eastern Pacific warm pool (EPWP). GDGT and alkenone concentrations showed higher values in marine isotope stage (MIS)-2 and MIS-6, which suggests the enhancement of primary production at glacial maxima. $ {\text{TEX}}_{86}^{\text{H}} $ - and $ U_{ 3 7^\prime }^{\text{K}} $ -derived temperature depicted different temperature evolutions. $ U_{ 3 7^\prime }^{\text{K}} $ -derived temperature was marked by small variation during the glacial–interglacial cycles, whereas $ {\text{TEX}}_{86}^{\text{H}} $ -derived temperature showed pronounced glacial–interglacial variation that was similar to Mg/Ca-derived temperature records from nearby cores in the EPWP. Given that enhanced primary production during glacial maxima suggests nutricline shoaling, unchanged $ U_{ 3 7^\prime }^{\text{K}} $ over glacial–interglacial cycles can be interpreted as the shift of alkenone production depth. $ {\text{TEX}}_{86}^{\text{H}} $ seems not to be influenced by glacial–interglacial changes in nutricline depths, recording an integrated temperature in surface and thermocline water. The shallow nutricline in the EPWP during glacial maxima most likely reflected the intense formation of Antarctic intermediate water.