基于卫星云图的轮廓法自动识别热带气旋研究

以卫星云图的热带气旋云系为研究对象,基于Canny边缘检测、Contour轮廓提取等方法,研究了通过图像的形状特征和范围大小对热带气旋云系进行自动识别,得到了较好的效果。这种方法有助于提高热带气旋云系识别的自动化程度,实现热带气旋的自动跟踪,从而为热带气旋的预报提供了便利。     

Availability of an equation to evaluate error by optical path discretization in radiative transfer computation based on the successive order of scattering method

The availability of an equation to evaluate the influence of multiple scattering in the single scattering process corresponding to a layer of arbitrary optical thickness was established. In order to confirm the validity of this equation, the radiance distribution in this layer was computed using a plane–parallel layer model based on the successive order of scattering method. The relative errors in a radiance distribution computed were evaluated as a function of optical thickness by the derived equation. It was shown that this equation provides a theoretical background for determining layer thickness using the plane–parallel layer model.     

Iceberg drifting and distribution in the Vilkitsky Strait studied by detailed satellite radar and optical images

This paper is devoted to the detection and identification of icebergs in the Russian Arctic Seas from the use of high- and medium-resolution radar and optical images from EROS-B, Radarsat-1, Radarsat-2, SPOT-4 and SPOT-5 Earth observation satellites. In July–September of 2011–2013, the SCANEX Research and Development Center, the Federal State Unitary Enterprise Atomflot, and other partner organizations provided operational satellite monitoring of icebergs in the Kara Sea and the Laptev Sea. More than 130 highly detailed optical and radar images were received and processed. The Vilkitsky Strait—one of the narrowest and most dangerous places within the Northern Sea Route—was chosen as an experimental polygon. As a result, iceberg location in the strait during the 2011–2013 navigation periods was analyzed, as were the iceberg size, area, drift direction, and height.     

Collocation resolution of the oceanic geoid based on satellite altimetry

