Analysis and comparison of heat flux of landfast ice during 2016 in the Prydz Bay,Antarctica

Long term in situ atmospheric observation of the landfast ice nearby Zhongshan Station in the Prydz Bay was performed from April to November 2016. The in situ observation, including the conventional meteorological elements and turbulent flux, enabled this study to evaluate the sea ice surface energy budget process. Using in situ observations, three different reanalysis datasets from the European Centre for Medium-Range Weather Forecasts Interim Re-analysis(ERA-Interim), National Centers for Environmental Prediction Reanalysis2(NCEP R2), and Japanese 55-year Reanalysis(JRA55), and the Los Alamos sea ice model, CICE, output for surface fluxes were evaluated. The observed sensible heat flux(SH) and net longwave radiation showed seasonal variation with increasing temperature. Air temperature rose from the middle of October as the solar elevation angle increased.The ice surface lost more energy by outgoing longwave radiation as temperature increased, while the shortwave radiation showed obvious increases from the middle of October. The oceanic heat flux demonstrated seasonal variation and decreased with time, where the average values were 21 W/m~2 and 11 W/m~2, before and after August,respectively. The comparisons with in situ observations show that, SH and LE(latent heat flux) of JRA55 dataset had the smallest bias and mean absolute error(MAE), and those of NCEP R2 data show the largest differences.The ERA-Interim dataset had the highest spatial resolution, but performance was modest with bias and MAE between JRA55 and NCEP R2 compare with in situ observation. The CICE results(SH and LE) were consistent with the observed data but did not demonstrate the amplitude of inner seasonal variation. The comparison revealed better shortwave and longwave radiation stimulation based on the ERA-Interim forcing in CICE than the radiation of ERA-Interim. The average sea ice temperature decreased in June and July and increased after September,which was similar to the temperature measured by buoys, with a bias and MAE of 0.9°C and 1.0°C, respectively.     

基于走航观测数据的ERA5南海北部夏季短波辐射收支评估

海表短波辐射收支是海–气界面能量交换的重要物理过程。本研究利用2019年南海北部夏季科考航次的走航观测数据,评估了ERA5再分析数据的海表短波辐射通量收支。结果表明,ERA5的向下短波辐射相比观测偏小,11时和15时(北京时间)的偏差最大,可达-100 W/m²。与此同时,ERA5的海表反照率整体偏低,其中高太阳高度角时段偏差较小,约为-0.03,低太阳高度角时段偏差较大,约为-0.15。向下短波辐射和反照率的偏差共同造成ERA5白天平均海表净短波辐射通量比观测偏小约25.4 W/m²;其中,反照率低估抵消了约50%向下短波辐射偏差的贡献。研究表明,在不同大气透射率情况下,ERA5的海表辐射收支偏差存在不同表现。ERA5海表反照率的低估可能与其采用的参数化方案在南海北部的适用性不足有关。基于观测本研究也给出了一个简单的参数优化方案。     

北极冰海耦合模式对两种不同大气再分析资料响应的分析

本文中我们比较了Climate Forecast System Reanalysis(CFSR)高分辨率的再分析数据集和低分辨率的Japanese 25-year Reanalysis Project(JRA25)再分析数据集在向下短波辐射、向下长波辐射、10m风场、近地面气温、降水、湿度上的不同,发现二者差异最大的为降水数据,其次为向下短波辐射数据、向下长波辐射数据。用这两个数据集驱动同一冰海耦合模式,CFSR强迫的海冰、北冰洋中层水和加拿大海盆温盐结构与实测相比有很大差距,等密度面上的地转流速在加拿大海盆和欧亚海盆比JRA25强迫的结果高20%,同时等密度面的深度偏深、位温偏高,在弗拉姆海峡的流通量也比海洋再分析数据Simple Ocean Data Assimilation(SODA)偏多。CFSR的向下辐射数据更加接近实测,采用此数据的敏感性实验模拟结果与实测符合的更好。对于海冰的模拟,云量起着至关重要的作用,降水带来的淡水通量通过影响大西洋入流水携带的热量进而影响到冰区。此外,CFSR过量的降水也是二者对于北冰洋温盐结构、弗拉姆海峡流通量以及地转流强度模拟产生偏差的主要原因。尽管风场的分辨率不同,在海盆尺度上对于海冰和海水温盐结构的影响并不大。     

