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.     

Size trends and population characteristics of striped marlin,Tetrapturus audax caught in the New Zealand recreational fishery

Size trends and population characteristics of striped marlin, Tetrapturus audax, caught in the New Zealand recreational fishery were evaluated using club records from the Bay of Islands Sword‐fish Club (BOISC) from 1925 to 2003 (n = 15 114) and biological samples from northern New Zealand collected from 1985 to 1994 (n = 684). Average weight (kg) ± SE of striped marlin declined 6–9 kg every 20 years between 1925 (124.2 ± 1.0 kg) and 2003 (100.1 kg ± 0.4). Weight (H14843 = 444.58, P < 0.001) and condition (H443 = 37.54, P < 0.001) increased significantly throughout the recreational fishing season (December‐May). No differences were observed in the lower jaw‐fork length (LJFL)‐weight (W) relationship between sexes or years but average weight of females (106.2 ±1.1 kg) was significantly greater than males (90.2 ±1.2 kg) from 1985 to 1994 (F684 = 88.37, P < 0.001). Negative allometric growth (W = 2E ‐08 LJFL^2.88) was recorded for all striped marlin. Individual growth was modelled by fitting backcalculated LJFL to eight age classes derived from dorsal spine measurements (n = 94) using the von Bertalanffy growth equation (r² = 0.96). The parameters estimated in the von Bertalanffy growth equation were asymptotic length (L ) = 3010 mm, annual growth rate (K) = 0.22, age (yr) at hypothetical length 0 (t0) = ‐0.04.     

Uptake of nitrate and ammonium,and distribution of related variables,in the upwelled plume of western Cook Strait/Taranaki Bight,New Zealand

To distinguish the manner of spread of upwelling effects in the Cook Strait/Taranaki Bight (New Zealand) region, nitrogen uptake and the distribution of nutrients and chlorophyll a are described. NO3 concentrations were closely linked to upwelled water, but this was distributed irregularly and not necessarily with the upwelling focus. NH4, SRP, and chlorophyll a showed varying degrees of association with recycling processes, but linear correlations were not present. NO3 uptake was strongly inhibited by ambient NH4, but not completely. There was no evidence of either N or P limitation of total planktonic biomass, and there was indirect evidence from C/N uptake and cellular ratios that a high growth rate prevailed. N uptake increased with irradiance (= decreasing depth) similar to photosynthesis, but without high light inhibition, and could be related to light by a simple equation. The irregularity of distribution of features suggests that future modelling of the phytoplankton biomass in the region could benefit from analysis of event occurrence rather than of evenly progressive changes.     

Tidal observations on the West Coast,South Island,New Zealand

Harmonic tidal constants, calculated from sea surface elevation observations at Jackson Bay on the West Coast of the South Island, are consistent with available semi‐diurnal and diurnal tidal phase distributions. Current observations taken over a 111 day period at mid‐depth in 1505 m of water on the southern flank of the Challenger Plateau and over a 240 day period in 1430 m of water on the South Island western coast continental slope, are subject to tidal analysis. At both sites there is a component of energy flux directed across the isobaths and only at the northern site for the M2 tide is the phase consistent with a dominant progressive barotropic tide. The successive 30 day harmonic constants at the southern continental slope site exhibit a trend in the M2 tidal ellipse speed and ellipticity suggesting the presence of a regular internal tide. Superposition of ‘internal tidal’ and barotropic tidal flows, as prescribed from progressive‐ and standing‐wave elevations, to fit the observations indicates that the ‘internal tide’ is probably associated with the first baroclinic mode. At the current‐meter depths the speeds of the ‘internal tide’ for the M2 tide are about the same as the barotropic speeds, whereas, the S2 ‘internal tide’ speeds are larger than those of the barotropic tide. The consistency of the trend in ellipse parameters lends support to the theoretical progressive trapped barotropic tidal flows being a good approximation to the actual barotropic tide. Some support for the hypothesis that the S2 tide on the West Coast of New Zealand has a substantial standing wave contribution is given by the northern observations, where the ratio of the S2: M2 internal tidal ellipse current amplitudes are substantially larger than the ratio of the elevations, the internal tide being generated by across‐isobath flows.     

Nitrate levels in subsurface waters of the Waimea Plain,Nelson

