An improved TES algorithm based on the corrected ALPHA difference spectrum

Different from visible signals, thermal infrared radiances depend on both temperature and emissivity. It is a key problem for us to separate temperature and emissivity in thermal infrared remote sensing re- search. Another difficulty encountered in the retrieval of surface temperature is the correction of downwelling sky irradiance, because it is closely related to surface emissivity. When emissivity is un- known, the downwelling sky irradiance is difficult to be removed. In this paper, we introduce a correc- tion term of downwelling sky irradiance developed by Li and Becker into Wien’s approximation, to de- rive an improved ALPHA difference spectrum which is independent of temperature, and furthermore develop a correction term to remove the error of Wien’s approximation. Under the support of the above work, attractive features of Alpha derived emissivity method and ASTER TES algorithm are combined together to acquire a new Improved TES algorithm based on Corrected ALPHA Difference Spectrum (ICADS TES). Because a multi-band inversion technique is applied, and the operations of band ratios and differences are included in the algorithm, it can partly remove the influence of atmosphere and noises. Numerical simulation experiments show that for various combinations of atmosphere, land covers and surface temperatures, the algorithm is applicable and stable. Its accuracy for temperature is 0―1.5 K, and that for emissivity is 0―0.015. Compared with current TES algorithms, our method has clear physical meaning, is easy to be implemented, and is applicable for a wide temperature range and surface types. The results are not influenced by the directional characteristic of emissivity. Because ICADS TES does not need the support of a priori information of surface types, it is also not influenced by the accuracy of classification and the problem of mixture pixels. Compared with our former TES algorithm based on corrected Alpha difference spectra (CADS TES), the new algorithm takes the effect of downwelling atmospheric radiation into account. When the quantity of atmosphere radiation can be estimated precisely, the performance of ICADS TES is much better.     

Rapid estimation of soil heavy metal nickel content based on optimized screening of near-infrared spectral bands

In near-infrared spectroscopy,the traditional feature band extraction method has certain limitations.Therefore,a band extraction method named the three-step extraction method was proposed.This method combines characteristic absorption bands and correlation coefficients to select characteristic bands corresponding to various spectral forms and then uses stepwise regression to eliminate meaningless variables.Partial least squares regression(PLSR)and extreme learning machine(ELM)models were used to verify the effect of the band extraction method.Results show that the differential transformation of the spectrum can effectively improve the correlation between the spectrum and nickel(Ni)content.Most correlation coefficients were above 0.7 and approximately 20%higher than those of other transformation methods.The model effect established by the feature variable selection method based on comprehensive spectral transformation is only slightly affected by the spectral transformation form.Infive types of spectral transformation,the RPD values of the proposed method were all within the same level.The RPD values of the PLSR model were concentrated between 1.6 and 1.8,and those of the ELM model were between 2.5 and2.9,indicating that this method is beneficial for extracting more complete spectral features.The combination of the three-step extraction method and ELM algorithm can effectively retain important bands associated with the Ni content of the soil.The model based on the spectral band selected by the three-step extraction method has better prediction ability than the other models.The ELM model of the first-order differential transformation has the best prediction accuracy(RP^2=0.923,RPD=3.634).The research results provide some technical support for monitoring heavy metal content spectrum in local soils.     

Improvement of dark object method in atmospheric correction of hyperspectral remotely sensed data

The existing methods in atmospheric correction of hyperspectral data usually focus on removing the effects of water vapor and other absorptive gases, while this paper mainly studies the method of re- moving the influence of the aerosol and the water vapor simultaneously. Because the hyperspectral data has a larger number of bands, the conventional dark object method cannot be applied to the at- mospheric correction of the hyperspectral data which can be improved, as described in this paper, by adequately making use of spectral characteristics of the hyperspectral data with an iterative correction during the whole process. The effects of the aerosol and water vapor are eliminated at the same time finally. The improved dark object method is used to do the atomospheric correction of the Hyperion data in Yanzhou, Shandong Province as an example. And the result indicates that it can correct the atmospheric influence of the hyperspectral data quickly and remarkably.     

Estimating instantaneous concentration of suspended sediment using acoustic backscatter from an ADV

The measurement of instantaneous sediment concentration remains a challenging task. In this study, a three-step procedure is proposed to estimate instantaneous sediment concentration using acoustic backscatter from Acoustic Doppler Velocimeter (ADV). The influences of acoustic noise and particle diameter on an ADV's performance was first tested in the laboratory, then the three-step procedure was verified based on field measurements in the Zhongxian and Fengjie reaches in the Three Gorges Reservoir. The first step involves reconstructing the backscatter signal time series. Due to contamination from both the noise floor and spurious spikes, the denoising-despiking method was applied instead of the traditional velocity-despiking methods, and this approach performed well based on spectrum ana-lysis. The second step involves calibrating the sediment concentration against the backscatter signal. A linear relation, whose slope and intercept were calibrated to be dependent on particle diameter, is proposed in double logarithmic coordinates. The third step involves calculating the instantaneous sediment concentration using reconstructed instantaneous backscatter based on the proposed relation. The accuracy of the proposed method was evaluated through consideration of the concentration spec-trum and sediment flux, indicating that the proposed three-step procedure is effective for the mea-surement of instantaneous sediment concentration.     

