Improving the accuracy of precipitation estimates in a typical inland arid area of China using a dynamic Bayesian model averaging approach

Xinjiang Uygur Autonomous Region is a typical inland arid area in China with a sparse and uneven distribution of meteorological stations,limited access to precipitation data,and significant water scarcity.Evaluating and integrating precipitation datasets from different sources to accurately characterize precipitation patterns has become a challenge to provide more accurate and alternative precipitation information for the region,which can even improve the performance of hydrological modelling.This study evaluated the applicability of widely used five satellite-based precipitation products(Climate Hazards Group InfraRed Precipitation with Station(CHIRPS),China Meteorological Forcing Dataset(CMFD),Climate Prediction Center morphing method(CMORPH),Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record(PERSIANN-CDR),and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis(TMPA))and a reanalysis precipitation dataset(ECMWF Reanalysis v5-Land Dataset(ERA5-Land))in Xinjiang using ground-based observational precipitation data from a limited number of meteorological stations.Based on this assessment,we proposed a framework that integrated different precipitation datasets with varying spatial resolutions using a dynamic Bayesian model averaging(DBMA)approach,the expectation-maximization method,and the ordinary Kriging interpolation method.The daily precipitation data merged using the DBMA approach exhibited distinct spatiotemporal variability,with an outstanding performance,as indicated by low root mean square error(RMSE=1.40 mm/d)and high Person's correlation coefficient(CC=0.67).Compared with the traditional simple model averaging(SMA)and individual product data,although the DBMA-fused precipitation data were slightly lower than the best precipitation product(CMFD),the overall performance of DBMA was more robust.The error analysis between DBMA-fused precipitation dataset and the more advanced Integrated Multi-satellite Retrievals for Global Precipitation Measurement Final(IMERG-F)precipitation product,as well as hydrological simulations in the Ebinur Lake Basin,further demonstrated the superior performance of DBMA-fused precipitation dataset in the entire Xinjiang region.The proposed framework for solving the fusion problem of multi-source precipitation data with different spatial resolutions is feasible for application in inland arid areas,and aids in obtaining more accurate regional hydrological information and improving regional water resources management capabilities and meteorological research in these regions.

Cross Validation Based Model Averaging for Varying-Coefficient Models with Response Missing at Random

In this paper, a model averaging method is proposed for varying-coefficient models with response missing at random by establishing a weight selection criterion based on cross-validation. Under certain regularity conditions, it is proved that the proposed method is asymptotically optimal in the sense of achieving the minimum squared error.

Improving microwave brightness temperature predictions based on Bayesian model averaging ensemble approach

The choices of the parameterizations for each component in a microwave emission model have significant effects on the quality of brightness temperature （Tb） sim- ulation. How to reduce the uncertainty in the Tb simulation is investigated by adopting a statistical post-processing procedure with the Bayesian model averaging （BMA） ensemble approach. The simulations by the community microwave emission model （CMEM） cou- pled with the community land model version 4.5 （CLM4.5） over China's Mainland are con- ducted by the 24 configurations from four vegetation opacity parameterizations （VOPs）, three soil dielectric constant parameterizations （SDCPs）, and two soil roughness param- eterizations （SRPs）. Compared with the simple arithmetical averaging （SAA） method, the BMA reconstructions have a higher spatial correlation coefficient （larger than 0.99） than the C-band satellite observations of the advanced microwave scanning radiometer on the Earth observing system （AMSR-E） at the vertical polarization. Moreover, the BMA product performs the best among the ensemble members for all vegetation classes, with a mean root-mean-square difference （RMSD） of 4 K and a temporal correlation coefficient of 0.64.

Response of Growing Season Gross Primary Production to El Nino in Different Phases of the Pacific Decadal Oscillation over Eastern China Based on Bayesian Model Averaging

