Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free-falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of I m x I m, the expected number of received free-failing raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re-detachment amount. The re-detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free-falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re-detachment amount were small parts of the total splash amount. Their proportions were 0.15% and 2.6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil-splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright (c) 2007 John Wiley & Sons, Ltd.

P wave velocity of the pumice from the Okinawa Trough at high temperature and high pressure.

P wave velocity of the pumice sample from the middle Okinawa Trough and andesite sample from vicinity Yingdao volcanic island, Kyushu Japan were measured at temperature (from room temperature to 1500 C) and pressure (from room pressure to 2.4GPa) using a multi-anvil pressure apparatus called the YJ-3000 press. The measured data shows that at low temperature and low pressure (<1GPa, <800degreesC), the P wave velocity of pumice is lower than that of andesite, while at high temperature and high pressure (>1GPa, >800degreesC) the P wave velocity of pumice and andesite. becomes consistent (5.9km/s). The paper points out that 1GPa/800degreesC is the point of thermodynamic phase transformation Okinawa Trough pumice and vicinity andesite, and the point is deeper than 18km.

On SAR wind speed ambiguities and related geophysical model functions

Recent investigations show that normalized radar cross sections for C-band microwave sensors decrease under high wind conditions with certain incident angles instead of increase, as is the case for low to moderate wind speeds. This creates the problem of ambiguities in high wind speed retrievals from synthetic aperture radar (SAR). In the present work, four geophysical model functions (GMFs) are studied, namely the high wind C-band model 4 (CMOD4HW), C-band model 5 (CMOD5), the high wind vertical polarized GMF (HWGMF_VV), and the high wind horizontal polarized GMF (HWGMF_HH). Our focus is on model behaviours relative to wind speed ambiguities. We show that, except for CMOD4HW, the other GMFs exhibit the wind speed ambiguity problem. To consider this problem in high wind speed retrievals from SAR, we focus on hurricanes and propose a method to remove the speed ambiguity using the dominant hurricane wind structure.

Speed ambiguity in hurricane wind retrieval from SAR imagery

Under strong ocean surface wind conditions, the normalized radar cross section of synthetic aperture radar (SAR) is dampened at certain incident angles, compared with the signals under moderate winds. This causes a wind speed ambiguity problem in wind speed retrievals from SAR, because two solutions may exist for each backscattered signal. This study shows that the problem is ubiquitous in the images acquired by operational space-borne SAR sensors. Moreover, the problem is more severe for the near range and range travelling winds. To remove this ambiguity, a method was developed based on characteristics of the hurricane wind structure. A SAR image of Hurricane Rita (2005) was analysed to demonstrate the wind speed ambiguity problem and the method to improve the wind speed retrievals. Our conclusions suggest that a speed ambiguity removal algorithm must be used for wind retrievals from SAR in intense storms and hurricanes.

A new wind vector algorithm for C-band SAR

Ocean wind speed and wind direction are estimated simultaneously using the normalized radar cross sections or' corresponding to two neighboring (25-km) blocks, within a given synthetic aperture radar (SAR) image, having slightly different incidence angles. This method is motivated by the methodology used for scatterometer data. The wind direction ambiguity is removed by using the direction closest to that given by a buoy or some other source of information. We demonstrate this method with 11 EN-VISAT Advanced SAR sensor images of the Gulf of Mexico and coastal waters of the North Atlantic. Estimated wind vectors are compared with wind measurements from buoys and scatterometer data. We show that this method can surpass other methods in some cases, even those with insufficient visible wind-induced streaks in the SAR images, to extract wind vectors.

Improvement of different source function expressions in SWAN model for the Bohai Sea

The simulating wave nearshore (SWAN) wave model has been widely used in coastal areas, lakes and estuaries. However, we found a poor agreement between modeling results and measurements in analyzing the chosen four typical cases when we used the default parameters of the source function formulas of the SWAN to make wave simulation for the Bohai Sea. Also, it was found that at the same wind process the simulated results of two wind generation expressions (Komen, Janssen) demonstrated a large difference. Further study showed that the proportionality coefficient alpha in linear growth term of wave growth source function plays an unperceived role in the process of wave development. Based on experiments and analysis, we thought that the coefficient alpha should change rather than be a constant. Therefore, the coefficient alpha changing with the variation of friction velocity U (*) was introduced into the linear growth term of wave growth source function. Four weather processes were adopted to validate the improvement in the linear growth term. The results from the improved coefficient alpha agree much better with the measurements than those from the default constant coefficient alpha. Furthermore, the large differences of results between Komen wind generation expression and Janssen wind generation expression were eliminated. We also experimented with the four weather processes to test the new white-capping mechanisms based on the cumulative steepness method. It was found that the parameters of the new white-capping mechanisms are not suitable for the Bohai Sea, but Alkyon's white-capping mechanisms can be applicable to the Bohai Sea after amendments, demonstrating that this improvement of parameter alpha can improve the simulated results of the Bohai Sea.

SUSPENDED SEDIMENT TRANSPORT IN THE OFFSHORE NEAR YANGTZE ESTUARY