依据密度分布频率特性的重力反演方法研究

Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.     

A cyclic UH model for sand

A simple cyclic elastoplastic constitutive model for sand is proposed based on the UH model for overconsolidated clay. The proposed model has the following features. First, in order to describe the stress-induced anisotropy in sand, a rotational hardening rule is introduced for the evolution of the yield surface axis with development of plastic deviatoric strain in the principal stress space. Second, the relationship between the rotational axis and stress-induced anisotropy is modeled by introducing the slope of rotational axis into the yield function. The fl atness of the yield surface can be determined by the slope of rotational axis. Finally, a revised unifi ed hardening parameter is proposed to incorporate the stress-induced anisotropy. The model capability in describing the cyclic response of sand is verifi ed by comparing the simulations with available test results.     

Multi—scale expression of spatial activity anomalies of earthquakes and its indicative significance on the space and time attributes of strong earthquakes

The noise model based on a trous wavelet algorithm produces a multi-scale expression of image through the combination of wavelet transform and a testing model of statistical significance.This kind of expression not only gives the formation and location of image structure on different scales,but also eliminates the influence of noise Since the algorithm does not need any priori hypotheses,it is suitable for the data with complex structure.The research line is employed in this paper to analyze the spatial activity of earthquake.The method of how to recognize and describe the multi-scale space activity of earthquake is emphatically discussed in this paper.Taking typical sequences in Southwest China as research cases,we systematically study the structure characters of spatial activity of earthquake on different scales.Results show that multi-scale space structure to some extent possesses indicative effect on strong epicenters.And the foreshock anomalies of Songpan seismic sequence also reveal interesting pattem during the spatial-temporal evolvement.     

华北时变重力场离散小波多尺度分解

On the basis of the absolute and relative gravity observations in North China, spatial dynamic variation of regional gravity fields is obtained. A multi-scale decomposition technique is used to separate anomalies at different depths, and give some explanation to gravity variation at different time space scales. Gravity variation trends in North China are improved. Based on this result and the analysis of wavelet power spectrum, the images of the depth of wavelet approximation and detail are obtained. The results obtained are of scientific significance for the deep understanding of potential seismic risk in North China from gravity variations in different time space scales.     

A simplified methd of evaluating the seismic performance of buildings

This paper presents a simplified method of evaluating the seismic performance of buildings. The proposed method is based on the transformation of a multiple degree of freedom (MDOF) system to an equivalent single degree of freedom (SDOF) system using a simple and intuitive process. The proposed method is intended for evaluating the seismic performance of the buildings at the intermediate stages in design, while a rigorous method would be applied to the final design. The performance of the method is evaluated using a series of buildings which are assumed to be located in Victoria in western Canada, and designed based on the upcoming version of the National Building Code of Canada which is due to be published in 2005. To resist lateral loads, some of these buildings contain reinforced concrete moment resisting frames,while others contain reinforced concrete shear walls. Each building model has been subjected to a set of site-specific seismic spectrum compatible ground motion records, and the response has been determined using the proposed method and the general method for MDOF systems. The results from the study indicate that the proposed method can serve as a useful tool for evaluation of seismic performance of buildings, and carrying out performance based design.     

A comparison of atmospheric CO2 concentration GOSAT-based observations and model simulations

Satellite observations of atmospheric CO2 are able to truly capture the variation of global and regional CO2 concentration.The model simulations based on atmospheric transport models can also assess variations of atmospheric CO2 concentrations in a continuous space and time,which is one of approaches for qualitatively and quantitatively studying the atmospheric transport mechanism and spatio-temporal variation of atmospheric CO2 in a global scale.Satellite observations and model simulations of CO2 offer us two different approaches to understand the atmospheric CO2.However,the difference between them has not been comprehensively compared and assessed for revealing the global and regional features of atmospheric CO2.In this study,we compared and assessed the spatio-temporal variation of atmospheric CO2 using two datasets of the column-averaged dry air mole fractions of atmospheric CO2(XCO2)in a year from June 2009 to May 2010,respectively from GOSAT retrievals(V02.xx)and from Goddard Earth Observing System-Chemistry(GEOS-Chem),which is a global 3-D chemistry transport model.In addition to the global comparison,we further compared and analyzed the difference of CO2 between the China land region and the United States(US)land region from two datasets,and demonstrated the reasonability and uncertainty of satellite observations and model simulations.The results show that the XCO2 retrieved from GOSAT is globally lower than GEOS-Chem model simulation by 2 ppm on average,which is close to the validation conclusion for GOSAT by ground measures.This difference of XCO2 between the two datasets,however,changes with the different regions.In China land region,the difference is large,from 0.6 to 5.6 ppm,whereas it is 1.6 to 3.7 ppm in the global land region and 1.4 to 2.7 ppm in the US land region.The goodness of fit test between the two datasets is 0.81 in the US land region,which is higher than that in the global land region(0.67)and China land region(0.68).The analysis results further indicate that the inconsistency of CO2concentration between satellite observations and model simulations in China is larger than that in the US and the globe.This inconsistency is related to the GOSAT retrieval error of CO2 caused by the interference among input parameters of satellite retrieval algorithm,and the uncertainty of driving parameters in GEOS-Chem model.     

