Sea surface temperature retrieval based on simulated microwave polarimetric measurements of a one-dimensional synthetic aperture microwave radiometer

Compared with traditional real aperture microwave radiometers, one-dimensional synthetic aperture microwave radiometers have higher spatial resolution. In this paper, we proposed to retrieve sea surface temperature using a one-dimensional synthetic aperture microwave radiometer that operates at frequencies of 6.9 GHz, 10.65 GHz,18.7 GHz and 23.8 GHz at multiple incidence angles. We used the ERA5 reanalysis data provided by the European Centre for Medium-Range Weather Forecasts and a radiation transmission forward model to calculate the model brightness temperature. The brightness temperature measured by the spaceborne one-dimensional synthetic aperture microwave radiometer was simulated by adding Gaussian noise to the model brightness temperature.Then, a backpropagation(BP) neural network algorithm, a random forest(RF) algorithm and two multiple linear regression algorithms(RE1 and RE2) were developed to retrieve sea surface temperature from the measured brightness temperature within the incidence angle range of 0°–65°. The results show that the retrieval errors of the four algorithms increase with the increasing Gaussian noise. The BP achieves the lowest retrieval errors at all incidence angles. The retrieval error of the RE1 and RE2 decrease first and then increase with the incidence angle and the retrieval error of the RF is contrary to that of RE1 and RE2.     

The brightness reversal of submarine sand waves in “HJ-1A/B” CCD sun glitter images

The brightness reversal of submarine sand waves appearing in the small satellite constellation for environment and disaster monitoring and forecasting("HJ-1A/B") CCD sun glitter images can affect the observation and depth inversion of sand wave topography. The simulations of the normalized sun glitter radiance on the submarine sand waves confirm that the reversal would happen at a specific sensor viewing angle, defined as the critical angle. The difference between the calculated critical angle position and the reversal position in the image is about 1, which is excellent in agreement. Both the simulation and actual image show that sand wave crests would be indistinct at the reversal position, which may cause problems when using these sun glitter images to analyze spatial characteristics and migration of sand waves. When using the sun glitter image to obtain the depth inversion, one should take the advantage of image properties of sand waves and choose the location in between the reversal position and the brightest position. It is also necessary to pay attention to the brightness reversal when using "HJ-1A/B" CCD images to analyze other oceanic features, such as internal waves, oil slicks, eddies, and ship wakes.     

利用2012年中国北极考察现场观测数据和卫星遥感资料评估海洋2号卫星扫描辐射计数据反演北极海冰密集度的研究

A retrieval algorithm of arctic sea ice concentration(SIC) based on the brightness temperature data of "HY-2" scanning microwave radiometer has been constructed. The tie points of the brightness temperature were selected based on the statistical analysis of a polarization gradient ratio and a spectral gradient ratio over open water(OW), first-year ice(FYI), and multiyear ice(MYI) in arctic. The thresholds from two weather filters were used to reduce atmospheric effects over the open ocean. SIC retrievals from the "HY-2" radiometer data for idealized OW, FYI, and MYI agreed well with theoretical values. The 2012 annual SIC was calculated and compared with two reference operational products from the National Snow and Ice Data Center(NSIDC) and the University of Bremen. The total ice-covered area yielded by the "HY-2" SIC was consistent with the results from the reference products. The assessment of SIC with the aerial photography from the fifth Chinese national arctic research expedition(CHINARE) and six synthetic aperture radar(SAR) images from the National Ice Service was carried out. The "HY-2" SIC product was 16% higher than the values derived from the aerial photography in the central arctic. The root-mean-square(RMS) values of SIC between "HY-2" and SAR were comparable with those between the reference products and SAR, varying from 8.57% to 12.34%. The "HY-2" SIC is a promising product that can be used for operational services.     

Detection of nighttime sea fog/stratus over the Huanghai Sea using MTSAT-1R IR data

A dual channel difference (DCD) method is applied to detect nighttime sea fog/stratus over the Huanghai Sea using the infrared (IR) data of shortwave (3.5–4.0 μm) and longwave (10.3–11.3 μm) channels from the Multi-functional Transport Satellite (MTSAT)-1R, i.e., shortwave minus longwave brightness temperature difference (SLTD). Twenty-four sea fog events over the Huanghai Sea during March to July of 2006 and 2007 are chosen to determine a suitable value of SLTD for nighttime sea fog/stratus detection, and ...     

Validation of the “HY-2” altimeter wet tropospheric path delay correction based on radiosonde data

