An approach for Baltic Dry Index analysis based on empirical mode decomposition

The bulk shipping market is seasonal, cyclical and highly volatile. Due to the nonstationary and nonlinear nature of price series and the complexity of influencing factors, it is difficult to analyse the fluctuations in the bulk shipping market. In this study, a method based on empirical mode decomposition (EMD) is proposed to investigate the volatility of the Baltic Dry Index (BDI). In this method, the original freight price series is decomposed into several independent intrinsic modes, using EMD first. Then, the intrinsic modes are composed into three components: short-term fluctuations caused by normal market activities, the effect of extreme events and a long-term trend. Numerical experiments indicate that the proposed method can effectively reveal the characteristics of bulk freight price series with different economic meanings and decrease error accumulation. Meanwhile, by decomposition of intrinsic modes, the complexity of the model formulation can be controlled and the operability of the model can be improved.     

EMD在海洋内波研究中的应用

采用南海内波的实测资料,利用EMD方法对内波波动进行水平和垂向的分解,并由低频模态拟合函数并叠加得到了内波引起的温度波动的函数,由垂向位移得到了内波垂向模态的一种获取方法,为非线性的海洋内波研究提供了依据。     

基于Hilbert-Huang变换的船舶结构冲击响应研究

针对水下爆炸时船舶结构遭受的冲击响应信号具有非线性非平稳的特点,提出将Hilbert-Huang变换用于船舶结构冲击响应信号处理.Hilbert-Huang变换是基于经验模态分解(EMD)和Hilbert谱的一种信号处理方法,用EMD分解把时间序列信号分解成不同特征时间尺度的固有模态函数(IMF),然后对IMF分量进行Hilbert变换,从本质上分析船舶结构冲击响应信号的组成成分及特点.与FFT变换和小波变换相比,Hilbert-Huang变换体现出自适应性和先进性,可以有效提取船舶冲击响应信号的时频特征,揭示船舶结构自身的动态特性对冲击响应的影响.     

Long-term sea surface temperature baselines—time series, spatial covariation and implications for biological processes

Coastal areas such as estuaries, bays and fjords usually have hydrographic characteristics (e.g., temperature, salinity) which differ from those at larger spatial scales and in offshore areas. The differences can arise if the areas are subject to different climatic forcing or if they are relatively isolated from each other due to topographic and ocean circulation features which inhibit advective inputs of water mass properties. Local differences in hydrographic conditions can therefore potentially limit the applicability of existing long time series of coastally monitored temperatures for addressing questions at large spatial scales, such as the response of species distributions and phenologies to climate change. In this study we investigate the spatial synchrony of long-term sea surface temperatures in the North Sea–Baltic Sea region as measured daily at four coastal sites (Marsdiep, Netherlands; Torungen, Norway; Skagens Reef, Denmark; and Christiansø, Denmark) and in several large offshore areas. All time series, including two series reconstructed and intercalibrated for this study (Skagens Reef and Christiansø, Denmark), began during 1861–1880 and continue until at least 2001. Temperatures at coastal sites co-varied strongly with each other and with opportunistically measured offshore temperatures despite separation distances between measuring locations of 20–1200 km. This covariance is probably due to the influence of large-scale atmospheric processes on regional temperatures and is consistent with the known correlation radius of atmospheric fluctuations (ca. 1000 km). Differences (e. g, long-term trends, amplitude of seasonal variations) between coastal temperatures and those measured in adjacent offshore areas varied nonrandomly over time and were often significantly autocorrelated up to 2 years. These differences suggest that spatial variations in physical oceanographic phenomena and sampling heterogeneities associated with opportunistic sampling could affect perceptions of biological responses to temperature fluctuations. The documentation that the coastally measured temperatures co-vary with those measured opportunistically in offshore areas suggests that the coastal data, which have been measured daily using standardized methods and instruments, contain much of the variability seen at larger spatial scales. We conclude that both types of time series can facilitate assessments of how species and ecosystems have responded to past temperature changes and how they may react to future temperature changes.     

发动机轴承故障诊断方法的研究

针对EMD（Empirical Mode Decomposition）方法中存在的端点效应和IMF（Intrinsic Mode Function）虚假分量过多的问题,提出了基于互相关的EMD方法。首先,对非平稳性信号进行互相关延拓消除端点效应;其次,对消除端点效应的信号进行EMD分解,并将分解后的IMF与原信号作互相关,保留与原信号最相关的IMF;最后,作出信号的Hilbert边际谱,识别信号的频域特征。仿真结果表明该方法能够有效地克服端点效应,分离出真实IMF。将其应用于船舶发动机滚动轴承故障诊断中,能有效地识别出故障特征。     

