Assessment of reprocessed sea surface height measurements derived from HY-2A radar altimeter and its application to the observation of 2015–2016 El Ni?o

Haiyang-2A(HY-2A) is China's first ocean dynamic environment satellite and the radar altimeter is one of its main payloads. One of the main purposes of the radar altimeter is to measure the sea surface height(SSH). The SSH determined from the altimeter range measurements includes some range and geophysical corrections. These corrections largely affect the accuracy of the SSH measurements. The range and the geophysical corrections are reprocessed and the altimeter waveforms in HY-2A sensor interim geophysical data set records(S-IGDR) are retracked from June 1, 2014 to June 14, 2014, and the accuracy of the reprocessed SSH measurements is evaluated.The methods of the range and geophysical corrections used to reprocess HY-2A altimeter data are validated by using these methods to reprocess the Jason-2 range and geophysical corrections and comparing the results with the range and geophysical corrections in Jason-2 geophysical dataset records(GDR) product. A crossover analysis is used to evaluate the accuracy of the reprocessed HY-2A SSH measurements. The standard deviation(STD) of the crossover SSH differences for HY-2A is around 4.53 cm while the STD of the SSH differences between HY-2A and Jason-2 is around 5.22 cm. The performance of the reprocessed HY-2A SSH measurements is significantly improved with respect to the SSH measurements derived from HY-2A interim geophysical dataset records(IGDR)product. The 2015–2016 El Ni?o has been the strongest El Ni?o event since 1997–1998. The range and the geophysical corrections in HY-2A IGDR are reprocessed and sea level anomalies are used to monitor the2015–2016 El Ni?o. The results show that the HY-2A altimeter can well observe the 2015–2016 El Ni?o.     

Trends of sea surface temperature and sea surface temperature fronts in the South China Sea during 2003–2017

The trends of the sea surface temperature(SST) and SST fronts in the South China Sea(SCS) are analyzed during2003–2017 using high-resolution satellite data. The linear trend of the basin averaged SST is 0.31°C per decade,with the strongest warming identified in southeastern Vietnam. Although the rate of warming is comparable in summer and winter for the entire basin, the corresponding spatial patterns of the linear trend are substantially different between them. The SST trend to the west of the Luzon Strait is characterized by rapid warming in summer, exceeding approximately 0.6°C per decade, but the trend is insignificant in winter. The strongest warming trend occurs in the southeast of Vietnam in winter, with much less pronounced warming in summer. A positive trend of SST fronts is identified for the coast of China and is associated with increasing wind stress. The increasing trend of SST fronts is also found in the east of Vietnam. Large-scale circulation, such as El Ni?o, can influence the trends of the SST and SST fronts. A significant correlation is found between the SST anomaly and Ni?o3.4 index, and the ENSO signal leads by eight months. The basin averaged SST linear trends increase after the El Ni?o event(2009–2010), which is, at least, due to the rapid warming rate causing by the enhanced northeasterly wind. Peaks of positive anomalous SST and negatively anomalous SST fronts are found to co-occur with the strong El Ni?o events.     

结合极化分解特征的SVM溢油检测研究

Marine oil spills have caused major threats to marine environment over the past few years.The early detection of the oil spill is of great significance for the prevention and control of marine disasters.At present,remote sensing is one of the major approaches for monitoring the oil spill.Full polarization synthetic aperture radarc SAR data are employed to extract polarization decomposition parameters including entropy(H) and reflection entropy(A).The characteristic spectrum of the entropy and reflection entropy combination has analyzed and the polarization characteristic spectrum of the oil spill has developed to support remote sensing of the oil spill.The findings show that the information extracted from(1–A)×(1–H) and(1–H)×A parameters is relatively evident effects.The results of extraction of the oil spill information based on H×A parameter are relatively not good.The combination of the two has something to do with H and A values.In general,when H0.7,A value is relatively small.Here,the extraction of the oil spill information using(1–A)×(1–H) and(1–H)×A parameters obtains evident effects.Whichever combined parameter is adopted,oil well data would cause certain false alarm to the extraction of the oil spill information.In particular the false alarm of the extracted oil spill information based on(1–A)×(1–H) is relatively high,while the false alarm based on(1–A)×H and(1–H)×A parameters is relatively small,but an image noise is relatively big.The oil spill detection employing polarization characteristic spectrum support vector machine can effectively identify the oil spill information with more accuracy than that of the detection method based on single polarization feature.     

孔压探杆贯入及潮汐作用下超孔压响应规律研究

