A numerical study of tidal asymmetry: preferable asymmetry of nonlinear mechanisms in Xiangshan Bay,East China Sea

In-situ measurements in Xiangshan Bay, the East China Sea, show that the duration of the rising tide is shorter than that of the falling tide around the bay mouth, while it becomes much longer in the inner bay. A finite volume coastal ocean model (FVCOM) with an unstructured mesh was applied to simulate the asymmetric tidal field of Xiangshan Bay. The model reproduced the observed tidal elevations and currents successfully. Several numerical experiments were conducted to clarify the roles of primary mechanisms underlying the asymmetric tidal field. According to the model results, the time-varying channel depth and nonlinear advection prefer shorter duration of the rising tide in Xiangshan Bay, while the time-varying bay width favors longer duration of the rising tide. The overtides generated by these two opposite types of nonlinear mechanisms are out of phase, resulting in smaller M₄ amplitude than the sumfold of each individual contribution. Although the bottom friction as a nonlinear mechanism contributes little to the generation of overtide M₄, it is regarded as a mechanism that could cause a shorter duration of the rising tide, for it can slow down the M₂ phase speed much more than it slows down the M₄ phase speed. The time-varying depth, nonlinear advection and bottom friction are dominating factors around the bay mouth, while the time-varying width dominates in the inner bay, causing the tidal elevation asymmetry to be inverted along the bay.     

ON THE BEHAVIORS OF OCEAN BOTTOM EDDY

On the basis of scaling analysis,this paper presents a frontal geostrophic dynamic system which canbe used to describe the dynamics of an isolated bottom eddy and shows that five possible dynamicsubregimes are included in the system.A particle-in-cell method is used to investigate the behavior of the eddy in the two dynamic subregimes,with 1)nonlinear effect(ε)relevant to the sloping bottom-inducedbeta effect(β_1)and 2)seamount-induced beta effect(β_2)dominant.In the dynamic subregime of β_1～ ε～10~(-2)and β_2=0,the eddy can keep as a whole for a long time and migrates for a long distance.Inthe dynamic subregime where the seamount-induced beta effect (β_2～10~(-1))is dominant,the seamounttraps the water particles,causing the eddy to break up into a discrete set of eddies moving around the seamount.     

Observation and numerical experiments for drag coefficient under typhoon wind forcing

This paper presents a study on drag coefficients under typhoon wind forcing based on observations and numerical experiments. The friction velocity and wind speed are measured at a marine observation platform in the South China Sea. Three typhoons: SOULIK(2013), TRAMI(2013) and FITOW(2013) are observed at a buoy station in the northeast sea area of Pingtan Island. A new parameterization is formulated for the wind drag coefficient as a function of wind speed. It is found that the drag coefficient(Cd) increases linearly with the slope of 0.083′10~(-3) for wind speed less than 24 m s~(-1). To investigate the drag coefficient under higher wind conditions, three numerical experiments are implemented for these three typhoons using SWAN wave model. The wind input data are objective reanalysis datasets, which are assimilated with many sources and provided every six hours with the resolution of 0.125?×0.125?. The numerical simulation results show a good agreement with wave observation data under typhoon wind forcing. The results indicate that the drag coefficient levels off with the linear slope of 0.012′10~(-3) for higher wind speeds(less than 34 m s~(-1)) and the new parameterization improvese the simulation accuracy compared with the Wu(1982) default used in SWAN.     

Study on the Determination of Trace Ni（Ⅱ） by the Catalytic Kinetic Spectrophotometric Method

A new kinetic spectrophotometric method has been developed for the determination of trace Ni (II) in natural water.The method is based on the catalytic effect of Ni (II) on the oxidation of weak acid brilliant blue dye (RAWL) by KIO4 in acid medium.The concentration of nickel (II) can be determined spectrophotometrically by measuring the decrease of absorbance of RAWL at λ= 626 nm using the fix-time method.The influencing factors are investigated by the orthogonal experimental design.The obtained optimum analytical conditions are:pH=2.00,c RAWL=5.00×10-5 mol L-1,cKIO 4= 2.00×10-5 mol L-1,the reaction time t=10 min and the temperature T=25℃..Under the optimum conditions,the developed method allows the measurement of Ni (II) in a range of 0-40.0 ng mL-1.The standard deviation of eleven independent measurements of blank reaction is S=3.08×10-3 and the limit of detection is 2.20 ng mL-1.The relative standard deviations (RSDs) in six replicate determinations of 5 ng mL-1 and 8 ng mL-1 Ni (II) are 2.87% and 1.11%,respectively.Moreover,the experiments show few cations and anions can interfere with the measurement of Ni (II).The recovery efficiencies of this method are in a range of 97.0%-102.5% in freshwater samples.But there is a decreasing effect,which is about 0.2 times the added Ni (II) in seawater medium.After reasonable calibration this processing method is used for the determination of Ni (II) in seawater samples successfully.The results show this developed method has high accuracy and precision,high sensitivity,large range of linearity and high speed.The method can,therefore,be employed at room temperature.     

