Runout prediction of potential landslides based on the multi-source data collaboration analysis on historical cases

Long runout landslides involve a massive amount of energy and can be extremely hazardous owing to their long movement distance,high mobility and strong destructive power.Numerical methods have been widely used to predict the landslide runout but a fundamental problem remained is how to determine the reliable numerical parameters.This study proposes a framework to predict the runout of potential landslides through multi-source data collaboration and numerical analysis of historical landslide events.Specifically,for the historical landslide cases,the landslide-induced seismic signal,geophysical surveys,and possible in-situ drone/phone videos(multi-source data collaboration)can validate the numerical results in terms of landslide dynamics and deposit features and help calibrate the numerical(rheological)parameters.Subsequently,the calibrated numerical parameters can be used to numerically predict the runout of potential landslides in the region with a similar geological setting to the recorded events.Application of the runout prediction approach to the 2020 Jiashanying landslide in Guizhou,China gives reasonable results in comparison to the field observations.The numerical parameters are determined from the multi-source data collaboration analysis of a historical case in the region(2019 Shuicheng landslide).The proposed framework for landslide runout prediction can be of great utility for landslide risk assessment and disaster reduction in mountainous regions worldwide.

Full-Scale Numerical Simulation of the Local Scour Under Combined Current and Wave Conditions Based on Field Data

The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of wind turbine foundations.In this paper,a full-scale numerical model was developed and validated based on field data from Rudong,China.The scour development around monopiles was investigated,and the effects of waves and the Reynolds number Re were analyzed.Several formulas for predicting the scour depth in the literature have been evaluated.It is found that waves can accelerate scour development even if the KC number is small(0.78<KC<1.57).The formula obtained from small-scale model tests may be unsafe or wasteful when it is applied in practical design due to the scale effect.A new equation for predicting the scour depth based on the average pile Reynolds number(Rea)is proposed and validated with field data.The equilibrium scour depth predicted using the proposed equation is evaluated and compared with those from nine equations in the literature.It is demonstrated that the values predicted from the proposed equation and from the S/M(Sheppard/Melville)equation are closer to the field data.

An Analysis of Observed Ocean Currents Data and Numerical Simulation in the Offshore Area of Jiaonan

Based on the diurnal consecutively observed data in the offshore area of Jiaonan in 2005, the paper tries to make a preliminary analysis of the specificity of ocean currents, tidal current property and residual current property in the area in observing dates. Then on the basis of observed data analysis and by employing the split-step method, the paper conducts a numerical simulation of the tidal current field, which can show the M2 tidal constituent tidal wave system, current ellipse distribution, maximum current velocity distribution and time-dependent current field. The calculated results agree well with the observed data, which can on the one hand reflect the basic specificities of temporal and spatial distribution of the M2 tidal constituent current field to some extent, and, on the other hand, offer more information about the hydrodynamic condition. So the paper would provide a scientific basis for the making of sea environment protection plans in the offshore area of Jiaonan under certain conditions.

Workflow-Based Platform Design and Implementation for Numerical Weather Prediction Models and Meteorological Data Service

In this paper, we present a set of best practices for workflow design and implementation for numerical weather prediction models and meteorological data service, which have been in operation in China Meteorological Administration (CMA) for years and have been proven effective in reliably managing the complexities of large-scale meteorological related workflows. Based on the previous work on the platforms, we argue that a minimum set of guidelines including workflow scheme, module design, implementation standards and maintenance consideration during the whole establishment of the platform are highly recommended, serving to reduce the need for future maintenance and adjustment. A significant gain in performance can be achieved through the workflow-based projects. We believe that a good workflow system plays an important role in the weather forecast service, providing a useful tool for monitoring the whole process, fixing the errors, repairing a workflow, or redesigning an equivalent workflow pattern with new components.

System of Multigrid Nonlinear Least-squares Four-dimensional Variational Data Assimilation for Numerical Weather Prediction(SNAP):System Formulation and Preliminary Evaluation