The approach presented is directed toward a specific adaptation of the least‐squares collocation with noise, yielding smooth predictions of geophysical quantities. The smoothing corresponds here to a truncated gravity field equivalent to an (n’, n') spherical‐harmonic expansion. This is reflected in the truncation, at the degree n‘, of the pertinent covariance and cross‐covariance functions in most (but not all) instances. The smooth predictions of geophysical quantities, made in an equilateral grid corresponding to the truncation degree n‘, serve in constructing contour maps after having been densified for the needs of a contour routine. Such a densification is carried out efficiently via errorless collocation with the degree truncation n‘ throughout. Consistent with this procedure, “residuals” at observation points (i.e., discrepancies between the contour map and the data) are computed using the same algorithm. The complete collocation approach is utilized for a 2° resolution of the earth's gravity field with emphasis on the oceanic geoid, based on the residuals from a global spherical‐harmonic adjustment of SEASAT altimetry. The presented results include contour maps of geoid undulations and gravity anomalies. They are compared to the results of a point‐mass adjustment, another technique based on the spherical‐harmonic adjustment. The agreement between these two techniques is found to be excellent.     

Wave height forecasting by the transfer function model

Three completing methods were introduced in this paper. The data were completed from 1983 to 1988 at the KL and PTC stations. The completed data were called original data. The completing methods were the transfer function with fixed-parameters then built the Transfer Function (TF) model between the two stations. The TF model was used to estimate the wave data of the future. The difference of the estimated data and the original data were conferred in this paper. Two types of forecasting differences were discussed here: (1) the input data year increases from 1 to 6 years, and the difference of the estimated data and the original data were then compared. (2) Use one-year data as the input series to forecast the data of the next five years, and then compare the differences of the estimated data and the original data. After comparing the differences we can know the results. The mean value and the standard deviation of the differences decrease and go more stable as the quantity of the input data year increases. The mean value and the standard deviation of the differences increase and go more unstable as the quantity of the forecasted data year increases. This model can be used to forecast the wave data.     

基于核函数NPSSB模型的HY-2A卫星雷达高度计海况偏差研究

利用机器学习的方法,对14个周期HY-2A卫星高度计数据:风速、有效波高和海面高度差值进行训练,探究海况偏差和风速、有效波高之间的关系,创建海况偏差核函数非参数模型(NPSSB),并与参数模型中具有代表性的BM3、BM4模型进行对比。研究表明:(1)核函数NPSSB模型能够很好的反映SSB与U、SWH之间的关系,SSB与U呈二次函数关系,SSB与SWH呈反比例函数关系;(2)核函数NPSSB模型对SSB的模拟能力与训练数据集相关,数据量越多,模拟能力越好;(3)核函数NPSSB模型与BM3、BM4模型都存在0～-0.03 m的差值,随着风速和有效波高的增加,差值的绝对值越大。     

基于遥感数据的温度剖面预报研究

为了提高温度剖面的预报精度,提出了一种基于遥感数据的温度剖面预报方法.文中使用27个实测温度剖面和遥感SST、SSHA数据实现了对温度剖面的预报,并用该点的ARGO数据进行了检验.实验结果表明,将遥感数据同化到温度剖面的预报中是可行的,并能有效的提高温度剖面的预报精度.     

Mapping urban heat islands of arctic cities using combined data on field measurements and satellite images based on the example of the city of Apatity (Murmansk Oblast)

This article presents the results of a study of the urban heat island (UHI) in the city of Apatity during winter that were obtained according to the data of field meteorological measurements and satellite images. Calculations of the surface layer temperature have been made based on the surface temperature data obtained from satellite images. The experimental data on air temperature were obtained as a result of expeditionary meteorological observations, and the experimental data on surface temperature were obtained based on the data of the space hyperspectral Moderate-Resolution Imaging Spectroradiometer (MODIS) system, channels 31 and 32 (10.78–11.28 and 11.77–12.27 micrometers, respectively). As a result of the analysis of temperature fields, an intensive heat island (up to 3.2°C) has been identified that was estimated based on the underlying surface temperature, and its mean intensity over the observation period significantly exceeds the representative data for European cities in winter. It has also been established that the air temperature calculated according to the MODIS data is systematically higher under winter conditions than the air temperature from direct measurement data.     

健跳港上游引水对海床冲淤影响的数值模拟分析

健跳港是一个狭长的淤泥质港汊,针对枯季和洪季健跳港海床不同的冲淤特点,利用一维潮流泥沙模型复演了该港海床的冲淤过程,通过调整挟沙力系数来模拟健跳港枯水期的淤积和洪水期的冲刷这两种不同的海床变化形态,并预测了健跳港上游引水后对该港河床的淤积影响。     

Distribution of suspended particulate matter determined by in-situ observations and satellite images in the NW Aegean Sea (Greece)

Suspended particulate matter (SPM) patterns in the surface waters of the NW Aegean Sea were studied by (1) determining SPM concentration by water filtration, (2) measuring light transmission, and (3) evaluating satellite images. The SPM signals of the three major rivers discharging into the study area were recorded by all three methods, thereby providing information about the sources, transport pathways, and regional dispersion patterns of the SPM. The filtration of water samples and light transmission measurements were found to be good indicators of SPM concentrations in surface waters. Most of the SPM is composed of terrigenous minerals, thus explaining the correlation between the beam attenuation coefficient and the SPM concentration. A Landsat image obtained for the study period was found to adequately reveal regions with high SPM concentrations. Low concentrations, on the other hand, remain obscured. Received: 3 August 1999 / Revision accepted: 15 March 2000     

Characteristics of radiative heat transfer in the atmospheric surface layer from the results of direct measurements

Results of direct measurements of the long-wavelength (LW) radiative heat influx (RHI) in the atmospheric surface layer (ASL) are presented. These measurements were performed in August 2003 at the IAP RAS base in Tsimlyansk under the conditions of unstable and stable stratification during a weak wind and a cloudless sky and under nonsteady conditions during cumulus cloudiness in the daytime. The underlying surface was dry steppe with spars grass. The in situ RHI measurements were performed with an original optoacoustic receiver having a quasi-spherical angle of view at heights from 0.15 to 4 m. It is shown that the radiative heating in the ASL was many times the actual heating, especially during near-noon hours. In the daytime, the radiative heating attained its maximum at the heights of measurements 0.15–1 m and decreased with height. The radiative heating at these heights in the near-noon hours was on average about 20 K/h, attaining 60 K/h under a cloudless sky and a weak wind. Under inversion stratification, the radiative cooling usually exceeded the actual cooling, amounting on average from 0 to ?8 K/h and changing with height only slightly. Periods with close (in phase) fluctuations of the radiative and actual cooling, sometimes changing to heating, were observed during the night. Regression equations, showing a high correlation between the RHI values at the heights of measurements 0.5 and 1 m and the soil-air temperature differences at the height of measurements, are obtained for different heights. The diurnal mean RHI profiles are characterized by a heating on the order of several K/h in the lower part of the layer of measurements, which decreases with height and changes to cooling at heights of up to 4 m. A change in the effective radiation with height in the layer of measurements, which was obtained through the summation of RHI values at several heights, was significant, attaining on average ?25 W/m² in the near-noon hours and +10 W/m² in the evening hours. The nonradiative (turbulent) heat influx, obtained as the difference between the rates of actual and radiative temperature variations measured in situ, decreased the radiative heating in the daytime many times. The main sources of error in direct RHI measurements are estimated.     

Asymptotic estimates of the mean probability of radiative transfer through an exponentially correlated stochastic medium

Various models of exponentially correlated random fields associated with Poisson point ensembles, as well as algorithms for simulating radiative transfer in stochastic media of that type, are considered. Asymptotic estimates for the mean probability of particle (quantum of radiation) passage, assuming that the flow of trajectory intersections with domains of constant random density is Poissonian and using the central limit theorem for the corresponding optical length, are made.     

Convective structures in the Lofoten Basin based on satellite and Argo data

We discuss the possibility of detecting deep convection in the Lofoten Basin of the Norwegian Sea based on the eddy structures revealed from the satellite data. Satellite altimetry, SAR imagery, and MODIS satellite spectral radiometer sea-surface temperature (SST) data are used in the analysis, along with the data of oceanographic Argo floats. It is shown that the eddies identified from the satellite data correspond to the convective cells in the same region according to the data of the Argo floats. We consider several examples of the summer eddy and one winter eddy and the corresponding structures in the ocean measured by the Argo floats when they were located close to the identified eddies. As this method develops and improves, it can be used for the analysis of the dynamic of oceanic eddies in the region of the Lofoten Basin, and possibly in other regions with active deep convection.     

Spectral characteristics of Rossby waves in the Northwestern Pacific based on satellite altimetry

Using satellite altimetry measurement data for 1993–2013, we study the spectral characteristics of Rossby waves in the Northwestern Pacific (25°–50° N, 140°–180° E). For each latitude degree, we draw integral plots of spectral power density calculated with a two-dimensional Fourier transform (2D-FFT). We compare the dispersion equations of Rossby waves calculated from the WKB-approximation and an approximation of a two-layer ocean model with the empirical velocities determined by the slope of isopleths by the Radon method; also, we compare the dispersion equations with the spectral distributions of level variations. It is shown that the main energy of Rossby waves in the Northwestern Pacific corresponds to the first baroclinic mode. At almost all latitudes, there is good agreement between the empirical phase velocities calculated by isopleths by the Radon method and the theoretical values; also, the spectral peaks correspond to graphs of the dispersion equations for the first baroclinic mode Rossby waves, except for the Kuroshio region, where some peaks correspond to the second mode.     

Analysis of horizontal and vertical ring structure based on analytical model and satellite data: Application to the North Brazil Current Rings

A theoretical analysis of the principal oceanic ring parameters (tangential and angular velocities, shape function, horizontal and vertical scale of the ring) is described. The theoretical model consists of a reduced gravity model of the lens-like vortices with solid body rotation. The application of this approach is tested by comparison with data from the North Brazil Current Rings Experiment and remote sensing data. Specifically, we used the data corresponding to the surfaceintensified North Brazil Current ring (R-3) surveyed in February–March 1999, using direct velocity and hydrographical measurements. The theoretical model was used for evaluating the geometrical structure of the surface-intensified rings that produce remarkable signals in satellite data. The principal ring parameters from the model were compared with those from satellite data (altimetry and drifter information), which were estimated by using the method of minimization of the multivariable objective functions. Although the proposed model is linear in its conception, a good agreement was observed between the model and the primary characteristics of the observed rings. The model, however, allows for improvement in its assumptions, since its application is rather limited to intense ocean rings.     

Modeling radiative forcing by background aerosol on the basis of measurement data

The solar radiation flux incident on the land surface and aerosol radiative forcing are calculated from measurements carried out under clear skies during the summers of 2004 and 2005 at the Zvenigorod Scientific Station (ZSS) of the A.M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences. In these calculations, the following five models of aerosol optical parameters are used: a standard model of continental aerosol, a semiempirical model based on measurement data, and three models obtained from calculations on the basis of the Mie scattering theory for the mean size distribution of aerosol particles (which is retrieved from aureole measurements). A study of the sensitivity of the flux and forcing to the choice of aerosol model has shown that the relative error related to the model choice is not large (<5%) for the incident radiation flux and reaches 120% for aerosol forcing at the atmospheric upper boundary. The aerosol radiative forcing at the atmospheric upper boundary is also estimated and, according to our calculations, varies from −15 to −2 W/m². The use (in calculations) of the values of the albedo of single-scattering and the factor of scattering-indicatrix asymmetry at a wavelength of 550 nm, which were obtained for each of the observation periods, has made it possible to significantly refine the value of aerosol radiative forcing (compared to that calculated for the means of these parameters over all the periods). An even more accurate determination of aerosol radiative forcing is possible only when data on the size-distribution of aerosol particles and their chemical composition for a certain situation under analysis are available. As a result of this study, two models of the optical parameters of background aerosol have been proposed for use in radiation calculations: a semiempirical model and a mean model obtained from the calculations based on the theory of scattering. Both models use the values of the albedo of single scattering and the factor of asymmetry (at a wavelength of 550 nm) retrieved from the measurements for each of the observation periods.     

Surface currents around Hokkaido in the late fall of 1981 obtained from analysis of satellite images

This paper describes the characteristics of currents around Hokkaido using a current vector map compiled by chasing the displacement of sea marks on a pair of successive thermal infrared images taken from a satellite, NOAA-6. The points of some sea surface patterns showing distinctive features which can be commonly identified in both images are called Sea Marks≓. This sea mark chase method≓ has a great advantage over velocity measurements by boats or buoys, in that it gives a synoptic view of the velocity distribution over a broad sea area extending for some hundreds of miles on a short time scale of half a day.In order to investigate the current in the late fall of 1981, we used the data taken at 19:05 JST on 30 October and at 7:20 JST on 31 October. With these data taken with a 12 hr difference, the measurement accuracy of the speed of sea marks reached ±0.1 knot through geometrical correction. The velocity vectors of sea marks agreed with results of GEK measurements performed on those days, and also with the mean current pattern obtained in the past based on sea surface data.The most distinct features recognized were some cyclonic and anticyclonic eddies of the order of 100 km in diameter which dominated in the Kuril Basin of the Okhotsk Sea. They appeared clearly in the original NOAA images and their physical parameters were determined quantitatively from the vector map. These eddies were located in the region of decay of the Soya Warm Current and were rotating at a speed nearly as large as that of the current. The problems of what feeds energy to these eddies and how long they live remain to be solved in future studies.     

Retrieval of sea surface specific humidity based on AMSR-E satellite data

A new algorithm using a multivariate regression technique for retrieving sea surface specific humidity (Q) from remote sensing data from the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) is proposed. Daily and monthly specific humidity data from the National Center for Environmental Prediction (NCEP) reanalysis dataset and data of sea surface temperature, atmospheric total water vapor, and wind speed from AMSR-E oceanographic products were used to derive the regression coefficients of the algorithm, and all the data for derivation are from the year 2003. An F-test was applied to the regression, and small P-values indicate that the regressions are significant to a high level of confidence. The derived coefficients have been validated using similar data from the year 2004. The root mean square (rms) error of the algorithm for daily retrieved Q over the global oceans is 1.05 g kg⁻¹, and the rms error for monthly retrieved Q is 0.61 g kg⁻¹.     

