Tower-based conditional sampling for measuring ecosystem-scale carbon dioxide exchange in coastal wetlands

Long-term measurements of CO₂ exchange between coastal wetlands and the atmosphere are necessary to improve our understanding of the role these ecosystems play in the global carbon cycle, and the response of these systems to environmental change. We conducted research to adapt and evaluate tower-based conditional sampling as a method for measuring net CO₂ exchange (NCE) at the ecosystem scale on a continuous basis. With conditional sampling, NCE is determined from the product of the standard deviation of vertical wind velocity, the difference in CO₂ concentration between updrafts and downdrafts in the constant flux portion of the boundary layer above the surface, and an empirical coefficient. We constructed a system that used a sonic anemometer to measure vertical wind velocity (w) and control a high-speed three-way valve that diverted air from updrafts and downdrafts into separate sample lines, depending on the direction ofw. an infrared gas analyzer was used to measure the concentration difference. The conditional sampling system was installed and tested in a marsh in the Nueces River Delta near Corpus Christi, Texas, as part of a long-term study of effects of freshwater inflow on CO₂ flux. System accuracy was evaluated by comparing conditional sampling measurements of water vapor flux with independent estimates obtained with the Bowen ratio method. Average daily flux estimates for the two methods agreed to within 13%. Measurements showed that freshwater inflow due to flooding of the Nueces River increased NCE by increasing CO₂ assimilation and decreasing CO₂ efflux. Over a 65-d period, daily NCE varied from a maximum gain of 0.16 mol CO₂ m⁻² d⁻¹ during flooding to a maximum loss of −0.14 mol CO₂ m⁻² d⁻¹ when the marsh dried. Our study showed that conditional sampling was well suited for quantifying CO₂ exchange in coastal wetlands on a diel, daily, and seasonal basis.     

The value ofC efor the Arabian Sea during summer monsoon

We estimate, from the moisture budget the bulk aerodynamic coefficient for latent heat flux (C _e) during the monsoon season over the central Arabian Sea. The average value ofC _eunder active monsoon conditions was found to be 2.25 × 10⁻³ which is nearly 60% higher than those previously used.     

A study of turbulent characteristics of atmospheric boundary layer over monsoon trough region using Kytoon and Doppler sodar

As a part of the MONTBLEX-90 observational programme, Kytoon and Doppler sodar observations were taken at Kharagpur. These data are analysed to study the turbulent characteristics of the atmospheric boundary layer in terms of stability, temperature structure function (C _T ² ) and velocity structure function (C _v ² ).C _T ² follows aZ ^−4/3 law on most of the days, whereas the variation ofC _V ² is not systematic.C _V ² andC _T ² values are found to vary between 10⁻⁵−10⁻¹ m^4/3s⁻² and 10⁻⁵−10⁻²°C² m^−2/3 respectively.     

Some observations from the data taken in and around Kharagpur during the onset of the monsoon, 1990

During MONTBLEX 1990, various observational platforms were operated at Kharagpur and the nearby Kalaikunda Air Base. Using the data from all the platforms, one can draw the following conclusions. The temperature and wind data obtained from various sensors have overall compatibility. Sodar wind data indicate the presence of a low level jet at around 300 m above ground. The inversion height may be evaluated from the vertical profile of the sodar back-scatter echo intensity. The sub-synoptic or synoptic scale convergence modulates the inversion height and the presence of cloud-base within the inversion height in turn modulates the sensible heat and momentum fluxes.     

Dust emission and transport in Mali, West Africa