北黄海春季海面太阳辐射船基观测与分析

利用"中国近海海洋综合调查与评价"2007年北黄海春季航次的资料,综合分析了短波、长波辐射及海面短波、长波反射辐射特征及与云量、气温、水汽压、海面红外温度等的对应关系。结果显示海面太阳入射辐射远大于海面反射辐射,二者相差约1个量级;海面长波辐射总体上大于大气逆辐射,二者相差约40 W/m2,平均净辐射为153.1 W/m2。观测期间大气透射率约为0.45,并据此分析探讨了云对太阳辐射传输的影响,建立了考虑总云量、低云量及水汽压等因素的海面太阳短波辐射经验估算公式,计算结果与实测结果具有较好一致性,均方根误差为11.3 W/m2。最后讨论了太阳辐射船基观测由于船体晃动可能造成的误差影响。     

Evaluation of Vector Winds Observed by NSCAT in the Seas around Japan

In order to validate wind vectors derived from the NASA Scatterometer (NSCAT), two NSCAT wind products of different spatial resolutions are compared with observations by buoys and research vessels in the seas around Japan. In general, the NSCAT winds agree well with the wind data from the buoys and vessels. It is shown that the root-mean-square (rms) difference between NSCAT-derived wind speeds and the buoy observations is 1.7 ms^–1, which satisfies the mission requirement of accuracy, 2 ms^–1. However, the rms difference of wind directions is slightly larger than the mission requirement, 20°. This result does not agree with those of previous studies on validation of the NSCAT-derived wind vectors using buoy observations, and is considered to be due to differences in the buoy observation systems. It is also shown that there are no significant systematic trends of the NSCAT wind speed and direction depending on the wind speed and incidence angle. Comparison with ship winds shows that the NSCAT wind speeds are lower than those observed by the research vessels by about 0.7 ms^–1 and this bias is twice as large for data observed by moving ships than by stationary ships. This result suggests that the ship winds may be influenced by errors caused by ship's motion, such as pitching and rolling.     

ERA-Interim和NCEP/NCAR再分析资料在我国东南近海适用性分析

利用我国东南近海5个浮标站观测资料,对2012—2016年ERA-Interim和NCEP/NCAR再分析资料10 m风、2 m气温、海平面气压的适用性进行了评估。结果表明:NCEP/NCAR的再分析10 m风适用性更好,ERA-Interim的2 m气温适用性更好,海平面气压两者差异不大。风速再分析值与观测值具有较好的一致性,相关系数达0.8～0.9,但再分析风速总体上有偏小的趋势,平均偏差在-1.3～0 m/s之间,均方根误差在1.5～3 m/s。再分析资料的平均风向有顺时针偏差的趋势,温州浮标偏右达14°以上,均方根误差大多在40°～50°。不管风速还是风向,5个浮标站中均以舟山浮标的再分析值与观测值最为接近;分析还表明,再分析资料的冬季风代表性相对较差,这是造成风速和风向系统性偏差的主要原因。再分析资料与观测2 m气温相关系数均在0.95以上,且有偏高的趋势,NCEP偏高更为明显,有4个浮标站平均偏差达1～2℃,而ERA-I仅1个浮标站偏差1～2℃,4个在1℃以内。春季和冬季气温偏高最为明显,春季升温过程存在异常偏高的可能,秋季气温与观测值最为接近。海平面气压适用性较好,总体优于10 m风和2 m气温,且季节间差异也不大。     

Characteristics of sea ice kinematics from the marginal ice zone to the packed ice zone observed by buoys deployed during the 9th Chinese Arctic Expedition

Sea ice growth and consolidation play a significant role in heat and momentum exchange between the atmosphere and the ocean. However, few in situ observations of sea ice kinematics have been reported owing to difficulties of deployment of buoys in the marginal ice zone (MIZ). To investigate the characteristics of sea ice kinematics from MIZ to packed ice zone (PIZ), eight drifting buoys designed by Taiyuan University of Technology were deployed in the open water at the ice edge of the Canadian Basin. Sea ice near the buoy constantly increased as the buoy drifted, and the kinematics of the buoy changed as the buoy was frozen into the ice. This process can be determined using sea ice concentration, sea skin temperature, and drift speed of buoy together. Sea ice concentration data showed that buoys entered the PIZ in mid-October as the ice grew and consolidated around the buoys, with high amplitude, high frequency buoy motions almost ceasing. Our results confirmed that good correlation coefficient in monthly scale between buoy drift and the wind only happened in the ice zone. The correlation coefficient between buoys and wind was below 0.3 while the buoys were in open water. As buoys entered the ice zone, the buoy speed was normally distributed at wind speeds above 6 m/s. The buoy drifted mainly to the right of the wind within 45° at wind speeds above 8 m/s. During further consolidation of the ice in MIZ, the direct forcing on the ice through winds will be lessened. The correlation coefficient value increased to 0.9 in November, and gradually decreased to 0.7 in April.     

