基于浮游植物吸收的海洋初级生产力模型的不确定性分析

Satellite-derived phytoplankton pigment absorption(a_(ph)) has been used as a key predictor of phytoplankton photosynthetic efficiency to estimate global ocean net primary production(NPP). In this study, an a_(ph)-based NPP model(Ab PM) with four input parameters including the photosynthetically available radiation(PAR), diffuse attenuation at 490 nm(K_d(490)), euphotic zone depth(Z_(eu)) and the phytoplankton pigment absorption coefficient(a_(ph)) is compared with the chlorophyll-based model and carbon-based model. It is found that the Ab PM has significant advantages on the ocean NPP estimation compared with the chlorophyll-based model and carbonbased model. For example, Ab PM greatly outperformed the other two models at most monitoring sites and had the best accuracy, including the smallest values of RMSD and bias for the NPP estimate, and the best correlation between the observations and the modeled NPPs. In order to ensure the robustness of the model, the uncertainty in NPP estimates of the Ab PM was assessed using a Monte Carlo simulation. At first, the frequency histograms of simple difference(δ), and logarithmic difference(δ~(LOG)) between model estimates and in situ data confirm that the two input parameters(Z_(eu) and PAR) approximate the Normal Distribution, and another two input parameters(a_(ph) and K_d(490)) approximate the logarithmic Normal Distribution. Second, the uncertainty in NPP estimates in the Ab PM was assessed by using the Monte Carlo simulation. Here both the PB(percentage bias), defined as the ratio of ΔNPP to the retrieved NPP, and the CV(coefficient of variation), defined as the ratio of the standard deviation to the mean are used to indicate the uncertainty in the NPP brought by input parameter to Ab PM model. The uncertainty related to magnitude is denoted by PB and the uncertainty related to scatter range is denoted by CV.Our investigations demonstrate that PB of NPP uncertainty brought by all parameters with an annual mean of5.5% covered a range of –5%–15% for the global ocean. The PB uncertainty of Ab PM model was mainly caused by a_(ph); the PB of NPP uncertainty brought by a_(ph) had an annual mean of 4.1% for the global ocean. The CV brought by all the parameters with an annual mean of 105% covered a range of 98%–134% for global ocean. For the coastal zone of Antarctica with higher productivity, the PB and CV of NPP uncertainty brought by all parameters had annual means of 7.1% and 121%, respectively, which are significantly larger than those obtained in the global ocean. This study suggests that the NPPs estimated by Ab PM model are more accurate than others, but the magnitude and scatter range of NPP errors brought by input parameter to Ab PM model could not be neglected,especially in the coastal area with high productivity. So the improving accuracy of satellite retrieval of input parameters should be necessary. The investigation also confirmed that the SST related correction is effective for improving the model accuracy in low temperature condition.     

Determining the influence of wind-wave breaking on the dissipation of the turbulent kinetic energy in the upper ocean and on the dependence of the turbulent kinetic energy on the stage of wind-wave development

New experimental data that make it possible to explain and predict the observed variability of turbulent-energy dissipation in the upper ocean are discussed. For this purpose, the dependence of the energy dissipation rate of breaking wind waves on their propagation velocity (see [1]) is used. The turbulent-energy dissipation values obtained earlier in [2, 3] by a direct method are compared to the results of radar measurements of individual breaking events presented in [1]. On the basis of this comparison, a strong dependence of the turbulent-energy dissipation value on the stage of wind-wave development, which is characterized by the ratio U _a /c _p (U _a is the wind speed and c _p is the phase speed of the peak of the wind-wave spectrum) is confirmed. This dependence was found earlier purely empirically. Moreover, it is shown that the theoretically obtained dependence (c _p /U _a )⁴, does not contradict the available empirical data. The results of this study opens possibilities for scientifically substantiated calculations of greenhouse-gas exchange (specifically, CO₂ exchange between the ocean and the atmosphere).     

Effect of the nonlinear interaction of tidal harmonics on their spatial structure as applied to the system of the white and Barents Seas