Vertical dust fluxes were measured in the Inland Delta region of the Niger River, Mali, West Africa, during April-June, 1989 and 1990. Measurements of dust flux represent, for the most part, non-dust storm conditions or ‘dust haze’ periods. The observed concentration versus height relationships are similar to data presented by other investigators. The relationship between wind shear velocity (u_*, m s^?1) and vertical dust flux (F, μg m^?2 s^?1) can be described by a relationship in which F is proportional to u_*⁴. However, there is considerable scatter within the data set which is attributable to textural controls and surface conditions. The vertical dust fluxes measured in Mali are compared to dust fluxes measured in Texas, USA, and Yukon Territories, Canada. The significantly different values for the constant of proportionality (a) in the F α a u_*⁴ relationship for these geographically diverse areas is a function of surficial controls on the release of sediments to the air stream. Dust concentrations measured in Mali were found to be uniformly high and in general exceed WHO health standards for acceptable total suspended particulate loadings. As a result background dust may be considered a long term stress on health for the people of this region of Mali.     

Electrical conductivity of hydrous basaltic melts: implications for partial melting in the upper mantle

The Earth’s uppermost asthenosphere is generally associated with low seismic wave velocity and high electrical conductivity. The electrical conductivity anomalies observed from magnetotelluric studies have been attributed to the hydration of mantle minerals, traces of carbonatite melt, or silicate melts. We report the electrical conductivity of both H₂O-bearing (0–6 wt% H₂O) and CO₂-bearing (0.5 wt% CO₂) basaltic melts at 2 GPa and 1,473–1,923 K measured using impedance spectroscopy in a piston-cylinder apparatus. CO₂ hardly affects conductivity at such a concentration level. The effect of water on the conductivity of basaltic melt is markedly larger than inferred from previous measurements on silicate melts of different composition. The conductivity of basaltic melts with more than 6 wt% of water approaches the values for carbonatites. Our data are reproduced within a factor of 1.1 by the equation log σ = 2.172 − (860.82 − 204.46 w ^0.5)/(T − 1146.8), where σ is the electrical conductivity in S/m, T is the temperature in K, and w is the H₂O content in wt%. We show that in a mantle with 125 ppm water and for a bulk water partition coefficient of 0.006 between minerals and melt, 2 vol% of melt will account for the observed electrical conductivity in the seismic low-velocity zone. However, for plausible higher water contents, stronger water partitioning into the melt or melt segregation in tube-like structures, even less than 1 vol% of hydrous melt, may be sufficient to produce the observed conductivity. We also show that ~1 vol% of hydrous melts are likely to be stable in the low-velocity zone, if the uncertainties in mantle water contents, in water partition coefficients, and in the effect of water on the melting point of peridotite are properly considered.     

Variability in the surface wind direction at a coastal site of complex terrain

Variability in the standard deviation of surface wind direction (σ_θ), under different Pasquill stability regimes on diurnal, seasonal and interannual scales has been investigated making use of a 10-year data set collected at Visakhapatnam (17°42′ N., 82° 18′ E) during January, April, August and October for winter, pre-monsoon, monsoon and post-monsoon seasons respectively. The diurnal scale variability in σ_θ is more pronounced during day time than in night. The seasonal variability in σ_θ is only moderate around noon while relatively large fluctuations are noticed on inter-annual scale only during day time in January and August. The seasonal dispersion in σ_θ decreased from most unstable regime to most stable regime.     

An analysis of MONTBLEX data on heat and momentum flux at Jodhpur

Parameterization of sensible heat and momentum fluxes as inferred from an analysis of tower observations archived during MONTBLEX-90 at Jodhpur is proposed, both in terms of standard exchange coefficientsC _H andC _D respectively and also according to free convection scaling. Both coefficients increase rapidly at low winds (the latter more strongly) and with increasing instability. All the sensible heat flux data at Jodhpur (wind speed at 10 m Ū₁₀ < 8 ms⁻¹) also obey free convection scaling, with the flux proportional to the ‘4/3’ power of an appropriate temperature difference such as that between 1 and 30 m. Furthermore, for Ū₁₀ < 4 ms⁻¹ the momentum flux displays a linear dependence on wind speed.     

Movement of seasonal eddies and its relation with cyclonic heat potential and cyclogenesis points in the Bay of Bengal

