Unsteadiness as a cause of non-equality of eddy diffusivities for heat and vapour at the base of an advective inversion

This paper examines the effect of non-stationarity of the wind on similarity of the eddy diffusivities for heat and vapour within a stable layer at the bottom of an internal boundary layer formed downwind of a dry-to-wet transition. First, we present some experimental data taken above a rice crop downwind of very extensive dry range lands at Warrawidgee, NSW, Australia. These data establish that periods of higher wind speed were associated with periods of higher saturation deficit in the canopy of the rice crop, and lower Bowen ratio. It is shown that Bowen ratios calculated for 30-second sub-intervals varied three-fold within a single 20-minute averaging period. Thus periods of higher wind speed corresponded to periods of higher moisture flux and smaller sensible heat flux.An idealized situation is then analysed theoretically. It is assumed that the time scale of the slow variations of the wind is long compared with the surface-layer time scale and that fetch is sufficient that the air near the ground is in continuous equilibrium with the surface. Using a two-scale Reynolds decomposition of the fluctuating wind and scalar variables into active and inactive components, it is shown that unsteadiness can lead to an eddy diffusivity for saturation deficit, calculated as the ratio of average flux to average gradient, that is larger than that for total energy calculated in a similar way. Using this ratio to calculate the ratio of diffusivities for temperature and humidity, K_T/K_q, it is found that the latter can be much larger than one if the Bowen ratio is small and negative. Despite this, assuming K_T = K_q and using the Bowen ratio method to calculate surface energy fluxes will usually incur only minor errors.     

SOUTHERN BOUNDARY FORCING RELATIVE TO GENESIS, MAINTENANCE AND OSCILLATION OF A SUBTROPICAL HIGH*

A nonlinear critical layer and a Kelvin cat's eye excited thereupon are simulated through four schemes in the context of a nonlinear quasi-geostrophic barotropic vorticity equation model with forced stationary wave acting along the southern boundary to investigate effects of tropical steady forcing on the genesis,maintenance and oscillation of a subtropical high(STH).Evidence suggests that the southern forcing is responsible for the planetary quasi-steady anticyclonic Kelvin cat's eye-form flow field inside the nonlinear critical layer,with the eye shifting,vigor and shape changing quite similar to the behaviors of a summer STH,in striking contrast to the northern stationary forcing.As such,the southern boundary-caused cat's eye is likely to be an even more important mechanism for STH genesis and evolution.In addition,a physical mechanism is introduced for quasi-steady planetary wave moving through the critical layer at subtropical latitudes.     

辽东湾东海岸地区热内边界层特性研究

利用辽东湾东海岸三个低空探测点的实测资料计算了各低空探测点的热内边界层平均高度,并将观测值与用经验公式得到的估算值进行了比较,给出了可用于这一地区的计算公式。     

大气动力学方程组的定性理论及其应用

基于完整的湿大气动力学方程组,利用无穷维动力系统的新理论和新方法,系统讨论了强迫耗散的非线性大气系统的定性理论及其应用。将完整的强迫耗散非线性湿大气动力学方程组化为Hilbert空间中一个等价的算子方程,研究了算子的性质及其物理意义,在此基础上得到湿大气系统全局吸引子的存在性定理,揭示出系统向外源的非线性适应特征,并把结果推广到有地形动力作用和非定常外源强迫的情形。同时探讨了大气方程组惯性流形的存在,大气多平衡态产生的根源以及强迫、耗散和非线性对系统解的渐近行为的影响。在理论结果的基础上,提出强迫耗散的非线性动力系统中存在三类时间边界层、方程组简化准则、分解算法的算子约束原则以及支撑吸引子基底的少数自由度的构造方法,探讨了理论在非线性发展方程差分格式的设计和计算稳定性分析、多平衡态的数值分析、数值模式延伸预报的改进、短期气候预测以及一类中尺度系统分析与预测中的应用,指出描述长期过程动力学模式的必备条件,给出初值与模式相协调的合理解释。最后,对今后的研究方向作了展望。     

区域气候模式侧边界的处理对东亚夏季风降水模拟的影响

