An environmentally based growth model that uses finite difference calculus with maximum likelihood method: its application to the brackish water bivalve Corbicula japonica in Lake Abashiri, Japan

We present a growth analysis model that combines large amounts of environmental data with limited amounts of biological data and apply it to Corbicula japonica. The model uses the maximum-likelihood method with the Akaike information criterion, which provides an objective criterion for model selection. An adequate distribution for describing a single cohort is selected from available probability density functions, which are expressed by location and scale parameters. Daily relative increase rates of the location parameter are expressed by a multivariate logistic function with environmental factors for each day and categorical variables indicating animal ages as independent variables. Daily relative increase rates of the scale parameter are expressed by an equation describing the relationship with the daily relative increase rate of the location parameter. Corbicula japonica grows to a modal shell length of 0.7 mm during the first year in Lake Abashiri. Compared with the attain-able maximum size of about 30 mm, the growth of juveniles is extremely slow because their growth is less susceptible to environmental factors until the second winter. The extremely slow growth in Lake Abashiri could be a geographical genetic variation within C. japonica.

Biology and population dynamics of cowcod (Sebastes levis) in the southern California Bight

Cowcod (Sebastes levis) is a large (100-cm-FL), long-lived (maximum observed age 55 yr) demersal rockfish taken in multispecies commercial and recreational fisheries off southern and central California. It lives at 20–500 m depth: adults (>44 cm TL) inhabit rocky areas at 90–300 m and juveniles inhabit fine sand and clay at 40–100 m. Both sexes have similar growth and maturity. Both sexes recruit to the fishery before reaching full maturity. Based on age and growth data, the natural mortality rate is about M =0.055/yr, but the estimate is uncertain. Biomass, recruitment, and mortality during 1951–98 were estimated in a delay-difference model with catch data and abundance indices. The same model gave less precise estimates for 1916–50 based on catch data and assumptions about virgin biomass and recruitment such as used in stock reduction analysis. Abundance indices, based on rare event data, included a habitat-area–weighted index of recreational catch per unit of fishing effort (CPUE index values were 0.003–0.07 fish per angler hour), a standardized index of proportion of positive tows in CalCOFI ichthyoplankton survey data (binomial errors, 0–13% positive tows/yr), and proportion of positive tows for juveniles in bottom trawl surveys (binomial errors, 0–30% positive tows/yr). Cowcod are overfished in the southern California Bight; biomass during the 1998 season was about 7% of the virgin level and recent catches have been near 20 metric tons (t)/yr. Projections based on recent recruitment levels indicate that biomass will decline at catch levels > 5 t/yr. Trend data indicate that recruitment will be poor in the near future. Recreational fishing effort in deep water has increased and has become more effective for catching cowcod. Areas with relatively high catch rates for cowcod are fewer and are farther offshore. Cowcod die after capture and cannot be released alive. Two areas recently closed to bottom fishing will help rebuild the cowcod stock.

Validation of back-calculation equations for juvenile bluefish (Pomatomus saltatrix) with the use of tetracycline-marked otoliths

In recent years, a decrease in the abundance of bluefish (Pomatomus saltatrix) has been observed (Fahay et al., 1999; Munch and Conover, 2000) that has led to increased interest in a better understanding the life history of the species. Estimates of several young-of-the-year (YOY) life history characteristics, including the importance and use of estuaries as nursery habitat (Kendall and Walford, 1979) and size-dependant mortality (Hare and Cowen, 1997), are reliant upon the accuracy of growth determination. By using otoliths, it is possible to use back-calculation formulae (BCFs) to estimate the length at certain ages and stages of development for many species of fishes. Use of otoliths to estimate growth in this way can provide the same information as long-term laboratory experiments and tagging studies without the time and expense of rearing or recapturing fish. The difficulty in using otoliths in this way lies in validating that 1) there is constancy in the periodicity of the increment formation, and 2) there is no uncoupling of the relationship between somatic and otolith growth. To date there are no validation studies demonstrating the relationship between otolith growth and somatic growth for bluefish. Daily increment formation in otoliths has been documented for larval (Hare and Cowen, 1994) and juvenile bluefish (Nyman and Conover, 1988). Hare and Cowen (1995) found ageindependent variability in the ratio of otolith size to body length in early age bluefish, although these differences varied between ontogenetic stages. Furthermore, there have been no studies where an evaluation of back-calculation methods has been combined with a validation of otolithderived lengths for juvenile bluefish.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.