Movement of seasonal eddies in the Bay of Bengal (BOB) and its relation with cyclonic heat potential (CHP) and cyclogenesis points have been investigated in this study using 6 years (2002–2007) of global ocean monthly analysis datasets based on the Simple Ocean Data Assimilation (SODA) package (SODA v2.0.4) of Carton et al. (2005) and Indian Meteorological Department cyclogenesis points. The region dominated by anticyclonic eddies with CHP greater than 70 × 10⁷ J/m² as well as good correlations (>0.9) with sea surface height (SSH) and 26°C isothermal depth (D ₂₆) can be a potential region of cyclogenesis. The region dominated by cyclonic eddies with CHP greater than 50 × 10⁷ J/m² and good correlation (>0.9) with both SSH and D ₂₆ can serve as a potential region of high-level depression. Potential cyclogenesis regions are the southern BOB (5°N–12°N) for the post-monsoon season and the head of BOB (north of 15°N) during southwest monsoon. Seven potential regions are identified for the eddy formation for different seasons, which are consistent with the cyclogenesis points. The CHP distributions alone are able to explain the cyclone tracks for the pre-monsoon and post-monsoon seasons but not for the monsoon season.     

Identification of a surface layer structure and analysis of humidity data in two weather situations at Jodhpur (26°18′N, 73°04′E), India,during MONTBLEX 1990

The Monsoon Trough Boundary Layer Experiment held in 1990 was a multi-institutional effort to probe the atmospheric boundary layer over the monsoon trough over northern India. For this experiment, four micrometeorological towers were set up at four different locations along the normal position of the trough. One such tower of 30m height was located at Jodhpur (26‡18′N, 73‡04′E), Rajasthan. The fast and slow response data available during the experiment have been used in the present study to determine a suitable layer-structure of the surface layer for evaluation of sensible heat flux using the multilayer hypothesis of Kramm (1989).     

Temporal variability in summertime bottom hypoxia in shallow areas of Mobile Bay,Alabama

This paper addresses temporal variability in bottom hypoxia in broad shallow areas of Mobile Bay, Alabama. Time-series data collected in the summer of 2004 from one station (mean depth of 4 m) exhibit bottom dissolved oxygen (DO) variations associated with various time scales of hours to days. Despite a large velocity shear, stratification was strong enough to suppress vertical mixing most of the time. Bottom DO was closely related to the vertical salinity gradient (ΔS). Hypoxia seldom occurred when ΔS (over 2.5 m) was <2 psu and occurred almost all the time when ΔS was >8 psu in the absence of extreme events like hurricanes. Oxygen balance between vertical mixing and total oxygen demand was considered for bottom water from which oxygen demand and diffusive oxygen flux were estimated. The estimated decay rates at 20°C ranging between 0.175–0.322 d⁻¹ and the corresponding oxygen consumption as large as 7.4 g O₂ m⁻² d⁻¹ fall at the upper limit of previously reported ranges. The diffusive oxygen flux and the corresponding vertical diffusivity estimated for well mixed conditions range between 8.6–9.5 g O₂ m⁻² d⁻¹ and 2.6–2.9 m² d⁻¹, respectively. Mobile Bay hypoxia is likely to be associated with a large oxygen demand, supported by both water column and sediment oxygen demands, so that oxygen supply from surface water during destratification events would be quickly exhausted to return to hypoxic conditions within a few hours to days after destratification events are terminated.     

Homogeneous Indian Monsoon rainfall: Variability and prediction