Wet tropospheric path delay （PD） is a highly variable term for the altimeter measurement of a sea surface height, caused by the refraction effect of atmospheric water vapor and cloud liquid water. In order to esti- mate PD values, the ＂HY-2＂ system includes a calibration microwave radiometer （CMR） operating at 18.7, 23.8 and 37 GHz. The PD data of the CMR were compared and validated by coincident radiosonde profiles from ten globally distributed radiosonde stations during October 2011 to August 2012. The temporal interval was 1 h. In order to avoid land contamination, different spatial intervals between these two data sets were tested. The empirical fit function of PD uncertainty and spatial interval was found and extrapolated to the ideal situation that the data of CMR and radiosonde were totally coincident. The stability of the brightness temperature of the CMR and its impact on the PD correction was also studied. Consequently, the uncertainty of the PD algorithm of the CMR was estimated to be 2.1 cm.     

基于一维综合孔径微波辐射计的海面温度反演研究

Due to the low spatial resolution of sea surface temperature(T_S) retrieval by real aperture microwave radiometers,in this study, an iterative retrieval method that minimizes the differences between brightness temperature(T_B)measured and modeled was used to retrieve sea surface temperature with a one-dimensional synthetic aperture microwave radiometer, temporarily named 1 D-SAMR. Regarding the configuration of the radiometer, an angular resolution of 0.43° was reached by theoretical calculation. Experiments on sea surface temperature retrieval were carried out with ideal parameters; the results show that the main factors affecting the retrieval accuracy of sea surface temperature are the accuracy of radiometer calibration and the precision of auxiliary geophysical parameters. In the case of no auxiliary parameter errors, the greatest error in retrieved sea surface temperature is obtained at low T_S scene(i.e., 0.710 6 K for the incidence angle of 35° under the radiometer calibration accuracy of0.5 K). While errors on auxiliary parameters are assumed to follow a Gaussian distribution, the greatest error on retrieved sea surface temperature was 1.330 5 K at an incidence angle of 65° in poorly known sea surface wind speed(W)(the error on W of 1.0 m/s) over high W scene, for the radiometer calibration accuracy of 0.5 K.     

基于Sentinel-3载荷OLCI和SRAL数据的内波同步探测研究

The ocean and land color instrument(OLCI) and synthetic aperture radar altimeter(SRAL) installed aboard the Sentinel-3 satellite have been in orbit for operational uses. In this study, data collected from Sentinel-3 are used to investigate internal waves in the South China Sea. An internal wave is detected using an OLCI image with a resolution of 300 m, and an analysis was performed with a quasi-synchronous moderate-resolution imaging spectroradiometer(MODIS) image. The opposite characteristics of OLCI and MODIS images of the same internal wave are explained by the critical angle in brightness reversals. The unique observational geometry of the OLCI image and its influence on observations of internal waves are discussed. The distribution of σ0 and sea surface height anomalies(SSHAs) induced by internal waves are studied using SRAL records. The σ0 records of SRAL occasionally show less sensitivity to the modulation of internal waves, which may be attributed to the observational geometry, while SSHAs show obvious variations. The synchronous pairing of OLCI images and SRAL records are analyzed to extract the three-dimensional sea surface signatures induced by internal waves. The analysis demonstrates that the profile of SSHAs in the surface shows an opposite phase to the profiles of internal waves in the ocean. The opposite phase relationship, observed in the remote sensing view, is also confirmed with a laboratory experiment.     

A new scheme for retrieving ocean surface salinity from simulated multi-angular SMOS brightness temperature

The European Space Agency will launch the first salinity satellite for remotely sensing the global soil moisture and ocean salinity (SMOS) at a sun-synchronous orbit in 2009. One of the payloads on the satellite is a synthetic aperture microwave radiometer (MIRAS), which is an innovative instrument designed as a two-dimensional (2D) interferometer for acquiring brightness temperature (TB) at L-band (1.4 GHz). MIRAS allows measuring TB at a series of incidences for full polarizations. As the satellite travels, a given location within the 2D field of view is observed from different incidence angles. The authors develop a new scheme to retrieve the sea-surface salinity (SSS) from SMOS’s TB at multi-incidence angles in a pixel, utilizing the properties of emissivity changing with incidence angles. All measurements of a given Stokes parameter in a pixel are first fitted to incidence angles in three order polynomial, and then the smoothed data are used for retrieving the SSS. The procedure will remove the random noise in TB greatly. Furthermore, the new method shows that the error in retrieved SSS is very sensitive to the system biases in the calibrated TB of the sensor, but the error in the retrieval is also a system bias, which can be corrected by post-launch validation. Therefore, this method may also serve as a means to evaluate the calibration precision in TB.     

An analysis on the error structure and mechanism of soil moisture and ocean salinity remotely sensed sea surface salinity products