Atmospheric forcing of the southeast Bering Sea Shelf during 1995–99 in the context of a 40-year historical record

The atmospheric forcing of the Bering Sea over its eastern shelf is estimated using the 40-year record of daily data from the NCEP/NCAR Reanalysis. This data set includes estimates of the processes responsible for the atmospheric forcing, namely the surface fluxes of momentum, sensible and latent heat, and longwave and shortwave radiation, and therefore permits quantifying effects that previously could be inferred only from the large-scale nature of the flow. The forcing in 1995–1999 is described in detail using daily time series; historical context for these results is provided with seasonal averages for the years 1959–1999.The analysis for winter concentrates on aspects related to the formation and advection of sea ice. Results indicate that the presence of sea ice is strongly related to the net surface-heat fluxes as well as the cross-shelf component of the wind. The 40-year record lacks any discernible long-term trend in the winter forcing and response. There was a notably cold period in the early to middle 1970s, and a warm period from the late 1970s into the early 1980s, but conditions during the 1990s are similar to those in the late 1950s and 1960s.The analysis for the warm season focuses on the mechanisms responsible for the variability in SST warming. Much of the intraseasonal and interannual variability in this warming can be attributed to variations in the downward shortwave radiation (solar heating). The 40-year record does indicate a long-term trend toward increased solar heating, and reduced surface latent-heat fluxes (evaporative cooling). These changes have led to August SSTs in the 1990s that are roughly 1°C warmer than in the 1960s.     

Sea surface temperature observation by Global Imager (GLI)/ADEOS-II: Algorithm and accuracy of the product

A sea surface temperature (SST) retrieval algorithm for Global Imager (GLI) aboard the ADEOS-II satellite has been developed. The algorithm is used to produce the standard SST product in the Japan Aerospace Exploration Agency (JAXA). The algorithm for cloud screening is formed by combinations of various types of tests to detect cloud-contaminated pixels. The combination is changed according to the solar zenith angle, which enables us to detect clouds even in the sun glitter region in daytime. The parameters in the cloud-detection tests have been tuned using the GLI global observations. SST is calculated by the Multi-Channel SST (MCSST) technique from the detected clear pixels. Using drifting buoy measurements, match-up data are produced to derive the coefficients of the MCSST equations and to examine their performance. The bias and RMSE of the GLI SST are 0.03 K and 0.66 K for daytime and, −0.01 K and 0.70 K for nighttime, respectively.     

The seasonal variability of an air-sea heat flux in the northern South China Sea

The seasonal variabilities of a latent-heat flux (LHF), a sensible-heat flux (SHF) and net surface heat flux are examined in the northern South China Sea (NSCS), including their spatial characteristics, using the in situ data collected by ship from 2006 to 2007. The spatial distribution of LHF in the NSCS is mostly controlled by wind in summer and autumn owing to the lower vertical gradient of air humidity, but is influenced by both wind and near-surface air humidity vertical gradient in spring and winter. The largest area-averaged LHF is in autumn, with the value of 197.25 W/m 2 , followed by that in winter; the third and the forth are in summer and spring, respectively. The net heat flux is positive in spring and summer, so the NSCS absorbs heat; and the solar shortwave radiation plays the most important role in the surface heat budget. In autumn and winter, the net heat flux is negative in most of the observation region, so the NSCS loses heat; and the LHF plays the most important role in the surface heat budget. The net heating is mainly a result of the offsetting between heating due to the shortwave radiation and cooling due to the LHF and the upward (outgoing) long wave radiation, since the role of SHF is negligible. The ratio of the magnitudes of the three terms (shortwave radiation to LHF to long-wave radiation) averaged over the entire year is roughly 3:2:1, and the role of SHF is the smallest.     