Using the three-dimensional nonlinear finite-element thermohydrodynamic model QUODDY-4, we obtain two solutions giving an idea of the role of the nonlinear interaction of tidal harmonics in the formation of their spatial structure. The first solution is induced by the total tide (M ₂ + S ₂ + K ₁ + O ₁) at the open boundary and by the total static tide inside the area under consideration; then, the solution obtained is subject to harmonic analysis. The second solution is obtained by specifying such tidal sea level elevations at the open boundary that meet individual tidal harmonics. These two solutions are compared. It is shown that the differences between the solutions for the S ₂, K ₁, and O ₁ waves can be significant, especially near the open boundary between the White and Barents seas. This conclusion remains valid also for the maximum velocities (major semiaxes of the ellipses) of the barotropic (depth-average) tidal current as well as for the average (over the tidal cycle) densities of the total tidal energy and components of the tidal-energy budget. The emergence of this feature indicates that there are resonance modes with frequencies that differ from those of the S ₂, K ₁, and O ₁ harmonics to a lesser extent than the M ₂ harmonic frequency. The same conclusion can be made by comparing the values of the amplification factor, which is defined as the ratio between the actual and static tides, for the system of the Barents and White seas as a whole and for Mezen Bay in the White Sea and Czech Bay in the Barents Sea, taken separately.     

Oceanic processes of upper ocean heat content associated with two types of ENSO

The spatiotemporal variability of equatorial Pacific upper ocean heat content (HC) and subsurface heat during two types of El Niño-Southern Oscillation (ENSO), namely eastern and central Pacific (EP and CP) types, is investigated using subsurface ocean heat budget analysis. Results show that HC tendencies during both types of ENSO are mainly controlled by oceanic heat advection beneath the mixed layer to the thermocline, and the role of net surface heat flux can be neglected. The most important three terms are the zonal and vertical advections of anomalous heat by climatological currents (QU ₀ T′, QW ₀ T′) and zonal advection of climatological heat by anomalous current (QU′T ₀). The large contribution of QU ₀ T′ extends from west to east along the equatorial Pacific. The considerable contribution of QU′T ₀ is confined to the east of 160°W, and that of the QW ₀ T′ is observed in the central Pacific between 180°E and 120°W. In particular, a major contribution of QW ₀ T′ is also observed in the far eastern Pacific east of 100°W during EP ENSO. There is also a small contribution from meridional advection of climatological heat by anomalous current (QV′T ₀). In contrast, the meridional advection of anomalous heat by climatological currents (QV ₀ T′) and vertical advection of climatological heat by anomalous current (QW′T ₀) are two damping factors in the HC tendency, with the former dominating. Differences in spatial distribution of the heat advection associated with the two types of ENSO are also presented. We define a warm water heat index (WWH) as integrated heat content above 26 kg m^?3 potential density (26σ _? ) isopycnal depth within 130°E–80°W and 5°S–5°N. Further examination suggests that the recharge–discharge of WWH is involved in both types of El Niño, though with some differences. First, it takes about 42 (55) months for the evolution of a recharge–discharge cycle during an EP (CP) ENSO. Second, the EP El Niño event peaks during the discharge phase, 7–8 months after the recharge time. The CP El Niño peaks during the recharge phase, 4–5 months before the recharge time. The locations of HC anomalies in the El Niño mature phase relative to those at recharged time explain why the EP and CP El Niño peak in different stages of the recharge–discharge process.     

2012年夏季海南岛东岸上升流区的混合观测