Application of new parameterizations of gas transfer velocity and their impact on regional and global marine CO2 budgets

One of the dominant sources of uncertainty in the calculation of air–sea flux of carbon dioxide on a global scale originates from the various parameterizations of the gas transfer velocity, k, that are in use. Whilst it is undisputed that most of these parameterizations have shortcomings and neglect processes which influence air–sea gas exchange and do not scale with wind speed alone, there is no general agreement about their relative accuracy.The most widely used parameterizations are based on non-linear functions of wind speed and, to a lesser extent, on sea surface temperature and salinity. Processes such as surface film damping and whitecapping are known to have an effect on air–sea exchange. More recently published parameterizations use friction velocity, sea surface roughness, and significant wave height. These new parameters can account to some extent for processes such as film damping and whitecapping and could potentially explain the spread of wind-speed based transfer velocities published in the literature.We combine some of the principles of two recently published k parameterizations [Glover, D.M., Frew, N.M., McCue, S.J. and Bock, E.J., 2002. A multiyear time series of global gas transfer velocity from the TOPEX dual frequency, normalized radar backscatter algorithm. In: Donelan, M.A., Drennan, W.M., Saltzman, E.S., and Wanninkhof, R. (Eds.), Gas Transfer at Water Surfaces, Geophys. Monograph 127. AGU,Washington, DC, 325–331; Woolf, D.K., 2005. Parameterization of gas transfer velocities and sea-state dependent wave breaking. Tellus, 57B: 87–94] to calculate k as the sum of a linear function of total mean square slope of the sea surface and a wave breaking parameter. This separates contributions from direct and bubble-mediated gas transfer as suggested by Woolf [Woolf, D.K., 2005. Parameterization of gas transfer velocities and sea-state dependent wave breaking. Tellus, 57B: 87–94] and allows us to quantify contributions from these two processes independently.We then apply our parameterization to a monthly TOPEX altimeter gridded 1.5° × 1.5° data set and compare our results to transfer velocities calculated using the popular wind-based k parameterizations by Wanninkhof [Wanninkhof, R., 1992. Relationship between wind speed and gas exchange over the ocean. J. Geophys. Res., 97: 7373–7382.] and Wanninkhof and McGillis [Wanninkhof, R. and McGillis, W., 1999. A cubic relationship between air−sea CO2 exchange and wind speed. Geophys. Res. Lett., 26(13): 1889–1892]. We show that despite good agreement of the globally averaged transfer velocities, global and regional fluxes differ by up to 100%. These discrepancies are a result of different spatio-temporal distributions of the processes involved in the parameterizations of k, indicating the importance of wave field parameters and a need for further validation.     

等效波的散货船舱口角疲劳评估应力组合方法(英文)

The stress combination method for the fatigue assessment of the hatch corner of a bulk carrier was investigated based on equivalent waves.The principles of the equivalent waves of ship structures were given,including the determination of the dominant load parameter,heading,frequency,and amplitude of the equivalent regular waves.The dominant load parameters of the hatch corner of a bulk carrier were identified by the structural stress response analysis,and then a series of equivalent regular waves were defined based on these parameters.A combination method of the structural stress ranges under the different equivalent waves was developed for the fatigue analysis.The combination factors were obtained by least square regression analysis with the stress ranges derived from spectral fatigue analysis as the target value.The proposed method was applied to the hatch corner of another bulk carrier as an example.This shows that the results from the equivalent wave approach agree well with those from the spectral fatigue analysis.The workload is reduced substantially.This method can be referenced in the fatigue assessment of the hatch corner of a bulk carrier.     

Microscale variability in the distributions of large fluorescent particles observed in situ with a planar laser imaging fluorometer

It has been known for decades that particle-size and biomass spectra show regular patterns in the ocean, and that these patterns often show systematic variations with other properties such as total biomass, nutrient concentration, season, and distance (both vertical and horizontal). The recent finding of the ubiquitous nature of layers of phytoplankton < 1 m thick prompted us to explore the fine- and microscale vertical variations of size- and fluorescence-abundance spectra in the ocean. Using a two-dimensional planar laser imaging system mounted on a free-falling platform, we quantified the properties of large fluorescent particles ( 20 μm–2 cm) through the water column, obtaining images every 10–30 cm. These images showed systematic relationships of the spectral properties to total chlorophyll: increased proportions of the smallest particles at high chlorophyll concentrations, and a lengthening of the spectral size range at high total chlorophyll concentrations (more large particles at high chlorophyll concentrations). Further, we observed significant variations of the spectral properties over scales of 1 m and less, and recorded the frequent occurrence of unusual layers of large particles. Our new instrument, which is sensitive to thin layers of enhanced phytoplankton biomass, shows the planktonic community to be highly structured vertically on scales of 1–2 m, particularly within the DCM.     