Japanese Ocean Flux Data Sets with Use of Remote Sensing Observations (J-OFURO)

We have constructed ocean surface data sets using mainly satellite data and called them Japanese Ocean Flux data sets with Use of Remote sensing Observations (J-OFURO). The data sets include shortwave radiation, longwave radiation, latent heat flux, sensible heat flux, and momentum flux etc. This article introduces J-OFURO and compares it with other global flux data sets such as European Centre for Medium Range Weather Forecasting (ECMWF) and National Center for Environmental Prediction (NCEP) reanalysis data and da Silva et al. (1994). The usual ECMWF data are used for comparison of zonal wind. The comparison is carried out for a meridional profile along the dateline for January and July 1993. Although the overall spatial variation is common for all the products, there is a large difference between them in places. J-OFURO shortwave radiation in July shows larger meridional contrast than other data sets. On the other hand, J-OFURO underestimates longwave radiation flux at low- and mid-latitudes in the Southern Hemisphere. J-OFURO latent heat flux in January overestimates at 10°N–20°N and underestimates at 25°N–40°N. Finally, J-OFURO shows a larger oceanic net heat loss at 10°N–20°N and a smaller loss north of 20°N in January. The data of da Silva et al. in July show small net heat loss around 20°S and large gain around 20°N, while the NCEP reanalysis (NRA) data show the opposite. The da Silva et al. zonal wind speed overestimates at low-latitudes in January, while ECMWF wind data seem to underestimate the easterlies. This revised version was published online in August 2006 with corrections to the Cover Date.     

Optical properties and surface energy flux of spring fast ice in the Arctic

Over the past decades, sea ice in the polar regions has been significantly affecting local and even hemispheric climate through a positive ice albedo feedback mechanism. The role of fast ice, as opposed to drift ice, has not been well-studied due to its relatively small coverage over the earth. In this paper, the optical properties and surface energy balance of land fast ice in spring are studied using in situ observations in Barrow, Alaska. The results show that the albedo of the fast ice varied between 0.57 and 0.85 while the transmittance increased from 1.3×10?3 to 4.1×10?3 during the observation period. Snowfall and air temperature affected the albedo and absorbance of sea ice, but the transmittance had no obvious relationship with precipitation or snow cover. Net solar shortwave radiation contributes to the surface energy balance with a positive 99.2% of the incident flux, with sensible heat flux for the remaining 0.8%. Meanwhile, the ice surface loses energy through the net longwave radiation by 18.7% of the total emission, while the latent heat flux accounts for only 0.1%. Heat conduction is also an important factor in the overall energy budget of sea ice, contributing 81.2% of the energy loss. Results of the radiative transfer model reveal that the spectral transmittance of the fast ice is determined by the thickness of snow and sea ice as well as the amount of inclusions. As major inclusions, the ice biota and particulates have a significant influence on the magnitude and distribution of the spectral transmittance. Based on the radiative transfer model, concentrations of chlorophyll and particulate in the fast ice are estimated at 5.51 mg/m2 and 95.79 g/m2, which are typical values in the spring in Barrow.     

Evaluation of reanalysis surface wind products with quality-assured buoy wind measurements along the north coast of the South China Sea

Three archived reanalysis wind vectors at 10 m height in the wind speed range of 2–15 m/s, namely, the second version of the National Centres for Environmental Prediction(NCEP) Climate Forecast System Reanalysis(CFSv2), European Centre for Medium-Range Weather Forecasting Interim Reanalysis(ERA-I) and NCEPDepartment of Energy(DOE) Reanalysis 2(NCEP-2) products, are evaluated by a comparison with the winds measured by moored buoys in coastal regions of the South China Sea(SCS). The buoy data are first quality controlled by extensive techniques that help eliminate degraded measurements. The evaluation results reveal that the CFSv2 wind vectors are most consistent with the buoy winds(with average biases of 0.01 m/s and 1.76°). The ERA-I winds significantly underestimate the buoy wind speed(with an average bias of –1.57 m/s), while the statistical errors in the NCEP-2 wind direction have the largest magnitude. The diagnosis of the reanalysis wind errors shows the residuals of all three reanalysis wind speeds(reanalysis-buoy) decrease with increasing buoy wind speed, suggesting a narrower wind speed range than that of the observations. Moreover, wind direction errors are examined to depend on the magnitude of the wind speed and the wind speed biases. In general, the evaluation of three reanalysis wind products demonstrates that CFSv2 wind vectors are the closest to the winds along the north coast of the SCS and are sufficiently accurate to be used in numerical models.     