A new forecasting system-the System of Multigrid Nonlinear Least-squares Four-dimensional Variational(NLS-4DVar)Data Assimilation for Numerical Weather Prediction(SNAP)-was established by building upon the multigrid NLS-4DVar data assimilation scheme,the operational Gridpoint Statistical Interpolation(GSI)−based data-processing and observation operators,and the widely used Weather Research and Forecasting numerical model.Drawing upon lessons learned from the superiority of the operational GSI analysis system,for its various observation operators and the ability to assimilate multiple-source observations,SNAP adopts GSI-based data-processing and observation operator modules to compute the observation innovations.The multigrid NLS-4DVar assimilation framework is used for the analysis,which can adequately correct errors from large to small scales and accelerate iteration solutions.The analysis variables are model state variables,rather than the control variables adopted in the conventional 4DVar system.Currently,we have achieved the assimilation of conventional observations,and we will continue to improve the assimilation of radar and satellite observations in the future.SNAP was evaluated by case evaluation experiments and one-week cycling assimilation experiments.In the case evaluation experiments,two six-hour time windows were established for assimilation experiments and precipitation forecasts were verified against hourly precipitation observations from more than 2400 national observation sites.This showed that SNAP can absorb observations and improve the initial field,thereby improving the precipitation forecast.In the one-week cycling assimilation experiments,six-hourly assimilation cycles were run in one week.SNAP produced slightly lower forecast RMSEs than the GSI 4DEnVar(Four-dimensional Ensemble Variational)as a whole and the threat scores of precipitation forecasts initialized from the analysis of SNAP were higher than those obtained from the analysis of GSI 4DEnVar.

Applying Hybrid Clustering in Pulsar Candidate Sifting with Multi-modality for FAST Survey

Pulsar search is always the basis of pulsar navigation,gravitational wave detection and other research topics.Currently,the volume of pulsar candidates collected by the Five-hundred-meter Aperture Spherical radio Telescope(FAST)shows an explosive growth rate that has brought challenges for its pulsar candidate filtering system.Particularly,the multi-view heterogeneous data and class imbalance between true pulsars and non-pulsar candidates have negative effects on traditional single-modal supervised classification methods.In this study,a multi-modal and semi-supervised learning based on a pulsar candidate sifting algorithm is presented,which adopts a hybrid ensemble clustering scheme of density-based and partition-based methods combined with a feature-level fusion strategy for input data and a data partition strategy for parallelization.Experiments on both High Time Resolution Universe SurveyⅡ(HTRU2)and actual FAST observation data demonstrate that the proposed algorithm could excellently identify pulsars:On HTRU2,the precision and recall rates of its parallel mode reach0.981 and 0.988 respectively.On FAST data,those of its parallel mode reach 0.891 and 0.961,meanwhile,the running time also significantly decreases with the increment of parallel nodes within limits.Thus,we can conclude that our algorithm could be a feasible idea for large scale pulsar candidate sifting for FAST drift scan observation.

Assimilating FY-4A AGRI Radiances with a Channel-Sensitive Cloud Detection Scheme for the Analysis and Forecasting of Multiple Typhoons

This paper presents an attempt at assimilating clear-sky FY-4A Advanced Geosynchronous Radiation Imager(AGRI)radiances from two water vapor channels for the prediction of three landfalling typhoon events over the West Pacific Ocean using the 3DVar data assimilation(DA)method along with the WRF model.A channel-sensitive cloud detection scheme based on the particle filter(PF)algorithm is developed and examined against a cloud detection scheme using the multivariate and minimum residual(MMR)algorithm and another traditional cloud mask–dependent cloud detection scheme.Results show that both channel-sensitive cloud detection schemes are effective,while the PF scheme is able to reserve more pixels than the MMR scheme for the same channel.In general,the added value of AGRI radiances is confirmed when comparing with the control experiment without AGRI radiances.Moreover,it is found that the analysis fields of the PF experiment are mostly improved in terms of better depicting the typhoon,including the temperature,moisture,and dynamical conditions.The typhoon track forecast skill is improved with AGRI radiance DA,which could be explained by better simulating the upper trough.The impact of assimilating AGRI radiances on typhoon intensity forecasts is small.On the other hand,improved rainfall forecasts from AGRI DA experiments are found along with reduced errors for both the thermodynamic and moisture fields,albeit the improvements are limited.

