Two Methods of Deriving Advanced Dynamic Site Equations Demonstrated on Inland Douglas-fir Site Curves

Authors

  • Chris J. Cieszewski Quantfor Co., Watkinsville, GA, USA

Keywords:

Site index, dynamic equations, integrability, base-age invariance, state-space models, UTADA, inland Douglas-fir

Abstract

Two complementary methods are used to derive base-age-invariant dynamic site equations. Method 1 tests the compatibility of separately proposed age-growth and site-response rates. Compatible rates are integrated to recover a common surface, from which the unobserved site coordinate is eliminated. An inland Douglas-fir example reconstructs the combined-base model of Cieszewski (2001) using Monserud's published coefficients; a second example demonstrates forward construction and the resulting parameter restrictions. Six solved routes relate this construction to invertible maps and scalar differential equations. Method 2 applies the Unified Theory of the Algebraic Differences Approaches (UTADA; Cieszewski 2021), starting with a rational yield–site relationship and introducing age through its parameters. This family is a subclass of compatible surfaces. The method recovers the three-parameter Cieszewski (2001) equation and its four-parameter branch with equal site exponents and nonzero coefficients. With three global parameters and estimated reference heights, the former approximates Monserud's five-parameter Douglas-fir-series curves with root mean square error (RMSE) of approximately 0.15 ft over ages 10–200. This is a deterministic curve approximation, not independent empirical validation. Mathematical compatibility, site inversion, and biological domain restrictions are distinguished from the statistical requirements of fitting and validating the equations.

References

Adame, P., Cañellas, I., Roig, S., and del Río, M. 2006. Modelling dominant height growth and site index curves for rebollo oak (Quercus pyrenaica Willd.). Ann. For. Sci. 63(8): 929–940.

Arnold, V.I. 1992. Ordinary Differential Equations, 3rd ed., transl. R. Cooke. Springer, Berlin/Heidelberg, 334 p.

Bailey, R.L., and Clutter, J.L. 1974. Base-age invariant polymorphic site curves. For. Sci. 20(2): 155–159.

Barrio-Anta, M., Castedo-Dorado, F., Diéguez-Aranda, U., Álvarez-González, J.G., Parresol, B.R., and Rodríguez-Soalleiro, R. 2006. Development of a basal area growth system for maritime pine in northwestern Spain using the generalized algebraic difference approach. Can. J. For. Res. 36(6): 1461–1474.

Barrio-Anta, M., Sixto-Blanco, H., Cañellas-Rey de Viñas, I., and Castedo-Dorado, F. 2008. Dynamic growth model for I-214 poplar plantations in the northern and central plateaux in Spain. For. Ecol. Manage. 255(3–4): 1167–1178.

Bhatia, H., Norgard, G., Pascucci, V., and Bremer, P.-T. 2013. The Helmholtz–Hodge decomposition—a survey. IEEE Trans. Vis. Comput. Graph. 19(8): 1386–1404.

Bravo-Oviedo, A., del Río, M., and Montero, G. 2007. Geographic variation and parameter assessment in generalized algebraic difference site index modelling. For. Ecol. Manage. 247(1–3): 107–119.

Bravo-Oviedo, A., Tomé, M., Bravo, F., Montero, G., and del Río, M. 2008. Dominant height growth equations including site attributes in the generalized algebraic difference approach. Can. J. For. Res. 38(9): 2348–2358.

Callen, H.B. 1985. Thermodynamics and an Introduction to Thermostatistics, 2nd ed. John Wiley & Sons, New York.

Castedo-Dorado, F., Diéguez-Aranda, U., Barrio-Anta, M., and Álvarez-González, J.G. 2007. Modelling stand basal area growth for radiata pine plantations in Northwestern Spain using the GADA. Ann. For. Sci. 64(6): 609–619.

