Ocean wave parameters retrieved directly from compact polarimetric SAR data

We aim to directly invert wave parameters by using the data of a compact polarimetric synthetic aperture radar(CP SAR) and validate the effectiveness of ocean wave parameter retrieval from the circular transmit/linear receive mode and π/4 compact polarimetric mode. Relevant data from the RADARSAT-2 fully polarimetric SAR on the C-band were used to obtain the compact polarimetric SAR images, and a polarimetric SAR wave retrieval algorithm was used to verify the sea surface wave measurements. Usin...     

北极各海域海冰覆盖范围的变化特征

Sea ice in the Arctic has been reducing rapidly in the past half century due to global warming.This study analyzes the variations of sea ice extent in the entire Arctic Ocean and its sub regions.The results indicate that sea ice extent reduction during 1979–2013 is most significant in summer,following by that in autumn,winter and spring.In years with rich sea ice,sea ice extent anomaly with seasonal cycle removed changes with a period of 4–6 years.The year of 2003–2006 is the ice-rich period with diverse regional difference in this century.In years with poor sea ice,sea ice margin retreats further north in the Arctic.Sea ice in the Fram Strait changes in an opposite way to that in the entire Arctic.Sea ice coverage index in melting-freezing period is an critical indicator for sea ice changes,which shows an coincident change in the Arctic and sub regions.Since 2002,Region C2 in north of the Pacific sector contributes most to sea ice changes in the central Aarctic,followed by C1 and C3.Sea ice changes in different regions show three relationships.The correlation coefficient between sea ice coverage index of the Chukchi Sea and that of the East Siberian Sea is high,suggesting good consistency of ice variation.In the Atlantic sector,sea ice changes are coincided with each other between the Kara Sea and the Barents Sea as a result of warm inflow into the Kara Sea from the Barents Sea.Sea ice changes in the central Arctic are affected by surrounding seas.     

Application of the HY-1 satellite to sea ice monitoring and forecasting

The HY-1A satellite is the first oceanic satellite of China. During the winter of 2002-2003, the data of the HY-1A were applied to the sea ice monitoring and forecasting for the Bohai Sea of China for the fhst time. The sea ice retrieval system of the HY-1A has been constructed. It receives 1B data from the satellite, outputs sea ice images and provides digital products of ice concentration, ice thickness and ice edge, which can be used as important information for sea ice monitoring and the initial fields of the numeric sea ice forecast and as one of the reference data for the sea ice forecasting verification. The sea ice retrieval system of the satellite is described, including its processes, methods and parameters. The retrieving results and their application to the sea ice monitoring and forecasting for the Bohai Sea are also discussed.     

基于高光谱遥感的渤海海冰厚度半经验模型

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness(SEMSIT) that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands(spectral dimension information) are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.     

渤海模型冰运动特征的现场观测研究

Sea ice drift is mainly controlled by ocean currents, local wind, and internal ice stress. Information on sea ice motion, especially in situ synchronous observation of an ice velocity, a current velocity, and a wind speed, is of great significance to identify ice drift characteristics. A sea ice substitute, the so-called "modelled ice", which is made by polypropylene material with a density similar to Bohai Sea ice, is used to complete a free drift experiment in the open sea. The trajectories of isolated modelled ice, currents and wind in the Bohai Sea during non-frozen and frozen periods are obtained. The results show that the currents play a major role while the wind plays a minor role in the free drift of isolated modelled ice when the wind is mild in the Bohai Sea. The modelled ice drift is significantly affected by the ocean current and wind based on the ice–current–wind relationship established by a multiple linear regression. The modelled ice velocity calculated by the multiple linear regression is close to that of the in situ observation, the magnitude of the error between the calculated and observed ice velocities is less than12.05%, and the velocity direction error is less than 6.21°. Thus, the ice velocity can be estimated based on the observed current velocity and wind speed when the in situ observed ice velocity is missing. And the modelled ice of same thickness with a smaller density is more sensitive to the current velocity and the wind speed changes. In addition, the modelled ice drift characteristics are shown to be close to those of the real sea ice, which indicates that the modelled ice can be used as a good substitute of real ice for in situ observation of the free ice drift in the open sea, which helps solve time availability, safety and logistics problems related to in situ observation on real ice.     