Parameterized Forward Operators for Simulation and Assimilation of Polarimetric Radar Data with Numerical Weather Predictions

Many weather radar networks in the world have now provided polarimetric radar data(PRD)that have the potential to improve our understanding of cloud and precipitation microphysics,and numerical weather prediction(NWP).To realize this potential,an accurate and efficient set of polarimetric observation operators are needed to simulate and assimilate the PRD with an NWP model for an accurate analysis of the model state variables.For this purpose,a set of parameterized observation operators are developed to simulate and assimilate polarimetric radar data from NWP model-predicted hydrometeor mixing ratios and number concentrations of rain,snow,hail,and graupel.The polarimetric radar variables are calculated based on the T-matrix calculation of wave scattering and integrations of the scattering weighted by the particle size distribution.The calculated polarimetric variables are then fitted to simple functions of water content and volumeweighted mean diameter of the hydrometeor particle size distribution.The parameterized PRD operators are applied to an ideal case and a real case predicted by the Weather Research and Forecasting(WRF)model to have simulated PRD,which are compared with existing operators and real observations to show their validity and applicability.The new PRD operators use less than one percent of the computing time of the old operators to complete the same simulations,making it efficient in PRD simulation and assimilation usage.

A statistical analysis of geomechanical data and its effect on rock mass numerical modeling:a case study

Geomechanical data are never sufficient in quantity or adequately precise and accurate for design purposes in mining and civil engineering.The objective of this paper is to show the variability of rock properties at the sampled point in the roadway's roof,and then,how the statistical processing of the available geomechanical data can affect the results of numerical modelling of the roadway's stability.Four cases were applied in the numerical analysis,using average values(the most common in geomechanical data analysis),average minus standard deviation,median,and average value minus statistical error.The study show that different approach to the same geomechanical data set can change the modelling results considerably.The case shows that average minus standard deviation is the most conservative and least risky.It gives the displacements and yielded elements zone in four times broader range comparing to the average values scenario,which is the least conservative option.The two other cases need to be studied further.The results obtained from them are placed between most favorable and most adverse values.Taking the average values corrected by statistical error for the numerical analysis seems to be the best solution.Moreover,the confidence level can be adjusted depending on the object importance and the assumed risk level.

A MATLAB-Based Numerical and GUI Implementation of Cross-Gradients Joint Inversion of Gravity and Magnetic Data

The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in the literature. We present a MATLAB-based three-dimensional cross-gradients joint inversion program with application to gravity and magnetic data. The input and output information was examined with care to create a rational, independent design of a graphical user interface (GUI) and computing kernel. For 3D visualization and data file operations, UBC-GIF tools are invoked using a series of I/O functions. Some key issues regarding the iterative joint inversion algorithm are also discussed: for instance, the forward difference of cross gradients, and matrix pseudo inverse computation. A synthetic example is employed to illustrate the whole process. Joint and separate inversions can be performed flexibly by switching the inversion mode. The resulting density model and susceptibility model demonstrate the correctness of the proposed program.

Current Status and Future Challenges of Weather Radar Polarimetry: Bridging the Gap between Radar Meteorology/Hydrology/Engineering and Numerical Weather Prediction