Structural studies on the oligomers from the polysaccharide ofGracilaria textorii (Rhodophyta) using β-agarase and13C-NMR spectroscopy

The 60% ethanolic extract fromGracilaria textorii (Rhodophyta) was degraded with β-agarase, and certain charged (sulfated) and neutral oligosaccharides were separated by using DEAE Sephadex A25 and Bio- gel P6, P2 chromatographic techniques. Some of the charged oligomers were verified to be neoagarotetraose-6³-sulfate (DP2), neoagarohexaose-6³, 6⁵-disulfate (DP3) and neoagarooctaose-6³, 6⁵, 6⁷-trisulfate (DP4) by using¹³C- and¹H-NMR spectroscopy. One neutral oligomer was assumed to be a mixture of methylated neoagarotetraoses (DP2) by¹H-NMR spectroscopy. These oligomers assigned by their chemical shifts may be used as the model compounds for the structural investigation of the agar-type sulfated polysaccharides using the β-agarase degradation method. Contribution No. 1749 from the Institute of Oceanology, Academia Sinica.     

Nonlinear interactions among internal tidal waves in the northeastern South China Sea

We used a set of 75-day long ADCP data from the northeastern South China Sea (SCS) to investigate nonlinear interactions among freely propagating internal tidal waves. The kinetic energy spectra displayed significant peaks at some higher tidal frequencies, such as O₁M₂ (O₁+M₂), and M₄ (M₂+M₂), where O₁ is the lunar diurnal internal tide, M₂ is the lunar semidiurnal internal tide, and M₄ is the first higher harmonic frequency of M₂. These higher tidal harmonic frequency peaks, as well as the fundamental tidal harmonic peaks, show a σ ^−2.3 spectral falloff rate with frequency. In addition, we explored the possible generation mechanism of higher tidal harmonics. Analysis on the rotary and bicoherence spectra suggests that strong forced non-resonant interaction induced by nonlinear advections was the dominant physical mechanism that induced these higher tidal harmonics. Moreover, the energetic, freely propagating semidiurnal (M₂) internal tidal wave played the most crucial role in these interactions. These results indicate that strong nonlinear forced non-resonant interactions among internal tides can be one of the processes responsible for the redistribution of energy in the internal wave spectrum.     

The inhibitory degree between Skeletonema costatum and dinoflagllate Prorocentrum donghaiense at different concentrations of phosphate and nitrate/phosphate ratios

Interactions between Skeletonema costatum (S. costatum) and Prorocentrum donghaiense (P. donghaiiense) were inves-tigated using bi-algal cultures at different concentrations of phosphate (PO4-P) and nitrate/phosphate (N/P) ratios. Experiments were conducted under P-limited conditions and the Lotka-Volterra mathematical model was used to simulate the growth of S. costatum and P. donghaiense in the bi-algal cultures. Both of these two species were inhibited significantly in bi-algal culture. The results of the simulation showed that the inhibitory degree of S. costaum by P. donghaiense was high when the concentration of PO4-P was low (0.1μmolL-1/2 d), but that of P. donghaiense by S. costaum was high with increased PO4-P supply (0.6μmolL-1/2 d). At low concen-tration of PO4-P (0.1μmolL-1/2 d), or high concentration of PO4-P (0.6μmolL-1/2 d) with high N/P ratio (160), the interactions be-tween S. costatum and P. donghaiense were dependent on the initial cell densities of both species. At high concentration of PO4-P (0.6μmolL-1/2 d) with low N/P ratio (25 or 80), S. costatum exhibited a survival strategy superior to that of P. donghaiense. The de-gree of inhibition of P. donghaiense by S. costaum increased with elevated N/P ratio when the medium was supplemented with con-centration 0.1μmolL-1/2 d of PO4-P. The degree of inhibition to P. donghaiense by S. costaum increased with elevated N/P ratio at low concentration of PO4-P (0.1 μmolL-1/2 d). This trend was conversed at high concentration of PO4-P (0.6μmolL-1/2 d). However, the degree of inhibition of S. costaum by P. donghaiense increased with the increased N/P ratio at different PO4-P concentrations (0.1μmolL-1/2 d and 0.6μmolL-1/2 d). These results suggested that both phosphate concentration and N/P ratio affected the competition between S. costaum and P. donghaiense: P. donghaiense is more competitive in environments with low phosphate or high N/P ratio and the influence of N/P ratio on the competition was more significant with lower phosphate concentration.     