Gross primary production(GPP) plays a crucial part in the carbon cycle of terrestrial ecosystems.A set of validated monthly GPP data from 1957 to 2010 in 0.5°× 0.5° grids of China was weighted from the Multi-scale Terrestrial Model Intercomparison Project using Bayesian model averaging(BMA).The spatial anomalies of detrended BMA GPP during the growing seasons of typical El Nino years indicated that GPP response to El Nino varies with Pacific Decadal Oscillation(PDO) phases: when the PDO was in the cool phase,it was likely that GPP was greater in northern China(32°–38°N,111°–122°E) and less in the Yangtze River valley(28°–32°N,111°–122°E);in contrast,when PDO was in the warm phase,the GPP anomalies were usually reversed in these two regions.The consistent spatiotemporal pattern and high partial correlation revealed that rainfall dominated this phenomenon.The previously published findings on how El Nino during different phases of PDO affecting rainfall in eastern China make the statistical relationship between GPP and El Nino in this study theoretically credible.This paper not only introduces an effective way to use BMA in grids that have mixed plant function types,but also makes it possible to evaluate the carbon cycle in eastern China based on the prediction of El Nino and PDO.

Optimal zero-crossing group selection method of the absolute gravimeter based on improved auto-regressive moving average model

An absolute gravimeter is a precision instrument for measuring gravitational acceleration, which plays an important role in earthquake monitoring, crustal deformation, national defense construction, etc. The frequency of laser interference fringes of an absolute gravimeter gradually increases with the fall time. Data are sparse in the early stage and dense in the late stage. The fitting accuracy of gravitational acceleration will be affected by least-squares fitting according to the fixed number of zero-crossing groups. In response to this problem, a method based on Fourier series fitting is proposed in this paper to calculate the zero-crossing point. The whole falling process is divided into five frequency bands using the Hilbert transformation. The multiplicative auto-regressive moving average model is then trained according to the number of optimal zero-crossing groups obtained by the honey badger algorithm. Through this model, the number of optimal zero-crossing groups determined in each segment is predicted by the least-squares fitting. The mean value of gravitational acceleration in each segment is then obtained. The method can improve the accuracy of gravitational measurement by more than 25% compared to the fixed zero-crossing groups method. It provides a new way to improve the measuring accuracy of an absolute gravimeter.

Climate change in the Tianshan and northern Kunlun Mountains based on GCM simulation ensemble with Bayesian model averaging

Climate change in mountainous regions has significant impacts on hydrological and ecological systems. This research studied the future temperature, precipitation and snowfall in the 21^（st） century for the Tianshan and northern Kunlun Mountains（TKM） based on the general circulation model（GCM） simulation ensemble from the coupled model intercomparison project phase 5（CMIP5） under the representative concentration pathway（RCP） lower emission scenario RCP4.5 and higher emission scenario RCP8.5 using the Bayesian model averaging（BMA） technique. Results show that（1） BMA significantly outperformed the simple ensemble analysis and BMA mean matches all the three observed climate variables;（2） at the end of the 21^（st） century（2070–2099） under RCP8.5, compared to the control period（1976–2005）, annual mean temperature and mean annual precipitation will rise considerably by 4.8°C and 5.2%, respectively, while mean annual snowfall will dramatically decrease by 26.5%;（3） precipitation will increase in the northern Tianshan region while decrease in the Amu Darya Basin. Snowfall will significantly decrease in the western TKM. Mean annual snowfall fraction will also decrease from 0.56 of 1976–2005 to 0.42 of 2070–2099 under RCP8.5; and（4） snowfall shows a high sensitivity to temperature in autumn and spring while a low sensitivity in winter, with the highest sensitivity values occurring at the edge areas of TKM. The projections mean that flood risk will increase and solid water storage will decrease.

Modeling continuous traffic flow with the average velocity effect of multiple vehicles ahead on gyroidal roads