After decades of research and development, the WSR-88 D(NEXRAD) network in the United States was upgraded with dual-polarization capability, providing polarimetric radar data(PRD) that have the potential to improve weather observations,quantification, forecasting, and warnings. The weather radar networks in China and other countries are also being upgraded with dual-polarization capability. Now, with radar polarimetry technology having matured, and PRD available both nationally and globally, it is important to understand the current status and future challenges and opportunities. The potential impact of PRD has been limited by their oftentimes subjective and empirical use. More importantly, the community has not begun to regularly derive from PRD the state parameters, such as water mixing ratios and number concentrations, used in numerical weather prediction(NWP) models.In this review, we summarize the current status of weather radar polarimetry, discuss the issues and limitations of PRD usage, and explore potential approaches to more efficiently use PRD for quantitative precipitation estimation and forecasting based on statistical retrieval with physical constraints where prior information is used and observation error is included. This approach aligns the observation-based retrievals favored by the radar meteorology community with the model-based analysis of the NWP community. We also examine the challenges and opportunities of polarimetric phased array radar research and development for future weather observation.

Initialization and Simulation of a Typhoon Using 4-Dimensional Variational Data Assimilation——Research on Typhoon Herb (1996)

The MM5 and its four dimensional variational data assimilation (4D-Var) system are used in this paper. Based on the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis data, the authors generate an optimal initial condition for a typhoon by using the bogus data assimilation (BDA) scheme. BDA is able to recover many of the structural features of typhoons including a warm-core vertex, the correct center position, and the strong circulation. As a result of BDA using a bogus surface low, dramatic improvement is achieved in the 72 h prediction of typhoon Herb. Through several cases, the initialization by BDA effectively generates the harmonious inner structure of the typhoon, but which is lacking in the original analysis field. Therefore the intensity forecast is improved greatly. Some improvements are made in the track forecast, but more work still needs to be done.

An approach to estimating and extrapolating model error based on inverse problem methods:towards accurate numerical weather prediction

Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.

The Features of Ozone Quasi-Biennial Oscillation in Tropical Stratosphere and Its Numerical Simulation

