洞庭湖镉迁移转化的马尔可夫模型评价

讨论了在现时泥沙、水体、鱼类、藻类的含镉状态下洞庭湖的输入输出通量；并以其为例，在简要介绍马尔可夫链数学方法的基础上，建立了湖泊生态系统中镉迁移转化的马尔可夫链数学模型，初步评价了当洞庭湖不断有镉加入、系统到达状态稳定时各种介质中的镉总量及镉的迁移转化趋势。     

Interpretation of a two-point flux stencil for skew parallelogram grids

Control-volume formulations for elliptic equations often use two-point flux stencils, even for skew grids. Any two-point flux stencil may be interpreted as a multipoint flux stencil. This yields a definition of the permeability (or conductivity) tensor. Formulas for calculating the permeability tensor, based on the user-specified quantities in the two-point flux stencil, are given. Numerical test examples demonstrate the validity of the derivation.      

内蒙古夏季干旱的水汽输送特征分析

对于大气中水汽的定性、定量测量是提高内蒙古干旱研究与预测的有效手段。利用NCEP再分析资料和内蒙古地区实测降水资料（1961-2000年），定性、定量地分析东亚水汽分布、水汽输送与内蒙古夏季干旱的关系。得出内蒙古地区上空低层大气的含水量偏少是内蒙古夏季干旱的一个重要因素。由于青藏高原大地形作用，孟加拉湾和南海的水汽北上受阻，不能达到内蒙古高原是形成内蒙古干旱的又一原因。受大气环流影响，内蒙古高原水汽辐合偏弱是内蒙古夏季干旱的重要原因。青藏高原的水汽通量与内蒙古东西部降水都存在明显负相关；孟加拉湾至内蒙古西部是内蒙古西部地区降水的水汽通道，南海至中国东北地区为内蒙古东部降水的水汽通道。     

Relationship between global seismicity and solar activities

The relations between sunspot numbers and earthquakes (M≧6), solar 10.7 cm radio flux and earthquakes, solar proton events and earthquakes have been analyzed in this paper. It has been found that: (1) Earthquakes occur frequently around the minimum years of solar activity. Generally, the earthquake activities are relatively less during the peak value years of solar activity, some say, around the period when magnetic polarity in the solar polar regions is reversed. (2) the earthquake frequency in the minimum period of solar activity is closely related to the maximum annual means of sunspot numbers, the maximum annual means of solar 10.7 cm radio flux and solar proton events of a whole solar cycle, and the relation between earthquake and solar proton events is closer than others. (3) As judged by above interrelationship, the period from 1995 to 1997 will be the years while earthquake activities are frequent. In the paper, the simple physical discussion has been carried out. These results supported the exploration and studies of some researchers to a certain extent. This work is supported by Foundation of the Chinese Academy of Sciences (major item).     

Estimation of Sensible Heat Flux from Surface Temperature Wave and One-Time-of-Day Air Temperature Observation

By means of the algorithm presented here, the temporal course H(t) and the daily mean H¯ of the sensible heat flux H can be estimated from measurements of the thermodynamic surface temperature (as a function of time) and from a one-time-of-day air temperature observation. In addition to these temperatures, one needs estimates for daily mean wind speed, for the roughness lengths of momentum and heat transfer, and for the displacement height. First, a quite general solution of the equation for heat conductance (equation for the vertical profile of potential temperature (z,t)) in the dynamic sublayer will be presented. The undetermined parameters in this solution will be defined with the aid of the above mentioned measurements. The influence of horizontal advection will be taken into account. After that, the sensible heat flux can be evaluated from the temperature difference between surface and air with the well known resistance formulae. In this paper the algorithm is derived for areas with homogeneous surfaces, i.e., with uniform surface temperatures. Finally, the method will be verified by measurements taken during the field campaigns HIBE 89 (Hildesheimer Börde in Germany) and EFEDA 91 (Spain). The root mean square errors (RMSE) for the comparison between measurement and model with regard to the temperature difference of surface and air amount to one or two degrees Kelvin, and the error of H¯ reaches 10 to 25 per cent. The method can be used to determine the sensible heat flux from measurements of surface temperatures by satellites (e.g., METEOSAT), but can also be applied to ground based measurements. For instance, horizontal temperature advection can be estimated from measurements at a single location, especially if more than one near-surface air temperature is available. The procedure can be generalized for larger areas, which consist of various surface types with different surface temperatures. This generalization of the algorithm is in progress and will be addressed in a subsequent paper. It will allow us to improve the estimates for H(t) by means of temperature measurements from, e.g., NOAA/AVHRR or LANDSAT/TM, taking into account the heterogeneity of the area that is contained in one METEOSAT pixel.     

Analysis of Astex-Stratocumulus Observational Data Using a Mass-Flux Approach

A mass-flux approach is applied to observational data obtained in a convective boundary layer topped with stratocumulus clouds. The observational data were obtained from aircraft measurements during the Atlantic Stratocumulus Transition Experiment (ASTEX). A conditional sampling method is used to calculate average updraft and downdraft values. The vertical fluxes calculated with the mass-flux approach are found to be proportional to the real (measured) fluxes, with a proportionality factor being about 0.6. This value is predicted by theory for two variables having a joint Gaussian distribution function; proportionality factor = 2^-1 0.637. The horizontal fractional entrainment and detrainment rates calculated from the data ( 1–2 × 10^-2 m^-1) are an order of magnitude higher than the rates obtained by large eddy simulations for cumulus convection ( 2–3 × 10^-3 m^-1) and two orders of magnitude higher than those used in modelling cumulus convection with a mass-flux scheme in an operational weather forecast model ( 3 × 10^-4 m^-1). A numerical mass-flux model for the thermodynamics was developed and showed that results are in good agreement when compared with measured profiles of the liquid water content.     

