机动式激光测风雷达三维风场反演模型仿真研究

激光测风雷达具有探测精度高、环境适应性好等优点,发展其机动式应用技术,在军事保障、抢险救灾、海洋调查等领域具有广阔应用前景.雷达在机动条件下工作时,会由于载体的运动及倾斜造成数据探测误差,且不同光束之间的实际探测仰角不一致,导致无法基于传统VAD算法反演大气三维风场.本文研究探索了激光测风雷达机动探测条件下的数据订正算法,及同一组VAD扫描中光束仰角不同条件下的三维风场反演算法,并通过数值仿真试验对算法进行验证.结果表明,数据订正算法能够有效订正载体运动带来的雷达探测误差,数据反演算法能够有效反演得到三维风场,但同一组VAD扫描中光束仰角差越大,三维风场反演误差越大.

     

Regional wind fields and dunefield migration,southern Brazil

This paper examines a seemingly anomalous situation in southern Brazil where the dunefields on Santa Catarina Island (e.g. Joaquina Beach) migrate to the NNW, almost completely the opposite direction (c. 160) to the dunefields immediately to the south (e.g. Pinheira Beach), and some much further to the north (e.g. Cabo Frio) which migrate to the SSW. A variety of mechanisms are examined to explain the differences in dunefield migration including grain size variations, topographic effects on local winds, shoreline orientation, and regional wind field changes. The mean grain sizes of the two beaches, Pinheira and Joaquina, are not sufficiently different to restrict aeolian sediment transport in either place, nor to account for a lack of transport from the NNE to the SSW in the case of Joaquina. Some topographic steering of the wind is likely but could not account for the long‐term average difference in migration trends of the island dunefields compared to the mainland dunefields. While the orientation of the shoreline to prevailing winds is an important control on beach and dune sediment transport, it is not the dominant controlling mechanism. An analysis of the regional wind patterns demonstrates that there is a major shift in the regional wind field near the island such that the dominant island winds blow from the SW/SSW while those further south blow from the NE. It is concluded that this is the predominant reason for the divergence in the direction of migration of the dunefields. Copyright © 2006 John Wiley & Sons, Ltd.     

Multivariate stationary non-Gaussian process simulation for wind pressure fields

Stochastic simulation is an important means of acquiring fluctuating wind pressures for wind induced response analyses in structural engineering. The wind pressure acting on a large-span space structure can be characterized as a stationary non-Gaussian field. This paper reviews several simulation algorithms related to the Spectral Representation Method (SRM) and the Static Transformation Method (STM). Polynomial and Exponential transformation functions (PSTM and ESTM) are discussed. Deficiencies in current algorithms, with respect to accuracy, stability and efficiency, are analyzed, and the algorithms are improved for better practical application. In order to verify the improved algorithm, wind pressure fields on a large-span roof are simulated and compared with wind tunnel data. The simulation results fit well with the wind tunnel data, and the algorithm accuracy, stability and efficiency are shown to be better than those of current algorithms.     

典型风场的分形特征及分形模拟

本文对于由风洞试验得出的不同风场类别的脉动风速进行了盒维数的分析,发现各类风场的盒维数明显不同,而采用传统数值模拟方法模拟的A、B、C三类风场的脉动风速的盒维数几乎没有变化。为了弥补脉动风速传统数值模拟方法的缺陷,采用包含分形维数的Weierstrass-Mandelbro函数模拟A、B、C三类风场的脉动风速,并与由风洞试验得到的各类风场的脉动风速的盒维数,湍流度,功率谱等特征进行对比。结果表明分形模拟方法所得的脉动风速的盒维数在各类风场中明显不同,且随着湍流度的增大而减小,并且分形模拟方法模拟的脉动风速谱和风洞试验的脉动风速谱更加接近。     

Dam breach simulations with multiple breach locations under wind and wave actions