The inlcrunnual variation of the vertical distribution of ozone in the tropical stratosphere and its quasi—biennial oscillation (QBO) is analyzed using HALOE data. The results are compared with the wind QBO. A numerical experiment is carried out to study the effects of wind QBO on the distribution, and variation of ozone in the stratosphere by using (he NCAR interactive chemical, dynamical, and radiative two—dimensional model (SOCRATES). Data analysis shows that the location of the maximum ozone mixing ratio in the stratosphere changes in the meridional and vertical directions, and assumes a quasi—biennial period. The meridional and vertical motion of the maximum mixing ratio leads to a QBO of column ozone and its hemispheric asymmetry. The QBO of the location of the maximum is closely connected with the zonal wind QBO. The data analysis also shows that in the tropical region, the phase of the QBO for ozone density changes many times with height. Numerical simulation shows that the meridional circulation induced by the wind QBO includes three pairs of cells in the stratosphere, which have hemispheric symmetry. The transport of ozone by the induced meridional circulation in various latitudes and heights is the main dynamic cause for the ozone QBO. Cells of the induced circulation in the middle stratosphere (25-35 km) play an important role in producing the ozone QBO.

A Forecast Error Correction Method in Numerical Weather Prediction by Using Recent Multiple-time Evolution Data

The initial value error and the imperfect numerical model are usually considered as error sources of numerical weather prediction （NWP）. By using past multi-time observations and model output, this study proposes a method to estimate imperfect numerical model error. This method can be inversely estimated through expressing the model error as a Lagrange interpolation polynomial, while the coefficients of polyno- mial are determined by past model performance. However, for practical application in the full NWP model, it is necessary to determine the following criteria： （1） the length of past data sufficient for estimation of the model errors, （2） a proper method of estimating the term ＂model integration with the exact solution＂ when solving the inverse problem, and （3） the extent to which this scheme is sensitive to the observational errors. In this study, such issues are resolved using a simple linear model, and an advection diffusion model is applied to discuss the sensitivity of the method to an artificial error source. The results indicate that the forecast errors can be largely reduced using the proposed method if the proper length of past data is chosen. To address the three problems, it is determined that （1） a few data limited by the order of the corrector can be used, （2） trapezoidal approximation can be employed to estimate the ＂term＂ in this study; however, a more accurate method should be explored for an operational NWP model, and （3） the correction is sensitive to observational error.

Comparison of satellite-estimated and model-forecasted rainfall data during a deadly debris-flow event in Zhouqu, Northwest China

The data of several rainfall products, including those estimated from satellite measurements and those forecasted via numerical weather modeling, for a severe debris-flow event in Zhouqu, Northwest China, are compared and analyzed in this paper. The satellite products, including CPC MORPHing technique（CMORPH）, TMPA-RT, and PERSIANN are all near-real-time retrieved with high temporal and spatial resolutions. The numerical weather model used in this paper for precipitation forecasting is WRF. The results show that all three satellite products can basically reproduce the rainfall pattern, distribution, timing, scale, and extreme values of the event, compared with gauge data. Their temporal and spatial correlation coefficients with gauge data are as high as about 0.6, which is statistically significant at 0.01 level. The performance of the forecasted results modeled with different spatial resolutions are not as good as the satellite-estimated results, although their correlation coefficients are still statistically significant at 0.05 level. From the total rainfall and extreme value time series for the domain, it is clear that, from the grid-to-grid perspective, the passive microwave-based CMORPH and TRMM products are more accurate than the infrared-based PERSIANN, while PERSIANN performs very well from the general point of view, especially when considering the whole domain or the whole convective precipitation system. The forecasted data — especially the highest resolution model domain data — are able to represent the total or mean precipitation very well in the research domain, while for extreme values the errors are large. This study suggests that satellite-retrieved and model-forecasted rainfall data are a useful complement to gauge data, especially for areas without gauge stations and areas not covered by weather radars.

ANALYSIS OF THE EFFECT OF 3DVAR AND ENSRF DIRECT ASSIMILATION OF RADAR DATA ON THE FORECAST OF A HEAVY RAINFALL EVENT

The present study designs experiments on the direct assimilation of radial velocity and reflectivity data collected by an S-band Doppler weather radar(CINRAD WSR-98D) at the Hefei Station and the reanalysis data produced by the United States National Centers for Environmental Prediction using the Weather Research and Forecasting(WRF) model,the WRF model with a three-dimensional variational(3DVAR) data assimilation system and the WRF model with an ensemble square root filter(EnSRF) data assimilation system.In addition,the present study analyzes a Meiyu front heavy rainfall process that occurred in the Yangtze-Huaihe River Basin from July 4 to July 5,2003,through numerical simulation.The results show the following.(1) The assimilation of the radar radial velocity data can increase the perturbations in the low-altitude atmosphere over the heavy rainfall region,enhance the convective activities and reduce excessive simulated precipitation.(2) The 3DVAR assimilation method significantly adjusts the horizontal wind field.The assimilation of the reflectivity data improves the microphysical quantities and dynamic fields in the model.In addition,the assimilation of the radial velocity and reflectivity data can better adjust the wind fields and improve the intensity and location of the simulated radar echo bands.(3) The EnSRF assimilation method can assimilate more small-scale wind field information into the model.The assimilation of the reflectivity data alone can relatively accurately forecast the rainfall centers.In addition,the assimilation of the radial velocity and reflectivity data can improve the location of the simulated radar echo bands.(4) The use of the 3DVAR and EnSRF assimilation methods to assimilate the radar radial velocity and reflectivity data can improve the forecast of precipitation,rain-band areal coverage and the center location and intensity of precipitation.

Insights into the Microwave Instruments Onboard the Fengyun 3D Satellite:Data Quality and Assimilation in the Met Office NWP System

This paper evaluates the microwave instruments onboard the latest Chinese polar-orbiting satellite, Fengyun 3D (FY- 3D). Comparing three months of observations from the Microwave Temperature Sounder 2 (MWTS-2), the Microwave Humidity Sounder 2 (MWHS-2), and the Microwave Radiation Imager (MWRI) to Met Office short-range forecasts, we characterize the instrumental biases, show how those biases have changed with respect to their predecessors onboard FY- 3C, and how they compare to the Advanced Technology Microwave Sounder (ATMS) onboard NOAA-20 and the Global Precipitation Measurement Microwave Imager (GMI). The MWTS-2 global bias is much reduced with respect to its predecessor and compares well to ATMS at equivalent channel frequencies, differing only by 0.36 ± 0.28 K (1σ) on average. A suboptimal averaging of raw digital counts is found to cause an increase in striping noise and an ascending- descending bias. MWHS-2 benefits from a new calibration method improving the 183-GHz humidity channels with respect to its predecessor and biases for these channels are within ± 1.9 K to ATMS. MWRI presents the largest improvements, with reduced global bias and standard deviation with respect to FY-3C;although, spurious, seemingly transient, brightness temperatures have been detected in the observations at 36.5 GHz (vertical polarization). The strong solar-dependent bias that affects the instrument on FY-3C has been reduced to less than 0.2 K on average for FY-3D MWRI. Experiments where radiances from these instruments were assimilated on top of a full global system demonstrated a neutral to positive impact on the forecasts, as well as on the fit to the background of independent instruments.

All-sky Data Assimilation of MWTS-2 and MWHS-2 in the Met Office Global NWP System.

Microwave radiances from passive polar-orbiting radiometers have been,until recently,assimilated in the Met Office global numerical weather prediction system after the scenes significantly affected by atmospheric scattering are discarded.Recent system upgrades have seen the introduction of a scattering-permitting observation operator and the development of a variable observation error using both liquid and ice water paths as proxies of scattering-induced bias.Applied to the Fengyun 3 Microwave Temperature Sounder 2(MWTS-2)and the Microwave Humidity Sounder 2(MWHS-2),this methodology increases the data usage by up to 8%at 183 GHz.It also allows for the investigation into the assimilation of MWHS-2118 GHz channels,sensitive to temperature and lower tropospheric humidity,but whose large sensitivity to ice cloud have prevented their use thus far.While the impact on the forecast is mostly neutral with small but significant short-range improvements,0.3%in terms of root mean square error,for southern winds and low-level temperature,balanced by 0.2%degradations of short-range northern and tropical low-level temperature,benefits are observed in the background fit of independent instruments used in the system.The lower tropospheric temperature sounding Infrared Atmospheric Sounding Interferometer(IASI)channels see a reduction of the standard deviation in the background departure of up to 1.2%.The Advanced Microwave Sounding Unit A(AMSU-A)stratospheric sounding channels improve by up to 0.5%and the Microwave Humidity Sounder(MHS)humidity sounding channels improve by up to 0.4%.

Use of Microwave Radiances from Metop-C and Fengyun-3 C/D Satellites for a Northern European Limited-area Data Assimilation System

MetCoOp is a Nordic collaboration on operational Numerical Weather Prediction based on a common limited-area km-scale ensemble system. The initial states are produced using a 3-dimensional variational data assimilation scheme utilizing a large amount of observations from conventional in-situ measurements, weather radars, global navigation satellite system, advanced scatterometer data and satellite radiances from various satellite platforms. A version of the forecasting system which is aimed for future operations has been prepared for an enhanced assimilation of microwave radiances. This enhanced data assimilation system will use radiances from the Microwave Humidity Sounder, the Advanced Microwave Sounding Unit-A and the Micro-Wave Humidity Sounder-2 instruments on-board the Metop-C and Fengyun-3 C/D polar orbiting satellites. The implementation process includes channel selection, set-up of an adaptive bias correction procedure, and careful monitoring of data usage and quality control of observations. The benefit of the additional microwave observations in terms of data coverage and impact on analyses, as derived using the degree of freedom of signal approach, is demonstrated. A positive impact on forecast quality is shown, and the effect on the precipitation for a case study is examined. Finally, the role of enhanced data assimilation techniques and adaptions towards nowcasting are discussed.