Castillo López, A., Vargas-Larreta, B., Corral Rivas, J.J., Nájera Luna, J.A., Cruz Cobos, F., and Hernández, F.J. 2013. Modelo compatible altura–índice de sitio para cuatro especies de pino en Santiago Papasquiaro, Durango. Rev. Mex. Cienc. For. 4(18): 86–103.

Cieszewski, C.J. 1994. Development of a variable density height-growth-model through defining multidimensional height growth spaces. Ph.D. thesis, Univ. of Alberta, Edmonton, AB, 72 p.

Cieszewski, C.J. 2000. Analytical site index solution for the generalized log-logistic height equation. For. Sci. 46(2): 291–296.

Cieszewski, C.J. 2001. Three methods of deriving advanced dynamic site equations demonstrated on inland Douglas-fir site curves. Can. J. For. Res. 31(1): 165–173.

Cieszewski, C.J. 2002. Comparing fixed- and variable-base-age site equations having single versus multiple asymptotes. For. Sci. 48(1): 7–23.

Cieszewski, C.J. 2003. Developing a well-behaved dynamic site equation using a modified Hossfeld IV function Y³ = (axᵐ)/(c + xᵐ⁻¹), a simplified mixed-model and scant subalpine fir data. For. Sci. 49(4): 539–554.

Cieszewski, C.J. 2004. GADA derivation of dynamic site equations with polymorphism and variable asymptotes from Richards, Weibull, and other exponential functions. Univ. of Georgia PMRC-TR 2004-5, Athens, GA, 16 p.

Cieszewski, C.J. 2021. UTADA: Unified Theory of the Algebraic Differences Approaches—Derivation of dynamic site equations from direct yield-site relationships. Math. Comput. For. Nat.-Res. Sci. 13(1): 36–43.

Cieszewski, C.J., and Bailey, R.L. 2000. Generalized algebraic difference approach: theory based derivation of dynamic site equations with polymorphism and variable asymptotes. For. Sci. 46(1): 116–126.

Cieszewski, C.J., and Bella, I.E. 1989. Polymorphic height and site index curves for lodgepole pine in Alberta. Can. J. For. Res. 19(9): 1151–1160.

Cieszewski, C.J., Bella, I.E., and Walker, D. 1999. An implementation of new dynamic site index models in a company's timber supply analysis. For. Chron. 75(6): 981–983.

Cieszewski, C.J., Harrison, M., and Martin, S.W. 2000. Practical methods for estimating non-biased parameters in self-referencing growth and yield models. Univ. of Georgia PMRC Tech. Rep. 2000-7, Athens, GA, 10 p.

Cieszewski, C.J., and Nigh, G.D. 2002. A dynamic equation for a published Sitka spruce site-dependent height-age model. For. Chron. 78(5): 690–694.

Cieszewski, C.J., and Strub, M. 2008. Generalized algebraic difference approach derivation of dynamic site equations with polymorphism and variable asymptotes from exponential and logarithmic functions. For. Sci. 54(3): 303–315.

Cieszewski, C.J., Strub, M., and Zasada, M. 2007. New dynamic site equation that fits best the Schwappach data for Scots pine (Pinus sylvestris L.) in Central Europe. For. Ecol. Manage. 243(1): 83–93.

Cieszewski, C.J., and Zasada, M. 2002. Dynamiczna forma anamorficznego modelu bonitacyjnego dla sosny pospolitej w Polsce. Sylwan 146(7): 17–24.

Cieszewski, C.J., Zasada, M., and Strub, M. 2006. Analysis of different base models and methods of site model derivation for Scots pine. For. Sci. 52(2): 187–197.

Clutter, J.L., Fortson, J.C., Pienaar, L.V., Brister, G.H., and Bailey, R.L. 1983. Timber Management: A Quantitative Approach. John Wiley & Sons, New York, NY, 333 p.

Diéguez-Aranda, U., Burkhart, H.E., and Amateis, R.L. 2006a. Dynamic site model for loblolly pine (Pinus taeda L.) plantations in the United States. For. Sci. 52(3): 262–272.