In 1969–72, nitrate determinations were made on waters from 122 wells on the Waimea Plain, Nelson, South Island, New Zealand. Four wells showed values above 30 mg NO3‐N/litre with a maximum of 50 mg/litre in one well, fifteen were in the 30–20 mg/litre range, forty‐six were between 20 mg and 10 mg/litre; fifty‐seven were below 10 mg/litre, and in only four wells was the nitrate concentration below 1.0 mg NO3‐N/litre. Where it entered the plain, the Wairoa River usually contained about 0.07 mg NO3‐N/litre, whereas the Waiiti River contained from 0.5 to 2.0 mg NO3‐N/litre.

Mapping of gradients of nitrate level indicated that these might contribute to aquifer definition and to sources of replenishing waters.     

Carbon dioxide fixation in aerial roots of the New Zealand mangrove Avicennia marina var. resinifera (Note)

The CO2 exchange characteristics and chlorophyll content of the pneumatophores of Avicennia marina (Forsk.) Vierh. var. resinifera (Forst.) Bakh. were examined. The results indicate that, in addition to their primary role in ventilation of the buried roots, the aerial roots have the ability to fix exogenous CO2 via photosynthesis. Consequently, during the day these organs are almost auto trophic in terms of carbon balance. When corrected for the respiration of non‐chlorophyllous tissue, they show photosynthetic rates comparable to those of many other woody plants.     

Ventilation and oxygen consumption in the spotted stargazer Genyagnus monopterygius †

As with other benthic fish which bury themselves in sand, in Genyagnus monopterygius it is the branchiostegal apparatus which makes the major contribution to ventilation of the gills. The prolonged abduction of the branchiostegal membranes draws water in through the mouth, and exhaled water moves away from the fish through mucus‐lined channels in the sand. The steady and continuous flow of water during ventilation may avoid detection of the stargazer by prey or predators. Associated with its sluggish habits and the low gill surface area available for gas exchange, Genyagnus has a low rate of oxygen consumption (15.8 ml kg^‐1 h^‐1 for a 100 g fish at 17°C). At partial pressures of oxygen (PO2) below 6 kPa (45 mm Hg) oxygen consumption is reduced, but the fish does not emerge from the sand until the PO2 falls below 2 kPa. There is evidence that the fish can sustain an oxygen debt. Respiration in the stargazer is discussed in relation to the fish's ability to withstand temporary hypoxia when it is stranded above the water level in the intertidal zone.     

Laboratory experiments on the effects of light and temperature on the uptake of nutrients by Lake Rotongaio phytoplankton

Laboratory experiments were carried out to investigate the effect of sustained light and temperature on the uptake of PO4 ³ ‐P, NH/‐N, and NO3 ^‐‐N and on the biomass of the phytoplankton of Lake Rotongaio. Two sets of experiments compared differences between the seasonally contrasting dominant phytoplankton species, which consist of a filamentous nitrogen‐fixing cyanophyte in summer and diatoms in winter. The extent to which nutrients became exhausted, cell N/P ratios changed, and nitrogen‐fixation was enhanced, depended on the amount of light supplied and strongly indicated that light is the limiting factor for phytoplankton growth under normal lake conditions. Re‐illumination after extended periods of darkness suggested that sedimentation below the euphotic zone could be tolerated by the phytoplankton, and that growth could be quickly reestablished upon re‐entry to better light conditions. Temperature could have some effect on metabolic rate but is probably of minor importance compared to light availability for the species selection and overall nutrient dynamics of the lake.     

Availability of soil and sediment phosphorus to a planktonic alga

Chlorophyll production by Chlorella vulgaris Beij. var. vulgaris was used to estimate alga‐available phosphorus in clays, soils, and lake sediments suspended in water at concentrations appropriate to lake inflows during floods (100–500 g/m³). Chlorella apparently used 24–81% of 0.5M H2SO4 extractable phosphorus in clays from topsoils, about 25% from lake sediments, and 0.3–1.0% from, subsoils low in phosphorus and with high phosphorus retention.

The presence of suspended soil material did not reduce the availability to Chlorella of inorganic phosphorus added to the cultures. Increasing the Chlorella population by adding inorganic phosphorus resulted in an apparent increase in availability of phosphorus from the soil, possibly as a result of enzymic mineralisation of organic soil phosphorus. The amount of available phosphorus in lake sediments was not a reliable guide to the trophic condition, of the lake.