The Indian summer monsoon rainfall is known to have considerable spatial variability, which imposes some limitations on the all-India mean widely used at present. To prepare a spatially coherent monsoon rainfall series for the largest possible area, fourteen subdivisions covering the northwestern and central parts of India (about 55% of the total area of the country), having similar rainfall characteristics and associations with regional/global circulation parameters are merged and their area-weighted means computed, to form monthly and seasonal Homogeneous Indian Monsoon (HIM) rainfall series for the period 1871–1990. This paper includes a listing of monthly and seasonal rainfall of HIM region. HIM rainfall series has been statistically analysed to understand its characteristics, variability and teleconnections for long-range prediction. HIM rainfall series isfound to be homogeneous, Gaussian distributed and free from persistence. The mean (R) rainfall is 757 mm (87% of annual) and standard deviation (S) 119 mm, with a Coefficient of Variation (CV) of 16%. There were 21 dry (K, -<R S) and 19 wet (R _i R + S) years during 1871–1990. There were clusters of frequent negative departures during 1899–1920 and 1965–1987 and positive departures during 1942–1961. The recent three decades show very high rainfall variability with 10 dry and 6 wet years. The decadal averages were alternatively positive and negative for three consecutive decades, viz., 1871–1900 (positive); 1901–1930 (negative); 1931–1960 (positive) and 1961–1990 (negative) respectively. Significant QBO and autocorrelation at 14th lag have been found in HIM rainfall series. To delineate the changes in the climatic regime of the Indian summer monsoon, sliding correlation coefficients (CCs) between HIM rainfall series and (i) Bombay msl pressure, (ii) Darwin msl pressure and (iii) Northern Hemisphere surface air temperature over the period 1871–1990 have been examined. The 31-year sliding CCs showed the systematic turning points of positive and negative CCs around the years, 1900 and 1940. In the light of other corroborative evidences, these turning points seem to delineate ‘meridional’ monsoon regime during 1871–1900 and 1940–1990 and ‘zonal’ monsoon regime during 1901–1940. The monsoon signal is particularly dominant in many regional and global circulation parameters, during 1951–1990. Using the teleconnections ofHIM series with 12 regional/global circulation parameters during the recent 36-year period 1951–86 regression models have been developed for long-range prediction. In the regression equations 3 to 4 parameters were entered, explaining upto 80% of the variance, depending upon the data period. The parameters that prominently enter the multiple regression equations are (i) Bombay msl pressure, (ii) April 500 mb Ridge at 75°E, (iii) NH temperature, (iv) Nouvelle minus Agalega msl pressure and (v) South American msl pressure. Eleven circulation parameters for the period 1951–80 were subjected to Principal Component Analysis (PCA) and the PC’s were used in the regression model to estimate HIM rainfall. The multiple regression with three PCs explain 72% of variance in HIM rainfall.     

Sensitivity of high-resolution tropical cyclone intensity forecasts to surface flux parameterization

Surface flux parameterization schemes used in current dynamic models are primarily based upon measurements at low and moderate wind speeds. Recent studies show that these parameterization schemes may be incorrect at high wind speeds (e.g., tropical cyclone forecasts). Five high-resolution numerical model experiments are designed to assess the sensitivity of tropical cyclone intensity forecasts to changes in the surface flux parameterization. The sensitivity experiments are conducted by running 48 h forecasts of the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) for six selected tropical cyclones with individual modifications to surface flux calculation that include: (1) limiting the surface stress for wind speeds greater than 33 m s⁻¹, or 64 knots (kt); (2) computing the stress at the top of the model bottom grid layer (MBGL) by averaging results from surface layer similarity and turbulence mixing parameterization for wind speeds greater than 33 m s⁻¹; (3) increasing the roughness lengths for heat and moisture transfer by a factor of ten; (4) setting the roughness lengths for heat and moisture transfer to 1/10 of the momentum roughness length; and (5) cooling the sea surface temperature (SST) by a prescribed rate at high winds. Averaged responses for the six storms to these sensitivity tests show that: (i) the limit on surface stress at high winds significantly increases the cyclone intensity in 48 h forecasts; (ii) the averaged surface layer stress at high winds increases the cyclone intensity but to a much lesser degree than limiting the surface stress; (iii) large increases in the roughness lengths for heat and moisture transfer are needed to significantly impact the intensity forecast; (iv) the different roughness length formula for surface transfer coefficients notably increases C _h/C _d ratio from 0.59 to 0.79 for 25 m s⁻¹ and 0.41 to 0.75 for 50 m s⁻¹ that significantly increases the predicted cyclone intensity; and (v) cooling of the SST by −5.8°C in 48 h reduces the maximum surface wind speed by −32 kt, or 16.5 m s⁻¹, at 48 h forecast. These results suggest that a surface flux parameterization scheme suitable for tropical cyclone intensity forecast must correctly model the leveling-off character of surface stress and C _h/C _d ratio at high winds. All modifications to surface flux calculation have little influence on 48 h track forecasts, even though they may significantly impact the intensity forecasts.