The turbulent mixing in the upwelling region east of Hainan Island in the South China Sea is analyzed based on in situ microstructure observations made in July 2012. During the observation, strong upwelling appears in the coastal waters, which are 3℃ cooler than the offshore waters and have a salinity 1.0 greater than that of the offshore waters. The magnitude of the dissipation rate of turbulent kinetic energy ε in the upwelling region is O(10–9 W/kg), which is comparable to the general oceanic dissipation. The inferred eddy diffusivity K_ρ is O(10–6 m~2/s), which is one order of magnitude lower than that in the open ocean. The values are elevated to K_ρ≈O(10–4 m~2/s) near the boundaries. Weak mixing in the upwelling region is consistent with weak instability as a result of moderate shears versus strong stratifications by the joint influence of surface heating and upwelling of cold water.The validity of two fine-scale structure mixing parameterization models are tested by comparison with the observed dissipation rates. The results indicate that the model developed by Mac Kinnon and Gregg in 2003 provides relatively better estimates with magnitudes close to the observations. Mixing parameterization models need to be further improved in the coastal upwelling region.     

Features of the atmospheric boundary layer block for three versions of the WAM wind wave model

The estimated characteristics of the atmospheric boundary layer, obtained by the simulation of wind wave fields using three versions of the WAM numerical model are compared with the well-known empirical dependences of drag coefficient C _d on wind speed U ₁₀ and wave age A, as well as with the dependence of dimensionless roughness height z _n on inverse wave age u*/с р. Calculations carried out for several years in the areas of the Pacific and Indian oceans, based on the ERA-interim and CFSR wind reanalyses have shown good agreement between the model and empirical dependences C _d (U ₁₀) and C _d (A). The range of estimated variability for z _n (u*/с _р ) has been found to be significantly less than empirical. It has been also found that estimated values of wind speed U _10W (t) are overestimated from 5 to 10% in all versions of WAM models compared with the input wind reanalysis U _10R (t) at the moments of appearance maximum values of wind U _10R (t). The reasons for the established features of the WAM model and their dependence on the model version are discussed.     

Modeling the underwater light field fluctuations in coastal oceanic waters: Validation with experimental data

Modeling of the wave-induced underwater light fluctuations at near-surface depths in coastal oceanic waters is challenging because of the surface roughness and strong anisotropic effects of the light field. In the present work, a simple and computationally efficient radiative transfer model is used for the wind-driven sea surface for simulating underwater light fields such as downwelling irradiance (E_d), upwelling irradiance (E_u), and upwelling radiance (L_u) in a spatial domain. It is an extension of our previous work that essentially combines the air–sea interface of the wind-driven sea surface with transmittance and reflectance along with the diffuse and direct components of the homogenous and inhomogeneous water column. The present model simulates underwater light fields for any possible values of absorption and backscattering coefficients. To assess the performance of the model, the E_d, E_u, and L_u profiles predicted by the model are compared with experimental data from relatively clear and turbid coastal waters. Statistical results show significantly low mean relative differences regardless of the wavelength. Comparison of the simulated and in-situ time series data measured over rough sea surfaces demonstrates that model-observation agreement is good for the present model. The Hydrolight model when implemented with the modified bottom reflectance and phase function provides significantly better results than the original Hydrolight model without consideration of the bottom slope and vertically varying phase function. However, these results are non-spatial and have errors fluctuating at different wavelengths. To further demonstrate the efficiency of the present model, spatial distribution patterns of the underwater light fields are simulated based on the measured data from a coastal station for different solar zenith angles (under sunny condition). Simulated wave-induced fluctuations of the underwater lights fields show a good consistency with in-situ data for a few near-surface depths. The present model also provides a reasonable approximation for simulating wave-induced effects on the downward irradiance field and its anisotropic conditions caused by the surface roughness, wavelength and angle of incidence.     

南极西北威德尔海上层海洋湍流混合研究

