• Animal Bioscience
• /
• v.35 no.5
• /
• pp.648-658
• /
• 2022

Objective: The identification of nonlinear mixed models that describe the growth trajectory of New Zealand rabbits was performed based on weight records and carcass measures obtained using ultrasonography. Methods: Phenotypic records of body weight (BW) and loin eye area (LEA) were collected from 66 animals raised in a didactic-productive module of cuniculture located in the southern Piaui state, Brazil. The following nonlinear models were tested considering fixed parameters: Brody, Gompertz, Logistic, Richards, Meloun 1, modified Michaelis-Menten, Santana, and von Bertalanffy. The coefficient of determination (R²), mean squared error, percentage of convergence of each model (%C), mean absolute deviation of residuals, Akaike information criterion (AIC), and Bayesian information criterion (BIC) were used to determine the best model. The model that best described the growth trajectory for each trait was also used under the context of mixed models, considering two parameters that admit biological interpretation (A and k) with random effects. Results: The von Bertalanffy model was the best fitting model for BW according to the highest value of R² (0.98) and lowest values of AIC (6,675.30) and BIC (6,691.90). For LEA, the Logistic model was the most appropriate due to the results of R² (0.52), AIC (783.90), and BIC (798.40) obtained using this model. The absolute growth rates estimated using the von Bertalanffy and Logistic models for BW and LEA were 21.51g/d and 3.16 cm², respectively. The relative growth rates at the inflection point were 0.028 for BW (von Bertalanffy) and 0.014 for LEA (Logistic). Conclusion: The von Bertalanffy and Logistic models with random effect at the asymptotic weight are recommended for analysis of ponderal and carcass growth trajectories in New Zealand rabbits. The inclusion of random effects in the asymptotic weight and maturity rate improves the quality of fit in comparison to fixed models.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Journal of the Korean GEO-environmental Society
• /
• v.11 no.1
• /
• pp.19-26
• /
• 2010

The batch experiments and response surface methodology (RSM) have been applied to the investigation of the cadmium (Cd) adsorption by coal fly ash (CFA). CFA having maximum Cd removal mass of 8.51 mg/g were calculated from Langmuir model. Cd removal reaction with different initial pH ranged from 4 to 9. When the initial pH was higher, Cd was removed more by adsorption and precipitation. These results suggest that the lower pH cause an increase of $H^+$ ion concentration which competed with Cd ions for exchange sites in CFA. Also, The Cd adsorption was mathematically described as a function of parameters initial Cd concentration ($X_1$), initial pH ($X_2$), and initial CFA mass ($X_3$) being modeled by use of the Box-Behnken methods. Empirical models were developed to describe relationship between the experimental variables and response. Statistical analysis indicates that tree factors ($X_1$, $X_2$, and $X_3$) on the linear term (main effects), and tree factors ($X_1X_2$, $X_1X_3$, and $X_2X_3$) on the non-linear term (Interaction effect; cross-product) had significant effects, respectively. In this case, the value of the adjusted determination coefficient (adjusted $R^2=0.9280$) was closed to 1, showing a high significance of the model. Statistical results showed the order of Cd removal at experimental factors to be initial initial pH > initial Cd concentration > initial CFA mass.

• Applied Biological Chemistry
• /
• v.43 no.2
• /
• pp.116-123
• /
• 2000

A numerical model of three-dimensional soil water distribution for drip irrigation management under cropped conditions was developed using Richards equation in Cartesian coordinates. The model accounts for both seasonal and diurnal changes in evaporation and transpiration, and the growth of plant root and the shape of root zone. Solutions were numerically approximated using the Crank-Nicolson implicit finite difference technique on the block-centered grid system and the Gauss-Seidel elimination in tandem. The model was tested under several conditions to allow the flow rates and configurations of drip emitters vary. In general, simulation results agreed well with experimental results and were as follows. The velocity of soil-water flow decreased drastically with distance from the drip source, and the rate of expansion of the wetted zone decreased rapidly during irrigation. The wetting front of wetted zone from a surface drip emitter traveled farther in vertical direction than in horizontal direction. Under this experimental weather condition, water use efficiency of a drip-irrigated apple field was greatest for 4-drip-emitter system buried at 25 cm, resulting from 10% increase in transpiration but 20% reduction in soil evaporation compared to those for surface 1-drip emitter system. Soil moisture retention curve obtained using disk tension infiltrometer showed significant difference from the curve obtained with pressure plate extractor.