In the future connected vehicle environment,the information of multiple vehicles ahead can be readily collected in real-time,such as the velocity or headway,which provides more opportunities for information exchange and cooperative control.Meanwhile,gyroidal roads are one of the fundamental road patterns prevalent in mountainous areas.To effectively control the system,it is therefore significant to explore the evolution mechanism of traffic flow on gyroidal roads under a connected vehicle environment.In this paper,we present a new continuum model with the average velocity of multiple vehicles ahead on gyroidal roads.The stability criterion and KdV-Burger equation are deduced via linear and nonlinear stability analysis,respectively.Solving the above KdV-Burger equation yields the density wave solution,which explores the formation and propagation property of traffic jams near the neutral stability curve.Simulation examples verify that the model can reproduce complex phenomena,such as shock waves and rarefaction waves.The analysis of the local cluster effect shows that the number of vehicles ahead and the radius information,and the slope information of gyroidal roads can exert a great influence on traffic jams.The effect of the first and second terms are positive,while the last term is negative.

An Improved Time Domain Approach for Analysis of Floating Bridges Based on Dynamic Finite Element Method and State-Space Model

The floating bridge bears the dead weight and live load with buoyancy,and has wide application prospect in deep-water transportation infrastructure.The structural analysis of floating bridge is challenging due to the complicated fluid-solid coupling effects of wind and wave.In this research,a novel time domain approach combining dynamic finite element method and state-space model(SSM)is established for the refined analysis of floating bridges.The dynamic coupled effects induced by wave excitation load,radiation load and buffeting load are carefully simulated.High-precision fitted SSMs for pontoons are established to enhance the calculation efficiency of hydrodynamic radiation forces in time domain.The dispersion relation is also introduced in the analysis model to appropriately consider the phase differences of wave loads on pontoons.The proposed approach is then employed to simulate the dynamic responses of a scaled floating bridge model which has been tested under real wind and wave loads in laboratory.The numerical results are found to agree well with the test data regarding the structural responses of floating bridge under the considered environmental conditions.The proposed time domain approach is considered to be accurate and effective in simulating the structural behaviors of floating bridge under typical environmental conditions.

A Review on High-Dimensional Frequentist Model Averaging

Model averaging has attracted increasing attention in recent years for the analysis of high-dimensional data. By weighting several competing statistical models suitably, model averaging attempts to achieve stable and improved prediction. To obtain a better understanding of the available model averaging methods, their properties and the relationships between them, this paper is devoted to make a review on some recent progresses in high-dimensional model averaging from the frequentist perspective. Some future research topics are also discussed.

Comparison between Different ESI Methods on Refractory Epilepsy Patients Shows a High Sensitivity for Bayesian Model Averaging

Electrical Source Imaging (ESI) is a non-invasive technique of reconstructing brain activities using EEG data. This technique has been applied to evaluate epilepsy patients being evaluated for epilepsy surgery, showing encouraging results for mapping interictal epileptiform discharges (IED). However, ESI is underused in planning epilepsy surgery. This is basically due to the wide availability of methods for solving the electromagnetism inverse problem (e-IP) associated to few studies using EEG setups similar to those most commonly used in clinical setting. In this study, we applied six different methods of solving the e-IP based on IEDs of 20 focal epilepsy patients that presented abnormalities in their MRI. We compared the ESI maps obtained by each method with the location of the abnormality, calculating the Euclidian distances from the center of the lesion to the closest border of the method solution (CL-BM) and also to the solution’s maxima (CL-MM). We also applied a score system in order to allow us to evaluate the sensitivity of each method for temporal and extra temporal patients. In our patients, the Bayesian Model Averaging method had a sensitivity of 86% and the shortest CL-MM. This method also had more restricted solutions that were more representative of epileptogenic activities than those obtained by the other methods.

Model Averaging by Stacking

The paper introduces a new Frequentist model averaging estimation procedure, based on a stacked OLS estimator across models, implementable on cross-sectional, panel, as well as time series data. The proposed estimator shows the same optimal properties of the OLS estimator under the usual set of assumptions concerning the population regression model. Relatively to available alternative approaches, it has the advantage of performing model averaging exante in a single step, optimally selecting models’ weight according to the MSE metric, i.e. by minimizing the squared Euclidean distance between actual and predicted value vectors. Moreover, it is straightforward to implement, only requiring the estimation of a single OLS augmented regression. By exploiting exante a broader information set and benefiting of more degrees of freedom, the proposed approach yields more accurate and (relatively) more efficient estimation than available expost methods.