Based on the Estuarine, Coastal and Ocean Modeling System with Sediments (ECOMSED) model, a 3-D hydrodynamic-transport numerical model was established for the offshore area near the Yangtze Estuary in the East China Sea. The hydrodynamic module was driven by tide and wind. Sediment module included sediment resuspension, transport and deposition of cohesive and non-cohesive sediment. The settling of cohesive sediment in the water column was modeled as a function of aggregation (flocculation) and deposition. The numerical results were compared with observation data for August, 2006. It shows that the sediment concentration reduces gradually from the seashore to the offshore area. Numerical results of concentration time series in the observation stations show two peaks and two valleys, according with the observation data. It is mainly affected by tidal current. The suspended sediment concentration is related to the tidal current during a tidal cycle, and the maximum concentration appears 1 h-4 h after the current maximum velocity has reached.

Comparison of nonlinear and linear PCA on surface wind, surface height, and SST in the South China Sea

We compared nonlinear principal component analysis (NLPCA) with linear principal component analysis (LPCA) with the data of sea surface wind anomalies (SWA), surface height anomalies (SSHA), and sea surface temperature anomalies (SSTA), taken in the South China Sea (SCS) between 1993 and 2003. The SCS monthly data for SWA, SSHA and SSTA (i.e., the anomalies with climatological seasonal cycle removed) were pre-filtered by LPCA, with only three leading modes retained. The first three modes of SWA, SSHA, and SSTA of LPCA explained 86%, 71%, and 94% of the total variance in the original data, respectively. Thus, the three associated time coefficient functions (TCFs) were used as the input data for NLPCA network. The NLPCA was made based on feed-forward neural network models. Compared with classical linear PCA, the first NLPCA mode could explain more variance than linear PCA for the above data. The nonlinearity of SWA and SSHA were stronger in most areas of the SCS. The first mode of the NLPCA on the SWA and SSHA accounted for 67.26% of the variance versus 54.7%, and 60.24% versus 50.43%, respectively for the first LPCA mode. Conversely, the nonlinear SSTA, localized in the northern SCS and southern continental shelf region, resulted in little improvement in the explanation of the variance for the first NLPCA.

Winter counter-wind currents off the southeastern China coast: A review

This review covers the discovery and studies of the year-round northeastward currents off the southeastern China coast, paying special attention to its upwind characteristic in winter, mainly focusing on work by Chinese oceanographers. This current system is a prominent and unique phenomenon in the shelf circulation of the world ocean. The general features of the current system are summarized. The evidence for the existence and the variation of the three parts of the currents-the South China Sea Warm Current, the Taiwan Strait Warm Current and the Taiwan Warm Current-are separately elucidated. The formation mechanisms of the current as a whole are explained using dynamic analysis and numerical simulation results. Some suggestions for further studies are also made.

Probability distribution of random wave-current forces

Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave-wave and wave-current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan-Pierson-Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. (c) 2006 Elsevier Ltd. All rights reserved.

Wind-vector estimation for RADARSAT-1 SAR images: Validation of wind-direction estimates based upon geometry diversity

In this letter, a new wind-vector algorithm is presented that uses radar backscatter sigma(0) measurements at two adjacent subscenes of RADARSAT-1 synthetic aperture radar (SAR) images, with each subscene having slightly different geometry. Resultant wind vectors are validated using in situ buoy measurements and compared with wind vectors determined from a hybrid wind-retrieval model using wind directions determined by spectral analysis of wind-induced image streaks and observed by colocated QuikSCAT measurements. The hybrid wind-retrieval model consists of CMOD-IFR2 [applicable to C-band vertical-vertical (W) polarization] and a C-band copolarization ratio according to Kirchhoff scattering. The new algorithm displays improved skill in wind-vector estimation for RADARSAT-1 SAR data when compared to conventional wind-retrieval methodology. In addition, unlike conventional methods, the present method is applicable to RADARSAT-1 images both with and without visible streaks. However, this method requires ancillary data such as buoy measurements to resolve the ambiguity in retrieved wind direction.

Numerical study on the velocity structure around tidal fronts in the Yellow Sea

The velocity components across tidal fronts are examined using the Blumberg and Mellor 3-D nonlinear numerical coastal circulation model incorporated with the Mellor and Yamada level 2.5 turbulent closure model based on the reasonable model output of the M-2 tide and density residual currents. In the numerical experiments, upwelling motion appears around all the fronts with different velocity structures, accounting for surface cold water around the fronts. The experiments also suggest that the location and formation of fronts are closely related to topography and tidal mixing, as is the velocity structure around the front.

Restructuring of a zooplankton community by perturbation from a wind-forced coastal jet

The impact of transient wind events on an established zooplankton community was observed during a, field survey in a, coastal region off northern Norway in May 2002. A transient wind event induced a coastal jet/filament intrusion of warm, saline water into our survey area where a semi-permanent eddy was present. There was an abrupt change in zooplankton community structure within 4-7 days of the wind event, with a change in the size structure, an increase in lower size classes less than 1 mm in equivalent spherical diameter (ESD) and a decrease in larger size classes greater than 1.5 mm in ESD. The slope of zooplankton biovolume spectra changed from -0.6 to -0.8, consistent with the size shifting towards smaller size classes. This study shows that even well established zooplankton communities are susceptible to restructuring during transient wind events, and in particular when wind forcing induces horizontal currents or filaments.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.