Diéguez-Aranda, U., Burkhart, H.E., and Rodríguez-Soalleiro, R. 2005. Modeling dominant height growth of radiata pine (Pinus radiata D. Don) plantations in north-western Spain. For. Ecol. Manage. 215(1–3): 271–284.

Diéguez-Aranda, U., Grandas-Arias, J.A., Álvarez-González, J.G., and von Gadow, K. 2006b. Site quality curves for birch stands in north-western Spain. Silva Fenn. 40(4): 631–644.

Elfving, B., and Kiviste, A. 1997. Construction of site index equations for Pinus sylvestris L. using permanent plot data in Sweden. For. Ecol. Manage. 98(2): 125–134.

Eriksson, H., Johansson, U., and Kiviste, A. 1997. A site-index model for pure and mixed stands of Betula pendula and Betula pubescens in Sweden. Scand. J. For. Res. 12(2): 149–156.

Frankot, R.T., and Chellappa, R. 1988. A method for enforcing integrability in shape from shading algorithms. IEEE Trans. Pattern Anal. Mach. Intell. 10(4): 439–451.

García, O. 1994. The state-space approach in growth modelling. Can. J. For. Res. 24(9): 1894–1903.

Hosoya, Y. 2017. The relationship between revealed preference and the Slutsky matrix. J. Math. Econ. 70: 127–146.

Hossfeld, J.W. 1822. Mathematik für Forstmänner, Ökonomen und Cameralisten. Gotha, Bd. 4, S. 310.

Johansson, T. 1999. Site index curves for common alder and grey alder growing on different types of forest soil in Sweden. Scand. J. For. Res. 14(5): 441–453.

Krumland, B., and Eng, H. 2005. Site Index Systems for Major Young-Growth Forest and Woodland Species in Northern California. Calif. Dept. of Forestry and Fire Protection, Calif. For. Rep. No. 4, Sacramento, CA, 220 p.

López-Senespleda, E., Bravo-Oviedo, A., Alonso, R., and Montero, G. 2014. Modeling dominant height growth including site attributes in the GADA approach for Quercus faginea Lam. in Spain. For. Syst. 23(3): 494–499.

Martín-Benito, D., Gea-Izquierdo, G., del Río, M., and Cañellas, I. 2008. Long-term trends in dominant-height growth of black pine using dynamic models. For. Ecol. Manage. 256(5): 1230–1238.

Monserud, R.A. 1984. Height growth and site-index curves for inland Douglas-fir based on stem analysis data and forest habitat type. For. Sci. 30(4): 943–965.

Monserud, R.A. 1985. Applying height growth and site index curves for inland Douglas-fir. USDA For. Serv. Res. Pap. INT-347, Intermountain Research Station, Ogden, UT.

Neyman, J., and Scott, E.L. 1948. Consistent estimates based on partially consistent observations. Econometrica 16(1): 1–32.

Nord-Larsen, T. 2006. Developing dynamic site index curves for European beech (Fagus sylvatica L.) in Denmark. For. Sci. 52(2): 173–181.

Northway, S.M. 1985. Fitting site-index equations and other self-referencing functions. For. Sci. 31(1): 233–235.

Nunes, L., Patrício, M., Tomé, J., and Tomé, M. 2011. Modeling dominant height growth of maritime pine in Portugal using GADA methodology with parameters depending on soil and climate variables. Ann. For. Sci. 68(2): 311–323.

Olver, P.J. 1993. Applications of Lie Groups to Differential Equations, 2nd ed. Graduate Texts in Mathematics 107. Springer, New York, NY, 513 p.

Quiñonez-Barraza, G., De los Santos-Posadas, H.M., Cruz-Cobos, F., Velázquez-Martínez, A., Ángeles-Pérez, G., and Ramírez-Valverde, G. 2015. Índice de sitio con polimorfismo complejo para masas forestales de Durango, México. Agrociencia 49(4): 439–454.

