台湾海峡西岸重力场变化梯度场图像分析

讨论了应用曲面拟合技术和利用离散的重力观测点的重力变化对特定研究区的重力场变化进行数值模拟以获取该区重力变化空间梯度场图像的具体方法;并应用该方法获得了2004年3月福建安溪3.6级地震前的重力变化空间梯度场图像,发现空间梯度场图像能清晰地反映出地震前重力场群体性变化的差异性;利用此方法对2005年6月~2006年3月间台湾海峡西岸重力场变化进行数值模拟,获得其空间梯度场图像。     

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

The nature of Europa's dark non-ice surface material: Spatially-resolved high spectral resolution spectroscopy from the Keck telescope

We present new 1.45-1.75 μm spectra of Europa's dark non-ice material with a spectral resolution (λ/δλ) of 1200, obtained by using adaptive optics on the Keck telescope to spatially separate the spectrum of the non-ice material from that of the surrounding ice-rich regions. Despite the great increase in spectral resolution over the previous best spectra of the non-ice material, taken with Galileo's near-infrared mapping spectrometer (NIMS) with λ/δλ=66, no new fine-scale spectral structure is revealed. The smoothness of the spectra is inconsistent with available laboratory spectra of crystalline hydrated salts at Europa temperatures, but is more consistent with various combinations of non-crystalline hydrated salts and/or hydrated sulfuric acid, as have been matched to the lower-resolution NIMS spectra.     

Atmospheric radiative transfer simulation for atmospheric correction of remote sensing data

IntroductionThe radiance leaving the earth-atmosphere sys-tem which can be sensed by a satellite borne ra-diometer is the sum of radiation emission fromtheearth surface and each atmospheric level that aretransmittedtothe top of the atmosphere.The radia-tion emissionfromthe earthsurface andthe radianceof each atmospheric level can be separated fromtheradiance at the top the atmospheric level a satellitemeasured.Thus,the earth surface parameters willbe retrieved from the surface radiance after a…     

SPECULATION ON DEVELOPING TOUR RESOURCES OF THE TOURNAMENT ATHLETICS IN HARBIN CITY

The exploitation and utilization of the tour resources of tournament athletics, including skiling, boat sail-ing, archery, ice engraving, snow engraving, has become a new trend of the development of Chinese tourism. Due tothe unique cold climate and superior geographic location, Harbin is a promising city for developing tour resources of tourna-ment athletics. Based on the analysis of the superiority and peculiarity of Harbin, the speculation on development of tour resources of tournament athletics in Harbin is proposed as follows : 1 ) Harbin should develop its special tour resources of tournament athletics associated with needs of market; 2) Harbin should take the advantages of rich resources and devel-op ice and snow entertainment in winter and travel for sight-seeing and spending summer;, 3) the adjustment of the layoutof ice and snow resources should be based on the idea of taking Harbin as the center and all-side opening at the largescale in the way of radiation; 4) tourism should be developed by the combined efforts of various departments to make feasi-ble plan, and the organizers should pay much attention to ensuring the safetv of tourists.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.