Application of the Fourier method to the numerical solution of moving boundary problem in heat conduction

An algorithm for the solution of a nonlinear problem of phase boundary movement and evolution of temperature distribution due to the perturbation in the basal heat flux has been discussed. The reduction of the problem to a system of nonlinear ordinary differential equations with the help of a Fourier series method leads to a stiff system. This stiffness is taken care of by the use of a modified Euler’s method. Various cases of basal heat flow variation have been considered to show the performance and stability of the technique for such a nonlinear system. The first case of step-wise function is taken to analyse the performance of the technique, and the study has been extended to other general cases of linear increase, periodic variation, and box and triangular function type variations in the heat flux. In the step-wise case the phase boundary attains a constant position rapidly if the supplied heat flux is sufficiently large. The effect of periodicity in the heat flow is clearly depicted in the phase boundary movement, where the phase boundary oscillates about the mean position at large times. The absence of any constant level in the case of linear increase in heat flux is due to a very large value of heat flux. In the cases of box car and triangular heat flux the boundary starts moving downward after the cessation of excess heat flux but does not immediately return to its original preperturbation state, instead approaches it at large times. This technique may be applied to more general cases of heat flow variation.     

Determination of paleoheat-flux using multiple thermal indicator tomography

Multiple thermal indicator tomography is a non-linear inversion procedure which permits simultaneous determination of geological, chemical kinetic, and thermal parameters in the burial and thermal histories of sediments. The procedure is based on inverting measurements with depth in a well of multiple thermal indicators, each of which must have experienced the same thermal history. For two wells from diverse geological settings, tomography is used to assess the paleoheat-flux variations with time most consistent with all the thermal indicator data. Resolution and sensitivity of the results are also addressed. For the case history of well X in the North Sea, with both vitrinite reflectance and sterane isomer thermal indicators, a generally decreasing heat flux from past to present is obtained roughly halving every 50 Ma. The activation energy for the sterane isomer is also determined simultaneously to be 42±20 kJ/mole. For the second case history of the Irma Mclean well in the onshore Gulf of Mexico, data from four thermal indicators are available: vitrinite reflectance, sterane and hopane isomers, and sterane aromaticity. A general decrease in paleoheat-flux from amount 2.8 heat flux units (HFU) at 100 MaBP to 1.25 HFU at the present day is obtained, the activation energy for the sterane isomer is simultaneously determined to be around 50±20 kJ/mole, while the corresponding values are around 80±20 kJ/mole for the hopane isomer, and around 65±20 kJ/mole for sterane aromatization.     

SEASONAL CHANGES OF SOIL MOISTURE INDUCED BY INCREASED CO2 AS SIMULATED BY THE ATMOSPHERIC GENERAL CIRCULATION AND MIXED LAYER OCEAN MODEL

The atmospheric general circulation model coupled to the mixed layer ocean model has been used to simulate thechanges of the global soil moisture.Comparing the simulated results with observations,it is shown that the model is ca-pable of doing sensitive experiments about the carbon dioxide change.The 2×CO_2/1×CO_2 comparison shows that there are the obvious changes of the soil moisture in the global forfour seasons.There are the wet soil moisture in the lower latitudes of both hemispheres and dry soil moisture in the mid-dle latitudes of both hemispheres for four seasons.The dry soil moisture in summer and wet in other seasons are foundin the northern higher latitudes.The analyses of the physical feedbacks responsible for the CO_2-induced changes of soil moisture show that the bud-gets of the surface water and heat are the important factors.     

INFLUENCE OF SOIL MOISTURE AND VEGETATION ON CLIMATE CHANGES INDUCED BY THERMAL FORCING*

A land-process scheme has been incorporated in a vertical one-dimensional time-dependent atmospheric model and numerical experiments have been performed with the coupled model to examine influences of soil wetness and vegetation on climate changes associated to thermal forcing.It is showed that response of land-surface temperature to the thermal forcing becomes small with increase of soil water content and vegetation cover.Furthermore,the response is more obvious in arid climate region than in humid one.The result also shows that there exist two patterns of corresponding relation between variations in air temperature and humidity on the land surface in response to hydrologic and thermal focing.     

若羌地区降水量变动的特征

采用若羌36年(1953——1988年)的降水资料讨论了该地区降水特征,通过分析得到:1.若羌降水量的均值、方差、相关系数和降水变率等气候统计量的特征;2.若羌降水量有增多的趋势。本世纪50年代年降水16.2mm,60年代18.0mm,70年代20.3mm,而80年代达们40.2mm,为50年代的2.5倍;3.1981及1988年7月份大暴雨,说明干旱区降水量年际变化大,并从中得到有关水汽来源及其水汽输送路径的认识。     