Improving the simulation of terrestrial water storage anomalies over China using a Bayesian model averaging ensemble approach

作为全球能量和水分循环的关键参量,陆地水储量包括土壤水、地表水、地下水、积雪和生物体水等,在水文、气候、农业、生态等众多领域起重要影响。与地面观测和遥感反演相比,陆面模式在刻画陆地水储量的时空变率等方面具有明显优势。然而不同模式参数化方案以及大气强迫驱动导致陆地水储量模拟存在不确定。为了减少陆地水储量模拟不确定性,本研究建立了基于贝叶斯模型平均(BMA)和多强迫多模式集合的陆地水储量模拟系统,获得了中国区域1979–2008年陆地水储量数据集。选取2004–08年的数据与GRACE重力卫星数据比较分析,结果显示BMA集合模拟的陆地水储量异常(Terrestrial water storage anomalies,TWSA)优于所有单个模拟结果,与GRACE观测的TWSA有更高的相关系数和更小的误差。

A Mixture-Based Bayesian Model Averaging Method

Bayesian model averaging (BMA) is a popular and powerful statistical method of taking account of uncertainty about model form or assumption. Usually the long run (frequentist) performances of the resulted estimator are hard to derive. This paper proposes a mixture of priors and sampling distributions as a basic of a Bayes estimator. The frequentist properties of the new Bayes estimator are automatically derived from Bayesian decision theory. It is shown that if all competing models have the same parametric form, the new Bayes estimator reduces to BMA estimator. The method is applied to the daily exchange rate Euro to US Dollar.

Frequentist Model Averaging and Applications to Bernoulli Trials

In several instances of statistical practice, it is not uncommon to use the same data for both model selection and inference, without taking account of the variability induced by model selection step. This is usually referred to as post-model selection inference. The shortcomings of such practice are widely recognized, finding a general solution is extremely challenging. We propose a model averaging alternative consisting on taking into account model selection probability and the like-lihood in assigning the weights. The approach is applied to Bernoulli trials and outperforms Akaike weights model averaging and post-model selection estimators.

Least Squares Matrix Algorithm for State-Space Modelling of Dynamic Systems

This work presents a novel least squares matrix algorithm (LSM) for the analysis of rapidly changing systems using state-space modelling. The LSM algorithm is based on the Hankel structured data matrix representation. The state transition matrix is updated without the use of any forgetting function. This yields a robust estimation of model parameters in the presence of noise. The computational complexity of the LSM algorithm is comparable to the speed of the conventional recursive least squares (RLS) algorithm. The knowledge of the state transition matrix enables feasible numerical operators such as interpolation, fractional differentiation and integration. The usefulness of the LSM algorithm was proved in the analysis of the neuroelectric signal waveforms.

Application of the Moving Averaging Technique in Surplus Production Models

Surplus production models are the simplest analytical methods effective for fish stock assessment and fisheries management. In this paper, eight surplus production estimators(three estimation procedures) were tested on Schaefer and Fox type simulated data in three simulated fisheries(declining, well-managed, and restoring fisheries) at two white noise levels. Monte Carlo simulation was conducted to verify the utility of moving averaging(MA), which was an important technique for reducing the effect of noise in data in these models. The relative estimation error(REE) of maximum sustainable yield(MSY) was used as an indicator for the analysis, and one-way ANOVA was applied to test the significance of the REE calculated at four levels of MA. Simulation results suggested that increasing the value of MA could significantly improve the performance of the surplus production model(low REE) in all cases when the white noise level was low(coefficient of variation(CV) = 0.02). However, when the white noise level increased(CV= 0.25), adding the value of MA could still significantly enhance the performance of most models. Our results indicated that the best model performance occurred frequently when MA was equal to 3; however, some exceptions were observed when MA was higher.

