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.

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.

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.

Robust inversion analysis of local gravity anomalies caused by geological dislocation model of faults

Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.

Inversion Methods of Optical Constants of Semitransparent Solid Materials from Transmittance Spectrograms

The direct calculation models of spectral transmittance of single and double slabs consisted of semitransparent solid materials were developed based on ray trace method,and a new inversion method of optical constants ( k is extinction coefficient and n is refractive index) of materials was proposed based on transmittance spectrograms of double slabs. Differences between the new method and two others currently used methods were studied,and application range of methods was also investigated. Optical constants of selenide glass attained in references were selected as true values,and spectral transmittances of glass simulated based on direct calculation model were regarded as experimental values. Optical constants of selenide glass were achieved by inverse models. Influences of measurement error on inverse results were also determined. The results showed that: (1) based on transmittance spectrograms of double slabs in which thickness of single slab is the same,the new proposed method can attain optical constants of materials; (2) the effect of optical constants n and k on three inversion methods are urgent larger,but inversed calculation precision of optical constants are higher in most application ranges; (3) the influence of measurement errors existed in experimental datum on the inverse precision of three methods are urgent

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.

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.

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.

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.

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.

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.

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.

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.

Inverse Calculation of Wave-Absorbing Structure Dimensions Based on Extended ANFIS Model

A new wave energy dissipation structure is proposed, aiming to optimize the dimensions of the structure and make the reflection of the structure maintain a low level within the scope of the known frequency band. An optimal extended ANFIS model combined with the wave reflection coefficient analysis for the estimation of the structure dimensions is established. In the premise of lower wave reflection coefficient, the specific sizes of the structure are obtained inversely, and the contribution of each related parameter on the structural reflection performance is analyzed. The main influencing factors are determined. It is found that the optimal dimensions of the proposed structure exist, which make the wave absorbing performance of the structure reach a perfect status under a wide wave frequency band.

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.

Nonlinear inverse modeling of sensor based on back-propagation fuzzy logical system

Objective To correct the nonlinear error of sensor output,a new approach to sensor inverse modeling based on Back-Propagation Fuzzy Logical System(BP FS) is presented.Methods The BP FS is a computationally efficient nonlinear universal approximator,which is capable of implementing complex nonlinear mapping from its input pattern space to the output with fast convergence speed.Results The neuro-fuzzy hybrid system,i.e.BP FS,is then applied to construct nonlinear inverse model of pressure sensor.The experimental results show that the proposed inverse modeling method automatically compensates the associated nonlinear error in pressure estimation,and thus the performance of pressure sensor is significantly improved.Conclusion The proposed method can be widely used in nonlinearity correction of various kinds of sensors to compensate the effects of nonlinearity and temperature on sensor output.

Current studies on air distributions in commercial airliner cabins

Abstract Air distribution in commercial airliner cabins is very important for the comfort and health of passengers and crew. Experimental measurements, computational fluid dynamics （CFD） simulations, and inverse modeling are state-of-the-art methods available for studying the air distri- bution. This paper gave an overview of the different experimental models, such as scale models, simplified models, full-scale mockups, and actual air cabins. Although experimental measurements were expensive and time consuming, the data were essential for validating CFD simulations. Different modeling strategies for CFD simulations were also discussed in this paper, including large eddy simulations and Reynolds averaged Navier-Stokes equation modeling. CFD simulations were main stream approaches for studying the air distribution but they could not easily lead to optimal design. Inverse modeling of air distribution has recently emerged into a very powerful and attractive tool for designing the air distribution in airliner cabins, although most of the studies were preliminary.