Application of swat model in the upstream watershed of the Luohe River

1INTRODUCTIONIntheHuanghe(Yellow) Riverbasin, soilerosionisaseriousproblem,whilerunoffandsedimentyieldsim-ulation hasnotbeenextensivelystudiedonthebasisofGIS(GeographicInformationSystem) and dis-tributedhydrologicalmodel.Inthisstudy,theLushiwatershed,whichislocatedattheupstreamoftheLushiHydrologicalStationintheLuoheRiver—thebiggesttributary oftheHuanghe Riveranddown-streamofXiaolangdiDam,isselectedasthestudyarea.ThelevelofsoilerosioninLushiwatershedismoderatein theHuangheRiverbas…     

Maximum Sustainable Yield Estimates of Ladypees,Sillago sihama(Forsskl),Fishery in Pakistan Using the ASPIC and CEDA Packages

Using surplus production model packages of ASPIC(a stock-production model incorporating covariates) and CEDA(Catch effort data analysis),we analyzed the catch and effort data of Sillago sihama fishery in Pakistan.ASPIC estimates the pa-rameters of MSY(maximum sustainable yield),Fmsy(fishing mortality),q(catchability coefficient),K(carrying capacity or unexploited biomass) and B1/K(maximum sustainable yield over initial biomass).The estimated non-bootstrapped value of MSY based on logistic was 598 t and that based on the Fox model was 415 t,which showed that the Fox model estimation was more conservative than that with the logistic model.The R2 with the logistic model(0.702) is larger than that with the Fox model(0.541),which indicates a better fit.The coefficient of variation(cv) of the estimated MSY was about 0.3,except for a larger value 88.87 and a smaller value of 0.173.In contrast to the ASPIC results,the R2 with the Fox model(0.651-0.692) was larger than that with the Schaefer model(0.435-0.567),indicating a better fit.The key parameters of CEDA are:MSY,K,q,and r(intrinsic growth),and the three error assumptions in using the models are normal,log normal and gamma.Parameter estimates from the Schaefer and Pella-Tomlinson models were similar.The MSY estimations from the above two models were 398 t,549 t and 398 t for normal,log-normal and gamma error distributions,re-spectively.The MSY estimates from the Fox model were 381 t,366 t and 366 t for the above three error assumptions,respectively.The Fox model estimates were smaller than those for the Schaefer and the Pella-Tomlinson models.In the light of the MSY estimations of 415 t from ASPIC for the Fox model and 381 t from CEDA for the Fox model,MSY for S.sihama is about 400 t.As the catch in 2003 was 401 t,we would suggest the fishery should be kept at the current level.Production models used here depend on the assumption that CPUE(catch per unit effort) data used in the study can reliably quantify temporal variability in population abundance,hence the mod-eling results would be wrong if such an assumption is not met.Because the reliability of this CPUE data in indexing fish population abundance is unknown,we should be cautious with the interpretation and use of the derived population and management parameters.     

Effect of Temporal Resolution of NDVI on Potential Evapotranspiration Estimation and Hydrological Model Performance

Normalized difference vegetation index (NDVI) data, obtained from remote sensing information, are essential in the Shuttleworth-Wallace (S-W) model for estimation of evapotranspiration. In order to study the effect of temporal resolution of NDVI on potential evapotranspiration (PET) estimation and hydrological model performance, monthly and 10-day NDVI data set were used to estimate potential evapotranspiration from January 1985 to December 1987 in Huangnizhuang catchment, Anhui Province, China. The differences of the two calculation results were analyzed and used to drive the block-wise use of the TOPMODEL with the Muskingum-Cunge routing (BTOPMC) model to test the effect on model performance. The results show that both annual and monthly PETs estimated by 10-day NDVI are lower than those estimated by monthly NDVI. Annual PET from the vegetation root zone (PETr) lowers 9.77%-13.64% and monthly PETr lowers 3.28%-17.44% in the whole basin. PET from the vegetation interception (PETi) shows the same trend as PETr. In addition, temporal resolution of NDVI has more effect on PETr in summer and on PETi in winter. The correlation between PETr as estimated by 10-day NDVI and pan measurement (R2= 0.835) is better than that between monthly NDVI and pan measurement (R2 = 0.775). The two potential evapotranspiration estimates were used to drive the BTOPMC model and calibrate parameters, and model performance was found to be similar. In summary, the effect of temporal resolution of NDVI on potential evapotranspiration estimation is significant, but trivial on hydrological model performance.     

