• Journal of Korean Society of Forest Science
• /
• v.112 no.2
• /
• pp.209-216
• /
• 2023

This study aimed to use national forest inventory data to develop a forest productivity index and yield prediction model of a Tilia amurensis stand. The site index displaying the forest productivity of the Tilia amurensis stand was developed as a Schumacher model, and the site index classification curve was generated from the model results; its distribution growth in Korea ranged from 8-16. The growth model using age as an independent variable for breast height and height diameter estimation was derived from the Chapman-Richards and Weibull model. The Fitness Indices of the estimation models were 0.32 and 0.11, respectively, which were generally low values, but the estimation-equation residuals were evenly distributed around 0, so we judged that there would be no issue in applying the equation. The stand basal area and site index of the Tilia amurensis stand had the greatest effect on the stand-volume change. These two factors were used to derive the Tilia amurensis stand yield model, and the model's determination coefficient was approximately 94%. After verifying the residual normality of the equation and autocorrelation of the growth factors in the yield model, no particular problems were observed. Finally, the growth and yield models of the Tilia amurensis stand were used to produce the makeshift stand yield table. According to this table, when the Tilia amurensis stand is 70 years old, the estimated stand-volume per hectare would be approximately 208 m³ . It is expected that these study results will be helpful for decision-making of Tilia amurensis stands management, which have high value as a forest resource for honey and timber.

• Journal of Korean Society of Forest Science
• /
• v.90 no.5
• /
• pp.650-662
• /
• 2001

The growth and yield models for five different kinds of natural forest types were systemically developed in the natural Broadleaved-Korean pine Forests in Northeast China. The data were collected from 359 temporary plots and 58 permanent plots with area ranged from 0.06 ha to 1.0 ha, ranging in stand age from 43 to 364 years. The Site Class Index (SCI) was introduced to evaluate site quality and the Crown Competition Factor (CCF) was selected as a measure of stand density for the mixed natural forest. The Chapman-Richards function was adopted to develop SCI equation and height-diameter curve. The Schumacher growth function was selected as base model to develop the DBH, basal area, and stand volume growth models by using re-parameterized method. In modeling mean DBH and basal area growth, it was found that the asymptotic parameter A of Schumacher function was exponentially related to site quality (SCI) and stand density (CCF). The rate parameter k was related to stand density and it was independent of SCI. Several validation measures for predicted stand variables were evaluated in the growth and yield models using independent data sets. The results indicated that relative mean errors (RME) in predicted stand attributes were less than ${\pm}5%$ and the estimated precision values of the stand variables were all greater than 95%.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• Wind and Structures
• /
• v.14 no.2
• /
• pp.167-186
• /
• 2011

This paper presents the results of measurements relating to the aerodynamic forces on flat square plates which were allowed to rotate at different speeds about their horizontal axis, by modifying the velocity of the incoming flow. A 1 m square test-sheet and a 0.3 m square test-sheet were fitted with a number of pressure sensors in order to obtain information relating to the instantaneous pressure distribution acting on the test-sheet; a compact gyroscope to record the angular velocity during the rotational motion was also implemented. Previous work on autorotation has illustrated that the angular velocity varies with respect to the torque induced by the wind, the thickness and aspect ratio of the test-sheet, any frictional effects present at the bearings, and the vorticity generated through the interaction between the plate and the wind flow. The current paper sets out a method based on the solution of the equation of motion of a rotating plate which enables the determination of angular velocities on autorotating elements to be predicted. This approach is then used in conjunction with the experimental data in order to evaluate the damping introduced by the frictional effects at the bearings during steady autorotation.

• Geomechanics and Engineering
• /
• v.10 no.2
• /
• pp.225-245
• /
• 2016

Settlement of foundations in permafrost regions primarily results from three physical and mechanical processes such as thaw consolidation of permafrost layer, creep of warm frozen soils and the additional deformation of seasonal active layer induced by freeze-thaw cycling. This paper firstly establishes theoretical models for the three sources of settlement including a statistical damage model for soils which experience cyclic freeze-thaw, a large strain thaw consolidation theory incorporating a modified Richards' equation and a Drucker-Prager yield criterion, as well as a simple rheological element based creep model for frozen soils. A novel numerical method was proposed for live computation of thaw consolidation, creep and freeze-thaw cycling in corresponding domains which vary with heat budget in frozen ground. It was then numerically implemented in the FISH language on the FLAC platform and verified by freeze-thaw tests on sandy clay. Results indicate that the calculated results agree well with the measured data. Finally a model test carried out on a half embankment in laboratory was modeled.