• Journal of Korean Society of Forest Science
• /
• v.100 no.1
• /
• pp.62-69
• /
• 2011

Six commonly used non-linear growth functions were fitted to individual tree height-dbh data of eight major tree species measured by the $5^{th}$ National Forest Inventory in Chungcheong province. A total of 2,681 trees were collected from permanent sample plots across Chungcheong province. The available data for each species were randomly splitted into two sets: the majority (90%) was used to estimate model parameters and the remaining data (10%) were reserved to validate the models. The performance of the models was compared and evaluated by $R^2$, RMSE, mean difference (MD), absolute mean difference (AMD) and mean difference(MD) for diameter classes. The combined data (100%) were used for final model fitting. The results showed that these six sigmoidal models were able to capture the height-diameter relationships and fit the data equally well, but produced different asymptote estimates. Sigmoidal growth models such as Chapman-Richards, Weibull functions provided the most satisfactory height predictions. The effect of model performance on stem volume estimation was also investigated. Tree volumes of different species were computed by the Forest Resources Evaluation and Prediction Program using observed range of diameter and the predicted tree total height from the six models. For trees with diameter less than 30 cm, the six height-dbh models produced very similar results for all species, while more differentiation among the models was observed for large-sized trees.

• Journal of Korean Society of Forest Science
• /
• v.111 no.2
• /
• pp.311-318
• /
• 2022

We conducted this study to derive the site index, which is a criterion for the planting of Robinia pseudoacacia, a honey plant, and to investigate the diameter distribution change by derived site index. We applied the Chapman-Richards equation model to estimate the site index of the Robinia pseudoacacia stand. The site index was distributed within the range of 16-22 when the base age was 30 years. The fitness index of the site index estimation model was low, but we judged that there was no problem in the application because the residual distribution of the equation had not shifted to one side. We used the Weibull diameter distribution function to determine the diameter distribution of the Robinia pseudoacacia stand by site index. We used the mean diameter and the dominant tree height as independent variables to present the diameter distribution, and our analysis procedure was to estimate and recover the parameters of the Weibull diameter distribution function. We used the mean diameter and the dominant tree height of the Robinia pseudoacacia stand to show distribution by diameter class, and the fitness index for dbh distribution estimation was about 80.5%. As a result of schematizing the diameter distribution by site indices as a 30-year-old, we found that the higher the site index, the more the curve of the diameter distribution moved to the right. This suggests that if the plantation were to be established in a high site index stand, considering the suitable trees on the site, the growth of Robinia pseudoacacia woul d become active, and not onl y the production of wood but al so the production of honey would increase. We therefore anticipate that the site index classification table and curve of this Robinia pseudoacacia stand will become the standard for decision making in the plantation and management of this tree.

• Journal of Korean Society of Forest Science
• /
• v.109 no.4
• /
• pp.477-483
• /
• 2020

We constructed a stand yield table for Chamaecyparis obtusa based on data from an actual forest. The previous stand yield table had a number of disadvantages because it was based on actual forest information. In the present study we used data from more than 200 sampling plots in a stand of Chamaecyparis obtusa. The analysis included theestimation, recovery and prediction of the distribution of values for diameter at breast height (DBH), and the result is a valuable process for the preparation ofstand yield tables. The DBH distribution model uses a Weibull function, and the site index (base age: 30 years), the standard for assessing forest productivity, was derived using the Chapman-Richards formula. Several estimation formulas for the preparation of the stand yield table were considered for the fitness index, and the optimal formula was chosen. The analysis shows that the site index is in the range of 10 to 18 in the Chamaecyparis obtusa stand. The estimated stand volume of each sample plot was found to have an accuracy of 62%. According to the residuals analysis, the stands showed even distribution around zero, which indicates that the results are useful in the field. Comparing the table constructed in this study to the existing stand yield table, we found that our table yielded comparatively higher values for growth. This is probably because the existing analysis data used a small amount of research data that did not properly reflect. We hope that the stand yield table of Chamaecyparis obtusa, a representative species of southern regions, will be widely used for forest management. As these forests stabilize and growth progresses, we plan to construct an additional yield table applicable to the production of developed stands.