Excess pore water pressure is an important parameter that can be used to analyze certain physical characteristics of sediment. In this paper, the excess pore water pressure of subseafloor sediment and its variation with tidal movement was measured following the installation of a wharf in Qingdao, China by using a fiber Bragg grating(FBG) piezometer. The results indicated that this FBG piezometer is effective in the field. The measured variation of excess pore water pressure after installation is largely explained by the dissipation of excess pore water pressure. The dissipation rate can be used to estimate the horizontal consolidation coefficient, which ranged from1.3×10~(–6) m~2/s to 8.1×10~(–6) m~2/s. The measured values during tidal phases are associated with the variability of tidal pressure on the seafloor and can be used to estimate the compressibility and the permeability of the sediment during tidal movement. The volume compression coefficient estimated from tidal oscillation was approximately2.0×10~(–11) Pa~(–1), which was consistent with the data from the laboratory test. The findings of this paper can provide useful information for in situ investigations of subseafloor sediment.     

Mechanistic Drifting Forecast Model for A Small Semi-Submersible Drifter Under Tide–Wind–Wave Conditions

Understanding the drifting motion of a small semi-submersible drifter is of vital importance regarding monitoring surface currents and the floating pollutants in coastal regions. This work addresses this issue by establishing a mechanistic drifting forecast model based on kinetic analysis. Taking tide–wind–wave into consideration, the forecast model is validated against in situ drifting experiment in the Radial Sand Ridges. Model results show good performance with respect to the measured drifting features, characterized by migrating back and forth twice a day with daily downwind displacements. Trajectory models are used to evaluate the influence of the individual hydrodynamic forcing. The tidal current is the fundamental dynamic condition in the Radial Sand Ridges and has the greatest impact on the drifting distance. However, it loses its leading position in the field of the daily displacement of the used drifter. The simulations reveal that different hydrodynamic forces dominate the daily displacement of the used drifter at different wind scales. The wave-induced mass transport has the greatest influence on the daily displacement at Beaufort wind scale 5–6; while wind drag contributes mostly at wind scale 2–4.     

上新世以来南海海水表层温度（SST）的变化

The reconstruction of sea surface temperature(SST) is a key issue in research on paleoceanography. The recently related studies are mainly concentrated on the Quaternary. The skeletons of a typical species of benthonic foraminifer(Amphistegina radiata) in the top 0–375.30 m interval of Well "Xike-1" reef core, Shidao Island, the Xisha Islands, are uniformly selected. The ratios of magnesium to calcium concentrations and other indicators are analyzed by an electron microprobe analysis(EMPA) with the purpose of estimating the paleo-SSTs since the Pliocene and further investigating the periodic change law of paleoclimate in the South China Sea(SCS).Meanwhile, the geologic significance of paleoclimatic events in the SCS is discussed with global perspectives. The results indicate that the paleo-SSTs reconstructed by the ratios of magnesium to calcium concentrations in the SCS show a general periodic trend of "high–low–high–low" in values changes since the Pliocene. By comparison,the fluctuations of reconstructed paleo-SSTs are much stronger in the Quaternary. Moreover, the variations of the ratios of magnesium to calcium concentrations in the A. radiata skeletons have recorded a series of major paleoclimatic events since the Pliocene, of which the Quaternary glaciation events and the Arctic ice cap forming events during the late Pliocene are more significant. Thus, using the changes of the ratios of magnesium to calcium concentrations in the A. radiata skeletons from Well "Xike-1" reef core to reflect the relative changes of paleo-SSTs is a relatively feasible and reliable way in the SCS, which is also proved by the correlation of drilling cores characteristics in this area.     

一种利用SMOS卫星观测的海表面盐度反演海水盐度廓线的方法

This paper proposes a new method to retrieve salinity profiles from the sea surface salinity(SSS) observed by the Soil Moisture and Ocean Salinity(SMOS) satellite. The main vertical patterns of the salinity profiles are firstly extracted from the salinity profiles measured by Argo using the empirical orthogonal function. To determine the time coefficients for each vertical pattern, two statistical models are developed. In the linear model, a transfer function is proposed to relate the SSS observed by SMOS(SMOS_SSS) with that measured by Argo, and then a linear relationship between the SMOS_SSS and the time coefficient is established. In the nonlinear model, the neural network is utilized to estimate the time coefficients from SMOS_SSS, months and positions of the salinity profiles. The two models are validated by comparing the salinity profiles retrieved from SMOS with those measured by Argo and the climatological salinities. The root-mean-square error(RMSE) of the linear and nonlinear model are 0.08–0.16 and 0.08–0.14 for the upper 400 m, which are 0.01–0.07 and 0.01–0.09 smaller than the RMSE of climatology. The error sources of the method are also discussed.     

