四川盆地东部再生稻高温伏旱区腋芽萌发期气象适宜度

以四川盆地东部再生稻高温伏旱区为研究区,基于气象条件对腋芽萌发期再生稻生长发育的影响机理,采用结构方程模型探明了腋芽萌发的影响因素、影响路径及影响强度。进一步结合隶属函数和层次分析法,构建了再生稻腋芽萌发气象适宜度模型,并分析了1981—2021年研究区腋芽萌发期气象影响因素及气象适宜度的变化特征。结果表明：气温、空气湿度、降水是四川盆地东部再生稻高温伏旱区腋芽萌发的关键影响因素,基于以上因素的气象适宜度模型能较好地评价再生稻腋芽萌发期的气象影响。1981—2021年研究区气温适宜度、空气湿度适宜度及综合气象适宜度均呈下降趋势,降水适宜度没有表现出明显的变化趋势。气温适宜度、空气湿度适宜度及综合气象适宜度总体呈“西高东低”的空间分布特征,降水适宜度则表现为“东西高、中部低”。1981—2021年研究区气温上升趋势与空气湿度下降趋势显著,导致致害高温（日平均气温≥32℃）积温和致害低湿（日平均相对湿度≤65％）日数明显增加,这是再生稻腋芽萌发期气象适宜度总体呈下降趋势的诱因。     

Enhancement of the channel routing module in SWAT

The objective of this study is to demonstrate that the channel flow routing techniques used in the Soil and Water Assessment Tool (SWAT) can be sometimes inappropriate for runoff simulation in small basins. Simulated hydrographs have a tendency to underestimate peak flows or may send a false signal during the recession periods when SWAT is applied to the Mihocheon Basin in South Korea. This was particularly evident for sub‐basins that had a short travel time of much less than a day. In order to enhance the channel routing module in SWAT, an alternative routing technique in which Manning relationship is combined with a simple channel reach continuity equation is proposed in the present study. The advantage of the proposed routing technique is that parameters are readily available from channel morphological data and that it is applicable to small basins. Through simulations of specific events, this routing module that is incorporated into SWAT is almost free from instability errors and produces realistic flow depths and velocities unlike the original version of SWAT with Muskingum method. Copyright © 2009 John Wiley & Sons, Ltd.     

A simple model to estimate monthly forest evapotranspiration in Japan from monthly temperature

The estimation of evapotranspiration (E) in forested areas is required for various practical purposes (e.g. evaluation of drought risks) in Japan. This study developed a model that estimates monthly forest E in Japan with the input of monthly temperature (T). The model is based on the assumptions that E equals the equilibrium evaporation rate (E_eq) and that E_eq is approximated by a function of T. The model formulates E as E (mm month⁻¹) = 3·48 T ( °C) + 32·3. The accuracy of the model was examined using monthly E data derived using short‐term water balance (WB) and micrometeorological (M) methods for 15 forest sites in Japan. The model estimated monthly E more accurately than did the Thornthwaite and Hamon equations according to regression analysis of the estimated E and E derived using the WB and M methods. Although the model tended to overestimate monthly E, the overestimation could be reduced by considering the effect of precipitation on E. As T data are commonly available all over Japan, the model would be a useful tool to estimate forest E in Japan. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Water waves generated by instantaneous disturbances at the bed of a sloping beach

The two-dimensional problem of the generation of water waves due to instantaneous disturbances prescribed at the bed of a beach sloping at an arbitrary angle is studied here. It is formulated in terms of an initial-boundary-value problem for the velocity potential describing the motion in the fluid region assuming the linear theory. Using the Laplace transform in time and the Mellin transform in distance, the problem is reduced to solving a difference equation whose method of solution is of considerable importance in the literature. The form of the free surface is obtained in terms of a multiple infinite integral that is evaluated by the method of steepest-descent. For some prescribed forms of the disturbance at the bed of the beach, the free surface is depicted in a number of figures for different beach angles. It is observed that as the beach angle decreases, the maximum wave height increases, which is plausible.     

Curve‐Number/Green–Ampt mixed procedure for streamflow predictions in ungauged basins: Parameter sensitivity analysis

A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.     

Falling sphere observations of anisotropic gravity wave motions in the upper stratosphere over Australia

A study of inertial scale gravity wave motions in the region of the atmosphere between 30 and 60 km has been undertaken, using wind and temperature data derived from rocket-borne falling sphere density experiments performed over Woomera, Australia between 1962 nad 1976. The gross features of the wave field compare favorably with those found in similar northern hemispheric studies. Wave propagation is found to be both vertically and horizontally anisotropic. A rotary spectral analysis indicates predominately upgoing wave energy, suggesting that the majority of sources of these waves lie below 30 km. A detailed statistical investigation of the waves, made using the Stokes parameters technique, reveals that phase progression is also highly directional in the horizontal, with a significant zonal component in summer, but with a strong meridional component in winter. Propagation towards the southeast is inferred in summer, with the waves possibly emanating from tropospheric sources in equatorial regions to the north of Australia. The technique also shows that, on average, the waves appear to have mean ellipse eccentricities (=f/) around 0.4–0.45. Indirect estimates of a number of important wave parameters are made. In particular,v andw flux estimates are made over several height intervals. The vertical gradient of density weighted flux implies wave-induced mean flow accelerations of the order 0.1–1 ms^–1day^–1. This suggests that dissipating gravity waves are a significant source of the momentum residuals that are encountered in studies of satellite data from this region.     

Evaluation of rainfall-runoff performance using the stochastic integral equation method

The stochastic integral equation method (S.I.E.M.) is used to evaluate the relative performance of a set of both calibrated and uncalibrated rainfall-runoff models with respect to prediction errors. The S.I.E.M. is also used to estimate confidence (prediction) interval values of a runoff criterion variable, given a prescribed rainfall-runoff model, and a similarity measure used to condition the storms that are utilized for model calibration purposes.Because of the increasing attention given to the issue of uncertainty in rainfall-runoff modeling estimates, the S.I.E.M. provides a promising tool for the hydrologist to consider in both research and design.     

位场DFT算法研究

本文推广了经典的抽样定理,并据此导出了函数有限离散傅里叶变换误差方程（简称DFT误差方程,下同）。该方程把有限离散傅里叶变换中固有的离散效应和有限效应表示为确切的数学形式。离散效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数;有限效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数的DFT。 基于DFT误差方程和位函数特点,作者提出了两种位场数值傅里叶变换新算法--移样法和等效源续尾叠样法。移样法可近百倍地提高位场数值傅里叶反变换的精度,等效源续尾叠样法可数十倍地提高正变换精度。两种算法都不需要增加资料长度和取样密度,因而基本不需要增加计算机时间和内存。文中给出了算例。