Characteristics of earthquake-triggered landslides and post-earthquake debris flows in Beichuan County

Field investigations and aerial photography after the earthquake of May 12,2008 show a large number of geo-hazards in the zone of extreme earthquake effects.In particular,landslides and debris flows,the geo-hazards that most threaten post-disaster reconstruction,are widely distributed.We describe the characteristics of these geo-hazards in Beichuan County using high-resolution remote sensing of landslide distribution,and the relationships between the area and volume of landslides and the peak-discharges of debris flows both pre-and post-earthquake.The results show:1) The concentration(defined as the number of landslide sources per unit area:Lc) of earthquaketriggered landslides is inversely correlated with distance from the earthquake(DF) fault.The relationship is described by the following equation:Lc = 3.2264exp(-0.0831DF)(R2 = 0.9246);2) 87 % of the earthquake-triggered landslides were less than 15×104 m2 in area,and these accounted only for 50% of the total area;84% of the landslide volumes were less than 60×104 m3,and these accounted only for 50% of the total volume.The probability densities of the area and volume distributions are correlated:landslide abundance increases with landslide area and volume up to maximum values of 5 × 104 m2 and 30 × 104 m3,respectively,and then decreases exponentially.3) The area(AL) and volume(VL) of earthquake-triggered landslides are correlated as described with the following equation:VL=6.5138AL1.0227(R2 = 0.9131);4) Characteristics of the debris flows changed after the earthquake because of the large amount of landslide material deposited in the gullies.Consequently,debris flow peak-discharge increased following the earthquake as described with the following equation:Vpost = 0.8421Vpre1.0972(R2 = 0.9821)(Vpre is the peak discharge of pre-earthquake flows and the Vpost is the peak discharge of post-earthquake flows).We obtained the distribution of the landslides based on the above analyses,as well as the magnitude of both the landslides and the post-earthquake debris flows.The results can be useful for guiding post-disaster reconstruction and recovery efforts,and for the future mitigation of these geo-hazards.However,the equations presented are not recommended for use in site-specific designs.Rather,we recommend their use for mapping regional seismic landslide hazards or for the preliminary,rapid screening of sites.     

Numerical study on the bifurcation of the North Equatorial Current

A 1.5-layer reduced-gravity model forced by wind stress is used to study the bifurcations of the North Equatorial Current(NEC).The authors found that after removing the Ekman drift,the modelled circulations can serve well as a proxy of the SODA circulations on the σθ=25.0 kg m～-3 potential density surface based on available long-term reanalysis wind stress data.The modelled results show that the location of the western boundary bifurcation of the NEC depends on both zonal averaged and local zero wind stress curl latitude.The effects of the anomalous wind stress curl added in different areas are also investigated and it is found that they can change the strength of the Mindanao Eddy(ME),and then influence the interior pathway.     

210Pb chronology and trace metal geochemistry in the intertidal sediment of Qinjiang River estuary, China

Historical records of metal inputs were studied by using a sediment core collected from a sand-rich mudflat in the Qinjiang River estuary, China. 210Pb chronology was used to reconstruct the fluxes of Hg, Cu, Pb, Zn, Cd, Cr and As to the core site during the last 86 years. Based on the constant initial concentration model, the sedimentation rates are 1.18 cm year-1 in the top 30 cm sandy layer and 0.92 cm year-1 in the muddy bottom layer. To compensate for grain-size and mineralogy effects on metal concentra-tions, aluminum was used as the normalizing element. The enrichment factors (EF) indicate that the natural inputs had prevailed up to the early 1980s. After this period, the intensity of human activities has resulted in continual increasing trend of metals towards the surface. Recent sediment samples from the Qinjiang River estuary are found moderately enriched by Cd (EF>1.5) and slightly enriched by other metals (EF<1.5). Considering that the drainage area of the Qinjiang River is mostly agricultural land, the increased Cd may be due to the usage of fertilizers and pesticides in agricultural activities and the combustion of fossil fuels.     

Acute Toxicity of Formalin to the Red-tide Dinoflageilate Prorocentrum minimum Schiller （Protozoa, Mastigophora）