Estimation of tropical cyclone parameters and wind fields from SAR images

The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.     

Study on the atmospheric heat engine efficiency and heat source characteristics of the Qinghai-Tibet Plateau in summer

There are many types of atmospheric heat engines in land-air systems. The accurate definition, calculation and interpretation of the efficiency of atmospheric heat engines are key to understanding energy transfer and transformation of landair systems. The atmosphere over the Qinghai-Tibet Plateau(QTP) in summer can be regarded as a positive heat engine. The study of the heat engine efficiency is helpful to better understand land-air interaction and thermal-dynamic processes on the QTP.It also pr...     

Estimation of clear-sky longwave downward radiation from HJ-1B thermal data

To satisfy the requirement of surface energy budget research on the meso-and micro-scale,a parameterization is developed to calculate high spatial resolution,clear-sky downward longwave radiation(DLR)from HJ-1B thermal data.The DLR algorithm is established based on extensive radiative transfer simulation and statistical analysis.To address the problem that HJ-1B has a single thermal channel and lacks atmospheric information,the brightness temperature of HJ-1B and water vapor content are used in the algorithm.An accuracy evaluation and error analysis for the algorithm is conducted using a simulated radiation dataset.The result shows that the algorithm performs well in most circumstances,but there is obvious underestimation when water vapor content is greater than 4g/cm2 .Error analysis indicates the accuracy of estimated DLRs is affected by uncertainties in input parameters,including water vapor content and top-of-atmosphere radiance.It is also affected by the difference between ground and near-surface air temperature.The algorithm is applied to actual HJ-1B data,and validated by ground data from six stations in the Heihe River and Haihe River basins.The estimated DLRs have good consistency with measured data except at Huazhaizi,and root mean square errors at most sites are around 20W/m2 ,which is slightly better than the result of MODIS. There is significant overestimation of DLR at Huazhaizi during summer,which is mainly produced by the large ground-air temperature difference.A correction process based on temperature difference is proposed and applied at Huazhaizi.The result shows that the positive bias is largely diminished after correction.     

New supervised learning classifiers for structural damage diagnosis using time series features from a new feature extraction technique

The motivation for this article is to propose new damage classifiers based on a supervised learning problem for locating and quantifying damage. A new feature extraction approach using time series analysis is introduced to extract damagesensitive features from auto-regressive models. This approach sets out to improve current feature extraction techniques in the context of time series modeling. The coefficients and residuals of the AR model obtained from the proposed approach are selected as the ...     

THE ACTIVITY OF A FAULT ZONE INFERRED FROM MICROSTRUCTURAL STUDY OF DEFORMED ROCKS

The neotectonic activity of a minor fault zone on the Southwest margin of Yangzi platform has been studied based on analysis and comparison of the microstructural features of deformed rocks within and outside the fault zone.According to thermoluminesccne dating it is proposed that after Pleistocene there would be more than one active period of the fault.The strees of the main period was up to 700 bars.The strain and stress of the rocks within the fault during the faulting process were significantly greater than those in the rocks outside the fault     

27.3-day and 13.6-day atmospheric tide

