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

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.     

Drag and Sediment Dispersion Over Sand Waves

Estimates of the drag coefficient over sand waves during calm weather in the southern North Sea have been obtained from measurements of the water slope and currents at different heights (z) above the sea-bed using the log profile and momentum balance methods. An observed phase difference between principal terms in the momentum balance equation is examined theoretically. Drag coefficient estimates are found to agree broadly with previous studies. Owing to bedform asymmetry, average drag coefficient values obtained atz=1 m (C₁₀₀) are found to be 0·0021 and 0·0029 for flood and ebb tides, respectively. Systematic changes in bed roughness are not detected. Using a momentum balance approach, the average drag coefficient value (C_d) atz=10 m is found to be 0·0056. Changes in 10-min averageC_dvalues over sand waves during the tidal cycle are found to be small with bedform asymmetry having no detectable effect. Correlation betweenC_dandC₁₀₀is found to be poor and separation of skin friction and form drag terms is not possible with existing measurements. The inclusion of form drag inC₁₀₀values at the present site leads to over-estimation of the bed shear stress ({q) available to mobilize and transport sediment. Mobile sediment, detected through the use of tracers and a transmissometer, was not found to have any measurable effect on eitherC_dorC₁₀₀in calm weather conditions.     

The roughness height and drag law over the water surface based on the hypothesis of local equilibrium

A series of measurements of winds and wind-waves were carried out in wind-wave flumes. A data analysis based on the hypothesis of local equilibrium yielded a new empirical formula on the controversial quantity of roughness heightz ₀ over the water surface: , where the nondimensional roughness height is defined bygz ₀/u _* ² and the wave-wind parameterũ byω _p u _*/g, g being the gravitational acceleration,u _* the friction velocity of air,ω _p the peak frequency of wind-wave spectra. The obtained formula is compared with Charnock's (1955) and Toba's (1979) proposals; is constant in the former and inversely proportional toũ in the latter. As in Toba's, this formula immediately leads to a practically important conclusion that the drag coefficientC _d depends not merely on the usual variableU ₁₀ (wind velocity at 10m height over the water surface), but also on the surface state represented by wind-waves. An explicit expression is provided for the drag coefficient incorporating the wave-wind parameter; it covers the range ofC _d calculated from most of the previous drag formulas, by varying the wave-wind parameter.     

WAVE-CURRENT FORCE ON VERTICAL PILE

The drag and the interia coefficients Cd,CM of wave-current force on vertical pile are well related to the redefined Keulegan-Carpenter number N' Kc by the test data of the authors. The relation could be also used for irregular wave-current force to calculate in time domain. A simplified method for the calculation of cumulative probabilistic distribution of peak value of irregular wave-current force is also recommended in this paper. These methods were justified by the model test of the authors.     

Ice forces acting on inclined wedges and cones

The increased off shore activity in arctic waters has occasioned investigations with the purpose of reducing the ice forces acting on structures. Korzhavin (Novosibirsk, U.S.S.R.) has shown that the ice forces will be appreciably reduced if wedges or cones are introduced instead of vertical faces. A theory is presented in the following, by which it is possible to calculate the forces based on the physical and geometrical parameters of the system. The theory has been verified by model tests with both artificial and natural ice.For an inclined wedge or cone the force will F = C_FF_max where C_F is a reduction factor and F_max = r_ced is the maximum force acting on a vertical faced structure, r_c being the compression strength, e the thickness of the ice and d the width of the structure. The value of C_F will normally vary from 0·1 to 0·4.     

强风条件下海面拖曳系数和粗糙长度研究

进一步研究强风条件下海-气湍流动量交换以及海浪特征,有助于提高数值天气模式对台风强度演变、移动路径以及恶劣海况的预报能力。依照前人的方法将台风分为风向与浪向(1)相同,(2)相反,和(3)交叉3个扇形区,并结合台风路径数据,得到了浮标数据相对于台风的方位。分别对3种类型的浮标数据进行分析,进而发现了波浪高度和相速度随风速增加而变化的规律。并利用GWW参数化方案计算出摩擦速度(u_*)、拖曳系数(C_DN)和粗糙长度(z₀)。将这些结果与前人代表性的研究论文中所用观测数据和所得研究结论进行比较,结果表明二者有较强的一致性。该研究证明GWW参数化方案在强风条件下依然有很好的适用性。     

Seismic stratigraphy of the quaternary Yellow River delta,Bohai Sea,eastern China

The upper 40 m of stratigraphy of the Yellow River (Huang He) subaqueous delta has been well documented, but the nature of the underlying strata is currently unknown at high-resolution. To address this deficiency we used a Geopulse seismic system to image shallow sedimentary deposits up to 120 m deep on the Yellow River delta. High-resolution seismic reflection images were processed with a series of specific techniques (e.g. swelling attenuation, dynamic s/n filter; f-x deconvolution, predictive deconvolution dipscan stack), and used with borehole data to investigate the Quaternary offshore sequences in the Yellow River (Huang He) delta. Repetitive sequences were observed and interpreted as containing layers of transgressive and regressive deposits. Six seismic transgressive and regressive cycles are identified. Unit M₆F–C₆F correlates with a relative sea-level rise (173–157 ka) and fall (231–173 ka), while Unit M₅F–C₅F is associated with a relative sea-level rise (124–100 ka) and fall (157–124 ka). Unit M₄F–C₄F spans a period of sea-level fall at 100–87 ka, followed by a rise at 87–76 ka. Unit M₃F–C₃F is a transgressive–regressive cycle dated as 76–58 ka. Unit M₂F–C₂F correlates with relative sea level fall at 58.2–36 ka and subsequent rise at 36–22 ka. Unit M₁F–C₁F was deposited during relative sea level fall (22–18 ka), followed by a rise, especially since 8.5 ka.     

Composition research of <Emphasis Type="Italic">n</Emphasis>-alkanes in the samples of hydrothermal deposits of the Mid-Atlantic Ridge by means of gas chromatography-mass spectrometry

The optimum conditions were selected for the chromatographic separation of model mixtures of C₁₂–C₄₀ n-alkanes. For one of the samples of hydrothermal deposits, the extraction conditions of the hydrocarbons were studied and a sample preparation procedure was selected. A procedure is proposed to determine the n-alkanes in samples of hydrothermal deposits by means of gas chromatography-mass spectrometry (GC-MS). The detection limit for the n-alkanes amounted to 3 × 10⁻⁹ to 10⁻⁸% depending on the components. On the basis of the procedure proposed, the composition of the n-alkanes was studied in the samples of hydrothermal deposits collected at the Mid-Atlantic Ridge (the Broken Spur, Lost City, and Rainbow fields). The analysis performed showed that the samples treated contained C₁₄-C₃₅ n-alkanes. The concentrations of the n-alkanes considered were rather low and varied from 0.002 to 0.038 μg/g. Hypotheses concerning the genesis of the n-alkanes identified were formed.     

Model test of the surface and submerged vehicles with the micro-bubble drag reduction

This study is based on the effective experiment observation and measuring technology to discuss the interaction influence between liquid turbulent boundary layer and a crowded group micro-bubbles. It is in order to understand and quantify the micro-bubbles clouds inside the turbulent boundary layer to eliminate the capacity of skin friction drag. Whenever the micro-bubbles are over supplied, pile up effect happened which makes micro-bubbles to integrate to each other as a large-size air film. Although they still have the drag reduction effect, the efficiency of drag reduction slowed down at this transition period. In the experiment of vertical type circulating water tunnel, when 1 μm porous medium is at 7 m/s flow speed, the C_v value at 0.056 has the best drag reduction efficiency of 26%. While 10 μm porous medium is at the same flow speed, the drag reduction efficiency is only around 23%.     

Contrasting chemical and isotopic compositions of organic matter in Changjiang (Yangtze River) estuarine and East China Sea shelf sediments

To examine the source and preservation of organic matter in the shelf sediments of the East China Sea (ECS), we measured bulk C/N and isotopes, organic biomarkers (n-alkanes and fatty acids) and compound-specific (fatty acids) stable carbon isotope ratios in three sediment cores collected from two sites near the Changjiang Estuary and one in the ECS shelf. Contrasting chemical and isotopic compositions of organic matter were observed between the estuarine and shelf sediments. The concentrations of total n-alkanes and fatty acids in the shelf surface sediments (0–2 cm) were 5–10 times higher than those in estuarine surface sediments but they all decreased rapidly to comparable levels below the surface layer. The compositions of n-alkanes in the estuarine sediments were dominated by C₂₆-C₃₃ long-chain n-alkanes with a strong odd-to-even carbon number predominance. In contrast, the composition of n-alkanes in the shelf sediment was dominated by nC₁₅ to nC₂₂ compounds. Long-chain (>C₂₀) fatty acids (terrestrial biomarkers) accounted for a significantly higher fraction in the estuarine sediments compared to that in the shelf sediment, while short-chain (<C₂₀) saturated and unsaturated fatty acids were more abundant in the shelf surface sediments than in the estuarine sediments. Stable carbon isotopic ratios of individual fatty acids showed a general positive shift from estuarine to shelf sediments, consistent with the variations in bulk δ ¹³C_TOCTOC. These contrasts between the estuarine and shelf sediments indicate that terrestrial organic matter was mainly deposited within the Changjiang Estuary and inner shelf of ECS. Post-depositional diagenetic processes in the surface sediments rapidly altered the chemical compositions and control the preservation of organic matter in the region.     

The rates of organic matter destruction in the central part of the Arctic Basin

The intensity of the processes of organic matter (OM) destruction in the deep-water part of the Arctic Basin is considered in terms of the activities of hydrolytic enzymes (alkaline phosphatase and proteases) and of redox enzymes of the electron transport system (ETS). High concentrations of mineral phosphorus are shown to cause low phosphatase activity (0.000–0.005 μM P/l h). The average time of the phosphate recycling amounted to about 240 h. The ranges of the total and specific activity of protease are 0.000–0.192 mg of azocasein/(l h) and 0.00–3.57 mg of azocasein/h mg of protein, respectively. It is found that the microplankton production of exo- and endopeptidase in the Arctic Basin is controlled by the concentration of nitrates in the water. The values of the total and the specific activities of the ETS enzymes vary within 0.27–8.96 μl O₂/(l h) and 8.4–583 μl O₂/mg of protein h, respectively. It is shown that the maximum values of the total activities of the hydrolytic and redox enzymes are characteristic for the photic layer. The vertical fluxes of organic carbon from the photic layer, being calculated by the activity of the ETS enzymes, amounted to 220–600 g of C_org/year m². The presence of ice intensifies the redox processes by factors of 2–15.     

The virtual mass of a rectangular flat plate of finite aspect ratio

The skin friction of a two-dimensional planing flat plate is made up of two opposing components; a drag force from the flow aft of the stagnation line and an opposing thrust force from the jet flow. This paper is concerned only with the drag term and the wake velocity defect which it causes in the water behind the transom.It is concluded that the skin friction is less than would be expected from a flat plate at ambient static pressure (D_fo say) and is approximately equal to D_fo (1 − C_R), where C_R is the normal force coefficient based on wetted area. The wake velocity decrement due to this drag is found to be significant, particularly for surface piercing propellers.     

Currents and turbulent diffusion in the bottom boundary layer of the Barents Sea

The numerical analysis of the stationary field of current velocity on the upper boundary of the bottom boundary layer in the Barents Sea is performed on the basis of a simplified model taking into account the fields of wind velocity and density of water for the principal periods of the seasonal cycle and the bottom topography. The analysis is based on the climatic BarKode database and the data on the wind velocity over the Barents Sea for the last 50 yr. The numerical results demonstrate that the field of bottom currents is fairly nonuniform and the current velocities vary from several fractions of 1 cm/sec to 5 cm/sec in the zones with noticeable slopes of the bottom. The estimates of the thickness of the bottom boundary layer are obtained for the constant coefficient of bottom friction C _f = 0.04. In the major part of the water area of the Barents Sea, the thickness of the bottom boundary layer is close to 1 m. In the regions with significant slopes of the bottom, it increases to 2–2.5 m and, in the two zones of intensification of the bottom currents, becomes as large as 5 m. The maximum estimate of the coefficient of turbulent viscosity is close to 5 cm²/sec. The mean value of the coefficient of vertical density diffusion K _S is equal to 2.34 cm²/sec and its standard deviation is equal to 1.52 cm²/sec. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 4, pp. 31–49, September–October, 2007.     

Reproduction of the large-scale state of water and sea ice in the Arctic Ocean from 1948 to 2002: Part II. The state of ice and snow cover

This paper presents the results of reconstructing the state of ice and snow covers on the Arctic Ocean from 1948 to 2002 obtained with a couplod model of ocean circulation and sea-ice evolution. The area of the North Atlantic and Arctic Ocean north of 65° N, excluding Hudson Bay, is considered. The monthly mean ice areas and extents are analyzed. The trends of these areas are calculated separately for the periods of 1970–1979, 1979–1990, and 1990–2002. A systematic slight underestimation by the model is observed for the ice extent. This error is estimated to fit the error of 100 km in determining the position of the ice edge (i.e., close to the model resolution). In summer the model fails to reproduce many observed polynias: by observational data, the ice concentration in the central part of the Arctic Ocean constitutes around 0.8, while the model yields around 0.99. The average trend for the area of ice propagation in 1960–2002 is 13931 km²/year (or approximately 2% per decade); the trend of the ice area is 17643 km²/year (or approximately 3% per decade). This is almost three times lower than satellite data. The calculated data for ice thickness in the late winter varies from 3.5 to 4.8 m, with a clear indication of periods of thick ice (the 1960s–1970s) and relatively thin ice (the 1980s); 1995 is the starting point for quick ice-area reduction. The maximum ice accumulation is in 1977 and 1988; here, the average trend is negative: −121 km³/year (or approximately 5.5% per decade). In 1996–2002, the average change in the ice thickness constituted +1.7 cm/year. This speaks to the relatively fast disappearance of thin-ice fractions. This model also slightly underestimates the snow mass with a trend of −2.5 km³/year (almost 0.35 mm of snow per year or 0.1 mm of liquid water per year). An analysis of the monthly mean ice-drift velocity indicates the good quality of the model. Data on the average drift velocity and the results of comparisons between the calculated and satellite data for individual months are presented. A comparison with observational data from 1990–1996 in the Fram Strait shows that the model yields 3.28 m for the average ice thickness against the observed value of approximately 3.26 m. For the same period, the model yields a monthly mean transport of 291.29 km³ as compared to the observed value of 237.17 km³. A comparison between the measured and calculated drift velocities in the Fram Strait indicates that the model value is around 9.78 cm/s, which is comparable to the measured value of 10.2 cm/s. The existing problems with describing the ice redistribution by thickness gradations are illustrated by comparing data on ice thickness in the Fram Strait.     

Dependence of sea surface drag coefficient on wind-wave parameters

The relationships between sea surface roughness z 0 and wind-wave parameters are analyzed,and spurious self-correlations are found in all of the parameterization schemes.Sea surface drag coefficient C D is fitted by four wind-wave parameters that are wave age,wave steepness,windsea Reynolds number R B and R H ,and the analyzed data are divided into laboratory,field and combined data sets respectively.Comparison and analysis of dependence of C D on wind-wave parameters show that R B can fit the C D most appropriately.Wave age and wave steepness are not suitable to fit C D with a narrow range data set.When the value of wave age has a board range,R H is not suitable to fit C D either.Three relationships between C D and R B are integrated into the bulk algorithm COARE to calculate the observational friction velocity,and the results show that the relationship between C D and R B which is fitted with field data set can describe the momentum transfer in the open ocean,under low-moderate wind speed condition,most appropriately.     

CFD simulations of the effects of fouling and antifouling

Biofouling is a global problem in the marine industry though its effects on lift and drag are rarely discussed. This paper seeks to employ Computational Fluid Dynamics (CFD) method to quantify the effects of this problem based on the obtained flow field information. The simulation is carried out on a NACA 4424 airfoil and Defence Research Establishment Atlantic (DREA) Submarine Hull. Different levels of fouling are studied with quantified fouling height and density. As an extension, the effects of two common antifouling methods, i.e. tin-free Self-Polishing Copolymer (SPC) and Foul Release on the drag of NACA 4424 foil and the submarine hull is investigated. For NACA 4424 airfoil, fouling reduces lift-drag ratio (C_L/C_D) by up to 80% in maximum and therefore result in the significant increase in fuel consumption. Predicted flow data shows this is related to the increased flow separation region caused by the fouling. It is found that pressure gradient gradually increases from the smallest fouling height to the largest but does not vary that much for fouling of varying densities. The general trend of C_L/C_D varies with angle of attack agrees well with others experimental data. Computed results also show good agreement with experimental data for the DREA bare hull. As to antifouling, Foul Release, despite being 30% more expensive than SPC, exhibits 10-40% higher C_L/C_D as compared to SPC for NACA foil and submarine hull.     

Hydrodynamic coefficients for inclined cylinders

The wave induced dynamic pressures around a circular cylinder of diameter 0.2 m due to regular waves were measured in a wave flume in a water depth of 1 m and in a wave basin in a water depth of 3 m. The experimental investigations were carried out with the cylinder inclined along and against the direction of wave propagation. The least-squares technique was employed to evaluate the coefficients of drag (C_D) and inertia (C_M) from the sectional force time histories obtained by integrating the measured circumferential pressure distribution. The variation of drag and inertia coefficients are presented as a function of Keulegan-Carpenter number (KC) for different inclinations of the cylinder. The comparison between the measured and the theoretical force derived from the evaluated hydrodynamic coefficients is found to be good.     

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

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.     

Studies of hydrocarbons in the waters and snow-ice cover of the southeast sector of the antarctic

Data are presented on the content of hydrocarbons (HC) relative to the concentrations of particulate matter, lipids, C_org, and chlorophyll a in the surface waters and snow-ice cover of the East Antarctic coastal areas. It was shown that the growth of the concentrations of aliphatic HC (AHC) to 30 μg/l in the surface waters takes place in the frontal zones and under the young ice formation. The AHC content in the snow increases with the growth of the aerosol content in the atmosphere. In the lower part of the ice at the boundary with the seawater, despite the low temperatures, the autochthonous processes may provide high AHC concentrations (up to 289 μg/l). Within the snow-ice cover on fast ice, synchronous content variations of all the compounds considered take place.     

Correlation of chlorophyll concentrations and suspended solids with near-surface upward irradiance within LANDSAT bands 4, 5 and 6

Near-surface upward irradiances within each band of the MSS (multi spectral scanner) of LANDSAT, were computed from observed spectral upward irradiances measured just beneath the sea surface. Computation of these irradiances is a useful first step in the analysis of remotely-sensed radiance because the data are not affected by the effects of the atmosphere and sea surface. The computed irradianceE in LANDSAT bands 4, 5 and 6 was normalized by forming ratios,i.e., R ₄₅=E₄/E₅, R₄₆=E₄/E₆, andR ₅₆=E₅/E₆, and then the relationship between these ratios and chlorophylla plus phaeopigment concentration (C) observed simultaneously at 59 stations was examined by linear regression analysis. The observedC ranged over nearly three orders of magnitude, and there is a close to linear relationship between log₁₀-R ₄₅ and log₁₀ C. This suggests that it may be possible to evaluateC from LANDSAT data, if an appropriate atmospheric correction is made. The relationship between suspended solids determined at 35 stations and the same ratios of upward irradiance is also presented and discussed.