Multiple breaches of a dam resulting from wind-generated waves and wave overtopping are studied for a hypothetical long non-cohesive earthen dam with an uneven crest. Both wind speed and direction affect breach locations and outflow for a particular reservoir surface geometry. Locations on the dam with longer fetches along the wind direction are more subject to wave overtopping and breaching than other locations. Higher wind speeds lead to wave overtopping and dam breaches under larger freeboards than lower wind speeds. For a specified inflow hydrograph and spillway configuration, there exists a location at which the smallest estimated peak outflow occurs among all possible breach locations and the pool drops too quickly for additional breaches to develop. Using this location for a fuse plug or a pilot channel could minimize downstream impact, perhaps as an interim or emergency measure for a dam with inadequate spillway capacity.     

淀山湖风浪场的数值模拟

对淀山湖的风场和风浪场进行了观测,根据所观测的风场利用SWAN模型计算了湖区的风浪场.通过将计算的波高和观测值进行比较,说明计算结果的变化趋势和观测结果的变化趋势相吻合;观测结果的波动幅度较大,而计算结果比较平滑.总体说来,基于SWAN模型所计算的波高可信、可靠.在此基础上,利用SWAN模型较为系统地计算了不同水位和不...     

Estimation of tropical cyclone parameters and wind fields from SAR images

The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.     

Ocean wind fields retrieved from the advanced synthetic aperture radar aboard ENVISAT

In this paper an algorithm is presented which enables high-resolution ocean surface wind fields to be retrieved from the advanced synthetic aperture radar (ASAR) data acquired by the European remote sensing satellite ENVISAT. Wind directions are extracted from wind-induced streaks that are visible in ASAR images at scales above 200 m and that are approximately in line with the mean surface wind direction. Wind speeds are derived from the normalized radar cross section (NRCS) and image geometry of the calibrated ASAR images, together with the local ASAR-retrieved wind direction. Therefore the empirical C-band model CMOD4, which describes the dependency of the NRCS on wind and image geometry, is used. CMOD4 is a semi-empirical model, which was originally developed for the scatterometer of the European remote sensing satellites ERS-1 and 2 operating at C-band with vertical polarization. Consequently, CMOD4 requires modification when applied to ASAR images that were acquired with horizontal polarization in transmitting and receiving. This is performed by considering the polarization ratio of the NRCS. To demonstrate the applicability of the algorithm, wind fields were computed from several ENVISAT ASAR images of the North Sea and compared to atmospheric model results of the German weather service.Acknowledgements The authors were supported by the German Bundesministerium für Bildung und Forschung (BMBF) in the framework of the project. A new perspective of the Ocean ENVISAT Oceanography (ENVOC). The ENVISAT ASAR data were kindly made available by the European Space Agency in the framework of the ENVISAT Project AO-ID 220, Biological and geophysical parameters from synthetic aperture radar over the ocean (BIGPASO).     

A stochastic model reduction method for nonlinear unconfined flow with multiple random input fields

In this paper we present a stochastic model reduction method for efficiently solving nonlinear unconfined flow problems in heterogeneous random porous media. The input random fields of flow model are parameterized in a stochastic space for simulation. This often results in high stochastic dimensionality due to small correlation length of the covariance functions of the input fields. To efficiently treat the high-dimensional stochastic problem, we extend a recently proposed hybrid high-dimensional model representation (HDMR) technique to high-dimensional problems with multiple random input fields and integrate it with a sparse grid stochastic collocation method (SGSCM). Hybrid HDMR can decompose the high-dimensional model into a moderate M-dimensional model and a few one-dimensional models. The moderate dimensional model only depends on the most M important random dimensions, which are identified from the full stochastic space by sensitivity analysis. To extend the hybrid HDMR, we consider two different criteria for sensitivity test. Each of the derived low-dimensional stochastic models is solved by the SGSCM. This leads to a set of uncoupled deterministic problems at the collocation points, which can be solved by a deterministic solver. To demonstrate the efficiency and accuracy of the proposed method, a few numerical experiments are carried out for the unconfined flow problems in heterogeneous porous media with different correlation lengths. The results show that a good trade-off between computational complexity and approximation accuracy can be achieved for stochastic unconfined flow problems by selecting a suitable number of the most important dimensions in the M-dimensional model of hybrid HDMR.     