Turbulent mixing in the upper ocean(30-200 m) of the northwestern Weddell Sea is investigated based on profiles of temperature,salinity and microstructure data obtained during February 2014.Vertical thermohaline structures are distinct due to geographic features and sea ice distribution,resulting in that turbulent dissipation rates(ε) and turbulent diffusivity(K) are vertically and spatially non-uniform.On the shelf north of Antarctic Peninsula and Philip Ridge,with a relatively homogeneous vertical structure of temperature and salinity through the entire water column in the upper 200 m,both ε and K show significantly enhanced values in the order of O(10~(-7))-O(10~(-6)) W/kg and O(10~(-3))-O(10~(-2)) m~2/s respectively,about two or three orders of magnitude higher than those in the open ocean.Mixing intensities tend to be mild due to strong stratification in the Powell Basin and South Orkney Plateau,where s decreases with depth from O(10~(-8)) to O(10~(-9)) W/kg,while K changes vertically in an inverse direction relative to s from O(10~(-6)) to O(10~(-5)) m~2/s.In the marginal ice zone,K is vertically stable with the order of10~(-4) m~2/s although both intense dissipation and strong stratification occur at depth of 50-100 m below a cold freshened mixed layer.Though previous studies indentify wind work and tides as the primary energy sources for turbulent mixing in coastal regions,our results indicate weak relationship between K and wind stress or tidal kinetic energy.Instead,intensified mixing occurs with large bottom roughness,demonstrating that only when internal waves generated by wind and tide impinge on steep topography can the energy dissipate to support mixing.In addition,geostrophic current flowing out of the Weddell Sea through the gap west of Philip Passage is another energy source contributing to the local intense mixing.     

Variability of Particulate Organic Phosphorus in the Northwestern Part of the Black Sea

Based on long-term (1985–1995) monitoring data, the paper considers the peculiarities of seasonal variability in the spatial and vertical distribution of particulate organic phosphorus (Р_POM) in the surface layer and in the photosynthetic zone in the northwestern Black Sea. Regression equations, experimental data, and satellite observations for the chlorophyll a concentration allowed us to evaluate the seasonal longterm (1979–1995) variability in Р_POM in the surface layer and photosynthesis zone. The ratios of the concentrations of particulate organic carbon, nitrogen, phosphorus, and chlorophyll a are calculated and statistical estimates of seasonal changes in the Р_POM in the areas with different degrees of influence of river runoff and water of open seas are obtained. The consistency of intra-annual changes in the concentrations of Р_POM, chlorophyll a, and phytoplankton biomass is shown, which indicates the role of phytoplankton in the formation of Р_POM and in its intra- and interannual variability in the northwestern part of the sea. It is shown that long-term seasonal variations in Р_POM and related changes in the concentration of chlorophyll a depend on the variability of bulk river runoff, the extent of its abundance in the northwestern shelf, and regional hydrometeorological conditions.     

Seasonal dynamics and estimation of the annual primary production of phytoplankton in the Black Sea

Seasonal variations in the surface chlorophyll a concentrations (Chl _s) and the integrated primary production (PP _inf) were investigated for ten regions of the Black Sea based on long term observations (1973–1997). Two or three maximums of both Chl _s and PP _inf were registered in most of the shelf regions (SR, <200 m), the continental slope (CS, 200–1500 m), and the deep regions (DSR >1500 m) in February–March, June–August, and October–November. Such a pattern suggests that the seasonal dynamics of PP _inf strongly depend on the Chl _s variability. The mean annual values of the PP _inf comprised 130–420, 130–150, and 140–150 g C m^?2 in the SR, CS, and DSR, respectively. These values are mainly typical of the eutrophic layer and the transition between the eutrophic and mesotrophic waters (SR) or for the upper boundary of the mesotrophic waters (CS and DSR). The maximal contribution of the wintertime (December–March) to the total PP _inf values (40–42%) was observed in the DSR. In the SR and the adjacent eastern CS areas, the proportion of the PP _inf summertime production (June–September) reaches 40–60% and is higher than the wintertime production. The lowest values of PP _inf (9–17%) were produced in the spring and autumn periods. The total annual values of PP _inf in the Black Sea are close to 50–70 Mt C.     

Seasonal and Interannual Variations of Heat Fluxes in the Barents Sea Region