降雨导致的海面粗糙度对Aquarius卫星盐度反演的影响研究

Rainfall has two significant effects on the sea surface, including salinity decreasing and surface becoming rougher,which have further influence on L-band sea surface emissivity. Investigations using the Aquarius and TRMM 3B42 matchup dataset indicate that the retrieved sea surface salinity(SSS) is underestimated by the present Aquarius algorithm compared to numerical model outputs, especially in cases of a high rain rate. For example, the bias between satellite-observed SSS and numerical model SSS is approximately 2 when the rain rate is 25 mm/h. The bias can be eliminated by accounting for rain-induced roughness, which is usually modeled by rain-generated ring-wave spectrum. The rain spectrum will be input into the Small Slope Approximation(SSA) model for the simulation of sea surface emissivity influenced by rain. The comparison with theoretical model indicated that the empirical model of rain spectrumis more suitable to be used in the simulation. Further, the coefficients of the rain spectrum are modified by fitting the simulations with the observations of the 2–year Aquarius and TRMM matchup dataset. The calculations confirm that the sea surface emissivity increases with the wind speed and rain rate. The increase induced by the rain rate is rapid in the case of low rain rate and low wind speed. Finally, a modified model of sea surface emissivity including the rain spectrum is proposed and validated by using the matchup dataset in May 2014. Compared with observations, the bias of the rain-induced sea surface emissivity simulated by the modified modelis approximately 1e–4, and the RMSE is slightly larger than 1e–3. With using more matchup data, thebias between model retrieved sea surface salinities and observationsmay be further corrected,and the RMSE may be reduced to less than 1 in the cases of low rain rate and low wind speed.     

Study on Sound-Speed Dispersion in A Sandy Sediment at Frequency Ranges of 0.5–3 kHz and 90–170 kHz

In order to study the properties of sound-speed dispersion in a sandy sediment, the sound speed was measured both at high frequency(90–170 k Hz) and low frequency(0.5–3 k Hz) in laboratory environments. At high frequency, a sampling measurement was conducted with boiled and uncooked sand samples collected from the bottom of a large water tank. The sound speed was directly obtained through transmission measurement using single source and single hydrophone. At low frequency, an in situ measurement was conducted in the water tank, where the sandy sediment had been homogeneously paved at the bottom for a long time. The sound speed was indirectly inverted according to the traveling time of signals received by three buried hydrophones in the sandy sediment and the geometry in experiment. The results show that the mean sound speed is approximate 1710–1713 m/s with a weak positive gradient in the sand sample after being boiled(as a method to eliminate bubbles as much as possible) at high frequency, which agrees well with the predictions of Biot theory, the effective density fluid model(EDFM) and Buckingham's theory. However, the sound speed in the uncooked sandy sediment obviously decreases(about 80%)both at high frequency and low frequency due to plenty of bubbles in existence. And the sound-speed dispersion performs a weak negative gradient at high frequency. Finally, a water-unsaturated Biot model is presented for trying to explain the decrease of sound speed in the sandy sediment with plenty of bubbles.     

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

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

Case study on the three-dimensional structure of meso-scale eddy in the South China Sea based on a high-resolution model

Meso-scale eddies are important features in the South China Sea(SCS). The eddies with diameters of 50–200 km can greatly impact the transport of heat, momentum, and tracers. A high-resolution wave-tide-circulation coupled model was developed to simulate the meso-scale eddy in the SCS in this study. The aim of this study is to examine the model ability to simulate the meso-scale eddy in the SCS without data assimilations The simulated Sea Surface Height(SSH) anomalies agree with the observed the AVISO SSH anomalies well. The simulated subsurface temperature profiles agree with the CTD observation data from the ROSE(Responses of Marine Hazards to climate change in the Western Pacific) project. The simulated upper-ocean currents also agree with the main circulation based on observations. A warm eddy is identified in winter in the northern SCS. The position and domain of the simulated eddy are confirmed by the observed sea surface height data from the AVISO. The result shows that the model has the ability to simulate the meso-scale eddy in the SCS without data assimilation.The three-dimensional structure of the meso-scale eddy in the SCS is analyzed using the model result. It is found that the eddy center is tilted vertically, which agrees with the observation. It is also found that the velocity center of the eddy does not coincide with the temperature center of the eddy. The result shows that the model has the ability to simulate the meso-scale eddy in the SCS without data assimilations. Further study on the forming mechanism and the three-dimensional structure of the meso-scale eddies will be carried out using the model result and cruise observation data in the near future.     