基于MODIS热红外数据的渤海海冰厚度反演

Level ice thickness distribution pattern in the Bohai Sea in the winter of 2009–2010 was investigated in this paper using MODIS night-time thermal infrared imagery.The cloud cover in the imagery was masked out manually.Level ice thickness was calculated using MODIS ice surface temperature and an ice surface heat balance equation.Weather forcing data was from the European Centre for Medium-Range Weather Forecasts(ECMWF) analyses.The retrieved ice thickness agreed reasonable well with in situ observations from two off-shore oil platforms.The overall bias and the root mean square error of the MODIS ice thickness are –1.4 cm and 3.9 cm,respectively.The MODIS results under cold conditions(air temperature –10°C) also agree with the estimated ice growth from Lebedev and Zubov models.The MODIS ice thickness is sensitive to the changes of the sea ice and air temperature,in particular when the sea ice is relatively thin.It is less sensitive to the wind speed.Our method is feasible for the Bohai Sea operational ice thickness analyses during cold freezing seasons.     

The estimate of sea ice resources quantity in the Bohai Sea based on NOAA/AVHRR data

The research on sea ice resources is the academic base of sea ice exploitation in the Bohai Sea. According to the ice-water spectrum differences and the correlation between ice thickness and albedo, this paper comes up with a sea ice thickness inversion model based on the NOAA/AVHRR data. And then a sea ice resources quantity (SIQ) time series of Bohai Sea is established from 1987 to 2009. The results indicate that the average error of inversion sea ice thickness is below 30%. The maximum sea ice resources quantity is about 6 × 10 9 m 3 and the minimum is 1.3 × 10 9 m 3 . And a preliminary analysis has been made on the errors of the estimate of sea ice resources quantity (SIQ).     

海洋二号微波散射计极地海冰识别算法研究

An Antarctic sea ice identification algorithm on the HY-2A scatterometer(HSCAT) employs backscattering coefficient(σ0) and active polarization ratio(APR) for a preliminary sea ice identification.Then standard deviation(STD) filtering and space filtering are carried out.Finally,it is used to identify sea ice.A process uses a σ0,STD threshold and an APR as sea ice indicators.The sea ice identification results are verified using the sea ice distribution data of the ASMR2 released by the National Snow and Ice Data Center as a reference.The results show very good consistence of sea ice development trends,seasonal changes,area distribution,and sea ice edge distribution of the sea ice identification results obtained by this algorithm relative to the ASMR2 sea ice results.The accuracy of a sea ice coverage is 90.8% versus the ASMR2 sea ice results.This indicates that this algorithm is reliable.     

溢油污染的渤海海冰反射光谱特征实测研究

Oil spilled on the sea ice surface in the Bohai Sea of China is studied through the field measurements of the reflectance of a simulated sea ice-oil film mixed pixel. The reflection characteristics of sea ice and oil film are also analyzed. It is found that the mixed pixel of sea ice and oil film is a linear mixed pixel. The means of extracting sea ice pixels containing oil film is presented using a double-band ratio oil-film sea-ice index(DROSI) and a normalized difference oil-film sea-ice index(NDOSI) through the analysis of the reflectance curves of the sea iceoil film pixel for different ratios of oil film. The area proportion of the oil film in the sea ice-oil film pixel can be accurately estimated by the average reflectance of the band of 1 610–1 630 nm, and the volume of the spilled oil can be further estimated. The method of the sea ice-oil film pixel extraction and the models to estimate the proportion of oil film area in the sea ice-oil film pixel can be applied to the oil spill monitoring of the ice-covered area in the Bohai Sea using multispectral or hyperspectral remote sensing images in the shortwave infrared band(1 500–1 780 nm).     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