Seasonal and interannual variations in adjective heat fluxes in the ocean (dQ_oc) and the convergence of advective heat fluxes in the atmosphere (dQ_atm) in the Barents Sea region have been investigated over the period of 1993–2012 using the results of the MIT regional eddy-permitting model and ERA-Interim atmospheric reanalysis. Wavelet analysis and singular spectrum analysis are used to reveal concealed periodicities. Seasonal 2- to 4- and 5- to 8-year cycles are revealed in the dQ_oc and dQ_atm data. It is also found that seasonal variations in dQ_oc are primarily determined by the integrated volume fluxes through the western boundary of the Barents Sea, whereas the 20-year trend is determined by the temperature variation of the transported water. A cross-wavelet analysis of dQ_oc and dQ_atm in the Barents Sea region shows that the seasonal variations in dQ_oc and dQ_atm are nearly in-phase, while their interannual variations are out-of-phase. It is concluded that the basin of the Barents Sea plays an important role in maintaining the feedback mechanism (the Bjerknes compensation) of the ocean–atmosphere system in the Arctic region.     

棉兰老流和棉兰老潜流在季节尺度上的关系研究

The ocean general circulation model for the earth simulator(OFES) products is applied to estimate the transports of the Mindanao Current(MC) and the Mindanao undercurrent(MUC) and explore the relation between them on seasonal scale. In general, the MUC is composed of the lower part of the Southern Pacific Tropical Water(SPTW)and Antarctic Intermediate Water(AAIW). While the deep northward core below 1 500 m is regarded as a portion of MUC. Both salinity and potential density restrictions become more reasonable to estimate the transports of MC/MUC as the properties of water mass having been taken into consideration. The climatological annual mean transport of MC is(37.4±5.81)×10~6 m~3/s while that of MUC is(23.92±6.47)×10~6 m~3/s integrated between 26.5 σ_θ and 27.7 σ_θ, and(17.53±5.45)×10~6 m~3/s integrated between 26.5 σ_θ and 27.5 σ_θ in the OFES. The variations of MC and MUC have good positive correlation with each other on the seasonal scale: The MC is stronger in spring and weaker in fall, which corresponds well with the MUC, and the correlation coefficient of them is 0.67 in the OFES.The same variations are also appeared in hybrid coordinate ocean model(HYCOM) results. Two sensitive experiments based on HYCOM are conducted to explore the relation between MC and MUC. The MUC(26.5σ_θ27.7) is strengthening as the MC increases with the enhancement of zonal wind field. It is shown,however, that the main part of the increasement is the deeper northward high potential density water(HPDW),while the AAIW almost remains stable, SPTW decreases, and vice versa.     

Evaluation of spatial distribution of turbulent mixing in the central Pacific

A long-term mean turbulent mixing in the depth range of 200–1000 m produced by breaking of internal waves across the middle and low latitudes (40°S–40°N) of the Pacific between 160°W and 140°W is examined by applying fine-scale parameterization depending on strain variance to 8-year (2005–2012) Argo float data. Results show that elevated turbulent dissipation rate (ε) is related to significant topographic regions, along the equator, and on the northern side of 20°N spanning to 24°N throughout the depth range. Two patterns of latitudinal variations of ε and the corresponding diffusivity (K_ρ) for different depth ranges are confirmed: One is for 200–450 m with significant larger ε and K_ρ, and the maximum values are obtained between 4°N and 6°N, where eddy kinetic energy also reaches its maximum; The other is for 350–1000 m with smaller ε and K_ρ, and the maximum values are obtained near the equator, and between 18°S and 12°S in the southern hemisphere, 20°N and 22°N in the northern hemisphere. Most elevated turbulent dissipation in the depth range of 350–1000 m relates to rough bottom roughness (correlation coefficient?=?0.63), excluding the equatorial area. In the temporal mean field, energy flux from surface wind stress to inertial motions is not significant enough to account for the relatively intensified turbulent mixing in the upper layer.     

Modeling of ocean dynamics with large variations in sea level