Numerical weight function method for the structural analysis of ships: a speedy direct calculation with condensed structural information

The weight function method was originally derived for crack problems to calculate stress intensity factors for arbitrary loading conditions. In this article, a numerical weight function method has been extended to formulate the structural response analyses of two-dimensional elasticity, plate-bending, and three-dimensional plate-structures by using the finite-element method. The solution procedure is based on the well-known Maxwell–Betti reciprocal theorem, which is applied to the original and properly defined auxiliary problems. The present numerical weight function may be considered as a finite-element version of a Green’s function in an integral equation solution scheme. Although ship structures are certainly analysed by the finite-element method in a practical design procedure, the weight function approach has not yet been realized. The method is very useful for the analysis of structures subjected to a vast range of loading conditions, because structural responses can simply be calculated by the inner product of the universal weight function and load vectors. The validity and convergence characteristics of the present method are investigated by two-dimensional elastic and plate-bending problems, respectively. Finally, the method is applied to the calculation of the response amplitude operator of a stress component at a critical structural detail of a double-hull tanker, and the speed and efficiency of the method are quantitatively discussed based on the practical results.     

船舶结构疲劳评估设计波法

以散货船船舯内底与底边舱斜板折角处节点为例,对船舶结构疲劳评估的设计波法进行研究.论述了设计波法的基本原理,给出了确定设计波的流程。通过分析该节点的受力情况,确定5个控制载荷参数及设计波参数。以谱分析法得到的参考应力范围为目标值,对5个控制载荷参数对应的设计波作用下的节点应力范围结果进行回归分析,从而得到了用于疲劳评估的设计波组合,并利用该方法对1艘散货船的船舯内底与底边舱斜板折角处的节点进行了计算。结果显示,设计波法与谱分析法的结果吻合较好,计算量大幅减少。     

The combined effect of rain and wind on air–water gas exchange: A feasibility study

A series of experiments were conducted at University of Delaware's Air–Sea Interaction Laboratory to examine the combined effects of rain and wind on air–water gas exchange. During this study, ASIL WRX I, a combination of 3 rain rates and 4 wind speeds were used, for a total of 12 different environmental conditions. The SF₆ evasion method was used to determine the bulk gas transfer velocities, and airside profiles of wind and CO₂ were used to estimate flux–profiles of momentum and carbon dioxide. In addition to measurements of fluxes with and without rain in a wind–wave boundary layer, measurements of wave properties were also obtained. Rain is shown to alter the wind profile in the flume, and dampen surface waves. Also, SF₆ evasion indicates that with the present experimental setup, for most of the experimental conditions, rain and wind combine linearly to influence air–water gas exchange. Flux–profile relationships for marine atmospheric boundary layers, which were performed to scale up to field measurements, were explored by a comparison between SF₆-derived bulk fluxes and airside CO₂ profile measurements.     

Tuning of transfer functions for analysis of wave–ship interactions

Transfer functions are often used together with a wave spectrum for analysis of wave–ship interactions, where one application addresses the prediction of wave-induced motions or other types of global responses. This paper presents a simple and practical method which can be used to tune the transfer function of such responses to facilitate improved prediction capability. The input to the method consists of a measured response, i.e. time series sequences from a given sensor, the 2D wave spectrum characterising the seaway in which the measurements are taken, and an initial estimate of the transfer function for the response in study. The paper presents results obtained using data from an in-service container ship. The 2D wave spectra are taken from the ERA5 database, while the transfer function is computed by a simple closed-form expression. The paper shows that the application of the tuned transfer function leads to predictions which are significantly improved compared to using the transfer function without tuning.     

大开口型船舶在不规则波中的波浪扭矩预报分析

综述了现有船舶波浪扭矩的计算分析方法，用STF切片理论方法计算了大开口船舶波浪扭矩的传递函数。根据不同波谱理论模型对不规则波中的扭矩作了统计预报，讨论了波浪扭矩予报对波谱的敏感性。研究结果表明，波浪扭矩对波谱谱型的依赖程度相当大；特征波浪周期对统计值的影响比有义波高大；对相同有义波高，PM谱相应的预报值比ITTC二参数谱的结果更偏于保守。     

