Three-dimensional numerical simulation of wind-induced barotropic circulation in the Gulf of Patras

The barotropic, wind-induced circulation, which develops in the Gulf of Patras in Western Greece during the winter, is studied using three-dimensional numerical simulations. The simulations are performed using the numerical code MIKE 3 FM (HD). The Gulf's basin is bracketed between two sills, one on the west at the opening with the Ionian Sea and the other on the east at the Straits of Rio-Antirio at the opening with the Gulf of Corinth. The simulations show that the wind-induced flow creates strong currents near the coasts, which determine the sense of rotation of the gyres that develop in the Gulf. Strong currents are also created at the Rio-Antirio Straits. The wind-induced, barotropic currents do not seem to contribute to the direct replenishment of bottom waters, which recirculate between the two sills. Depending on the wind-speed forcing of the flow, the residence time of the waters in the Gulf of Patras is estimated to range from one week to one month.     

基于GIS的陕西省旬阳地区滑坡灾害危险性区划

西部地区是我国地质灾害的重灾区。随着西部大开发战略的实施，该地区即将开展大规模的基础建设、能源开发等。区域内的经济发展与地质灾害的矛盾将不可避免地暴露出来。为解决这一问题，论文选取中国滑坡重灾区的江汉流域开展灾害危险性区划应用研究。研究区选在旬阳地区的县城近郊，通过MAPGIS软件平台及其二次开发的滑坡灾害分析系统，采用规则网格单元划分方法，运用信息量模型对该区斜坡稳定性进行了．空间定量预测，并依信息量法的结果编制了该区的危险性预测分区图。为政府部门进行土地规划决策、避免在地质灾害易发区进行大规模土地开发和工程建设提供了科学依据。     

Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF

The dispersal of volcanic ash from the May 18, 1980 eruption of Mount St. Helens (MSH) has been simulated using the Lagrangian ash-tracking model PUFF. Previous applications of the model were limited to smaller, short-lived eruptions with ash dispersal occurring mainly within the troposphere. Two high-resolution atmospheric reanalysis datasets (ERA-40 and NCEP/NCAR-40) allowed MSH ash cloud dispersal to be simulated up to 30 km elevation. The 1980 eruption was divided into two distinct eruptive phases, (1) an initial, relatively short-lived blast/surge phase that injected ash up to 30 km and (2) a subsequent nine-hour plinian phase that maintained an average eruption column height of 16 km. Using PUFF, the two phases of the MSH eruption were modeled separately based on a range of individual input parameters and then combined to produce an integrated simulation of the entire eruption. The trajectory and areal extent of the modeled atmospheric ash cloud best match the actual distribution of MSH ash when input parameters are set to values inferred from satellite and radar data collected on May 18, 1980. The prevailing wind field exerts the strongest control on the advection and ultimate position of the modeled ash cloud, making the maximum column height and the vertical distribution of ash the most sensitive of the PUFF input parameters for this event. The results indicate that the PUFF model works well at simulating the dispersal of ash injected well into the lower stratosphere from a moderate, relatively long-lived eruption, such as MSH. However, attempts to use PUFF to recreate some granulometric aspects of the MSH fallout deposit, such as the maximum particle size as a function of distance from source, were not successful. PUFF consistently predicts much greater fallout distances for small ash particles (< 500 µm) than actually observed in the MSH deposit. The effective settling velocities used by the PUFF model appear to be too slow to accurately predict fallout distances of small ash particles. As a consequence the PUFF model may overestimate the duration of ash loading in the atmosphere associated with the distal fine ash component of explosive eruptions.     

Direct Climatic Effect of Dust Aerosol in the NCAR Community Atmosphere Model Version 3 （CAM3）

Direct climate responses to dust shortwave and longwave radiative forcing (RF) are studied using the NCAR Community Atmosphere Model Version 3 (CAM3). The simulated RF at the top of the atmosphere (TOA) is-0.45 W m-2 in the solar spectrum and +0.09 W m-2 in the thermal spectrum on a global average. The magnitude of surface RF is larger than the TOA forcing, with global mean shortwave forcing of-1.76 W m-2 and longwave forcing of +0.31 W m-2 . As a result, dust aerosol causes the absorption of 1.1 W m-2 in t...     

Flood frequency analysis of southern Africa: II. Identification of regional distributions

Abstract