Simulation of the effect of wind speedup in the formation of transverse dune fields

A computer simulation model for transverse‐dune‐field dynamics, corresponding to a uni‐directional wind regime, is developed. In a previous formulation, two distinct problems were found regarding the cross‐sectional dune shape, namely the erosion in the lee of dunes and the steepness of the windward slopes. The first problem is solved by introducing no erosion in shadow zones. The second issue is overcome by introducing a wind speedup (shear velocity increase) factor, which can be accounted for by adding a term to the original transport length, which is proportional to the surface height. By incorporating these features we are able to model dunes whose individual shape and collective patterns are similar to those observed in nature. Moreover we show how the introduction of a non‐linear shear‐velocity‐increase term leads to the reduction of dune height, and this may result in an equilibrium dune field configuration. This is thought to be because the non‐linear increase of the transport length makes the sand trapping efficiency lower than unity, even for higher dunes, so that the incoming and the outgoing sand flux are in balance. To fully describe the inter‐dune morphology more precise dynamics in the lee of the dune must be incorporated. Copyright © 2000 John Wiley & Sons, Ltd.     

A stochastic approach to nonlinear unconfined flow subject to multiple random fields

In this study, the KLME approach, a moment-equation approach based on the Karhunen–Loeve decomposition developed by Zhang and Lu (Comput Phys 194(2):773–794, 2004), is applied to unconfined flow with multiple random inputs. The log-transformed hydraulic conductivity F, the recharge R, the Dirichlet boundary condition H, and the Neumann boundary condition Q are assumed to be Gaussian random fields with known means and covariance functions. The F, R, H and Q are first decomposed into finite series in terms of Gaussian standard random variables by the Karhunen–Loeve expansion. The hydraulic head h is then represented by a perturbation expansion, and each term in the perturbation expansion is written as the products of unknown coefficients and Gaussian standard random variables obtained from the Karhunen–Loeve expansions. A series of deterministic partial differential equations are derived from the stochastic partial differential equations. The resulting equations for uncorrelated and perfectly correlated cases are developed. The equations can be solved sequentially from low to high order by the finite element method. We examine the accuracy of the KLME approach for the groundwater flow subject to uncorrelated or perfectly correlated random inputs and study the capability of the KLME method for predicting the head variance in the presence of various spatially variable parameters. It is shown that the proposed numerical model gives accurate results at a much smaller computational cost than the Monte Carlo simulation.     

Turbulence fields in the lee of two‐dimensional transverse dunes simulated in a wind tunnel

The turbulence field of airflow in the lee of a dune has significant impacts on dune dynamics and related processes. We used particle image velocimetry in a wind tunnel simulation to obtain detailed velocity measurements in the lee of two‐dimensional transverse dune models, then used the results to analyse their turbulence fields. The dune models used in this study had a single lee angle of 30°, and a total of six stoss angles: 3°, 5°, 10°, 15°, 20° and 25°. We used vorticity, turbulence intensity, Reynolds stress and turbulent kinetic energy to characterize the turbulence fields. These parameters were functions of stoss angle, wind velocity, distance from the dune crest and height above the ground surface. The stoss angles could generally be divided into two groups based on the profiles of mean velocity, turbulence and Reynolds stress. Stoss angles of 3° and 5° usually had similar profiles, and angles of 15°, 20° and 25° formed a second group with similar profiles. The profiles for the stoss angle of 10° were usually transitional and were intermediate between the two groups. Vorticity, Reynolds stress and turbulent kinetic energy increased monotonically with increasing free‐stream wind velocity, but their variations with respect to the stoss angle were complex. The stoss angles of 15° and 20° had the maximum values of these three parameters, thus these angles may have special significance in dune development given the characteristics of the mean velocity fields and turbulence fields they produce within the lee airflow. It is the streamwise velocity component and its turbulence that determine the surface shear stress. Copyright © 2008 John Wiley and Sons, Ltd.     