Using OMP analysis to observe temporal variability in water mass distribution

Previous optimum multiparameter analyses have investigated spatial distribution of water masses and reported successful results. This study tests the method further with application to a time series to examine temporal variability in water mass distribution. Observations collected from the Sargasso Sea are utilised, mainly to detect well-documented property changes that have occurred in Labrador Sea Water, but also possible fluctuations in water mass contributions as a whole.The results show much structure and variability in the contributions of Central Water and upper deep waters in the depth range 900–2000 m. It is proposed that the passage of warm and cold core eddies is indicated by these changes in relative contributions of water masses, although quantitative comparison of satellite altimetry data with the model results is not conclusive.Significant irregularities in the distribution of Labrador Sea Water are also revealed by the model. Through-flow of modified (since the early 1990s) Labrador Sea Water is denoted in the results by an apparent absence from mid-1995 to early 1998. This suggests a transit time of 5–6 years between the Labrador Sea and Bermuda, which is consistent with observational findings described in the literature.     

Air–sea gas transfer velocity estimates from the Jason-1 and TOPEX altimeters: Prospects for a long-term global time series

Estimation of global and regional air–sea fluxes of climatically important gases is a key goal of current climate research programs. Gas transfer velocities needed to compute these fluxes can be estimated by combining altimeter-derived mean square slope with an empirical relation between transfer velocity and mean square slope derived from field measurements of gas fluxes and small-scale wave spectra [Frew, N.M., Bock, E.J., Schimpf, U., Hara, T., Hauβecker, H., Edson, J.B., McGillis, W.R., Nelson, R.K., McKenna, S.P., Uz, B.M., Jähne, B., 2004. Air–sea gas transfer: Its dependence on wind stress, small-scale roughness and surface films, J. Geophys. Res., 109, C08S17, doi: 10.1029/2003JC002131.]. We previously reported initial results from a dual-frequency (Ku- and C-band) altimeter algorithm [Glover, D.M., Frew, N.M., McCue, S.J., Bock, E.J., 2002. A Multi-year Time Series of Global Gas Transfer Velocity from the TOPEX Dual Frequency, Normalized Radar Backscatter Algorithm, In: Gas Transfer at Water Surfaces, editors: Donelan, M., Drennan, W., Saltzman, E., and Wanninkhof, R., Geophysical Monograph 127, American Geophysical Union, Washington, DC, 325–331.] for estimating the air–sea gas transfer velocity (k) from the mean square slope of short wind waves (40–100 rad/m) and derived a 6-year time series of global transfer velocities based on TOPEX observations. Since the launch of the follow-on altimeter Jason-1 in December 2001 and commencement of the TOPEX/Jason-1 Tandem Mission, we have extended this time series to 12 years, with improvements to the model parameters used in our algorithm and using the latest corrected data releases. The prospect of deriving multi-year and interdecadal time series of gas transfer velocity from TOPEX, Jason-1 and follow-on altimeter missions depends on precise intercalibration of the normalized backscatter. During the Tandem Mission collinear phase, both satellites followed identical orbits with a mere 73-s time separation. The resulting collocated, near-coincident normalized radar backscatter (σ°) data from both altimeters present a unique opportunity to intercalibrate the two instruments, compare derived fields of transfer velocity and estimate the precision of the algorithm. Initial results suggest that the monthly gas transfer velocity fields generated from the two altimeters are very similar. Comparison of along-track Ku-band and C-band σ° during the collinear phase indicates that observed discrepancies are due primarily to small offsets between TOPEX and Jason-1 σ°. The Jason-1 k values have an apparent bias of + 4% relative to TOPEX, while the precision estimated from the two observation sets is 5–7% and scales with k. The resultant long-term, global, mean k is 16 cm/h.     

Whitecap coverage in coastal environment for steady and unsteady wave field conditions

Breaking waves represent a “key” parameter for many applications involved with a large number of environmental phenomena. In particular, it is well recognized that the whitecap cover induced by breaking waves allows substantial enhancement of heat, momentum, gas and particle transfer at the air–sea interface. A large number of studies were conducted during the last decades on the variation of the whitecap fraction, commonly noted W. The results presented in this paper deal with the evolution of the whitecap coverage in coastal zone. In such areas, the wave field is often unsteady with an important variety of sea state developments. The present analysis is based on an extensive series of data obtained during an experimental campaign which took place on the Mediterranean coast in 2001. The results allow observation of the influence of the sea state conditions of the wave field on the whitecap coverage. In addition, this paper confirms the occurrence of a peak in the variations of the whitecap fraction with the wave age for coastal areas as suggested by Lafon et al. [Lafon, C., Piazzola, J., Forget, P., Le Calvé, O. and Despiau, S., 2004. Analysis of the variations of the whitecap fraction as measured in a coastal zone. Boundary-Layer Meteorol., 111: 339-360.]. A wave age dependent model for the whitecap fraction is then proposed, which takes into account both the wind and the wave influence, and hence, is characteristics of the different sea state conditions.     

