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.     

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.     

Age determination and growth rate of <Emphasis Type="Italic">Mactra chinensis</Emphasis> (Bivalvia: Mactridae) by external rings and chondrophore growth bands

Age, growth and mortality of Mactra chinensis were investigated during the period from October 2012 to September 2013 in Busan, South Korea. The monthly variation of the marginal index (MI) of the shell and chondrophore showed that the ring of this species was formed once a year during July. We estimated the age of M. chinensis by reading the external rings on the shell and the growth bands of the chondrophore to compare growth parameters between the two growth characters. The age of this species ranged from 0 to 8 years (shell-based age reading) and from 0 to 10 years (chondrophore-based age reading). Based on external rings and growth bands of chondrophore for the same period, the von Bertalanffy growth functions were expressed by the equation, L _t = 101.53[1-exp {-0.15(t + 0.75)}] and L _t = 90.03[1-exp {-0.20(t + 0.50)}], respectively. The likelihood test showed that there was a significant difference in L _∞ (P < 0.001), K (P < 0.001), to (P < 0.001) estimated from non-linear regression between the two growth characters.     

Sources and distribution of allochthonous organic matter in surface sediment from the Seomjin River to the southern inner shelf of Korea

The spatial distributions of δ¹³C, δ¹⁵N, and n-alkanes were investigated to determine the source and transportation of allochthonous organic matter from the mouth of the Seomjin River to the southern inner shelf break of Korea. Total organic carbon (%) ranged from 0.3% to 1.6% (average = 0.80%, n = 81), and the C/N ratio varied from 2.4 to 12.4 (average = 6.76, n = 81). The δ¹³C values ranged from ?25.86 to ?20.26‰ (average = ?21.47‰, n = 81), and δ¹⁵N values ranged from 4.37‰ to 8.57‰ (average = 6.72‰, n = 81). The contribution of the terrestrial fraction of organic matter to the total ranged from 4.4% to 97.7% (average = 24.4%, n = 81), suggesting higher amounts around the catchment area and lower amounts in the offshore area. The concentration of total n-alkanes (nC₂₅ ? nC₃₅) was higher at the boundary between the outer bay and inner shelf break (BOBIS). Average chain length and the carbon preference index both indicated that major leaf wax n-alkanes accounted for the observed distribution of terrestrial organic matter, and were dominant in the inner shelf break (around BOBIS) and outer shelf break. Based on the spatial distribution of the total n-alkanes and the sum of nC₂₇, nC₂₉, and nC₃₁, the terrestrial organic matter distribution was considered to be controlled by local oceanographic conditions, especially at the center of the BOBIS. In addition to enabling the distribution and source of terrestrial organic matter to be identified, the n-alkanes indicated that minor anthropogenic allochthonous organic materials were superimposed on the total organic materials in the central part of Yeosu Bay and the catchment area. The n-alkane indices revealed weathered petroleum contamination, with contamination levels being relatively low at the present time.     

2008年夏季加拿大海盆净二氧化碳吸收量评估

The third Chinese National Arctic Research Expedition(CHINARE) was conducted in the summer of 2008.During the survey,the surface seawater partial pressure of CO_2(pCO_2) was measured,and sea water samples were collected for CO_2 measurement in the Canada Basin.The distribution of pCO_2 in the Canada Basin was determined,the influencing factors were addressed,and the air-sea CO_2 flux in the Canada Basin was evaluated.The Canada Basin was divided into three regions:the ice-free zone(south of 77°N),the partially ice-covered zone(77°–80°N),and the heavily ice-covered zone(north of 80°N).In the ice-free zone,pCO_2 was high(320 to 368μatm,1 μatm=0.101 325 Pa),primarily due to rapid equilibration with atmospheric CO_2 over a short time.In the partially ice-covered zone,the surface pCO_2 was relatively low(250 to 270 μatm) due to ice-edge blooms and icemelt water dilution.In the heavily ice-covered zone,the seawater pCO_2 varied between 270 and 300 μatm due to biological CO_2 removal,the transportation of low pCO_2 water northward,and heavy ice cover.The surface seawater pCO_2 during the survey was undersaturated with respect to the atmosphere in the Canada Basin,and it was a net sink for atmospheric CO_2.The summertime net CO_2 uptake of the ice-free zone,the partially ice-covered zone and the heavily ice-covered zone was(4.14±1.08),(1.79±0.19),and(0.57±0.03) Tg/a(calculated by carbon,1Tg=10~(12) g),respectively.Overall,the net CO_2 sink of the Canada Basin in the summer of 2008 was(6.5±1.3) Tg/a,which accounted for 4%–10% of the Arctic Ocean CO_2 sink.     

