Improving the Accuracy of Vegetation Index Retrieval for Biomass by Combining Ground-UAV Hyperspectral Data-A New Method for Inner Mongolia Typical Grasslands

Grassland biomass is an important parameter of grassland ecosystems.The complexity of the grassland canopy vegetation spectrum makes the long-term assessment of grassland growth a challenge.Few studies have explored the original spectral information of typical grasslands in Inner Mongolia and examined the influence of spectral information on aboveground biomass(AGB)estimation.In order to improve the accuracy of vegetation index inversion of grassland AGB,this study combined ground and Unmanned Aerial Vehicle(UAV)remote sensing technology and screened sensitive bands through ground hyperspectral data transformation and correlation analysis.The narrow band vegetation indices were calculated,and ground and airborne hyperspectral inversion models were established.Finally,the accuracy of the model was verified.The results showed that:(1)The vegetation indices constructed based on the ASD FieldSpec 4 and the UAV were significantly correlated with the dry and fresh weight of AGB.(2)The comparison between measured R^(2) with the prediction R^(2) indicated that the accuracy of the model was the best when using the Soil-Adjusted Vegetation Index(SAVI)as the independent variable in the analysis of AGB(fresh weight/dry weight)and four narrow-band vegetation indices.The SAVI vegetation index showed better applicability for biomass monitoring in typical grassland areas of Inner Mongolia.(3)The obtained ground and airborne hyperspectral data with the optimal vegetation index suggested that the dry weight of AGB has the best fitting effect with airborne hyperspectral data,where y=17.962e^(4.672x),the fitting R^(2) was 0.542,the prediction R^(2)was 0.424,and RMSE and REE were 57.03 and 0.65,respectively.Therefore,established vegetation indices by screening sensitive bands through hyperspectral feature analysis can significantly improve the inversion accuracy of typical grassland biomass in Inner Mongolia.Compared with ground monitoring,airborne hyperspectral monitoring better reflects the inversion of actual surface biomass.It provides a reliable modeling framework for grassland AGB monitoring and scientific and technological support for grazing management.

A Hybrid Compensation Scheme for the Input Rate-Dependent Hysteresis of the Piezoelectric Ceramic Actuators

A hybrid compensation scheme for piezoelectric ceramic actuators(PEAs)is proposed.In the hybrid compensation scheme,the input rate-dependent hysteresis characteristics of the PEAs are compensated.The feedforward controller is a novel input rate-dependent neural network hysteresis inverse model,while the feedback controller is a proportion integration differentiation(PID)controller.In the proposed inverse model,an input ratedependent auxiliary inverse operator(RAIO)and output of the hysteresis construct the expanded input space(EIS)of the inverse model which transforms the hysteresis inverse with multi-valued mapping into single-valued mapping,and the wiping-out,rate-dependent and continuous properties of the RAIO are analyzed in theories.Based on the EIS method,a hysteresis neural network inverse model,namely the dynamic back propagation neural network(DBPNN)model,is established.Moreover,a hybrid compensation scheme for the PEAs is designed to compensate for the hysteresis.Finally,the proposed method,the conventional PID controller and the hybrid controller with the modified input rate-dependent Prandtl-Ishlinskii(MRPI)model are all applied in the experimental platform.Experimental results show that the proposed method has obvious superiorities in the performance of the system.

The sensitivity of the array resistivity log to mud filtrate invasion and its primary five-parameter inversion for improved oil water recognition

In order to improve reservoir fluid recognition, the sensitivity of array resistivity response to the difference of the invasion properties in both oil-bearing layers and water layers is analyzed. Then the primary inversion is carried out based on the array resistivity log. The mud invasion process is numerically simulated based on the oil-water flow equation and water convection diffusion equation. The results show that the radial resistivity of a fresh mud-invaded oil-bearing layer presents complex distribution characteristics, such as nonlinear increase, increasing to decreasing and low resistivity annulus, and the resistive invasion profile of a water layer is monotonic. Under specific conditions, array resistivity log can reflect these changes and the array induction log is more sensitive. Nevertheless, due to the effect of factors like large invasion depth, reservoir physical and oil-bearing properties, the measured apparent resistivity may differ greatly from the actual mud filtrate invasion profile in an oil-bearing layer. We proposed a five-parameter formation model to simulate the complex resistivity distribution of fresh mud-invaded formation. Then, based on the principle of non-linear least squares, the measured array resistivity log is used for inversion with the Marquardt method. It is demonstrated that the inverted resistivity is typically non-monotonic in oil-bearing layers and is monotonic in water layers. Processing of some field data shows that this is helpful in achieving efficient reservoir fluid recognition.