Excess heat capacity and entropy of mixing along the chlorapatite–fluorapatite binary join

The heat capacity at constant pressure, C _p, of chlorapatite [Ca₅(PO₄)₃Cl – ClAp], and fluorapatite [Ca₅(PO₄)₃F – FAp], as well as of 12 compositions along the chlorapatite–fluorapatite join have been measured using relaxation calorimetry [heat capacity option of the physical properties measurement system (PPMS)] and differential scanning calorimetry (DSC) in the temperature range 5–764 K. The chlor-fluorapatites were synthesized at 1,375–1,220°C from Ca₃(PO₄)₂ using the CaF₂–CaCl₂ flux method. Most of the chlor-fluorapatite compositions could be measured directly as single crystals using the PPMS such that they were attached to the sample platform of the calorimeter by a crystal face. However, the crystals were too small for the crystal face to be polished. In such cases, where the sample coupling was not optimal, an empirical procedure was developed to smoothly connect the PPMS to the DSC heat capacities around ambient T. The heat capacity of the end-members above 298 K can be represented by the polynomials: C _p^ClAp = 613.21 − 2,313.90T ^−0.5 − 1.87964 × 10⁷ T ⁻² + 2.79925 × 10⁹ T ⁻³ and C _p^FAp = 681.24 − 4,621.73 × T ^−0.5 − 6.38134 × 10⁶ T ⁻² + 7.38088 × 10⁸ T ⁻³ (units, J mol⁻¹ K⁻¹). Their standard third-law entropy, derived from the low-temperature heat capacity measurements, is S° = 400.6 ± 1.6 J mol⁻¹ K⁻¹ for chlorapatite and S° = 383.2 ± 1.5 J mol⁻¹ K⁻¹ for fluorapatite. Positive excess heat capacities of mixing, ΔC _p^ex, occur in the chlorapatite–fluorapatite solid solution around 80 K (and to a lesser degree at 200 K) and are asymmetrically distributed over the join reaching a maximum of 1.3 ± 0.3 J mol⁻¹ K⁻¹ for F-rich compositions. They are significant at these conditions exceeding the 2σ-uncertainty of the data. The excess entropy of mixing, ΔS ^ex, at 298 K reaches positive values of 3–4 J mol⁻¹ K⁻¹ in the F-rich portion of the binary, is, however, not significantly different from zero across the join within its 2σ-uncertainty.     

Can shells of freshwater mussels (Unionidae) be used to estimate low summer discharge of rivers and associated droughts?