The acute toxic levels of formalin to a marine red-tide dinoflagellate Prorocentrum minimum Schiller (Mastigophora, Dinoflagellida) were determined by linear regression analysis. The data obtained revealed that there is a close regression relationship between logarithmic concentrations of formalin and mortality probit scales of the organism when exposure times are 1, 2, 6, 12, 24 h. The median lethal concentrations (LC50 values) obtained from correlation analysis for these exposure times (with 95% confidence intervals), were 11.83, 6.76, 4.37, 4.27 and 3.98mgL^-1 respectively. A toxicity curve was obtained to connect the exposure time and LC50 value from the correlation between them. The results indicated that P. minimum could be killed or fixed by formalin in the concentration range of 4-12mgL^-1 within 1-24h.     

Numerical analysis of the nonlinear parameterization of waves in currents over a submerged sill with a non-hydrostatic model

An investigation of the effects of a uniform current strength direction (following or opposing wave propagation) on the nonlinear transformation of irregular waves over a submerged trapezoidal sill is carried out using SWASH, a non-hydrostatic numerical wave model. The nonlinear parameters (i.e., asymmetry, skewness, and kurtosis) are calculated, and the empirical formulas for these parameters are presented as a function of the local Ursell number based on the present numerical data measured. In the shoaling area of the submerged sill, the nonlinear characteristics of waves are more obvious when waves propagate in the same direction as the currents than when waves propagate in the opposite direction. Whereas nonlinear parameters grow with the strengthening of the following currents over the crest, they tend to decrease as the adverse current velocity increases over the crest area of the submerged sill.     

THE DYNAMIC EFFECTS OF A FINITE DEPTH AMBIENT FLUID ON BOTTOM EDDY

Scaling analysis shows that if o(ε~2, β_1ε, γε) ~o (δ), frontal geostrophic dynamics governs the behav-ior of an isolated bottom eddy in a finite depth ambient fluid; and that the ambient flow induced bybottom eddy migration satisfies quasi-geostrophic dynamics. This two layer model includes the impor-tant processes of advection of bottom eddy due to ambient flow, baroclinic instability, and form dragintroduced by Rossby waves. The numerical results show that the three processes enhance theinstability and alter the migration speed of the bottom eddy, and that the form drag induces asignificant meridional drift of the eddy.     

Particle image velocimetry and numerical simulations of the hydrodynamic characteristics of an artificial reef

This article reports a particle image velocimetry study and the comparative results of a numerical simulation into the hydrodynamic characteristics around an artificial reef.We reveal the process of flow separation and vortex evolution,and compare the force terms generated by our artificial reef model.The numerical simulation agrees well with experimental results,showing the applicability of computational fluid dynamics to the hydrodynamics of an artificial reef.Furthermore,we numerically simulate the hydrodynamics of the reef model for seven velocities.The results show that the drag coefficient is approximately 1.21 in a self-modeling region for Reynolds numbers between 2.123×104and 9×104.Therefore,the upwelling height and current width of the flow field do not change significantly when the inflow velocity increases.Our study indicates that computational fluid dynamics can be applied to study the hydrodynamics of an artificial reef and offer clues to its construction.     

Episodes of volcanic activity and their environmental effects in the Okinawa Trough during the last 150 ka

1 INTRODUCTIONA large-scale volcanic eruption can injectenormous amount of ash and dust into the atmos-phere, leading to local and short cooling due toshielding incoming radiation. During the past 2 to 3decades, accumulated evidence indicated thatlarge-scale volcano eruptions or series of volcanicactivity appears to coincide with climatologicaltransitions. For example, Robock (1979) indicatedthat volcanic eruptions played a role in forcing "Lit-tle Ice Age". Well-dated high-resolution co…     

Dynamic model for oil slick dispersion into a water column—A wind-driven wave tank experiment

A dynamic experiment for oil dispersion into a water column was performed with a 21 m long, 0.5 m wide, and 1 m high wind-driven wave tank. At wind velocity between 6–12 m/s and with the oil slick kept constant (about 1 μm), the rate of the oil content increase in the water column could be approximated from the difference between the dispersion rate (R) of the oil slick and the coagulation rate (R′) of the dispersed oil slick. Assuming the coagulation rate is directly proportional to the concentration of the water dispersed oil slick (i. e.R′-K ^*C), the integral form of the dynamic model can be expressed asC=R ^* [1−exp(−K ^*t)]/K and parametersR andK can be regressed with a computer. The relative deviation of model results from the experimental data was mainly less than 10%. The oil slick dispersion rate (R) had exponential relationship with the wind velocity (U), and can be fitted with a formulaR=A ^* (U+1) ^B . The fitted constant of the coagulation rate,K(0.8–3.0^* 10⁻³ min⁻¹) did not have significant relationship with the wind velocity.