Regime shift of the dominant factor for halocline depth in the Canada Basin during 1990–2008

The World Ocean Database(WOD) is used to evaluate the halocline depth simulated by an ice-ocean coupled model in the Canada Basin during 1990–2008. Statistical results show that the simulated halocline is reliable.Comparing of the September sea ice extent between simulation and SSM/I dataset, a consistent interannual variability is found between them. Moreover, both the simulated and observed September sea ice extent show staircase declines in 2000–2008 compared to 1990–1999. That supports that the abrupt variations of the ocean surface stress curl anomaly in 2000–2008 are caused by rapid sea ice melting and also in favor of the realistic existence of the simulated variations. Responses to these changes can be found in the upper ocean circulation and the intermediate current variations in these two phases as well. The analysis shows that seasonal variations of the halocline are regulated by the seasonal variations of the Ekman pumping. On interannual time scale, the variations of the halocline have an inverse relationship with the ocean surface stress curl anomaly after 2000,while this relationship no longer applies in the 1990 s. It is pointed out that the regime shift in the Canada Basin can be derived to illustrate this phenomenon. Specifically, the halocline variations are dominated by advection in the 1990 s and Ekman pumping in the 2000 s respectively. Furthermore, the regime shift is caused by changing Transpolar Drift pathway and Ekman pumping area due to spatial deformation of the center Beaufort high(BH)relative to climatology.     

Estimation of annual variation of water vapor in the Arctic Ocean between 80°–87°N using shipborne GPS data based on kinematic precise point positioning

The measurement of atmospheric water vapor(WV) content and variability is important for meteorological and climatological research. A technique for the remote sensing of atmospheric WV content using ground-based Global Positioning System(GPS) has become available, which can routinely achieve accuracies for integrated WV content of 1–2 kg/m2. Some experimental work has shown that the accuracy of WV measurements from a moving platform is comparable to that of(static) land-based receivers. Extending this technique into the marine environment on a moving platform would be greatly beneficial for many aspects of meteorological research, such as the calibration of satellite data, investigation of the air-sea interface, as well as forecasting and climatological studies. In this study, kinematic precise point positioning has been developed to investigate WV in the Arctic Ocean(80°–87°N) and annual variations are obtained for 2008 and 2012 that are identical to those related to the enhanced greenhouse effect.     

Sea-surface acoustic backscattering measurement at 6–25 kHz in the Yellow Sea

Sea-surface acoustic backscattering measurements at moderate to high frequencies were performed in the shallow water of the south Yellow Sea, using omnidirectional spherical sources and omnidirectional hydrophones. Sea-surface backscattering data for frequencies in the 6–25 k Hz range and wind speeds of(3.0±0.5)and(4.5±1.0) m/s were obtained from two adjacent experimental sites, respectively. Computation of sea-surface backscattering strength using bistatic transducer is described. Finally, we calculated sea-surface backscattering strengths at grazing angles in the range of 16°–85°. We find that the measured backscattering strengths agree reasonably well with those predicted by using second order small-roughness perturbation approximation method with "PM" roughness spectrum for all frequencies at grazing angles ranged from 40° to 80°. The backscattering strengths varied slightly at grazing angles of 16°–40°, and were much stronger than roughness scattering. It is speculated that scattering from bubbles dominates the backscattering strengths at high wind speeds and small grazing angles. At the same frequencies and moderate to high grazing angles, the results show that the backscattering strengths at a wind speed of(4.5±1.0) m/s were approximately 5 d B higher than those at a wind speed of(3.0±0.5) m/s. However, the discrepancies of backscattering strength at low grazing angles were more than 10 d B. Furthermore the backscattering strengths exhibited no significant frequency dependence at 3 m/s wind speed. At a wind speed of 4.5 m/s, the scattering strengths increased at low grazing angles but decreased at high grazing angles with increasing grazing angle.     

Sediment Transport Capacity Under the River–Tide Interaction in the Changjiang Estuary

Sediment transport capacity is a fundamental parameter in sediment transport theory and its accurate calculation is important from both theoretical and engineering viewpoints. The capacity of sediment transport has been studied extensively by many researchers in the last decades. Nevertheless, the underlying mechanism behind sediment transport capacity in estuaries remains poorly understood. The current study aims to explore the impact of the river–tide interaction on sediment transport and establish a formula of sediment transport capacity under the river–tide interaction. The impact of the river–tide interaction on the hydrodynamics and sediment dynamics in the Changjiang Estuary was analyzed, a practical method for describing the variation in tide-runoff ratio was established,and a formula of sediment transport capacity considering the impact of river–tide interaction was proposed by introducing the tide-runoff ratio. The new method bridged the gap between two well-known sediment transport capacity methods by considering the variation in the index a for the gravitational term and overcomes the drawback of distinguishing flood/dry season or spring/ebb tide in the calculation of estuarine sediment transport. A large amount of flow and sediment data obtained from the Changjiang Estuary were collected to verify the proposed formula. The effect of salt-fresh water mixture and the morphological evolution on sediment transport capacity of the Changjiang Estuary were discussed.     