Rodríguez-Carrillo, A., Cruz-Cobos, F., Vargas-Larreta, B., and Hernández, F.J. 2015. Compatible dominant height–site index model for juniper (Juniperus deppeana Steud.). Rev. Chapingo Ser. Cienc. For. Ambient. 21(1): 97–108.

Rojo-Alboreca, A., Cabanillas-Saldaña, A.M., Barrio-Anta, M., Notivol-Paíno, E., and Gorgoso-Varela, J.J. 2017. Site index curves for natural Aleppo pine forests in the central Ebro valley (Spain). Madera y Bosques 23(1): 143–159.

Seki, M., and Sakici, O.E. 2017. Dominant height growth and dynamic site index models for Crimean pine in the Kastamonu–Taşköprü region of Turkey. Can. J. For. Res. 47(11): 1441–1449.

Sghaier, T., Palahí, M., Garchi, S., Bonet, J.A., Ammari, Y., and Piqué, M. 2012. Modeling dominant height growth in planted Pinus pinea stands in Northwest of Tunisia. Int. J. For. Res. 2012: 902381.

Sharma, R.P., Brunner, A., Eid, T., and Øyen, B.-H. 2011. Modelling dominant height growth from national forest inventory individual tree data with short time series and large age errors. For. Ecol. Manage. 262(12): 2162–2175.

Sharma, R.P., Vacek, Z., Vacek, S., Jansa, V., and Kučera, M. 2017. Modelling individual tree diameter growth for Norway spruce in the Czech Republic using a generalized algebraic difference approach. J. For. Sci. 63(5): 227–238.

Sideris, T.C. 2013. Ordinary Differential Equations and Dynamical Systems. Atlantis Studies in Differential Equations 2. Atlantis Press, Paris, 225 p.

Socha, J., Ochał, W., Grabczyński, S., and Maj, M. 2015. Site index models for forest-forming tree species in Poland developed basing on the yield tables. Sylwan 159(8): 639–649.

Socha, J., and Orzeł, S. 2011. Dynamic site index curves for Scots pine stands in Niepołomice Primeval Forest. Sylwan 155(5): 301–312.

Socha, J., Pierzchalski, M., Bałaży, R., and Ciesielski, M. 2017. Modelling top height growth and site index using repeated laser scanning data. For. Ecol. Manage. 406: 307–317.

Socha, J., Tymińska-Czabańska, L., Grabska, E., and Orzeł, S. 2020. Site index models for main forest-forming tree species in Poland. Forests 11(3): 301.

Strub, M., and Cieszewski, C.J. 2006. Base-age invariance properties of two techniques for estimating the parameters of site index models. For. Sci. 52(2): 182–186.

Swesty, F.D. 1996. Thermodynamically consistent interpolation for equation of state tables. J. Comput. Phys. 127(1): 118–127.

Tait, D.E., Cieszewski, C.J., and Bella, I.E. 1988. The stand dynamics of lodgepole pine. Can. J. For. Res. 18(10): 1255–1260.

Vargas-Larreta, B., Aguirre-Calderón, O.A., Corral-Rivas, J.J., Crecente-Campo, F., and Diéguez-Aranda, U. 2013. A dominant height growth and site index model for Pinus pseudostrobus Lindl. in northeastern Mexico. Agrociencia 47(1): 91–106.

Vargas-Larreta, B., Álvarez-González, J.G., Corral-Rivas, J.J., and Aguirre-Calderón, O.A. 2010. Construcción de curvas dinámicas de índice de sitio para Pinus cooperi Blanco. Rev. Fitotec. Mex. 33(4): 343–351.

Published

2026-10-25

Issue

Section

Mathematical Modeling

How to Cite

Two Methods of Deriving Advanced Dynamic Site Equations Demonstrated on Inland Douglas-fir Site Curves. (2026). Mathematical and Computational Forestry & Natural-Resource Sciences (MCFNS), 0-15 (16). https://mcfns.com/index.php/Journal/article/view/341

Similar Articles

1-10 of 127

You may also start an advanced similarity search for this article.

Most read articles by the same author(s)