Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but t...Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.展开更多
Background: The Chapman-Richards distribution is developed as a special case of the equilibrium solution to the McKendrick-Von Foerster equation. The Chapman-Richards distribution incorporates the vital rate assumptio...Background: The Chapman-Richards distribution is developed as a special case of the equilibrium solution to the McKendrick-Von Foerster equation. The Chapman-Richards distribution incorporates the vital rate assumptions of the Chapman-Richards growth function, constant mortality and recruitment into the mathematical form of the distribution. Therefore, unlike 'assumed' distribution models, it is intrinsically linked with the underlying vital rates for the forest area under consideration. Methods: It is shown that the Chapman-Richards distribution can be recast as a subset of the generalized beta distribution of the first kind, a rich family of assumed probability distribution models with known properties. These known properties for the generalized beta are then immediately available for the Chapman-Richards distribution, such as the form of the compatible basal area-size distribution. A simple two-stage procedure is proposed for the estimation of the model parameters and simulation experiments are conducted to validate the procedure for four different possible distribution shapes. Results: The simulations explore the efficacy of the two-stage estimation procedure;these cover the estimation of the growth equation and mortality-recruitment derives from the equilibrium assumption. The parameter estimates are shown to depend on both the sample size and the amount of noise imparted to the synthetic measurements. The results vary somewhat by distribution shape, with the smaller, noisier samples providing less reliable estimates of the vital rates and final distribution forms. Conclusions: The Chapman-Richards distribution in its original form, or recast as a generalized beta form, presents a potentially useful model integrating vital rates and stand diameters into a flexible family of resultant distributions shapes. The data requirements are modest, and parameter estimation is straightforward provided the minimal recommended sample sizes are obtained.展开更多
Several equations were selected using nonlinear regression analysis for setting up growth and yield modehe of Dahurian larch (Laris gmelinii Rupr.) plantations. Data of 405 stem analysis trees were collected from 336 ...Several equations were selected using nonlinear regression analysis for setting up growth and yield modehe of Dahurian larch (Laris gmelinii Rupr.) plantations. Data of 405 stem analysis trees were collected from 336 temporary plots throughout the Daxing’an Mountains. Results showed that the Richards equation was the best model for estimating tree height, stand mean height and stand dominant height by age, the Power equation was the fdiest model for predicting tree volume by DBH and tree height, and the Logarithmic stand vofume equation was good for predicting stand volume from age, mean height. basal area and other stand variables. These models can be used to construct volume tabIes, site index table and other forestry tables for Dahurian ghantations.展开更多
Growth and yield models were developed for individual tress and stands based on 336 temporary plots with 405 stem analysis trees of dahurian larch ( Larix gmelinii( Rupr. )Rupr.) plantations throughout Daxing'anli...Growth and yield models were developed for individual tress and stands based on 336 temporary plots with 405 stem analysis trees of dahurian larch ( Larix gmelinii( Rupr. )Rupr.) plantations throughout Daxing'anling mountains. Several equations were selected using nonlinear regression analysis. Results showed that the Richards equation was the best model for estimating tree height, stand mean height and stand dominant height from age; The Power equation was the best model for prediction tree volume from DBH and tree height; The logarithmic stand volume equation was good for predicting stand volume from age, mean height, basal area and other stand variables. These models can be used to construct the volume table, the site index table and other forestry tables for dahurian larch plantations.展开更多
基金funding support from the science and technology innovation Program of Hunan Province(Grant No.2023RC1017)Hunan Provincial Postgraduate Research and Innovation Project(Grant No.CX20220109)National Natural Science Foundation of China Youth Fund(Grant No.52208378).
文摘Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.
基金partially supported by the USDA National Institute of Food and Agriculture,Mc Intire Stennis Project OKL0 3063the Division of Agricultural Sciences and Natural Resources at Oklahoma State Universityprovided by the USDA Forest Service,Research Joint Venture 17-JV-11242306045,Old-Growth Forest Dynamics and Structure,to Mark Ducey
文摘Background: The Chapman-Richards distribution is developed as a special case of the equilibrium solution to the McKendrick-Von Foerster equation. The Chapman-Richards distribution incorporates the vital rate assumptions of the Chapman-Richards growth function, constant mortality and recruitment into the mathematical form of the distribution. Therefore, unlike 'assumed' distribution models, it is intrinsically linked with the underlying vital rates for the forest area under consideration. Methods: It is shown that the Chapman-Richards distribution can be recast as a subset of the generalized beta distribution of the first kind, a rich family of assumed probability distribution models with known properties. These known properties for the generalized beta are then immediately available for the Chapman-Richards distribution, such as the form of the compatible basal area-size distribution. A simple two-stage procedure is proposed for the estimation of the model parameters and simulation experiments are conducted to validate the procedure for four different possible distribution shapes. Results: The simulations explore the efficacy of the two-stage estimation procedure;these cover the estimation of the growth equation and mortality-recruitment derives from the equilibrium assumption. The parameter estimates are shown to depend on both the sample size and the amount of noise imparted to the synthetic measurements. The results vary somewhat by distribution shape, with the smaller, noisier samples providing less reliable estimates of the vital rates and final distribution forms. Conclusions: The Chapman-Richards distribution in its original form, or recast as a generalized beta form, presents a potentially useful model integrating vital rates and stand diameters into a flexible family of resultant distributions shapes. The data requirements are modest, and parameter estimation is straightforward provided the minimal recommended sample sizes are obtained.
文摘Several equations were selected using nonlinear regression analysis for setting up growth and yield modehe of Dahurian larch (Laris gmelinii Rupr.) plantations. Data of 405 stem analysis trees were collected from 336 temporary plots throughout the Daxing’an Mountains. Results showed that the Richards equation was the best model for estimating tree height, stand mean height and stand dominant height by age, the Power equation was the fdiest model for predicting tree volume by DBH and tree height, and the Logarithmic stand vofume equation was good for predicting stand volume from age, mean height. basal area and other stand variables. These models can be used to construct volume tabIes, site index table and other forestry tables for Dahurian ghantations.
文摘Growth and yield models were developed for individual tress and stands based on 336 temporary plots with 405 stem analysis trees of dahurian larch ( Larix gmelinii( Rupr. )Rupr.) plantations throughout Daxing'anling mountains. Several equations were selected using nonlinear regression analysis. Results showed that the Richards equation was the best model for estimating tree height, stand mean height and stand dominant height from age; The Power equation was the best model for prediction tree volume from DBH and tree height; The logarithmic stand volume equation was good for predicting stand volume from age, mean height, basal area and other stand variables. These models can be used to construct the volume table, the site index table and other forestry tables for dahurian larch plantations.