This paper examines if shell oxygen isotope ratios (δ¹⁸O_ar) of Unio sp. can be used as a proxy of past discharge of the river Meuse. The proxy was developed from a modern dataset for the reference time interval 1997–2007, which showed a logarithmic relationship between discharge and measured water oxygen isotope ratios (δ¹⁸O_w). To test this relationship for past time intervals, δ¹⁸O_ar values were measured in the aragonite of the growth increments of four Unio sp. shells; two from a relatively wet period and two from a very dry time interval (1910–1918 and 1969–1977, respectively). Shell δ¹⁸O_ar records were converted into δ¹⁸O_w values using existing water temperature records. Summer δ¹⁸O_w values, reconstructed from δ¹⁸O_ar of 1910–1918, showed a similar range as the summer δ¹⁸O_w values for the reference time interval 1997–2007, whilst summer reconstructed δ¹⁸O_w values for the time interval 1969–1977 were anomalously high. These high δ¹⁸O_w values suggest that the river Meuse experienced severe summer droughts during the latter time interval. δ¹⁸O_w values were then applied to calculate discharge values. It was attempted to estimate discharge from the reconstructed δ¹⁸O_w values using the logarithmic relationship between δ¹⁸O_w and discharge. A comparison of the calculated summer discharge results with observed discharge data showed that Meuse low-discharge events below a threshold value of 6 m³/s can be detected in the reconstructed δ¹⁸O_w records, but true quantification remains problematic.     

A New Approach to the Modeling of the Pressure Dependency of the Strength of Rocks

Based on the hypothesis that shearing stresses and normal tensile stresses both play a decisive role in fracturing and brittle failure of rock material, a novel strength criterion was developed in one of the earlier works of the author. In the criterion, a certain parameter ν′ occurs which depends on the structure of the material. Originally, the parameter was treated as a constant, which resulted in a linear form of the strength function ^F σ ₁ = f(σ ₃), where σ ₃ = σ ₂. Although the linear strength criterion is sometimes found to be applicable to various particular rock materials, it is not, in general, of a universal character. Analysis of the triaxial test results for 75 different sandstones revealed that parameter ν′ usually increases in an exponential or linear manner as confining pressure increases, and only in isolated cases does it seem to be independent of the confining pressure. For these three types of function ν′ = f(p) appropriate strength criteria ^F σ ₁ = f(σ ₃) are given in the present paper. These criteria were used to fit all of the collected empirical data sets. In general, a very good fit to the data was obtained.     

贵州新民铝土矿矿床Li的地球化学特征与富集机制探究

新民铝土矿床位于黔北务川-正安-道真地区(简称务正道地区)大塘向斜东翼,铝土矿(岩)型Li资源丰富.含矿岩系大竹园组(P1d)不同岩性的w(Li)有差别:土状-半土状铝土矿平均w(Li)为16.34×10-6,致密块状铝土矿平均w(Li)是803.84×10-6,铝土岩平均w(Li)是1436.22×10-6,黏土岩的平均w(Li)是786.62×10-6,梁山组泥岩的平均w(Li)是51.82×10-,韩家店群泥岩(页岩)的平均w(Li)是48.52×10-6,黄龙组灰岩的平均w(Li)是 11.99×10-6.由此可知,研究区铝土矿(岩)型Li资源主要富集在含矿岩系中上部的致密块状铝土矿、铝土岩和黏土岩中,顶板、底板以及土状-半土状铝土矿w(Li)较低.Li主要赋存于高岭石中,伊利石(水云母)和残存的三水铝石、勃姆石也可富集少量Li,当样品中赋锂矿物(高岭石、伊利石(水云母)、三水铝石和勃姆石)都存在时,高w(Li),主量元素w(Al3O2)、w(SiO2)、w(MgO)、w(K2O)、w(TiO2)和w(TFe2O3)与w(Li)的相关性也证实了上述结果.研究区的地球化学比值 CIA、w(Sr)/w(Cu)、u(CaO)/w(MgO)、w(Sr)/w(Ba)、w(V)/w(V+Ni)和 w(La)/w(Y)综合显示,炎热潮湿的古气候下,有机质腐烂形成酸性环境,半封闭海湾环境下,黏土化阶段中等强烈的化学风化程度和适宜的风化暴露剥蚀时间利于铝硅酸盐矿物和硅酸盐矿物化学键断裂,富集Al3+形成高岭石;风化程度过于强烈和长时间的风化剥蚀,继续脱Si富Al形成以硬水铝石为主的铝土矿;弱酸性-弱碱性以及氧化-还原过渡条件下,使得高岭石能最大程度的吸附Li元素.     