The formulation and the algorithm of solving an ocean model for the prediction and assimilation of the observed data which makes it possible to reconstruct the circulation in the deep-water parts of the sea and at a shallow water shelf, as well as to describe the large time–space variability in the surface level, are considered. The model uses a vertical hybrid σ–z coordinate system: the several upper tens of meters of the ocean are described in the σ-coordinate system and the rest of the water column is described in the z coordinates. Such hybridization extends the possibilities of models for reconstructing thermo-hydrodynamic processes in different sea basins and the World Ocean. The differential formulation of the model in the σ–z coordinate system is presented; the simplified records of several operators that are allowable in the case of a small thickness of the ocean σ-layer are described. The construction of a computational grid, approximation of the bottom topography on it, and discretization of equations and boundary conditions of the models are considered; an approach to describing the bottom friction at shallow waters is offered. The results of the comparative experiments in the z and σ–z coordinate models are analyzed.     

Mesoscale wind stress–SST coupling in the Kuroshio extension and its effect on the ocean

Mesoscale perturbations (with a size of 100–1000 km) of wind stress magnitude, divergence and curl in the Kuroshio Extension (KE) are observed to tightly link to those of sea surface temperature (SST), and downwind and crosswind SST gradients, respectively. Based on long-term satellite observational data, their empirical relationships are established, which are further used to represent mesoscale wind stress–SST coupling in an ocean model that is based on the Regional Oceanic Modelling Systems (ROMS). The strength of mesoscale perturbations of wind stress and SST is observed to display a consistent seasonal variability, with the maximum appeared in winter while the minimum appeared in summer. This seasonal variability characteristic is also successfully simulated by ROMS with high resolution. Through comparing two experiments with and without the mesoscale wind stress–SST coupling, it is found that the mesoscale wind stress perturbation (τ _MS) has a negative feedback on SST perturbation (SST_MS). Analyses of sensitivity experiments suggest that the τ _MS acts to inhibit SST_MS mainly by means of surface heat flux. The τ _MS –SST_MS coupling also exerts influences on the ocean mean state and seasonal variability of SST in the KE. The effect of τ _MS on the SST is distinct in autumn and winter when the mesoscale perturbations are most active. Analyses of sensitivity experiments demonstrate that the τ _MS can affect the long term mean SST through either way of surface heat flux or momentum flux.     

东海春季赤潮过程中优势种的演替:生物-物理耦合数值模拟