• Journal of Korean Society of Forest Science
• /
• v.88 no.3
• /
• pp.285-291
• /
• 1999

Loblolly(Pinus taeda L.) and slash(Pinus elliottii Engelm.) pines are the most important timber producing species in the Southern United States. Site index equations to estimate site index curves(base age 25 years) for loblolly pine and slash pine plantations have been developed based on long-term repeated measurement data sets. To check magnitude of errors in estimating site index, each cumulative measurement cycle data sets and all combined data sets were used to recalculate site index values. The Chapman-Richards' growth function was selected for stand height prediction. Anamorphic base age invariant site index curves were presented based on this height prediction equation. Statistics used in the evaluation were mean of the differences and mean of the absolute differences. For plantation ages less than 5 years, site index values showed very sensitive fur both species based on the evaluation test.

• 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.

• Journal of Korean Society of Forest Science
• /
• v.112 no.1
• /
• pp.83-92
• /
• 2023

This study was established to investigate the site environment of mixed forests in Korea and to estimate the growth and yield of stands using national forest resources inventory data. The growth of mixed forests was derived by applying the Chapman-Richards model with diameter at breast height (DBH), height, and cross-sectional area at breast height (BA), and the yield of mixed forests was derived by applying stepwise regression analysis with factors such as cross-sectional area at breast height, site index (SI), age, and standing tree density per ha. Mixed forests were found to be growing in various locations. By climate zone, more than half of them were distributed in the temperate central region. By altitude, about 62% were distributed at 101-400 m. The fitness indexes (FI) for the growth model of mixed forests, which is the independent variable of stand age, were 0.32 for the DBH estimation, 0.22 for the height estimation, and 0.18 for the basal area at breast height estimation, which were somewhat low. However, considering the graph and residual between the estimated and measured values of the estimation equation, the use of this estimation model is not expected to cause any particular problems. The yield prediction model of mixed forests was derived as follows: Stand volume =-162.6859+6.3434 ∙ BA+9.9214 ∙ SI+0.7271 ∙ Age, which is a step- by-step input of basal area at breast height (BA), site index (SI), and age among several growth factors, and the determination coefficient (R²) of the equation was about 96%. Using our optimal growth and yield prediction model, a makeshift stand yield table was created. This table of mixed forests was also used to derive the rotation of the highest production in volume.

• Journal of Korean Society of Forest Science
• /
• v.86 no.2
• /
• pp.135-145
• /
• 1997

This study aimed at evaluating one curvilinear equation and nine non-linear equations for estimating stand growth characteristics(mean dbh, mean height and volume per ha) for the plantations of Pinus koraiensis and the natural stands of Quercus mongolica. The data were collected from 92 plots in Pines koraiensis stands and 83 plots in Quercus mongolica stands, and the site index of all the stands is 14. The curvilinear equation, $Y=at^be^{-c/t}$, used in preparing the yield tables was well fitted within the range of data, but was likely to give overestimates when extrapolating in old stage due to the tendency of linear increase. Among the non-linear equations, logistic equation and Sloboda equation gave overestimates in young stands and reached the asymptotic status early which means underestimates in old stage. Extrapolating in old stage, Hossfeld equation generally gave larger values than others due to its large estimates of parameter a, the maximum value. On the other hand, Bertalanffy equation gave underestimates in young and old stands and overestimates in middle-aged stands. The estimates with Korf equation was relatively low for Pinus koraiensis stands, and this tendency was more obvious in dbh growth of Quercus mongolica stands. Ueno-Ohzaki equation was liable to give over or underestimates depending on the value of parameter b when extrapolating in old stands. Considering the accuracy of estimates and the biological base of the growth equations, Gompertz equation, Chapman-Richards equation and Weibull equation were generally applicable for estimating the stand growth characteristics of both species in the whole range of stand ages including extrapolated range. To get more accurate and precise parameter estimates, more data, especially in old stands, should be required in further study.

• 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.