For the application of soil moisture and ocean salinity（SMOS） remotely sensed sea surface salinity（SSS） products,SMOS SSS global maps and error characteristics have been investigated based on quality control information.The results show that the errors of SMOS SSS products are distributed zonally,i.e.,relatively small in the tropical oceans,but much greater in the southern oceans in the Southern Hemisphere（negative bias） and along the southern,northern and some other oceanic margins（positive or negative bias）.The physical elements responsible for these errors include wind,temperature,and coastal terrain and so on.Errors in the southern oceans are due to the bias in an SSS retrieval algorithm caused by the coexisting high wind speed and low temperature; errors along the oceanic margins are due to the bias in a brightness temperature（TB） reconstruction caused by the high contrast between L-band emissivities from ice or land and from ocean; in addition,some other systematic errors are due to the bias in TB observation caused by a radio frequency interference and a radiometer receivers drift,etc.The findings will contribute to the scientific correction and appropriate application of the SMOS SSS products.     

基于AMSR-E数据的多年冰密集度反演算法研究

In recent years, the rapid decline of Arctic sea ice area(SIA) and sea ice extent(SIE), especially for the multiyear(MY) ice, has led to significant effect on climate change. The accurate retrieval of MY ice concentration retrieval is very important and challenging to understand the ongoing changes. Three MY ice concentration retrieval algorithms were systematically evaluated. A similar total ice concentration was yielded by these algorithms, while the retrieved MY sea ice concentrations differs from each other. The MY SIA derived from NASA TEAM algorithm is relatively stable. Other two algorithms created seasonal fluctuations of MY SIA, particularly in autumn and winter. In this paper, we proposed an ice concentration retrieval algorithm, which developed the NASA TEAM algorithm by adding to use AMSR-E 6.9 GHz brightness temperature data and sea ice concentration using 89.0GHz data. Comparison with the reference MY SIA from reference MY ice, indicates that the mean difference and root mean square(rms) difference of MY SIA derived from the algorithm of this study are 0.65×106 km2 and0.69×106 km2 during January to March, –0.06×106 km2 and 0.14×106 km2 during September to December respectively. Comparison with MY SIE obtained from weekly ice age data provided by University of Colorado show that, the mean difference and rms difference are 0.69×106 km2 and 0.84×106 km2, respectively. The developed algorithm proposed in this study has smaller difference compared with the reference MY ice and MY SIE from ice age data than the Wang's, Lomax' and NASA TEAM algorithms.     

绿潮Lansat影像滑动窗口自适应阈值全自动检测方法

Since the atmospheric correction is a necessary preprocessing step of remote sensing image before detecting green tide, the introduced error directly affects the detection precision. Therefore, the detection method of green tide is presented from Landsat TM/ETM plus image which needs not the atmospheric correction. In order to achieve an automatic detection of green tide, a linear relationship(y =0.723 x+0.504) between detection threshold y and subtraction x(x=λ_(nir)–λ_(red)) is found from the comparing Landsat TM/ETM plus image with the field surveys.Using this relationship, green tide patches can be detected automatically from Landsat TM/ETM plus image.Considering there is brightness difference between different regions in an image, the image will be divided into a plurality of windows(sub-images) with a same size firstly, and then each window will be detected using an adaptive detection threshold determined according to the discovered linear relationship. It is found that big errors will appear in some windows, such as those covered by clouds seriously. To solve this problem, the moving step k of windows is proposed to be less than the window width n. Using this mechanism, most pixels will be detected[n/k]×[n/k] times except the boundary pixels, then every pixel will be assigned the final class(green tide or sea water) according to majority rule voting strategy. It can be seen from the experiments, the proposed detection method using multi-windows and their adaptive thresholds can detect green tide from Landsat TM/ETM plus image automatically. Meanwhile, it avoids the reliance on the accurate atmospheric correction.     

A study on retrieval algorithms in oceanic passive microwave remote sensing using 1.4 GHz and 10.7 GHz bands

To retrieve sea-surface salinity (SSS) from radiometer data at 1.4 GHz, auxiliary data of sea-surface temperature (SST), surface roughness and meteorological variables are needed. The authors study oceanic passive polarimetric microwave remote sensing using 1.4 GHz and 10.7 GHz bands. A set of algorithms are developed for 1.4 GHz and 10.7 GHz microwave polarimetric radiometer at 50° incidence angle to retrieve wind vector, as well as other geophysical parameters, such as SSS, SST, atmospheric volumes of water vapor and liquid water. Idealized retrievals are conducted using 2324 simulated brightness temperatures of full Stokes parameters at 1.4 GHz and 10.7 GHz. Results indicate that SSS, SST, sea-surface wind speed, direction, atmospheric volumes of water vapor and liquid water can be inversed at the same time. This suggests an alternative way for SSS remote sensing.     