引入拖曳系数参数化的海冰自由漂流模拟研究

Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient(C_a) and oceanic drag coefficient(C_w).Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions.In the present study,the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model,and the wind factor α and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of C_a to C_w are studied to investigate their dependence on the impact factors such as local drag coefficients,floe and ridge geometry.The results reveal that in an idealized steady ocean,C_a/C_w increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone,while it remains at a steady level(0.2-0.25) for large floes in the central ice zone.The wind factor α increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%.And the deflection angle θ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than20%without a steady level like α.The values of these parameters agree well with the previously reported observations in Arctic.The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient(C_s' /C_s) and the ratio of ridge form drag coefficient(C_r'/C_r) to the value of C_a/C_w,α,and θ,because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces.Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.     

Effects of temperature on growth,photophysiology, Rubisco gene expression in <Emphasis Type="Italic">Prorocentrum donghaiense</Emphasis> and <Emphasis Type="Italic">Karenia mikimotoi</Emphasis>

To explore the effects of temperature changes on dinoflagellate bloom succession in the coastal waters of the East China Sea, changes in the growth, photophysiology, and Rubisco gene expression of Prorocentrum donghaiense and Karenia mikimotoi, two harmful algal species, were investigated at different temperatures (16 to 28°C). The maximal specific growth rate and the maximal mRNA expression of Rubisco gene in P. donghaiense and K. mikimotoi occurred at 20 and 24°C, respectively. The photosynthetic activity of P. donghaiense was generally stable, but K. mikimotoi photosynthesis increased when temperatures rose from 16 to 28°C. The effective photochemical efficiency (F _q ^′ /F _m ^′ ) and the maximal relative electron transfer rate (rETR_max) of K. mikimotoi increased significantly with increasing temperature, and the lowest and highest values occurred at 16 and 28°C, respectively. It seems that P. donghaiense has higher photosynthetic capacity than K. mikimotoi due to its higher F _q ^′ /F _m ^′ , rETR_max, and photosynthetic efficiency (α). However, K. mikimotoi has a higher growth rate than P. donghaiense. These results suggest that the photosynthetic activity and genetic responses of dinoflagellates are species-dependent. It is likely that temperature changes affect species composition during blooms, leading to the observed patterns of bloom succession.     

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.     

The Interactive Effects of Elevated CO<Subscript>2</Subscript> and Ammonium Enrichment on the Physiological Performances of <Emphasis Type="Italic">Saccharina japonica</Emphasis> (Laminariales,Phaeophyta)

Environmental challenges such as ocean acidification and eutrophication influence the physiology of kelp species. We investigated their interactive effects on Saccharina japonica (Laminariales, Phaeophyta) under two pH conditions [Low, 7.50; High (control), 8.10] and three NH ₄ ⁺ concentrations (Low, 4; Medium, 60; High, 120 μM). The degree of variation of pH values in the culture medium and inhibition rate of photosynthetic oxygen evolution by acetazolamide were affected by pH treatments. Relative growth rates, carbon, nitrogen, and the C:N ratio in tissue samples were influenced by higher concentrations of NH ₄ ⁺ . Rates of photosynthetic oxygen evolution were enhanced under elevated CO₂ or NH ₄ ⁺ conditions, independently, but these two factors did not show an interactive effect. However, rates of NH ₄ ⁺ uptake were influenced by the interactive effect of increased CO₂ under elevated NH ₄ ⁺ treatment. Although ocean acidification and eutrophication states had an impact on physiological performance, chlorophyll fluorescence was not affected by those conditions. Our results indicated that the physiological reactions by this alga were influenced to some extent by a rise in the levels of CO₂ and NH ₄ ⁺ . Therefore, we expect that the biomass accumulation of S. japonica may well increase under future scenarios of ocean acidification and eutrophication.     