Experimental Research on Quantitative Inversion Models of Suspended Sediment Concentration Using Remote Sensing Technology

Research on quantitative models of suspended sediment concentration （SSC） using remote sensing technology is very important to understand the scouting and siltation variation in harbors and water channels. Based onlaboratory study of the relationship between different suspended sediment concentrations and reflectance spectra measured synchronously, quantitative inversion models of SSC based on single factor, band ratio and sediment parameter were developed, which provides an effective method to retrieve the SSC from satellite images. Results show that the bl （430-500nm） and b3 （670-735nm） are the optimal wavelengths for the estimation of lower SSC and the b4 （780-835nm） is the optimal wavelength to estimate the higher SSC. Furthermore the band ratio B2/B3 can be used to simulate the variation of lower SSC better and the B4/B1 to estimate the higher SSC accurately. Also the inversion models developed by sediment parameters of higher and lower SSCs can get a relatively higher accuracy than the single factor and band ratio models.

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.

3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions

It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.

Model fitting of the seasonal population dynamics of the soybean aphid, Aphis glycines Matsumura, in the field

The soybean aphid, Aphis glycines Matsumura(Hemiptera: Aphididae), is one of the greatest threats to soybean production, and both trend analysis and periodic analysis of its population dynamics are important for integrated pest management(IPM). Based on systematically investigating soybean aphid populations in the field from 2018 to 2020, this study adopted the inverse logistic model for the first time, and combined it with the classical logistic model to describe the changes in seasonal population abundance from colonization to extinction in the field. Then, the increasing and decreasing phases of the population fluctuation were divided by calculating the inflection points of the models, which exhibited distinct seasonal trends of the soybean aphid populations in each year. In addition, multifactor logistic models were then established for the first time, in which the abundance of soybean aphids in the field changed with time and relevant environmental conditions. This model enabled the prediction of instantaneous aphid abundance at a given time based on relevant meteorological data. Taken as a whole, the successful approaches implemented in this study could be used to build a theoretical framework for practical IPM strategies for controlling soybean aphids.

Machine Learning-based Inverse Model for Few-Mode Fiber Designs

The medium for next-generation communication is considered as fiber for fast,secure communication and switching capability.Mode division and space division multiplexing provide an excellent switching capability with high data transmission rate.In this work,the authors have approached an inverse modeling technique using regression-based machine learning to design a weakly coupled few-mode fiber for facilitating mode division multiplexing.The technique is adapted to predict the accurate profile parameters for the proposed few-mode fiber to obtain the maximum number of modes.It is for a three-ring-core few-mode fiber for guiding five,ten,fifteen,and twenty modes.Three types of regression models namely ordinary least-square linear multi-output regression,k-nearest neighbors of multi-output regression,and ID3 algorithm-based decision trees for multi-output regression are used for predicting the multiple profile parameters.It is observed that the ID3-based decision tree for multioutput regression is the robust,highly-accurate machine learning model for fast modeling of FMFs.The proposed fiber claims to be an efficient candidate for the next-generation 5G and 6G backhaul networks using mode division multiplexing.

Nonlinear direct data-driven control for UAV formation flight system

This paper proposes the nonlinear direct data-driven control from theoretical analysis and practical engineering,i.e.,unmanned aerial vehicle(UAV)formation flight system.Firstly,from the theoretical point of view,consider one nonlinear closedloop system with a nonlinear plant and nonlinear feed-forward controller simultaneously.To avoid the complex identification process for that nonlinear plant,a nonlinear direct data-driven control strategy is proposed to design that nonlinear feed-forward controller only through the input-output measured data sequence directly,whose detailed explicit forms are model inverse method and approximated analysis method.Secondly,from the practical point of view,after reviewing the UAV formation flight system,nonlinear direct data-driven control is applied in designing the formation controller,so that the followers can track the leader’s desired trajectory during one small time instant only through solving one data fitting problem.Since most natural phenomena have nonlinear properties,the direct method must be the better one.Corresponding system identification and control algorithms are required to be proposed for those nonlinear systems,and the direct nonlinear controller design is the purpose of this paper.