海浪成长状态对TOPEX/Poseidon高度计风速反演的影响

以墨西哥湾同步高度计、浮标资料为例,研究了海浪成长状态对高度计风速反演的影响。同步的高度计风速和浮标风速比较显示,在墨西哥湾地区,海浪成长状态对高度计风速反演有较大影响。在考虑海浪成长状态影响的条件下,利用谱模型反演高度计风速,取得了较好的效果。与目前TOPEX/Poseidon高度计风速反演业务化算法相比,在海浪未充分成长条件下,考虑海浪成长状态影响后,根据谱模型反演获得的风速与浮标风速之间的均方根误差减小了30%,平均误差减小了83%。在利用谱模型算法反演高度计风速时,谱模型中的波龄因子(表示海浪成长状态)可以根据高度计测得的有效波高和风速获得,因此该方法具有广泛的适用性。     

An evaluation of the Arctic clouds and surface radiative fluxes in CMIP6 models

To assess the performances of state-of-the-art global climate models on simulating the Arctic clouds and surface radiation balance, the 2001–2014 Arctic Basin surface radiation budget, clouds, and the cloud radiative effects(CREs) in 22 coupled model intercomparison project 6(CMIP6) models are evaluated against satellite observations. For the results from CMIP6 multi-model mean, cloud fraction(CF) peaks in autumn and is lowest in winter and spring, consistent with that from three satellite observation products(Cloud Sat-CALIPSO, CERESMODIS, and APP-x). Simulated CF also shows consistent spatial patterns with those in observations. However,almost all models overestimate the CF amount throughout the year when compared to CERES-MODIS and APP-x.On average, clouds warm the surface of the Arctic Basin mainly via the longwave(LW) radiation cloud warming effect in winter. Simulated surface energy loss of LW is less than that in CERES-EBAF observation, while the net surface shortwave(SW) flux is underestimated. The biases may result from the stronger cloud LW warming effect and SW cooling effect from the overestimated CF by the models. These two biases compensate each other,yielding similar net surface radiation flux between model output(3.0 W/m~2) and CERES-EBAF observation(6.1 W/m~2). During 2001–2014, significant increasing trend of spring CF is found in the multi-model mean,consistent with previous studies based on surface and satellite observations. Although most of the 22 CMIP6 models show common seasonal cycles of CF and liquid water path/ice water path(LWP/IWP), large inter-model spreads exist in the amounts of CF and LWP/IWP throughout the year, indicating the influences of different cloud parameterization schemes used in different models. Cloud Feedback Model Intercomparison Project(CFMIP)observation simulator package(COSP) is a great tool to accurately assess the performance of climate models on simulating clouds. More intuitive and credible evaluation results can be obtained based on the COSP model output. In the future, with the release of more COSP output of CMIP6 models, it is expected that those inter-model spreads and the model-observation biases can be substantially reduced. Longer term active satellite observations are also necessary to evaluate models' cloud simulations and to further explore the role of clouds in the rapid Arctic climate changes.     

Performance of the CCIW wave direction buoy at ARSLOE

The performance at the Atlantic Remote Sensing Land Ocean Experiment (ARSLOE) of the Canada Centre for Inland Waters (CCIW) wave direction buoy is examined. The waveheight and period data show good agreement with data from neighboring buoys. Nondimensional energy and fetch for wind sea cases agree well with data collected from fixed sensors in Lake Ontario. Direction measurements for a particularly well-defined swell agree well with the Experimental Environmental Research Buoy (XERB) measurements. Meteorological data are compared with that collected at the XERB. Windspeed estimates are reasonably well-correlated and much of the scatter is due to the 24-kin separation between the buoys. There is some systematic bias in the wind direction estimates.     

First six months quality assessment of HY-2A SCAT wind products using in situ measurements