Acoustic backscattering measurement from sandy seafloor at 6–24 kHz in the South Yellow Sea

The acoustic bottom backscattering strength was measured at the frequency range of 6–24 kHz on a typical sandy bottom in the South Yellow Sea by using omnidirectional sources and omnidirectional receiving hydrophones. In the experiment, by avoiding disturbances due to scattering off the sea surface and satisfying the far-field condition, we obtained values of acoustic bottom backscattering strength ranging from –41.1 to –24.4 dB within a grazing angle range of 18°–80°. In the effective range of grazing angles, the acoustic scattering strength generally increases with an increase in the grazing angles, but trends of the variation were distinct in different ranges of frequency, which reflect different scattering mechanisms. The frequency dependence of bottom backscattering strength is generally characterized by a positive correlation in the entire frequency range of 6–24 kHz at the grazing angles of 20°, 40° and 60° with the linear regression slopes of 0.222 9 dB/kHz, 0.513 0 dB/kHz and 0.174 6 dB/kHz, respectively. At the largest grazing angle of 80°, the acoustic backscattering strength exhibits no evident frequency dependence.     

Lateral Global Buckling of Submarine Pipelines Based on the Model of Nonlinear Pipe–Soil Interaction

With the increasing development and utilization of offshore oil and gas resources, global buckling failures of pipelines subjected to high temperature and high pressure are becoming increasingly important. For unburied or semi-buried submarine pipelines, lateral global buckling represents the main form of global buckling. The pipe–soil interaction determines the deformation and stress distribution of buckling pipelines. In this paper, the nonlinear pipe–soil interaction model is introduced into the analysis of pipeline lateral global buckling, a coupling method of PSI elements and the modified RIKS algorithm is proposed to study the lateral global buckling of a pipeline, and the buckling characteristics of submarine pipeline with a single arch symmetric initial imperfection under different pipe–soil interaction models are studied. Research shows that, compared with the ideal elastic–plastic pipe–soil interaction model, when the DNV-RP-F109 model is adopted to simulate the lateral pipe–soil interactions in the lateral global buckling of a pipeline, the buckling amplitude increases, however, the critical buckling force and the initial buckling temperature difference decreases. In the DNV-RP-F109 pipe–soil interaction model, the maximum soil resistance, the residual soil resistance, and the displacement to reach the maximum soil resistance have significant effects on the analysis results of pipeline global buckling.     

Upper ocean high resolution regional modeling of the Arabian Sea and Bay of Bengal

In this paper, effort is made to demonstrate the quality of high-resolution regional ocean circulation model in realistically simulating the circulation and variability properties of the northern Indian Ocean(10°S–25°N,45°–100°E) covering the Arabian Sea(AS) and Bay of Bengal(BoB). The model run using the open boundary conditions is carried out at 10 km horizontal resolution and highest vertical resolution of 2 m in the upper ocean.The surface and sub-surface structure of hydrographic variables(temperature and salinity) and currents is compared against the observations during 1998–2014(17 years). In particular, the seasonal variability of the sea surface temperature, sea surface salinity, and surface currents over the model domain is studied. The highresolution model's ability in correct estimation of the spatio-temporal mixed layer depth(MLD) variability of the AS and BoB is also shown. The lowest MLD values are observed during spring(March-April-May) and highest during winter(December-January-February) seasons. The maximum MLD in the AS(BoB) during December to February reaches 150 m (67 m). On the other hand, the minimum MLD in these regions during March-April-May becomes as low as 11–12 m. The influence of wind stress, net heat flux and freshwater flux on the seasonal variability of the MLD is discussed. The physical processes controlling the seasonal cycle of sea surface temperature are investigated by carrying out mixed layer heat budget analysis. It is found that air-sea fluxes play a dominant role in the seasonal evolution of sea surface temperature of the northern Indian Ocean and the contribution of horizontal advection, vertical entrainment and diffusion processes is small. The upper ocean zonal and meridional volume transport across different sections in the AS and BoB is also computed. The seasonal variability of the transports is studied in the context of monsoonal currents.     

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.     

基于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.