Genetic Algorithm-Based Estimation of Nonlinear Transducer

This paper describes an innovative, genetic algorithm based inverse model of nonlinear transducer. In the inverse modeling, using a genetic algorithm, the unknown coefficients of the model are estimated accurately. The simulation results indicate that this technique provides greater flexibility and suitability than the existing methods. It is very easy to modify the nonlinear transducer on line. Thus the method improves the transducer's accuracy. With the help of genetic algorithm (GA), the model coefficients' training are less likely to be trapped in local minima than traditional gradient based search algorithms.

Similarity measure of sedimentary successions and its application in inverse stratigraphic modeling

This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.

Spectrum Characteristics of Cotton Canopy Infected with Verticillium Wilt and Applications

Hyper spectrum remote sensing with fine spectrum information is an efficient method to estimate the verticillium wilt of cotton. The research was conducted in Xinjiang, the largest cotton plant region of China, by using the data which were collected both by canopy spectrum infected with verticillium wilt and severity level （SL） in the year 2005-2006. The quantitative correlation was analyzed between SL and canopy of reflectance spectrum or derivative spectrum reflectance. The results indicated that spectrum characteristics of cotton canopy infected with verticillium wilt changed regularly with the increase of SL in different periods and varieties, Spectrum reflectance increased in the visible light region （620-700 nm） with the increase of the SL, which inverted in near-infrared region and was extremely significant in the region of （780-1 300 nm）. When SL attained b2 （DI = 25）, cotton canopy infected with verticillium wilt was used as a watershed and diagnosed index in the beginning stages of the disease. The results also indicated that there were marked different characteristics of the first derivative spectrum in these SL, it changed significantly in the red edge ranges （680-760 nm） with different SL, i.e., red edge swing decreased, and red edge position equally moved to the blue. In this study 1 001-1 110 nm and 1 205- 1 320 nm were selected out as sensitive bands for SL of canopy. Inversion models established for estimating cotton canopy infected with verticillium wilt reached the most significant level. Finally, the different spectrum characteristics of cotton canopy infected with verticillium wilt were marked, some inversion models were established, which could estimate SL of canopy infected with verticillium wilt. The best recognized model was the first derivative spectra at （FD 731 nm- FD 1317 nm）, and it might be used to forecast the position of cotton canopy infected with verticillium wilt quantitatively.

Retrieval of leaf biochemical properties by inversed PROSPECT model and hyperspectral indices:an application to Populus euphratica polymorphic leaves

Leaf biochemical properties have been widely assessed using hyperspectral reflectance information by inversion of PROSPECT model or by using hyperspectral indices, but few studies have focused on arid ecosystems. As a dominant species of riparian ecosystems in arid lands, Populus euphratica Oliv. is an unusual tree species with polymorphic leaves along the vertical profile of canopy corresponding to different growth stages. In this study, we evaluated both the inversed PROSPECT model and hyperspectral indices for estimating biochemical properties of P. euphratica leaves. Both the shapes and biochemical properties of P. euphratica leaves were found to change with the heights from ground surface. The results indicated that the model inversion calibrated for each leaf shape performed much better than the model calibrated for all leaf shapes, and also better than hyperspectral indices. Similar results were obtained for estimations of equivalent water thickness （EWT） and leaf mass per area （LMA）. Hyperspectral indices identified in this study for estimating these leaf properties had root mean square error （RMSE） and R2 values between those obtained with the two calibration strategies using the inversed PROSPECT model. Hence, the inversed PROSPECT model can be applied to estimate leaf biochemical properties in arid ecosystems, but the calibration to the model requires special attention.

Algebraic Method‑Based Point‑to‑Point Trajectory Planning of an Under‑Constrained Cable‑Suspended Parallel Robot with Variable Angle and Height Cable Mast