Evaluation of two atmospheric models for wind–wave modelling in the NW Mediterranean

The NW Mediterranean experiences, as illustrated by the last decade, strong and rapidly varying storms with severe waves and winds. This has motivated a continuous validation of models and the efforts to improve wave and wind predictions. In this paper we use two atmospherics models, MASS (from SMC-Meteorological Office of Catalunya) and ARPEGE (from Météo-France), to force two third generation wave models: WAM and SWAN. The evaluation and comparison has been carried out for two severe storms registered in November 2001 and March–April 2002.The ARPEGE and MASS models predicted higher 10 m wind speeds than coastal meteorological stations, a fact attributed to local land influences. Regarding the 10 m wind direction, models do not present large differences, although considerable deviations from recorded data were found during some dates. ARPEGE presents less scatter and lower errors than MASS when compared with QuikSCAT data.The 10m wind fields from both atmospheric models were used to force the two selected wave models and analyse the errors and sensitivities when predicting severe wave storms. The wave model simulations show some interesting results; during the storm, the spatial wave pattern using ARPEGE showed a higher maximum, although the values of significant wave height at the buoys were lower than the ones forced by MASS (with both WAM and SWAN). The SWAN simulations show a better agreement in predicting the growing and waning of the storm peaks. The prediction of mean period was improved when using the ARPEGE wind field. However the underestimation by SWAN due to the large energy at high frequencies was evident. Validation of spectral shape predictions showed that it still has considerable error when predicting the full frequency spectra. The storms showed bimodal spectral features which were not always reproduced by wave models and are likely to be responsible for part of the discrepancies.     

SWT方法在内部充水弹性球壳散射的应用

采用Sommerfeld-Watson变换(SWT)方法研究了内部充水弹性球壳的散射,应用绕数积分求根方法在频率域复平面和波数域复平面求解环绕波的频散方程。比较了弹性球壳在内部真空和内部充水两种情况下环绕波相速度曲线和衰减曲线,结果表明内部流体负载使得相对于内部真空情况有大量流体环绕波产生,相邻流体环绕波相速度曲线的"排斥"现象是明显的,流体波向周围媒质的再辐射作用在其相速度接近内部真空球壳弹性环绕波相速度的时候达到局部极值。用镜反射波和环绕波的叠加合成反向散射形态函数,与简正级数解符合良好。     

组合波阻技术的波动力响应矩阵分析法及其特性研究

基于有限元思想,综合运用波分析法和阻抗法,提出了波阻元件阻抑结构声传递的波动力响应矩阵分析法。将结构离散为多个波导单元和波阻单元,根据连接节点的位移连续,力与力矩平衡,建立附加波阻元件的结构连接广义波动力响应平衡方程,推导出波单元波动力响应特征矩阵及波阻元件附加波动响应特征矩阵,并代入波动力响应平衡方程求解得到波单元的振动幅值,从而求得传递效率与损失。运用该方法对阻振质量,粘弹性夹层以及动力吸振器的波阻特性进行了数值分析。最后,重点分析了组合波阻技术的波阻特性。研究表明,根据不同类型波阻元件的波阻特性,进行科学的组合与优化布置,并选择合理的设计参数,能显著提高组合波阻元件的阻抑效果。以上研究为组合波阻元件的声学设计提供了分析方法及新的控制策略,在舰船等结构减振降噪中具有重要的理论意义与工程应用价值。     

Acoustic excitation of hull surfaces by propeller sources

We present preliminary results from our new ship motion model that includes both strong and weak three-dimensional interactions between environmental surface waves and ship bodies in arbitrary water depth. The linear solutions of steady flow using the new model agree well with those obtained using the Green function methods. When the Froude number Fn is large, the fully nonlinear solutions of our model are significantly different from linear solutions, even in calm water. The interactions between the ship and incident gravity waves are completely different from those in linear solutions even with small Fn and moderate-amplitude surface waves (e.g., Fn = 0.25 and a significant wave height of about 1–3m).