Application of chemical tracers to an estimate of benthic denitrification in the Okhotsk Sea

To estimate benthic denitrification in a marginal sea, we assessed the usefulness of \({\text{N}}_{2}^{*}\) , a new tracer to measure the excess nitrogen gas (N₂) using dissolved N₂ and argon (Ar) with N* in the intermediate layer (26.6–27.4σ _θ ) of the Okhotsk Sea. The examined parameters capable of affecting \({\text{N}}_{2}^{*}\) are denitrification, air injection and rapid cooling. We investigated the relative proportions of these effects on \({\text{N}}_{2}^{*}\) using multiple linear regression analysis. The best model included two examined parameters of denitrification and air injection based on the Akaike information criterion as a measure of the model fit to data. More than 80 % of \({\text{N}}_{2}^{*}\) was derived from the denitrification, followed by air injection. Denitrification over the Okhotsk Sea shelf region was estimated to be 5.6 ± 2.4 μmol kg^?1. The distribution of \({\text{N}}_{2}^{*}\) was correlated with potential temperature (θ) between 26.6 and 27.4σ _θ (r = ?0.55). Therefore, we concluded that \({\text{N}}_{2}^{*}\) and N* can act complementarily as a quasi-conservative tracer of benthic denitrification in the Okhotsk Sea. Our findings suggest that \({\text{N}}_{2}^{*}\) in combination with N* is a useful chemical tracer to estimate benthic denitrification in a marginal sea.     

A new formula for saturated water steam pressure within the temperature range −25 to 220°C

Instead of approximation formula ln(E(t)/E(0)) = [(a ? bt)t/(c + T)] commonly used at present for representing dependence of pressure of saturated streams of liquid water E upon temperature we suggested new approximation formula of greater accuracy in the form ln(E(t)/E(0)) = [(A ? Bt + Ct ²)t/T], where t and T are temperature in °C and K respectively. For this formula with parameters A = 19.846, B = 8.97 × 10^?3, C = 1.248 × 10^?5 and E(0) = 6.1121 GPa with ITS-90 temperature scale and for temperature range from 0°C to 110°C relative difference of approximation applying six parameter formula by W. Wagner and A. Pruß 2002, developed for positive temperatures, is less than 0.005%, that is approximately 15 times less than accuracy obtained with the firs formula. Increase of temperature range results in relative difference increasing, but for even temperature range from 0°C to 220°C it does not higher than 0.1%. For negative temperatures relative difference between our formula and a formula of D. M. Murphy and T. Koop, 2005, is less than 0.1% for temperatures higher than ?25°C. This paper also presents values of coefficients for approximation of Goff and Grach formula recommended by IMO. The procedure of finding dew point T _d for known water steam pressure e _n based on our formula adds up to solving an algebraic equation of a third degree, which coefficients are presented in this paper. For simplifying this procedure this paper also includes approximation ratio applying a coefficient A noted above, in the form T _d (e _n ) = \(\frac{{AT_0 }}{{A - \varepsilon }}\) + 0.0866?² + 0.0116?^10/3, where ? = ln(e _n /E(T ₀)). Error of dew point recovery in this ratio is less than 0.005 K within the range from 0 to 50°C.     

Operating experience with a BM-04 wave-tide gauge

A BM-04 microprocessor wave-tide gauge developed at the Shirshov Institute of Oceanology is used for the acquisition of data on the parameters of surface waves (the mean height \(\bar h\), the mean period \(\bar \tau \), the frequency spectrum S(ω), and so on) and the mean sea level ζ. As a wave gauge, the instrument can be deployed on the sea floor, on piles, and so on at depths down to 10 m in the near-shore zone and can be mounted on subsurface buoys offshore. As a tide gauge, the instrument can be used at depths down to 6000 m. The instruments were successfully tested and intercalibrated in the Black and Baltic seas and the Sea of Okhotsk.     

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

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.     

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

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.