To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an under-constrained cable-suspended parallel robot(UCPR)with variable angle and height cable mast as described in this paper.The end-effector of the UCPR with three cables can achieve three translational degrees of freedom(DOFs).The inverse kinematic and dynamic modeling of the UCPR considering the angle and height of cable mast are completed.The motion trajectory of the end-effector comprising six segments is given.The connection points of the trajectory segments(except for point P3 in the X direction)are devised to have zero instantaneous velocities,which ensure that the acceleration has continuity and the planned acceleration curve achieves smooth transition.The trajectory is respectively planned using three algebraic methods,including fifth degree polynomial,cycloid trajectory,and double-S velocity curve.The results indicate that the trajectory planned by fifth degree polynomial method is much closer to the given trajectory of the end-effector.Numerical simulation and experiments are accomplished for the given trajectory based on fifth degree polynomial planning.At the points where the velocity suddenly changes,the length and tension variation curves of the planned and unplanned three cables are compared and analyzed.The OptiTrack motion capture system is adopted to track the end-effector of the UCPR during the experiment.The effectiveness and feasibility of fifth degree polynomial planning are validated.

Coseismic and postseismic slip ruptures for 2015Mw 6.4 Pishan earthquake constrained by static GPS solutions

On 3 July 2015, a Mw 6.4 earthquake occurred on a blind fault struck Pishan, Xinjiang,China. By combining Crustal Movement Observation Network of China（CMONOC） and other Static Global Positioning System（GPS） sites surrounding Pishan region, it provides a rare chance for us to constrain the slip rupture for such a moderate event. The maximum displacement is up to 12 cm, 2 cm for coseismic and postseismic deformation, respectively,and both the deformation patterns show a same direction moving northeastward. With rectangular dislocation model, a magnitude of Mw6.48, Mw6.3 is calculated based on coseismic, postseismic deformation respectively. Our result indicates the western Kunlun range is still moving toward Tarim Basin followed by an obvious postseismic slip associated with this earthquake. To determine a more reasonable model for postseismic deformation, a longer GPS dataset will be needed.

Direct shear tests on cemented paste backfill-rock wall and cemented paste backfill-backfill interfaces

Even though a large number of large-scale arch dams with height larger than 200 m have been built in the world, the transient groundwater flow behaviors and the seepage control effects in the dam foundations under difficult geological conditions are rarely reported. This paper presents a case study on the transient groundwater flow behaviors in the rock foundation of Jinping I double-curvature arch dam, the world＇s highest dam of this type to date that has been completed. Taking into account the geological settings at the site, an inverse modeling technique utilizing the time series measurements of both hydraulic head and discharge was adopted to back-calculate the permeability of the foundation rocks,which effectively improves the uniqueness and reliability of the inverse modeling results. The transient seepage flow in the dam foundation during the reservoir impounding was then modeled with a parabolic variational inequality（PVI） method. The distribution of pore water pressure, the amount of leakage, and the performance of the seepage control system in the dam foundation during the entire impounding process were finally illustrated with the numerical results.

A new method for multi-point pollution source identification

Rapid and accurate identification of the characteristics(source location,number,and intensity)of pollution sources is essential for emergency assessment of contamination events.Compared with single-point source iden-tification,the reconstruction of multiple sources is more challenging.In this study,a two-step inversion method is proposed for multi-point pollution source reconstruction from limited measurements with the number of sources unknown.The applicability of the proposed method is validated with a set of synthetic experiments correspond-ing to one-,two-,and three-point pollution sources.The results show that the number and locations of pollution sources are retrieved exactly the same as prescribed,and the source intensities are estimated with negligible errors.The algorithm exhibits good performance in single-and multi-point pollution source identification,and its accuracy and efficiency of identification do not deteriorate with the increase in the number of sources.Some limitations of the algorithm,together with its capabilities,are also discussed in this paper.

Quantitative description of infrared radiation characteristics of preflawed sandstone during fracturing process