THE CHARACTERISTICS OF TURBULENT TRANSPORT IN THE SURFACE LAYER OVER EASTERN TIBET DURING SUMMER 1986

The turbulent fluxes for sensible and latent heat and momentum are computed and analyses are carried outabout the factors in terms of the fluxes,with the profile methods,based on the data from the Tibetan(Xi-zang)Plateau Meteorological Experiment in 1986(TIPMEX-86).It is shown that the fluxes of various kindshave evident diurnal variation,and each decade mean diurnal variation is quite different from others.Thesensible heat flux is about 2/3 less in July than in June.The results indicate that the averaged drag coefficient,C_d,and the averaged bulk transfer coefficient of sensible heat,C_h,are 0.0052 and 0.0075 respectively,for theperiod 13 to 28 in June for Nagqu.But for Lhasa,the mean C_d is 0.0056,and the mean C_h is 0.0085,for the period of June 11 to July 20.It is found that C_d and C_h are not only the function of wind velocity,but also influenced by stability conditions and wind and temperature gradients.     

Entrainment of circumpolar water in the Indian Ocean region of the Antarctic

The net influx of the circumpolar water on the western (approximately along 10°E) and eastern (approximately 115°E) boundaries of the Indian Ocean, adopting the method of Montgomery and Stroup is computed on bivariate distribution of potential thermosteric anomaly and salinity to identify the characteristics of the flux. The zonal flux at both the boundaries indicates an alternate strong easterly and westerly flow between 36°S and 45°S, south of which the flow is mainly easterly but weak up to 56°S. At the western boundary the easterly flow is 146 Sv and westerly is 98.07 Sv, while at the eastern boundary (115°E) the corresponding fluxes are 123.46 Sv and 27.20 Sv respectively, indicating a net outflux of 48.33 Sv. This water should have been accounted by the melting of ice and influx of the Equatorial Pacific Ocean Water.     

东亚夏季风水汽输送带及其对中国大暴雨与洪涝灾害的影响

东亚夏季风每年给中国东部地区带来充沛的降水，是中国水资源的主要来源，同时也常常给中国造成严重的洪涝灾害。东亚夏季风水汽输送的强度、影响范围和持续性在极端暴雨过程中起着关键的作用。这支夏季风气流的水汽输送带可称为东亚季风水汽输送带，与国际上近期提出的"大气河"概念相近，但又不完全相同。东亚夏季风水汽输送带是东亚夏季风最具地区性的特征，也是东亚地区夏季大暴雨和洪涝的制造者。本文根据近百年来的资料，综合评述了东亚夏季风水汽输送带的特征和形成原因，并以海河、黄河、淮河与长江近百年最强的5次持续大暴雨过程为例，分析了季风水汽输送带的重要作用。最后，提出气候变暖可以通过4个方面影响全球水循环，包括气候变暖后大气可容纳更多的水汽、大气环流发生变化、辐射强迫改变以及气溶胶影响的区域性等，这些变化都会对季风水汽输送带产生重要影响。     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA*

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urban area in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered structure of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heights and morning boundary layer development are combined with surface eddy correlation measurements of kinematic heat and moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is presented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed during the transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

Transformation of Irradiance Measurements into Spectral Actinic Flux for Photolysis Rates Determination

The actinic flux is the only radiometric quantity suitable for photolysisfrequency determination. It is derived from solar spectral irradiancemeasurements performed by a portable spectroradiometer in the 300–850nm wavelength range. The spectral irradiance is first divided into a directpart and a diffuse part, according to the atmospheric conditions, and thenconverted into the corresponding actinic flux quantity. As an intermediary,the ratio of diffuse actinic flux to diffuseirradiance is calculated by the spherical harmonics radiative code withrespect to wavelength, solar zenith angle, surface albedo, and aerosolproperties.The results of extensive sensitivity studies of asa function of the main atmospheric parameters are discussed and lead to theconclusion that aerosol optical depth is the major critical value for aprecise determination. The global algorithm totransform spectral irradiance into actinic flux is finally applied forphotodissociation rate calculations by convolution of the obtained actinicflux spectra with the absorption cross-sections and quantum yields of themolecule studied. Photolysis rates of different atmospheric photooxidants havebeen measured with this spectroradiometric method during the summers of 1993and 1994 in Brittany and in Portugal. The ozone and nitrogen dioxidephotodissociation rates obtained present a good agreement with thej(O₃) and j(NO₂) actinometerresults, for the same experimental conditions in Brittany.