Derivation of internal solitary wave amplitude in the South China Sea deep basin from satellite images

In the present study, theories based on the Korteweg–de Vries equation are extended to the Benjamin–Ono equation to allow the determination of internal solitary wave (ISW) amplitude from satellite images. The free surface flow induced by an ISW is derived for deep water. As a coherent structure, the amplitude of the ISW has a unique relation to the convergence/divergence of surface flow, such that the flow convergence/divergence will increase/decrease the backscattering cross section and generate bright/dark bands in satellite images. The distance between bright and dark bands can be related to the amplitude of ISW. To validate the theory, a multi-ship measurement made on 9–11 May 2005 during the spring tide period is used. A systematic approach to determine the thickness and density of the upper and lower layers is also included so that the free surface flow can be determined with a relatively high accuracy.     

Simulation of wave overtopping by an incompressible SPH model

The paper presents an incompressible Smoothed Particle Hydrodynamics (SPH) model to investigate the wave overtopping of coastal structures. The SPH method is a grid-less Lagrangian approach which is capable of tracking the large deformations of the free surface with good accuracy. The incompressible algorithm of the model is implemented by enforcing the constant particle density in the pressure projection. The SPH model is employed to reproduce a transient wave overtopping over a fixed horizontal deck and the regular/irregular waves overtopping of a sloping seawall. The computations are validated against the experimental and numerical data and a good agreement is observed. The SPH modelling is shown to provide a promising tool to predict the overtopping characteristics of different waves. The present model is expected to be of practical purpose if further improvement in the spatial resolution and CPU time can be adequately made.