Previous studies show that infrared radiation temperature(IRT)abnormalities are always accompanied by the crack development in rocks under external loads.In this context,experiments were conducted on preflawed sandstone to investigate the infrared radiation characteristics during failure process.Two indicators were defined herein,i.e.coefficient of variation of IRT(CVIRT)and skewness of IRT(SIRT).The regression analysis shows that the IRT probability distributions during loading process fit the Gaussian model.The variations in the CVIRT are characterized by four stages:primary stage,steady stage,accelerating stage and post-peak stage.Besides,the variations in the SIRT are divided into three stages:primary stage,steady stage and failure and post-peak stage.The precursor point for preflawed rock failure is identified based on the CVIRTetime curve,with average precursor point of 83%of the peak stress.Compared with other IRT indicators,the proposed two IRT indicators have higher sensitivity to IRT abnormalities during failure process.Furthermore,the connection between the IRT indicators and the rock fracturing was investigated to interpret the IRT indicator abnormalities.Based on the Verhulst inverse function,a new quantitative model was presented to describe the primary stage,steady stage and accelerating stage of the CVIRTetime curve.The results obtained in this study can provide early-warning information for rock failure prediction.

Geophysical and Production Data History Matching Based on Ensemble Smoother with Multiple Data Assimilation

The Ensemble Kalman Filter(EnKF),as the most popular sequential data assimilation algorithm for history matching,has the intrinsic problem of high computational cost and the potential inconsistency of state variables updated at each loop of data assimilation and its corresponding reservoir simulated result.This problem forbids the reservoir engineers to make the best use of the 4D seismic data,which provides valuable information about the fluid change inside the reservoir.Moreover,only matching the production data in the past is not enough to accurately forecast the future,and the development plan based on the false forecast is very likely to be suboptimal.To solve this problem,we developed a workflow for geophysical and production data history matching by modifying ensemble smoother with multiple data assimilation(ESMDA).In this work,we derived the mathematical expressions of ESMDA and discussed its scope of applications.The geophysical data we used is P-wave impedance,which is typically included in a basic seismic interpretation,and it directly reflects the saturation change in the reservoir.Full resolution of the seismic data is not necessary,we subsampled the P-wave impedance data to further reduce the computational cost.With our case studies on a benchmark synthetic reservoir model,we also showed the supremacy of matching both geophysical and production data,than the traditional reservoir history matching merely on the production data:the overall percentage error of the observed data is halved,and the variances of the updated forecasts are reduced by two orders of the magnitude.

Study of Spectral Response Characteristics of Oilseed Rape (Brassica napus) to Particulate Matters Based on Hyper-Spectral Technique

Haze is mainly caused by the suspended particulate matters in the air,of which the particulate matters pollution harms leaf vegetables.In this paper,oilseed rapes at four different growing periods were investigated in a simulated particulate pollution environment.In combination of hyper-spectral technology and micro examination,the response of hyper-spectral characteristics of the leaf to particulate matters was investigated in-depth.The hyperspectral,chlorophyll content,net photosynthetic rate and stomatal conductance of leaf were obtained.The deposition and adsorption of particulate matters on the leaf were observed by Environmental Scanning Electron Microscope(ESEM).Normalized difference vegetation index(NDVI),modified red edge normalized(mNDVI705)and modified red edge simple ratio index(mSR705)were selected as characteristic parameters and the range of 510 nm~620 nm as the sensitive band.16 methods were used to establish the physiological information inversion model.The main results were as follows:Under the influence of particulate matters,the spectral reflectance decreased as a whole.With the increase of leaf age,the phenomenon of blue shift aggravated.The amplitude of yellow and blue edge decreased with overall decreasing vegetation indices.The furrows and irregular band protrusions in leaves were favorable for keeping particulate matters.With longer affecting time and more deposition of particle matters on the leaf,the stomatal opening became smaller.After comparing,principal component regression(PCR)+multiple scatter correction(MSC)+second derivative(SD)+Savitzky-Golay smooth(SG),and partial least square(PLS)+multiple scatter correction(MSC)+first derivative(FD)+Savitzky-Golay smooth(SG)were determined the best method to establish the inversion model of chlorophyll content and net photosynthetic rate respectively.This study may bring novel ideas for the diagnosis and analysis of the physiological response of leaf vegetables under particulate matters pollution using hyper-spectral technology.