The Shear Behavior of Bedding Planes of Weakness Between Two Different Rock Types with High Strength Difference

In this article, the shear behavior of discontinuities caused by bedding planes of weakness between two different rock types with high strength difference is investigated. The effect of roughness and compressive strength of joint wall in such discontinuities are studied. The designed profiles consist of two regular and three irregular artificial joints molded by three types of plaster mortars with different uniaxial compressive strengths. Firstly, it is demonstrated that the shear behavior of discontinuities with different joint wall compressive strengths (JCS) is different from rock joints with identical wall compressive strengths by showing that Barton’s empirical criterion is not appropriate for the former discontinuities. After that, some correlation equations are proposed between the joint roughness coefficient (JRC) parameter and some surface statistical/fractal parameters, and the normal stress range of Barton’s strength criterion is also modified to be used for such discontinuities. Then, a new empirical criterion is proposed for these discontinuities in such a way that a rational function is used instead of JRC log₁₀(JCS/σ _n) as i ₀(σ _c/σ _n)^a/[b + (σ _c/σ _n) ^a ] by satisfying the peak dilation angle boundary conditions under zero and very high normal stress (physical infinite normal stress causing zero peak dilation angle). The proposed criterion has three surface parameters: i ₀, a, and b. The reason for separation of i ₀ from JRC is indicated and the method of its calculation is mentioned based on the literature. The two remaining coefficients (a and b) are discussed in detail and it is shown that a shows a power-law relationship with b, introducing the coefficient c through b = c ^a . Then, it is expressed that a is directly related to discontinuity surface topography. Finally, it is shown that the coefficient c has higher values in irregular profiles in comparison with regular profiles and is dominated by intensity of peak dilation angle reduction (majorly related to the surface irregularity and minorly related to roughness). The coefficient c is to be determined by performing regression analysis on experimental data.     

Evaporation and energy balance estimates over a large inland lake in the Tibet-Himalaya

The process of evaporation from the lake surface is one of the main mechanisms in the energy and water budgets of the lake hydrologic cycle, and an essential component of the water balance especially for inland lakes. In this study, using routine meteorological data as input, a one-layer potential evaporation model was employed to simulate evaporation and energy fluxes over Lake Yamdrok Yumco, the largest high-elevation inland lake in the mountain area of the Tibet-Himalaya in China. Then, the calculation results were compared with the measured values from a big pan evaporator of 20 m² near the lake. The results show that the average annual input radiation flux R ^↓ is 128.2 W m⁻², the lake storage heat flux G is 19.4 W m⁻², the sensible heat flux H is 20.4 W m⁻² and the latent heat flux lE is 107.8 W m⁻². The R ^↓ and G exhibit similar seasonal variations. The lE reaches a maximum in October, lagging nearly 4 months behind the R ^↓ and the G, which indicates the large heat capacity of the lake. The simulated annual evaporation ranges from 1,113.2 to 1,429.1 mm and its mean value is 1,252.5 mm during 1961–2005. The simulated annual evaporation is in good agreement with the measured value, and the measured average lake temperature is as expected when compared with the measured lake surface temperature.     

Wind profiles in the boundary layer over Kharagpur associated with synoptic scale systems

Doppler sodar wind data for the boundary layer over Kharagpur obtained during MONTBLEX-1990 at a height interval of 30 m from surface up to 1500 m have been analysed for the periods when intense synoptic scale disturbances from north Bay of Bengal moved along the eastern end of the monsoon trough. The variation in the vertical wind profile in the lower boundary layer over Kharagpur during the passage of synoptic scale disturbances has been discussed in the paper. The analysis indicates that the mean winds over Kharagpur veered with height in the lower boundary layer near the surface suggesting divergence over Kharagpur when the system lay south/southwest of the station. No such veering has been noticed when the centre of the system lay very close to the station.