Regional frequency analysis of annual maximum flood data comprising 407 stations from 11 countries of southern Africa is presented. Forty-one homogeneous regions are identified. The L-moments of the observed data indicate that the possible underlying frequency distributions are Pearson type 3 (P3), lognormal 3-parameter (LN3), General Pareto (GPA) or General Extreme Value (GEV). Simulation experiments for the selection of the most suitable flood frequency procedure indicate that Pearson type 3/Probability Weighted Moments (P3/PWM) and log-Pearson type 3/Method of Moments (LP3/MOM) are suitable procedures for the region.     

A three-dimensional finite-element model of wind effects upon higher harmonics of the internal tide

A non-linear three-dimensional unstructured grid model of the M₂ tide in the shelf edge area off the west coast of Scotland is used to examine the spatial distribution of the M₂ internal tide and its higher harmonics in the region. In addition, the spatial variability of the tidally induced turbulent kinetic energy and associated mixing in the area are considered. Initial calculations involve only tidal forcing, although subsequent calculations are performed with up-welling and down-welling favourable winds to examine how these influence the tidal distribution (particularly the higher harmonics) and mixing in the region. Both short- and long-duration winds are used in these calculations. Tidal calculations show that there is significant small-scale spatial variability particularly in the higher harmonics of the internal tide in the region. In addition, turbulence energy and mixing exhibit appreciable spatial variability in regions of rapidly changing topography, with increased mixing occurring above seamounts. Wind effects significantly change the distribution of the M₂ internal tide and its higher harmonics, with appreciable differences found between up- and down-welling winds and long- and short-duration winds because of differences in mixing and the presence of wind-induced flows. The implications for model validation, particularly in terms of energy transfer to higher harmonics, and mixing are briefly discussed.     

Broadband synthetic aperture borehole radar interferometry

Trials in mines have established that wideband VHF borehole radars (BHR), working in the 10–100-MHz band, can be used to probe the rockmass between boreholes over ranges from <5 m to as much as 150 m with submeter resolution. There is evidence that ore bodies reflect these radar signals both specularly and diffusely, much as the ground/air interface does when overflown by synthetic aperture radar (SAR). In both SAR and BHR, multiple flight lines, together with diffuse reflections admit the possibility of developing interferometric 3D images of the object. This paper examines the possibility of imaging buried objects in three dimensions by interferometrically combining broadband VHF borehole radar profiles shot in adjacent pairs of boreholes. Broadbanding in BHR has the advantage of releasing the image from 2nπ phase ambiguities, but practically, interferometric borehole radar (InBHR) needs high signal-to-noise ratios (SNR) to avoid noise capture. This means that 3D InBHR is limited to ranges in wavelengths which are less than the rock's attenuation factor Q. Interferometric methods are developed which are capable of mapping ore bodies and other structures in three dimensions. Tangent plane migration methods are developed here in order to reconstruct surfaces that lie in the near-field of sparse interferometric arrays.     

The Interannual Variability of Climate in a Coupled Ocean-Atmosphere Model

In this paper, the interannual variability simulated by the coupled ocean-atmosphere general circulation model of the Institute of Atmospheric Physics (IAP CGCM) in 40 year integrations is analyzed, and compared with that by the corresponding IAP AGCM which uses the climatic sea surface temperature as the boundary condition in 25 year integrations.The mean climatic states of January and July simulated by IAP CGCM are in good agreement with that by IAP AGCM, i.e., no serious ‘climate drift’ occurs in the CGCM simulation. A comparison of the results from AGCM and CGCM indicates that the standard deviation of the monthly averaged sea level pressure simulated by IAP CGCM is much greater than that by IAP AGCM in tropical region. In addition, both Southern Oscillation (SO) and North Atlantic Oscillation (NAO) can be found in the CGCM simulation for January, but these two oscillations do not exist in the AGCM simulation.The interannual variability of climate may be classified into two types: one is the variation of the annual mean, another is the variation of the annual amplitude. The ocean-atmosphere interaction mainly increases the first type of variability. By means of the rotated EOF, the most important patterns corresponding to the two types of interannual variability are found to have different spatial and temporal characteristics.     

附加系统参数平差法在导线网中的应用

本文探讨了附加系统参数平差法在带有测距系统误差的导线网中的应用，说明了当系统误差比较显著时．采用附加系统参数平差法的可行性、必要性及其优越性。     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。