“海参-对虾”混养海水池塘的食物网结构及营养层级研究

The food sources of aquacultured Apostichopus japonicus and the trophic levels of organisms in a sea cucumber(A. japonicus) and prawn(Penaeus japonica) polyculture system in a saltwater pond in Zhuanghe, Liaoning Province were examined using carbon and nitrogen stable isotopes. Across organisms, δ13C ranged from(–25.47±0.20)‰ to(–16.48±0.17)‰(mean±SD), and δ15N ranged from(4.23±0.49)‰ to(12.44±0.09)‰. The δ13C and δ15N contents of A. japonicus, P. japonica and Fenneropenaeus chinensis were comparatively higher than those of other organisms. Values of δ13C and δ15N revealed that P. japonica, Hemigrapsus sanguineus and Neomysis japonica comprised the largest component of the diet of A. japonicus. The mean trophic level of the organisms in this saltwater pond polyculture system was(2.75±0.08). P. japonica, A. japonicus, F. chinensis,Synechogobius hasta and Neomysis japonica were in the 3rd trophic level(2–3); jellyfish, H. sanguineus and zooplankton were in the 2nd trophic level(1–2); and Enteromorpha prolifera, benthic microalgae, periphyton and suspended matter primarily consisting of phytoplankton, bacteria and humus were in the primary trophic level(0–1).     

斜向波浪作用下双层水平板式防波堤波浪荷载试验研究

In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate width B/L,wave height Hs/D and incident angle θ0 on the wave forces were analyzed and discussed. The results showed that:(1) The envelopes of the total wave pressure were generally symmetrical along the direction of plate width under the incident angles(θ0) being 0°, 15°, 30°, 45° and 60°. In particular, the envelopes of wave pressure atθ0=30° were larger than all other cases.(2) The synchronous pressure distribution of the breakwater under oblique wave action was more complicated comparing to the normal incident waves.(3) Based on data analysis, an empirical formula was obtained to estimate the total vertical force of the twin-plate breakwater.This empirical formula can be a good reference for the design basis of engineering applications under specified wave conditions.     

High mortality and poor growth of green mussels, <Emphasis Type="Italic">Perna viridis</Emphasis>, in high chlorophyll-<Emphasis Type="Italic">a</Emphasis> environment

The current study was carried out from May 2014 to April 2015 to estimate the stock status of P. viridis in Marudu Bay. The gonad development was monitored by histological examination, while the population parameters including asymptotic length (L_∞), growth coefficient (K), mortality rate (Z, F and M), exploitation level (E) and recruitment of P. viridis were estimated using the lengthfrequency data. Results of the current study demonstrated that P. viridis in Marudu Bay spawned throughout the year with two major peaks, one in April to May and another one in October to December. The recruitment pattern was continuous with the peak in May to June 2014, which corresponded to the first spawning peak in April. However, no significant recruitment was observed from the second spawning peak due to the difference in spawning timing between male and female populations. The estimated asymptotic length (L_∞), growth coefficient (K), total mortality (Z), natural mortality (M), fishing mortality (F) and growth performance (ф) of P. viridis in Marudu Bay were estimate to be 117 mm, 0.97 yr^-1, 4.39 yr^-1, 1.23 yr^-1, 3.16 yr^-1 and 4.123, respectively. The exponent b of the lengthweight relationship was 2.4 and exploitation level (E) was 0.72. The high mortality, low condition indices and negative allometric of P. viridis in Marudu Bay is caused by a lack of suitable food in the surrounding water.     

Age and growth of filefish, <Emphasis Type="Italic">Thamnaconus modestus</Emphasis> (Günther, 1877) off the Jeju Island of Korea