Origin and shallow-layer structure of the sediments in the Yellow Sea trough

On the basis of the analysis and study of the records of shallow-layer profiles and information from systematic analyses of some cores and their 14C datings, the sediments of Yellow Sea trough since the final slage the Late Pleistocene are found to have the following features. On the south slope of the trough, it is probably composed of island-shelf de posits. On the mouth of the trough, it is made up of eolian dunes and its deposits. In the centra! bottom region of the trough where intense desertization occurred in the early stage of the regression, there are parallel oblique beddings in the strata known as the "angle of repose" texture, which demonstrates that the deposits there are composed of eolian sands. And in the late stage, thin derivative deposits developed in the northern part of the trough. On the north slope of the trough, thick-layer derivative deposits developed. This primary sedimentary pattern still remains unchanged since the occurrence of the Holocene transgression.     

The Effect of Internal Fluid on the Response of Vortex-Induced Vibration of Marine Risers

Based on Iwan‘s wake oseillator model developed with the classical van der Pol equalion, the differential equation for the response of the vortex-induced vihration (VIV) of the riser considering the effect of the internal flowing fluid and the external marine environmental condition is derived. The effect of the intermal flowing fluid on the response of VIV of the riser is studied by means of the Finite Element Method. The results show that the effect of the internal fluid velocity on the VIV of the riser is strong when the natural frequency of the riser is close to the vortex shedding frequency. In addition, the increase of the top tension can decrease the sensitivity of the riser to the internal fluid velocity.     

A stability analysis for the sheared oil film contained by the boom

According to the observation in experiment of stability of the oil film, the assumption of velocity distribution for both the water flow and the oil film is introduced. On the basis of the assumption, Orr-Sommerfeld stability equation is applied to develop the method of determining the critical velocity of the oil film, and the criterion for stability of the oil film is obtained. Meanwhile, a formula describing the relation between the thickness of the oil film and the velocity of the water flow is also given and examined by the laboratory experiment.     

Some Problems in the Design of Breakwaters in Dayaowan Port Area

- In this paper, some problems in the design of the breakwater for Dayaowan port area are described including the optimal selection of the layout plan of the breakwater; the selection of structural type of the breakwater, the concrete block type for bank protection, the reversed L-shape parapet; the elevation of the parapet; the stability test of the whole revetment; as well as wave absorbing effects of the perforated     

Study on the Improvement of Sea-Entering Waterway in the Changjiang Estuary

The Changjiang Estuary is the doorway of the Shanghai Harbour and the throat for the six provinces in the Changjiang River drainage area. However the minimum water depth of the waterway in the estuary is only about six meters, so that the development of the Shanghai Harbour and the utilization of water transportation on the Changjiang River have been greatly restrained. Since 1975, the depth of the sea-entering waterway of the Changjiang Estuary has been successfully increased by one meter. This article has made a comprehensive summing-up about the selection of waterways, the determination of the line for dredging and disposal of the dredged material, thus providing reference data for future work in increasing the depth of the waterway of the Changjiang Estuary and for controlling similar waterways of other estuaries.     

Indicators and impact analysis of sediment from the Changjiang Estuary and East China Sea to

Based on the historical evolution of the Hangzhou Bay, by making use of the conclusions made by the previous research workers and the integration of concrete data, five distinct impact indicators of the sediment from the Changjiang Estuary and the East China Sea to the Hangzhou Bay are summarized. Numerical calculation and analysis indicate that the scouring and deposition of seabed in the Hangzhou Bay are subject to the direct impact of the evolution of the Changjiang Estuary, and the growth and decline and the direction of the sandy bar at Nanhuizui give traces to the sediment transport between the Changjiang Estuary and the Hangzhou Bay. The transport of sediment from the Changjiang Estuary to the Hangzhou Bay occurs mainly in winter and spring seasons and the increase of the Changjiang River runoff and the decrease of sediment charge have caused scouring in the northern coast of the Hangzhou Bay and the seabed erosion along the frontal margin of the Changjiang River Delta.     

Joint Distribution of Wave Periods and Rate of Change of Wave Surface Elevation

The rate of change of wave surface elevation is of much importance in ocean engineering,es-pecially for the determination of the limitation of wave breaking.This paper gives a kind of joint distribu-tion of wave periods and the rate of change of wave surface elevation by means of calculation of the two-or-der to four-order moment of the frequency spectrum based on the linear wave theory.For the first time,the distribution density function of wave periods determined by peaks is provided,and the conclusion isdrawn that the rate of change of wave surface elevation obeys the Rayleigh distribution.