在区域气候模式模拟中,侧边界的作用是引入大尺度强迫场。如何处理好侧边界,即大尺度强迫场和区域气候模式本身之间的关系问题,对于区域气候模式模拟和预报东亚夏季风降水具有重要意义。本文利用美国纽约州立大学Albany分校的区域气候模式（SUNYA-ReCM）,设计了两种不同的侧边界处理方法,来探讨侧边界对东亚夏季风降水模拟的影响。驱动区域模式的大尺度强迫场来自欧洲中期天气预报中心（ECMWF）及热带海洋大气研究计划（TOGA）的分析资料场。试验结果表明:（1）当模式的区域较大时,采用较小的侧边界缓冲区会在缓冲区与模式内部的交界处产生不连续;扩大缓冲区并且考虑不同尺度强迫在垂直方向上的不同作用,可以避免这一缺陷。（2）更重要的是采用后一种方案,即减少了区域气候模式在模拟大尺度环流场方面所起的作用,使得模式更多地依赖侧边界来得到更真实的、对东亚夏季风降水起重大影响的一些气流,如副高、西南季风和南海季风,对东亚夏季风降水无论是在大小上还是在雨带位置的演变上都能进行更好的模拟。     

Experimental and computational validation and verification of the Stokes–Darcy and continuum pipe flow models for karst aquifers with dual porosity structure

In our previous study, we developed the Stokes–Darcy (SD) model was developed for flow in a karst aquifer with a conduit bedded in matrix, and the Beavers–Joseph (BJ) condition was used to describe the matrix–conduit interface. We also studied the mathematical well‐posedness of a coupled continuum pipe flow (CCPF) model as well as convergence rates of its finite element approximation. In this study, to compare the SD model with the CCPF model, we used numerical analyses to validate finite element discretisation methods for the two models. Using computational experiments, simulation codes implementing the finite element discretisations are then verified. Further model validation studies are based on the results of laboratory experiments. Comparing the results of computer simulations and experiments, we concluded that the SD model with the Beavers–Joseph interface condition is a valid model for conduit–matrix systems. On the other hand, the CCPF model with the value of the exchange parameter chosen within the range suggested in the literature perhaps does not result in good agreement with experimental observations. We then examined the sensitivity of the CCPF model with respect to the exchange parameter, concluding that, as has previously been noted, the model is highly sensitive for small values of the exchange parameter. However, for larger values, the model becomes less sensitive and, more important, also produces results that are in better agreement with experimental observations. This suggests that the CCPF model may also produce accurate simulation results, if one chooses larger values of the exchange parameter than those suggested in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Assessing uncertainties in dust emission in the Aral Sea region caused by meteorological fields predicted with a mesoscale model

This study investigates how the choice of the planetary boundary layer (PBL) parameterization and dust emission scheme affects the prediction of dust entrainment simulated with a regional mesoscale model. For this analysis we consider a representative dust episode which occurred on April 2001 in the Aral Sea region. The meteorological fields were simulated using the PSU/NCAR MM5 modeling system considering two different boundary layer parameterizations. In each case, emitted dust fluxes were computed off-line by incorporating MM5 meteorological fields into the dust module DuMo. Several dust emission schemes with a prescribed erodible fraction and fixed threshold wind speed were the subject of our analysis. Implications to assessment of the anthropogenic fraction of dust emitted in the Aral region were investigated by conducting the full, half, and no lake modeling experiments.Our results show that the discrepancies in dust fluxes between the two different PBLs are much higher compared to the discrepancy associated with the use of considered dust production schemes. Furthermore, the choice of the PBL affects the timing and duration of a modeled dust event. We demonstrate that different combinations of the PBL parameterization and wind- or friction velocity-driven dust emission schemes can result in up to about a 50% difference in predicted dust mass caused by the Aral Sea desiccation. We found that the drying-up of the Aral cannot only affect the dust emission by expanding the source area, but also by affecting atmospheric characteristics, especially winds. These competitive factors add further complexity to quantification of the anthropogenic dust fraction in the region.     

Tests and applications of an approach to absorbing reflected waves towards incident boundary

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.     

Importance of footing–soil separation on dynamic stiffness of piled embedded footings

Different phenomena such as soil consolidation, erosion, and scour beneath an embedded footing supported on piles may lead to loss of contact between soil and the pile cap underside. The importance of this separation on the dynamic stiffness and damping of the foundation is assessed in this work. To this end, a numerical parametric analysis in the frequency domain is performed using a rigorous three‐dimensional elastodynamic boundary element–finite element coupling scheme. Dimensionless plots relating dynamic stiffness functions computed with and without separation effects are presented for different pile–soil configurations. Vertical, horizontal and rocking modes of oscillation are analyzed for a wide range of dimensionless frequencies. It is shown that the importance of separation is negligible for frequencies below those for which dynamic pile group effects start to become apparent. Redistribution of stiffness contributions between piles and footing is also addressed. Copyright © 2009 John Wiley & Sons, Ltd.