An analysis of time variations in the earth's length of day (LOD) for 25 years (1973-1998) versus at- mospheric circulation changes and lunar phase is presented. It is found that, on the average, there is a 27.3-day and 13.6-day period oscillation in global zonal wind speed, atmospheric geopotential height, and LOD following alternating changes in lunar phase. Every 5-9 days (6.8 days on average), the fields of global atmospheric zonal wind and geopotential height and LOD undergo a sudden change in rela- tion to a change in lunar declination. The observed atmospheric oscillation with this time period may be viewed as a type of atmospheric tide. Ten atmospheric tidal cases have been analyzed by comparing changes in LOD, global zonal wind speed and atmospheric geopotential height versus change in lunar declination. Taken together these cases reveal prominent 27.3-day and 13.6-day tides. The lunar forcing on the earth's atmosphere is great and obvious changes occur in global fields of zonal wind speed and atmospheric geopotential height over the equatorial and low latitude areas. The driving force for the 27.3-day and 13.6-day atmospheric tides is the periodic change in lunar forcing during the moon's revolution around the earth. When the moon is located on the celestial equator the lunar declination equals zero and the lunar tidal forcing on the atmosphere reaches its maximum, at this time the global zonal wind speed increases and the earth's rotation rate decreases and LOD increases. Conversely, when the moon reaches its most northern or southern positions the lunar declination is maximized, lunar tidal forcing decreases, global zonal wind speed decreases, earth's rotation rate increases and LOD decreases. 27.3-day and 13.6-day period atmospheric tides deserve deeper study. Lunar tidal forcing should be considered in models of atmospheric circulation and in short and medium range weather forecasting.     

Early Identification of the Jiangdingya Landslide of Zhouqu Based on SBAS-InSAR Technology

SBAS-InSAR technology is characterized by the advantages of reducing the influence of terrain-simulation error, time-space decorrelation, atmospheric error, thereby improving the reliability of surface-deformation monitoring. This paper studies the early landslide identification method based on SBAS-InSAR technology. Selecting the Jiangdingya landslide area in Zhouqu County, Gansu Province as the research area, 84 ascending-orbit Sentinel-1A SAR images from 2015 to 2019 are collected. In addition, using SBAS-InSAR technology, the rate and time-series results of surface deformation of the landslide area in Jiangdingya during this period are extracted, and potential landslides are identified. The results show that the early landslide identification method based on SBAS-InSAR technology is highly feasible and is a better tool for identifying potential landslides in large areas.     

Structural modal parameter identification and damage diagnosis based on Hilbert-Huang transform

Traditional modal parameter identifi cation methods have many disadvantages,especially when used for processing nonlinear and non-stationary signals.In addition,they are usually not able to accurately identify the damping ratio and damage.In this study,methods based on the Hilbert-Huang transform(HHT) are investigated for structural modal parameter identifi cation and damage diagnosis.First,mirror extension and prediction via a radial basis function(RBF) neural network are used to restrain the troublesome end-effect issue in empirical mode decomposition(EMD),which is a crucial part of HHT.Then,the approaches based on HHT combined with other techniques,such as the random decrement technique(RDT),natural excitation technique(NExT) and stochastic subspace identifi cation(SSI),are proposed to identify modal parameters of structures.Furthermore,a damage diagnosis method based on the HHT is also proposed.Time-varying instantaneous frequency and instantaneous energy are used to identify the damage evolution of the structure.The relative amplitude of the Hilbert marginal spectrum is used to identify the damage location of the structure.Finally,acceleration records at gauge points from shaking table testing of a 12-story reinforced concrete frame model are taken to validate the proposed approaches.The results show that the proposed approaches based on HHT for modal parameter identifi cation and damage diagnosis are reliable and practical.     

Three-dimensional decomposition method of global atmospheric circulation

By adopting the idea of three-dimensional Walker, Hadley and Rossby stream functions, the global atmospheric circulation can be considered as the sum of three stream functions from a global per- spective. Therefore, a mathematical model of three-dimensional decomposition of global atmospheric circulation is proposed and the existence and uniqueness of the model are proved. Besides, the model includes a numerical method leading to no truncation error in the discrete three-dimensional grid points. Results also show that the three-dimensional stream functions exist and are unique for a given velocity field. The mathematical model shows the generalized form of three-dimensional stream func- tions equal to the velocity field in representing the features of atmospheric motion. Besides, the vertical velocity calculated through the model can represent the main characteristics of the vertical motion. In sum, the three-dimensional decomposition of atmospheric circulation is convenient for the further in- vestigation of the features of global atmospheric motions.     

Earthquake accelerogram simulation with statistical law of evolutionary power spectrum

By using the technique for evolutionary power spectrum proposed by Nakayama and with reference to the Kameda formula, an evolutionary spectrum prediction model for given earthquake magnitude and distance is established based on the 80 near-source acceleration records at rock surface with large magnitude from the ground motion database of western U.S.. Then a new iteration method is developed for generation of random accelerograms non-stationary both in amplitude and frequency which are compatible with target evolutionary spectrum. The phase spectra of those simulated accelerograms are also non-stationary in time and frequency domains since the interaction between amplitude and phase angle has been considered during the generation. Furthermore, the sign of the phase spectrum increment is identified to accelerate the iteration. With the proposed statistical model for predicting evolutionary power spectra and the new method for generating compatible time history, the artificial random earthquake accelerograms non-stationary both in amplitude and frequency for certain magnitude and distance can be provided.