CONCENTRATION AND FLUX OF CO_2, TURBULENCE FLUXES AND RADIATION BALANCE IN THE NEAR SURFACE LAYER OVER WHEAT FIELD

The gradient of CO_2 concentration, photosynthetically active radiation (PAR). net radiation.soil heat flux, profiles of wind speed, and air temperature and humidity were measured above awheat field during May and June 1985 at Beijing Agro-Ecosystems Experimental Station. Beijing,China. Fluxes of carbon dioxide, sensible heat, latent heat and momentum were calculated byusing the aerodynamic method. The observation site. equipment, calibration techniques, theerrors associated with the measurement, and the computational procedures are described. Theresults show that the diurnal variations of amplitude of CO_2 concentrations were 103.8 to 27. 0. 86.3 to 22.8 and 69.8 to 11.6 ppm: the average CO_2 concentrations were 331.5. 339.9 and 364.6 ppmfor the photosynthesis type, and 369.6. 364.0 and 375.2 ppm for the respiration type at 1. 2 and10 m above surface, respectively, from May 14 to June 15. In the daytime, transfer direction ofthe CO_2 fluxes and gradients is from air to crop canopy, and at noon (1100 to 1300 BT (BeijingTime)) the transfer rate reaches negative maximum value. At night, transfer of CO_2 fluxes andgradients is in the reversed direction and reaches positive maximum in the early morning (0400 to0600 BT). There are strong negative correlations between CO_2 flux and the net radiation (Rn),available energy (H LE). photosynthetically active radiation (PAR) and momentum flux (τ).     

Chandra observations of neutron stars: an overview

We present an overview of Chandra X-ray Observatory observations of neutron stars. The outstanding spatial and spectral resolution of this great observatory have allowed for observations of unprecedented clarity and accuracy. Many of these observations have provided new insights into neutron star physics. We present an admittedly biased and overly brief review of these observations, highlighting some new discoveries made possible by the Observatory’s unique capabilities. This includes our analysis of recent multiwavelength observations of the putative pulsar and its pulsar-wind nebula in the IC443 SNR.      

A new soil freeze/thaw discriminant algorithm using AMSR‐E passive microwave imagery

The soil freeze–thaw controls the hydrological and carbon cycling and thus affects water and energy exchanges at land surface. This article reported a newly developed algorithm for distinguishing the freeze/thaw status of surface soil. The algorithm was based on information from Advanced Microwave Scanning Radiometer Enhanced (AMSR‐E) which records brightness temperature (Tb) in the afternoon and after midnight. The criteria and discriminant functions were obtained from both radiometer observations and model simulations. First of all, the microwave radiation from freeze–thaw soil was examined by carrying out experimental measurements at 18·7 and 36·5 GHz using a Truck‐mounted Multi‐frequency Microwave Radiometer (TMMR) in the Heihe River of China. The experimental results showed that the soil moisture is a key component that differentiates the microwave radiation behaviours during the freeze–thaw process, and the differences in soil temperature and emissivity between frozen and thawed soils were found to be the most important criteria. Secondly, a combined model was developed to consider the impacts of complex ground surface conditions on the discrimination. The model simulations quite followed the trend of in situ observations with an overall relation coefficient (R) of approximately 0·88. Finally, the ratio of Tb18·7H (horizontally polarized Tb at 18·7 GHz) to Tb36·5V was considered primarily as the quasi‐emissivity, which is more reasonable and explicit in measuring the microwave radiation changes in soil freezing and thawing than the spectral gradient. By combining Tb36·5V to indicate the soil temperature variety, a Fisher linear discrimination analysis was used to establish the discriminant functions. After being corrected by TMMR measurements, the new discriminant algorithm had an overall accuracy of 86% when validated by 4‐cm soil temperature. The multi‐year discriminant results also provided a good agreement with the classification map of frozen ground in China. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of seasonal crop water consumption in a vineyard using Bowen ratio‐energy balance method

Information about seasonal crop water consumption is useful to develop the appropriate irrigation scheme. Measurements of energy balance components using the Bowen ratio method were made for a complete growing season at a vineyard in the arid region of northwest China. Vine in the experiment was furrow‐irrigated using a trellis system. The measured evapotranspiration was compared with estimates using the soil water balance method. It is shown that the Bowen ratio method provided accurate estimates of evapotranspiration from the vineyard and this requires that the Bowen ratio system is appropriately installed. The energy balance components showed typical diurnal pattern with peaks that occurred around the midday, except for the ground heat flux which delayed its peak by 2–3 h. The sensible heat flux was greater than the latent heat flux and followed the net radiation closely. The ratio of the latent heat flux to net radiation was low in the early growing season and increased over time. Under the limited irrigation experienced in the vineyard, the latent heat flux was controlled by available soil moisture and the total evapotranspiration in the growing season was 253 mm. The seasonal progression of the crop coefficient is similar to that reported in the literature, with the maximum occurring during the month of September. The crop coefficient can be estimated as a non‐linear function of day of year (DOY) and used to estimate evapotranspiration from vineyards in the region. Copyright © 2007 John Wiley & Sons, Ltd.