High-order kernels for Riemannian wavefield extrapolation

Riemannian wavefield extrapolation is a technique for one‐way extrapolation of acoustic waves. Riemannian wavefield extrapolation generalizes wavefield extrapolation by downward continuation by considering coordinate systems different from conventional Cartesian ones. Coordinate systems can conform with the extrapolated wavefield, with the velocity model or with the acquisition geometry. When coordinate systems conform with the propagated wavefield, extrapolation can be done accurately using low‐order kernels. However, in complex media or in cases where the coordinate systems do not conform with the propagating wavefields, low order kernels are not accurate enough and need to be replaced by more accurate, higher‐order kernels. Since Riemannian wavefield extrapolation is based on factorization of an acoustic wave‐equation, higher‐order kernels can be constructed using methods analogous to the one employed for factorization of the acoustic wave‐equation in Cartesian coordinates. Thus, we can construct space‐domain finite‐differences as well as mixed‐domain techniques for extrapolation. High‐order Riemannian wavefield extrapolation kernels improve the accuracy of extrapolation, particularly when the Riemannian coordinate systems does not closely match the general direction of wave propagation.     

Simulated precipitation fields with variance-consistent interpolation

Abstract

Gridded meteorological data are available for all of Norway as time series dating from 1961. A new way of interpolating precipitation in space from observed values is proposed. Based on the criteria that interpolated precipitation fields in space should be consistent with observed spatial statistics, such as spatial mean, variance and intermittency, spatial fields of precipitation are simulated from a gamma distribution with parameters determined from observed data, adjusted for intermittency. The simulated data are distributed in space, using the spatial pattern derived from kriging. The proposed method is compared to indicator kriging and to the current methodology used for producing gridded precipitation data. Cross-validation gave similar results for the three methods with respect to RMSE, temporal mean and standard deviation, whereas a comparison on estimated spatial variance showed that the new method has a near perfect agreement with observations. Indicator kriging underestimated the spatial variance by 60–80% and the current method produced a significant scatter in its estimates.

Citation Skaugen, T. & Andersen, J. (2010) Simulated precipitation fields with variance-consistent interpolation. Hydrol. Sci. J. 55(5), 676–686.     

The influence of local topography for wind direction on Mars: two examples of dune fields in crater basins

In this work, we perform an analysis of large dark dunes within Moreux Crater and Herschel Crater on Mars using High Resolution Imaging Science Experiment (HiRISE) and Context Camera (CTX) data sets. These data allow us to conduct a detailed analysis of dune morphology and slip faces, concluding that the studied dune fields are influenced by topographically‐controlled complex wind directions. Our morphological analysis reveals that inside Moreux Crater in particular, the topographic setting dominates the wind flow direction, leading to the development of a sand transport pathway encircling the central peak of the crater. The dune fields in Herschel Crater are also affected by winds controlled by variable topography as suggested by the presence of complex dunes and dune fields. Our analysis indicate that the studied dune systems is not the result of paleo‐wind regimes. Furthermore, we perform thermal inertia measurements using thermal emission spectrometer (TES) data, which indicate that the studied dune fields consist of medium sand 250–500 µm in diameter. Copyright © 2012 John Wiley & Sons, Ltd.     

Computing three-dimensional gravitational fields with equivalent sources