The first Chinese microwave ocean environment satellite HY-2A, carrying a Ku-band scatteromenter （SCAT）, was successfully launched in August 2011. The first quality assessment of HY-2A SCAT wind products is presented through the comparison of the first 6 months operationally released SCAT products with in situ data. The in situ winds from the National Data Buoy Center （NDBC） buoys, R/V Polarstern, Aurora Australis, Roger Revelle and PY30-1 oil platform, were converted to the 10 m equivalent neutral winds. The temporal and spatial differences between the HY-2A SCAT and the in situ observations were limited to less than 5 min and 12.5 km. For HY-2A SCAT wind speed products, the comparison and analysis using the NDBC buoys yield a bias of-0.49 m/s, a root mean square error （RMSE） of 1.3 m/s and an increase negative bias with increasing wind speed observation above 3 m/s. Although less accurate of HY-2A SCAT wind direction at low winds, the RMSE of 19.19° with a bias of 0.92° is found for wind speeds higher than 3 m/s. These results are found consistent with those from R/Vs and oil platform comparisons. Moreover, the NDBC buoy comparison results also suggest that the accuracy of HY-2A SCAT winds is consistent over the first half year of 2012. The encouraging assessment results over the first 6 months show that wind products from HY-2A SCAT will be useful for scientific community.     

The seasonal cycle of diabatic heat storage in the Pacific Ocean

This study quantifies uncertainties in closing the seasonal cycle of diabatic heat storage (DHS) over the Pacific Ocean from 20°S to 60°N through the synthesis of World Ocean Circulation Experiment (WOCE) reanalysis products from 1993 to 1999. These products are DHS from Scripps Institution of Oceanography (SIO); near-surface geostrophic and Ekman currents from Earth and Space Research (ESR); and air-sea heat fluxes from Comprehensive Ocean-Atmosphere Data Set (COADS), National Centers for Environmental Prediction (NCEP), and European Center for Mid-Range Weather Forecasts (ECMWF). With these products, we compute residual heat budget components by differencing long-term monthly means from the long-term annual mean. This allows the seasonal cycle of the DHS tendency to be modeled. Everywhere latent heat flux residuals dominate sensible heat flux residuals, shortwave heat flux residuals dominate longwave heat flux residuals, and residual Ekman heat advection dominates residual geostrophic heat advection, with residual dissipation significant only in the Kuroshio-Oyashio current extension. The root-mean-square (RMS) of the differences between observed and model residual DHS tendencies (averaged over 10° latitude-by-20° longitude boxes) is <20 W m⁻² in the interior ocean and <100 W m⁻² in the Kuroshio-Oyashio current extension. This reveals that the residual DHS tendency is driven everywhere by some mix of residual latent heat flux, shortwave heat flux, and Ekman heat advection. Suppressing bias errors in residual air-sea turbulent heat fluxes and Ekman heat advection through minimization of the RMS differences reduces the latter to <10 W m⁻² over the interior ocean and <25 W m⁻² in the Kuroshio-Oyashio current extension. This reveals air-sea temperature and specific humidity differences from in situ surface marine weather observations to be a principal source of bias error, overestimated over most of ocean but underestimated near the Intertropical Convergence Zone.     

The winds of the comparison data set for the Seasat Gulf of Alaska experiment

The ship and data buoy winds used for comparison in the validation of Seasat-derived winds are described in terms of the time series of hourly wind observations from the buoys and in terms of the technique used to produce 20- and 30-min average winds from the ships. Sources of scatter in the comparison data are briefly reviewed.     

海洋观测浮标体水中平衡性分析

海洋观测浮标系统应用广泛,浮标体的平衡性能是影响海洋观测浮标使用的重要因素之一(着重观测波浪等水文参数除外)。以一种主体直径为2.3 m的浮标体为例,运用动力学平衡理论,通过数学分析法进行推演计算,分析后得到了此浮标体的重心和浮心的位置。其中,重心位于主浮体的中轴线偏下位置;按照浮标体倾斜一定角度的情况下,推算出浮标体浮心位置,并由此得出其稳心的位置。依据计算,得出浮标体倾斜20°的情况下,此浮标体依然处于稳态平衡,理论结果表明此浮标体设计具有一定的可行性。通过海上实际测试,验证了理论计算分析的正确性。文中对浮标体平衡性的研究采用了理论分析与实验印证相结合的方法,其分析计算方法具有普适性,适用于其他型号的海洋观测浮标体,乃至所有水面浮体的平衡性分析,可为水面浮体设计提供借鉴。