• Proceedings of the Korean Society of Crop Science Conference
• /
• 2017.06a
• /
• pp.315-315
• /
• 2017

In Japan, more than 80 % of soybean growing area is converted fields and excess water is one of the major problems in soybean production. For example, recent study (Yoshifuji et al., 2016) suggested that in the fields of shallow groundwater level (GWL) (< 1m depth), rising GWL even in a short period (e.g. 1 day) causes inhibition of soybean growth. Thus it becomes more and more important to predict GWL and soil moisture in detail. In addition to conventional surface drainage and underdrain, FOEAS (Farm Oriented Enhancing Aquatic System), which is expected to control GWL in fields adequately, has been developed recently. In this study we attempted to predict GWL and soil moisture condition at the converted field with FOEAS in Biwa lake reclamation area, Shiga prefecture, near the center of the main island of Japan. Two dimensional HYDRUS model (Simuinek et al., 1999) based on common Richards' equation, was used for the calculation of soil water movement. The calculation domain was considered to be 10 and 5 meter in horizontal and vertical direction, respectively, with two layers, i.e. 20cm-thick of plowed layer and underlying subsoil layer. The center of main underdrain (10 cm in diameter) was assumed to be 5 meter from the both ends of the domain and 10-60cm depth from the surface in accordance with the field experiment. The hydraulic parameters of the soil was estimated with the digital soil map in "Soil information web viewer" and Agricultural soil-profile physical properties database, Japan (SolphyJ) (Kato and Nishimura, 2016). Hourly rainfall depth and daily potential evapo-transpiration rate data were given as the upper boundary condition (B.C.). For the bottom B.C., constant upward flux, which meant the inflow flux to the field from outside, was given. Seepage face condition was employed for the surrounding of the underdrain. Initial condition was employed as GWL=60cm. Then we compared the simulated and observed results of volumetric water content at depth of 15cm and GWL. While the model described the variation of GWL well, it tended to overestimate the soil moisture through the growing period. Judging from the field condition, and observed data of soil moisture and GWL, consideration of soil structure (e.g. cracks and clods) in determination of soil hydraulic parameters at the plowed layer may improve the simulation results of soil moisture.

• Journal of Environmental Science International
• /
• v.17 no.8
• /
• pp.891-898
• /
• 2008

This study was conducted to develop the effects of site environmental factors on estimation of site index function for Chamaecyparis obtusa Endlicher stands. We derived nonlinear growth equation and the draw site index curves by applying this estimated equation. This study with Chapman-Richards function showed significant P-value which was less then 0.0001 and $R^2$ value 0.5947. This study was conducted to develop the feasible site index equation of Chamaecyparis obtusa Endlicher. For the table, the data of 82 sample areas that were thought to be without errors among the data of Chamaecyparis obtusa Endlicher sample area located on the value-oriented forest location chart were used and estimated. After analyzing the quantification method I based on 13 environmental factors to develop the score table for the site-index estimation of Chamaecyparis obtusa Endlicher, $R^2$ value of the model was 0.7555. It has been analyzed that the scope value of Soil moisture in horizon A was 7.5045, that of total soil depth was 6.3896, that of topography was 5.3471, that of slope was 4.7000 and that of aspect was 3.2038. After analyzing the partial correlation to examine the factors that affected most the site-index of Chamaecyparis obtusa Endlicher, it has been noted that the partial correlation of climatic zone was 0.4987, which was highest, and it was followed by Soil moisture in horizon A (0.4592), slope (0.4537), topography (0.3299) and total soil depth (0.1035). As a result of conducting the significance test for partial correlation, it has been found that topography, climatic zone, parent rock, slope, altitude, aspect, Soil moisture in horizon A, soil hardness in horizon A and total soil depth were recognized significant with 1% of significance level and sedimentary type and soil texture in horizon A were recognized significant with 5% of significance level.