Existing techniques for computing the gravitational field due to a homogeneous polyhedron all transform the required volume integral, expressing the field due to a volume distribution of mass, into a surface integral, expressing the potential due to a surface mass distribution over the boundary of the source body. An alternative representation is also possible and results in a surface integral expressing the potential due to a variable-strength double layer located on the polyhedral source boundary. Manipulation of this integral ultimately allows the gravitational field component in an arbitrary direction to be expressed as a weighted sum of the potentials due to two basic source distributions. These are a uniform-strength double layer located on all faces and a uniform-strength line source located along all edges. The derivatives of the gravitational field components can also be expressed in a similar form as can the magnetic field components due to a homogeneous magnetic polyhedron. It follows that the present approach can be used to generate a universal program capable of modelling all the commonly used potential field responses due to 3D bodies of arbitrary shape.     

Convective dynamos with intermediate and strong fields

Abstract

The weak-field Benard-type dynamo treated by Soward is considered here at higher levels of the induced magnetic field. Two sources of instability are found to occur in the intermediate field regime M ～ T ^1/12, where M and T are the Hartmann and Taylor numbers. On the time scale of magnetic diffusion, solutions may blow up in finite time owing to destabilization of the convection by the magnetic field. On a faster time scale a dynamic instability related to MAC-wave instability can also occur. It is therefore concluded that the asymptotic structure of this dynamo is unstable to virtual increases in the magnetic field energy.

In an attempt to model stabilization of the dynamo in a strong-field regime we consider two approximations. In the first, a truncated expansion in three-dimensional plane waves is studied numerically. A second approach utilizes an ad hoc set of ordinary differential equations which contains many of the features of convection dynamos at all field energies. Both of these models exhibit temporal intermittency of the dynamo effect.     

Second-order smoothing of spatial point patterns with small sample sizes: a family of kernels

We consider kernel-based non-parametric estimation of second-order product densities of spatial point patterns. We present a new family of optimal and positive kernels that shows less variance and more flexibility than optimal kernels. This family generalises most of the classical and widely used kernel functions, such as Box or Epanechnikov kernels. We propose an alternative asymptotically unbiased estimator for the second-order product density function, and compare the performance of the estimator for several members of the family of optimal and positive kernels through MISE and relative efficiency. We present a simulation study to analyse the behaviour of such kernel functions, for three different spatial structures, for which we know the exact analytical form of the product density, and under small sample sizes. Some known datasets are revisited, and we also analyse the IMD dataset in the Rhineland Regional Council in Germany.     

The effect of wind averaging time on wind erosivity estimation

The Wind Erosion Prediction System (WEPS) and Revised Wind Erosion Equation (RWEQ) are widely used for estimating wind‐induced soil erosion at a field scale. Wind is the principal erosion driver in the two models. Wind erosivity, which describes the capacity of wind to cause soil erosion, is defined as erosive wind power density (WPD) in WEPS, and wind value (W) in RWEQ. In this study, the daily average WPD (AWPD) and the daily average W (Wf) were chosen to investigate the effect of averaging time on wind erosivity estimation based on observed wind data. We compare the daily AWPD and Wf calculated from 1, 5, 10, 15, 30, and 60 minute average wind speed data. The results of comparisons indicate that averaging wind speed can significantly influence estimates of wind erosivity. Compared with the daily AWPD and Wf calculated from one minute average wind speed data, all daily AWPD and Wf values calculated from 5, 10, 15, 30, and 60 minute averaged wind speeds tend to be significantly lower than values calculated from one minute values. In general, longer averaging times tend to produce smaller values of daily AWPD or Wf, which may lead to an under‐estimation of wind erosion. Further studies are needed to extend and apply the findings obtained in this study to actual wind erosion predictions. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluating effectiveness of multiple tuned mass dampers for vibration control of jacket offshore wind turbines under onshore and seafloor earthquakes

The dynamic characteristics and structural responses of operation and grid loss offshore wind turbines(OWTs) under onshore and seafloor earthquakes are analyzed based on the established coupled seismic analysis model. In addition to the remarkable influence of the rotor system on the responses of the operation OWT under earthquakes, interactions among the natural modes of the grid loss OWT in the fore-aft and side-to-side directions are revealed. By comparing with the onshore earthquakes, the mo...