Survey of plankton blooms and red tides and their interaction with environmental factors through monitoring satellite images

One of the most important marine ecological phenomena is red tide which is created by increasing of phytoplankton population, influenced by different factors such as climate condition changes, utrification hydrological factors and can leave sever and undesired ecological and economical effects behind itself in the case of durability. Coast line of Hormozgan is about 900km from east to west, within the range of geographical coordinates of 56 16 23.8, 26 58 8.8 to 54 34 5.33 and 26 34 32 eastern longitude and northern latitude, seven sampling stations were considered and sampled for a period of one year from October 2008 to October 2009. after the analysis of Satellite images, monthly, during the best time. In several stages, samplings were performed. In each station, three samples were collected for identification and determination of Bloom- creating species abundance. Cochlodinium polykrikoides was the species responsible for the discoloration which occurred at October 2008 in Hormozgan marine water. Environmental parameters such as sea surface temperature, pH, salinity, Dissolved Oxygen concentration, Total Dissolved Solids (T.D.S.), conductivity, nitrate, nitrite and phosphate and also chlorophyll a were measured and calculated. Kruscal Wallis test was used to compare the densities between different months, seasons and the studied stations. Mann-whitney test from Nonparametric Tests was used for couple comparison. Pearson correlation coefficient was used to determine the relationship between physical and chemical data set and the abundance of Cochlodinium polykrikoides. Multivariate Regression and analysis of variance (ANOVA) also were used to obtain the models and equations of red tide occurrence relationship, environmental parameters and nutrient data. The highest density was 26 million cells per liter in Qeshm station. A meaningful difference was observed between sampling months and seasons but there was no between sampling stations which indicates that in favorable conditions, the occurrence of this phenomenon by the studied species is probable. Regarding to β coefficients of nitrate, temperature, phosphate, Total Dissolvable Solutions (T.D.S) and pH these parameters are effective on the abundance of this species and red tide occurrence. Increase in these factors can represent the effects and outcomes of human activities and increase in marine pollution.

Three dimensional circulation of water with variable wind in the Caspian Sea

This research is based on a numerical model for forecasting the three-dimensional behavior of (sea) water motion due to the effect of a variable wind velocity. The results obtained are then analyzed and compared with observation. This model is based on the equations that overcome the current and distribution of temperature by applying the method of finite difference with assuming Δx, Δy as constant and Δz, variable. The model is based on the momentum equation, continuity equation and thermodynamic energy equation and tension at the surface and middle layers and bottom stress. The horizontal and vertical eddy viscosity and thermal diffusivity coefficients we used in accordance with that of the Bennet on Outario Lake (1977). Considering the Caspian Sea dimension in numerical model the Coriolis parameter used with β effects and the approximation Boussines have been used. For the program controlling some simple experiment with boundary condition similar to that of the Caspian Sea have been done. For modeling the Caspian Sea the grid of the field was done as follows: At horizontal surface grid size is 10×10km extension and at vertical in 10 layers with varying thickness from surface to bed respectively as: 5, 10, 20, 3, 50, 100, 150, 200, 25, 500 and higher. The data of wind as velocity، direction and temperature of water related to 15th September 1995 at 6،12 and 18 o’clock were obtained from synoptic station at the Caspian Sea shore and the research marine of Haji Alief. The information concerning shore wind was measured and by the method of SPM (shore protection manual) was transferred to far shore winds through interpolation and by use of inverse square distance of position distribution of the wind velocity at the Caspian surface field was obtained. The model has been evaluated according to the reports and observations. Through studying the position of the current in different layers، the velocity in the cross section in the northern، southern and the middle layers، will be discussed. The results reveal the presence of the circulation cells in the three above mentioned areas. The circulation with depth is reduced too. The results obtained through the numerical solution of the temperature equation have been compared with the observation. The temperature change in different layers in cross section illustrates the relative accordance of the model mentioned.