In the East China Sea(ECS), the succession of causative species responsible for blooms is a recurrent phenomenon during the spring, which changes from diatoms to dinoflagellates. Observations from space and in situ cruises captured this pattern of succession during spring of 2005. In this study, we coupled two biological models, which were developed previously for Skeletonema costatum and Prorocentrum donghaiense,into a circulation model tailored for the ECS. The coupled biophysical model was used to hindcast the blooms and to test the hypothesis proposed in earlier studies that phosphate(PO4 3–) is the first-order decider of the succession. The coupled model successfully reproduced the hydrodynamics(as described in a companion paper by Sun et al.(1), the spatiotemporal distribution of the chlorophyll a(Chl a) concentration, and the species succession reasonably well. By analyzing the effects of different factors on the surface Chl a distribution, we confirmed that the offshore boundaries of the blooms were confined by PO4 3–. In addition, we suggest that surface wind fields may modulate the horizontal distribution of blooms. Thus, during the dispersal of blooms, surface winds coupled with PO4 3– may control the succession of blooms in the ECS. The proposed coupled model provides a benchmark to facilitate future improvements by including more size classes for organisms, multiple nutrient schemes, and additional processes.     

Assimilation of sea surface temperature in the Northwest Pacific Ocean and its marginal seas using the ensemble Kalman filter

Satellite-borne sea surface temperature (SST) data were assimilated with the ensemble Kalman filter (EnKF) in a Northwest Pacific Ocean circulation model to examine the effect of data assimilation. The model domain included the northwestern part of the Pacific Ocean and its marginal seas, such as the Yellow Sea and East/Japan Sea. The performance of the data assimilation was evaluated by comparing the simulated ocean state with that observed. Spatially averaged root-mean-squared errors in the SST and sea surface height (SSH) decreased by 0.44 °C and 4 cm, respectively, by the assimilation. The results of the numerical experiments substantiated the effectiveness of the SST assimilation via the EnKF for all marginal seas, as well as the Kuroshio region. The benefit of the data assimilation depended on the characteristics of each marginal sea. The variation of the SST in the East/Japan Sea and the Kuroshio extension (KE) region were improved 34% and those in the Yellow Sea 12.5%. The variation of the SSH was improved approximately 36% in the KE region. This large improvement was achieved in the deep-water regions because assimilation of SST data corrected the separation point of the western boundary currents, such as the Kuroshio and the East Korea Warm Current, and the associated horizontal surface currents. The SST assimilation via the EnKF also improved the subsurface temperature profiles. The effectiveness of SST assimilation was seasonally dependent, with the improvement being relatively larger in winter than in summer, which was related to the seasonal variation of the vertical mixing and stratification in the ocean surface layer.     

Features of the spatial distribution of phytoplankton in Nhatrang Bay of the South China Sea during the rainy season

The species composition, phytoplankton abundance, and relative yield of the variable fluorescence (F _v /F _m ) were determined in the mesotrophic Nhatrang Bay in October–November of 2004. The species diversity (250 taxonomic units) and heterogeneity of the phytoplankton structure were high. With respect to the number of species and their abundance, diatoms prevailed. In selected parts of the bay, dinoflagellates dominated. The mean biomass in the water column under 1 m² (B _t ) varied from 2.3 to 64.4 mg C/m³ being 31.0 mg C/m³ on average. The values of B _t were the lowest at the stations nearest to the river mouth. Seaward, B _t increased. The values of B _t increased with depth at some stations and decreased at others. In the surface sea layers, the biomass was lower than that in the underlying waters. The values of F _v /F _m ranged from 0.10 to 0.64 (at a mean value of 0.49). The lowest values of F _v /F _m were observed in the area close to the seaport. Over the greater part of the bay, the values of F _v /F _m were higher than 0.47. Such values are indicative of the relatively high potential photosynthetic activity of the phytoplankton. The abundance and species diversity were higher than those in the dry season (March–April).     

台湾海峡铀的来源及保守性混合

Seawater samples are collected in the spring of 2013 from the Taiwan Strait for the analysis of uranium(U)concentrations and isotopic compositions using MC-ICP-MS, and the geochemical behavior patterns of U in the Taiwan Strait are then investigated. Average concentrations of individual U isotopes are(3.23±0.14) μg/kg for 238 U,(2.34±0.09)×10~(–2) μg/kg for ~(235)U and(2.05±0.07)×10~(–4) μg/kg for 234 U. Correspondingly, the U isotopic compositions are 155±18 for δ234U and 138±2 for 238U:235U. The U concentrations and isotopic ratios in the Taiwan Strait are similar to those of open ocean seawater, suggesting the dominance of the open ocean input to the strait's U pool.However, river input, as suggested by the slightly lower salinity than that of the open ocean, also affected the U concentrations and isotopic compositions in the strait. From a compilation of U concentrations in the Taiwan Strait and adjacent areas, including the Jiulong Estuary and Zhujiang Estuary, the Xiamen Bay and the northern South China Sea, a strong and significant relationship between U concentration and salinity [U:S; U=(0.093 4±0.002 4)S+(0.092 0±0.061 5)] is revealed, suggesting conservative mixing of U in the Taiwan Strait. To better understand the U geochemistry in the Taiwan Strait, a multiple endmembers mixing model is applied to estimate the contributions of potential sources. The open ocean seawater contributed 69%–95% of U in the Taiwan Strait, with river water approximately 2%, and dust deposition only around 0.13%. Therefore, the model results supported the open ocean input source and the conservative mixing behavior of U derived from the observation of U concentrations and isotopic ratios and U:S ratios. The sediment interstitial water may be an important source of U to the Taiwan Strait with a possible contribution of 3%–29%, consistent with previous investigations based on radium isotopes.However, further investigations are warranted to examine the U concentration in the sediment interstitial water and its input to the overlying seawater in the Taiwan Strait.