The age and growth of filefish, Thamnaconus modestus (Günther 1877) in the southern waters of Korea were investigated. Samples were collected with commercial trawl catches during the period from May 2009 to December 2011. Of the 2,626 specimens collected, the sex ratio was not significantly different from 1:1 (P > 0.05). The total length ranged from 11.3 to 42.1 cm. The gonadosomatic index for both sexes was the highest in May to June, indicating that May to June is the main spawning period. The length of females at sexual maturity was 25.92 cm. The length-weight relationship of the filefish was TW = 0.0121TL^3.0536 (n = 1,692, r² = 0.9034, P < 0.001). The age of the sampled individuals was estimated by counting growth rings recorded on the 5^th vertebrae; ages ranged from 0 to 9 years. The filefish of the same age displayed a high individual variation in total length. Length-at-age data were fitted by using the Von Bertalanffy growth model. The estimated Von Bertalanffy growth parameters were L_∞ = 42.04 cm, k = 0.21 year^?1 and t₀ = ?1.56 for females, L_∞ = 41.20 cm, k = 0.18 year^?1 and t₀ = ?2.36 for males, and L_∞ = 43.16 cm, k = 0.17 year^?1 and t₀ = ?2.18 for the combination of both male and female. These data can be used as useful biological information for the future fishery management of filefish resources in Korean waters.     

C波段紧缩极化SAR海冰探测能力评估

The C-band synthetic aperture radar(SAR) data from the Bohai Sea of China, the Labrador Sea in the Arctic and the Weddell Sea in the Antarctic are used to analyze and discuss the sea ice full polarimetric information reconstruction ability under compact polarimetric modes. The type of compact polarimetric mode which has the highest reconstructed accuracy is analyzed, along with the performance impact of the reconstructed pseudo quad-pol SAR data on the sea ice detection and sea ice classification. According to the assessment and analysis, it is recommended to adopt the CTLR mode for reconstructing the polarimetric parameters σ_(HH)~0,σ_(VV)~0, H and α,while for reconstructing the polarimetric parameters σ_(HV)~0, ρ_(H-V), λ_1 and λ_2, it is recommended to use the π/4 mode.Moreover, it is recommended to use the π/4 mode in studying the action effects between the electromagnetic waves and sea ice, but it is recommended to use the CTLR mode for studying the sea ice classification.     

Blooms of <Emphasis Type="Italic">Trichodesmium erythraeum</Emphasis> in the South Eastern Arabian Sea during the onset of 2009 summer monsoon

This study presents in situ evidence for the blooms of Trichodesmium erythraeum observed in the shelf waters of the South Eastern Arabian Sea (SEAS) during the onset of the southwest monsoon in June 2009. Evidence showed that water surface discoloration was caused by the accumulation of T. erythraeum, and that the water column contained a colony of T. thiebautii. The surface water color in the bloom region varied from pale brown to pinkish red. Pale brown indicated healthy algae at the peak of its photosynthetic activity, while pinkish red indicated the presence of photosynthetically less active filaments. Zooplankton abundance, especially copepodites, in the bloom area substantiated the theory that Trichodesmium filaments are excellent epiphytes to which the copepodites cling. The bloom area was very fertile with copious quantities of dissolved oxygen (6.85 ml L^?1), PO₄-P (0.108 μmol L^?1) and SiO₄ (1.29 μmol L^?1). Lower NO₃-N (0.028 μmol L^?1) values in the bloom area did not appear to affect Trichodesmium growth from molecular nitrogen fixation. However, lower NO₃-N values altered the normal phytoplankton composition of this area.     

Temporal variations in the current structure and volume transport of the Coastal Oyashio revealed by direct current measurement