Suspended material from sediments, soils, and especially clay eroded from fertilised topsoils may provide phosphorus for algal growth in lakes. If allophanic clays are applied to lakes to sorb phosphorus and hence control eutrophication, the particles must settle out before planktonic algae in the photic zone can use the adsorbed phosphorus.     

McMurdo sound hydrological observations, 1972–73

Hydrological data, including current observations, from a line of five stations across seasonal ice in McMurdo Sound, Antarctica, are presented. The observations, which were made during the 1972–73 summer, gave evidence for a net flux of approximately 0.4 × 10⁶m³.s^‐1 flowing towards the south to pass under the Ross Ice Shelf. There was a small northward flux (0.04 × 10⁶m³.s^‐1) on the western side of the sound, but over most of the width the flow was to the south.

Diurnal tidal currents, which were analysed in terms of K1, P1, and O1 components, were superposed on the drift currents.

Temperature and salinity measurements revealed signs of surface water dilution and warming down to 200 m in the region of southerly drift. Temperatures associated with the northerly drift on the western side of the section were close to freezing at all depths. This water can be identified as coming from under the Ross Ice Shelf.     

Subtidal Galeolaria hystrix (Polychaeta: Serpulidae) reefs in Paterson Inlet,Stewart Island,New Zealand

Serpulid patch reefs of Galeolaria hystrix Mörch, 1863 were found in water depths of 9–16 m in Big Glory Bay, Paterson Inlet, Stewart Island, and here we report preliminary studies of these important habitat‐formers. This is the first observation of this species in subtidal patch reefs; 114 reefs were noted in a survey of 28 000 m². Most reefs were 1–5 m in diameter, and up to 1.5 m high. Up to 65% of the serpulid tubes were occupied by living G. hystrix during a mid‐winter diving survey; 64% of reefs observed were whole, whereas 36% were broken or dead. Radiometric dating of a basal specimen of reef carbonate showed it to be less than 50 years old. Production of high‐Mg calcite, ranging from 9 to 11 wt% MgCO3, by G. hystrix may be as much as 11 kg CaCO3 m^‐2y^‐1, but was not reflected in surrounding sediments, which were dominantly terrigenous muds. A rich reef fauna, both sessile and motile, was associated with the reefs. Further study of these unusual temperate reefs is strongly recommended.     

Factors affecting clarity of New Zealand lakes

Average values of inherent optical properties for the 400–700 nm waveband were estimated from quantum irradiance measurements on 27 New Zealand lakes of diverse optical character, using published nomograms. Secchi disc depths, turbidity, algal pigment, non‐volatile suspended solids, and absorption by membrane‐filtered samples at 440 nm (g 440 ) were also measured. Turbidity (NTU) correlated closely with the scattering coefficient (m^‐1) and these quantities were almost numerically equal, as found in other studies. The data were found to conform to an expression in the oceanographic literature relating Secchi disc depth to the beam attenuation and diffuse attenuation coefficients. Specific beam attenuation, scattering, and absorption coefficients were estimated from the coefficients of linear multiple regressions of the measured total coefficients on the variables: total pigment (chlorophyll a + phaeopigment), non‐volatile suspended solids, and g 440. The estimated values were in reasonable agreement with similar specific coefficients reported in the oceanographic literature. The coefficients provide a basis for predicting clarity in new impoundments or for predicting the effects of loading changes (e.g., of nutrients) on the clarity of existing lakes. The coefficients can be used for classifying lake waters into different optical types.     

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.     

Heterotrophic potentials and trophic status of ten New Zealand lakes

A comparative study of the heterotrophic potentials of ten New Zealand lakes was carried out. Two sampling periods were chosen, one during late August and early September 1972 (late winter, minimal thermal stratification) and another during February 1973 (mid summer, strong thermal stratification in all but two lakes). A wide range of heterotrophic potentials occurred which enabled the lakes to be placed in the following basic groups according to descending trophic levels:

Group I: Rotowhero (thermal, acid lake) >Rotorua>Okaro

Group II: Rotokakahi> Tikitapu >Okareka

Group III: Rotoma = Rotoroa (S.I.) = Rotoiti (S.I.) = Okataina

Maximum values for heterotrophic potential (Vm a x ) in micrograms of glucose per litre per hour were: Maximum heterotrophic activities were found within the metalimnion, except in Lakes Tikitapu and Rotokakahi, where maximum activities occurred within the hypolimnion.     

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.     

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.     

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.     

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.     

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.     

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⁵.