Partial Linear Model Averaging Prediction for Longitudinal Data

Prediction plays an important role in data analysis.Model averaging method generally provides better prediction than using any of its components.Even though model averaging has been extensively investigated under independent errors,few authors have considered model averaging for semiparametric models with correlated errors.In this paper,the authors offer an optimal model averaging method to improve the prediction in partially linear model for longitudinal data.The model averaging weights are obtained by minimizing criterion,which is an unbiased estimator of the expected in-sample squared error loss plus a constant.Asymptotic properties,including asymptotic optimality and consistency of averaging weights,are established under two scenarios:(i)All candidate models are misspecified;(ii)Correct models are available in the candidate set.Simulation studies and an empirical example show that the promise of the proposed procedure over other competitive methods.

Improved State-space Modelling for Microgrids Without Virtual Resistances

Power converters and their interfacing networks are often treated as modular state-space blocks for small-signal stability studies in microgrids;they are interconnected by matching the input and output states of the network and converters.Virtual resistors have been widely used in existing models to generate a voltage for state-space models of the network that require voltage inputs.This paper accurately quantifies the adverse impacts of adding the virtual resistance and proposes an alternative method for network modelling that eliminates the requirement of the virtual resistor when interfacing converters with microgrids.The proposed nonlinear method allows initialization,time-domain simulations of the nonlinear model,and linearization and eigenvalue generation.A numerically linearized small-signal model is used to generate eigenvalues and is compared with the eigenvalues generated using the existing modelling method with virtual resistances.Deficiencies of the existing method and improvements offered by the proposed modelling method are clearly quantified.Electromagnetic transient(EMT)simulations using detailed switching models are used for validation of the proposed modelling method.

Application of Seasonal Auto-regressive Integrated Moving Average Model in Forecasting the Incidence of Hand-foot-mouth Disease in Wuhan,China

Outbreaks of hand-foot-mouth disease（HFMD） have occurred many times and caused serious health burden in China since 2008. Application of modern information technology to prediction and early response can be helpful for efficient HFMD prevention and control. A seasonal auto-regressive integrated moving average（ARIMA） model for time series analysis was designed in this study. Eighty-four-month（from January 2009 to December 2015） retrospective data obtained from the Chinese Information System for Disease Prevention and Control were subjected to ARIMA modeling. The coefficient of determination（R^2）, normalized Bayesian Information Criterion（BIC） and Q-test P value were used to evaluate the goodness-of-fit of constructed models. Subsequently, the best-fitted ARIMA model was applied to predict the expected incidence of HFMD from January 2016 to December 2016. The best-fitted seasonal ARIMA model was identified as（1,0,1）（0,1,1）12, with the largest coefficient of determination（R^2=0.743） and lowest normalized BIC（BIC=3.645） value. The residuals of the model also showed non-significant autocorrelations（P_（Box-Ljung（Q））=0.299）. The predictions by the optimum ARIMA model adequately captured the pattern in the data and exhibited two peaks of activity over the forecast interval, including a major peak during April to June, and again a light peak for September to November. The ARIMA model proposed in this study can forecast HFMD incidence trend effectively, which could provide useful support for future HFMD prevention and control in the study area. Besides, further observations should be added continually into the modeling data set, and parameters of the models should be adjusted accordingly.

Average modeling of Single Stage Flyback PFC ＋ Flyback DC／DC converter

With the use of this novel average model for Single Stage Flyback PFC + Flyback DC/DC converter, voltage control mode, peak current control mode and average current control mode can be simulated easily by changing the model's parmneters. It can be used to do various analysis not only for small signal and static behavior but also for large signal and dynamic behavior of the converter. By using this average model the simulation speed can be improved by 2 orders of rrmgnitude above that obtained by using the conventional switched model. It can be applied to optimize the trade-off between high power factor, voltage stress, current stress and good output performance while designing this kind of single stage PFC converter. A 60W single stage power factor corrector was built to verify the proposed model. The modeling principle can be applied to other Single Stage PFC topologies.