Direct current measurements by a shipboard and bottom-mounted acoustic Doppler current profiler and concurrent hydrographic observations with a CTD were conducted off southeastern Hokkaido, Japan, between January and May 2005 to reveal temporal variations in the current structure and volume transport of the Coastal Oyashio (CO). The CO, which has a baroclinic jet structure with southwestward speeds exceeding 90 cm s^?1 and a width of 7–8 km, was associated with a surface-to-bottom density front and was formed on the offshore side of the shelf break. The volume transport of CO (T _CO) was estimated by integrating the fluxes of lower-density water that was trapped against the coast along the density front represented by the 26.2 σ _θ isopycnal line. This transport decreased monotonously from 0.79 Sv (1 Sv = 10⁶ m³ s^?1) in January to 0.21 Sv in March and subsequently to 0.12 Sv in May, possibly due to the decay of the East Sakhalin Current Water in the Okhotsk Sea. Accompanied by a decrease in T _CO, the location of the jet structure associated with the density front moved toward the coast while the maximum speed of the jet decreased and the tilt of the front became more horizontal. Consequently, more saline offshore Oyashio water flowed into the deep part of the shelf area, and the current structure altered from relatively barotropic in winter to baroclinic in spring. This study is the first to estimate the observed volume transport of the CO from direct current measurements.     

Vertical turbulent exchange in the Black Sea pycnocline and its relation to water dynamics

Estimates of vertical turbulent diffusion coefficient (K _t ) in the Black Sea pycnohalocline have been obtained from data of simultaneous observations of seawater temperature, salinity, density, and horizontal current velocity, obtained in the northeastern part of the Black Sea during 2013–2014 with a moored Aqualog profiler. A Munk and Andersson (1948) type parameterization, adapted for the Black Sea environment, is proposed for calculating K _t . Strong short-period (several days) variability of turbulent exchange is revealed, induced by vertical shear variations of the current velocity.     

Retrieval of total ozone in the mesosphere with a new model of electronic-vibrational kinetics of O<Subscript>3</Subscript> and O<Subscript>2</Subscript> photolysis products

The paper presents a new model of electronic-vibrational kinetics of the products of ozone and molecular oxygen photodissociation in the terrestrial middle atmosphere. The model includes 45 excited states of the oxygen molecules O₂(b ¹, Σ _g ⁺ ,v= 0−2), O₂ (a ¹Δ _g , v= 0−5), and O₂(X ³Σ _g ⁻ , v= 1−35) and of the metastable atom O (¹ D) and over 100 aeronomic reactions. The model takes into account the dependence of quantum yields of the production of O₂(a ¹Δ _g , v= 0−5) in a singlet channel of ozone photolysis in the Hartley band on the wavelength of photolytic emission. Taking account of the electronic-vibrational kinetics is important in retrieval of the vertical profiles of ozone concentration from measured intensities of the Atm and IR Atm emissions of the oxygen bands above 65 km and leads to an increase in the ozone concentration retrieved from the 1.27-μm emission, in contrast to the previous model of pure electronic kinetics. Sensitivity analysis of the new model is made for variations in the concentrations of atmospheric constituents ([O₂], [N₂], [O(³P)], [O₃], [CO₂]), the gas temperature, rate constants of the reactions, and quantum yields of the reaction products. A group of reactions that most strongly affect the uncertainty of ozone retrieval from measured intensities of atmospheric emissions of molecular oxygen O₂(b ¹Σ _g ⁺ , v) and O₂(a ¹Δ _g , v) has been determined. Original Russian Text ? V.A. Yankovsky, V.A. Kuleshov, R.O. Manuilova, A.O. Semenov, 2007, published in Izvestiya AN. Fizika Atmosfery i Okeana, 2007, Vol. 43, No. 4, pp. 557–569.     

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.     

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.     

Influence of Feedbacks in the Climate–Energetics System on the Intensity of an Urban Heat Island

Both horizontal and vertical heat exchanges and feedbacks between air temperature and anthropogenic heat fluxes significantly affect the characteristics of the urban heat island (UHI). The UHI intensity depends, in particular, on the ratio between the scales L_A (area of anthropogenic forcing) and L_γ (distance passed by an air particle of the oncoming stably stratified flow before its temperature approaches air temperature within the UHI). Both advection and feedback effects may be estimated based on the equation for the local heat balance of the underlying surface. In this case, heat advection is taken into account by calculating temperatures individually for the atmospheric boundary layer and the surface of the urban canopy layer. The estimates show that the asymptotics of strong advection is more characteristic of a typical city. However, under weak winds, with consideration for the feedback between air temperature and anthropogenic heat flux, some deviations from this asymptotics are probable.     

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

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.     

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.