A Copula Approach for Predicting Tree Sap Flow Based on Vapor Pressure Deficit
Abstract
:1. Introduction
2. Materials and Methods
2.1. Data Acquisition
2.2. Copula Analysis
- i.
- VPD data adjustment. Figure 3a shows the percentages of maximum diurnal peak values for the measured hourly sap flow and VPD. About 95% of the diurnal sap flow peaks occurred at 12 h, while about 95% of the diurnal VPD peaks took place at 16 h, indicating a 4 h lag with 95% maximum values in the sap flow and VPD data. To accurately predict hourly sap flows based on hourly VPDs, the VPDs need to be adjusted by shifting 4 h backward to align with the diurnal peaks of the sap flow (Figure 3b). This shift will obtain a better correlation between sap flow and VPD when randomly generating the correlation in copula analysis.
- ii.
- Histogram plot. A histogram plot showed a Gamma distribution of the sap flows and VPDs (Figure 4). Thus, the Gamma type of distribution was used as a marginal distribution function when building the bivariate distribution for the five copulas. In statistical analysis, the probability distribution of all its random variables is defined as a joint distribution, whereas the probability distribution of one random variable is called a marginal distribution. The marginal distribution functions play a vital role in determining dependence among random variables. Two random variables are dependent (or correlated) if and only if their joint distribution function is not equal to the product of their marginal distribution functions (https://www.statlect.com/glossary/marginal-distribution-function, accessed 10 April 2024).
- iii.
- Dependence of sap flow on VPD. In a copula analysis, the correlation (or dependence) of two or more variables is normally measured (or estimated) by Mann–Kendall’s tau (τ) and Spearman’s rho methods. In this study, Mann–Kendall’s τ value was used to measure the correlation (or dependence) between the sap flow and the VPD. The τ value ranges from −1 to 1. If τ = 0, no relationship exists, and if τ = 1 (or −1), a perfect relationship exists (with positive τ for an increasing trend and negative τ for a decreasing trend). The best copula function (with the highest τ value) was selected for further analysis.
- iv.
- Validation of the selected copula. The selected copula function was used to predict sap flows based on VPDs. The predicted sap flows were then compared with our field measurements. The goodness-of-validation was estimated with Mann–Kendall’s τ at p < 0.01.
- v.
- Randomly generate sap flow and VPD data. Once the selected copula was validated, it was used to randomly generate sap flow and VPD data.
- vi.
- Establish a copula-based regression equation. The randomly generated sap flow and VPD data (Step 5) were used to establish a copula-based regression equation in an Excel spreadsheet. This equation was applied to predict sap flows when the VPDs were given. It should be noted that the sap flows predicted by the copula-based regression equation based on VPDs do not include the time-series component. In the real world, however, sap flow varies with times in hourly, daily, monthly, and annual manners. In other words, the sap flows predicted by the copula-based regression equation cannot be directly used because they do not tell when the sap flows occurred. Therefore, the time-series component must be included in the copula analysis. Fortunately, VPD data are always associated with time series and therefore the sap flows predicted by the copula-based regression equation based on VPDs are time-series data.
3. Results
3.1. Copula Selection and Validation
3.2. Growing Season Sap Flow Prediction
3.3. Multiple Years Sap Flow Prediction
4. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Wullschleger, S.D.; Hanson, P.J.; Tschaplinski, T.J. Whole-plant water flux in understory red maple exposed to altered precipitation regimes. Tree Physiol. 1998, 18, 71–79. [Google Scholar] [CrossRef]
- Merlin, M.; Solarik, K.A.; Landhausser, S.M. Quantification of uncertainties introduced by data-processing procedures of sap flow measurements using the cut-tree method on a large mature tree. Agric. Forest. Meteorol. 2020, 287, 107926. [Google Scholar] [CrossRef]
- Oki, T.; Kanae, S. Global hydrological cycles and world water resources. Science 2006, 313, 1068–1072. [Google Scholar] [CrossRef] [PubMed]
- Ouyang, Y.; Leininger, T.D.; Moran, M. Impacts of reforestation upon sediment load and water outflow in the Lower Yazoo River Watershed, Mississippi. Ecol. Eng. 2013, 61, 394–406. [Google Scholar] [CrossRef]
- Phillips, N.; Oren, R. Intra- and inter-annual variation in transpiration of a pine forest. Ecol. Appl. 2001, 11, 385–396. [Google Scholar] [CrossRef]
- Wilcox, B.P.; Breshears, D.D.; Allen, C.D. Ecohydrology of a resource-conserving semiarid woodland: Effects of scale and disturbance. Ecol. Monogr. 2003, 73, 223–239. [Google Scholar] [CrossRef]
- Liu, X.; Zhang, B.; Zhuang, J.Y.; Han, C.; Zhai, L.; Zhao, W.R.; Zhang, J.C. The Relationship between Sap Flow Density and Environmental Factors in the Yangtze River Delta Region of China. Forests 2017, 8, 74. [Google Scholar] [CrossRef]
- Nadezhdina, N.; Cermák, J.; Ceulemans, R. Radial patterns of sap flow in woody stems of dominant and understory species: Scaling errors associated with positioning of sensors. Tree Physiol. 2002, 22, 907–918. [Google Scholar] [CrossRef]
- Wieser, G.; Gruber, A.; Oberhuber, W. Sap flow characteristics and whole-tree water use of across the treeline ecotone of the central Tyrolean Alps. Eur. J. Forest Res. 2014, 133, 287–295. [Google Scholar] [CrossRef]
- Ouyang, Y.; Leininger, T.D.; Renninger, H.; Gardiner, E.S.; Samuelson, L. A Model to Assess Eastern Cottonwood Water Flow Using Adjusted Vapor Pressure Deficit Associated with a Climate Change Impact Application. Climate 2021, 9, 22. [Google Scholar] [CrossRef]
- Chen, X.Y.; Miller, G.R.; Rubin, Y.; Baldocchi, D.D. A statistical method for estimating wood thermal diffusivity and probe geometry using in situ heat response curves from sap flow measurements. Tree Physiol. 2012, 32, 1458–1470. [Google Scholar] [CrossRef] [PubMed]
- Li, Y.; Guo, L.; Wang, J.; Wang, Y.; Xu, D.; Wen, J. An Improved Sap Flow Prediction Model Based on CNN-GRU-BiLSTM and Factor Analysis of Historical Environmental Variables. Forests 2023, 14, 1310. [Google Scholar] [CrossRef]
- Ouyang, Y. A Hybrid of Copula Prediction and Time Series Computation to Estimate Stream Discharge Based on Precipitation Data. J. Am. Water Resour. Assoc. 2022, 58, 471–484. [Google Scholar] [CrossRef]
- Sklar, A. Fonctions de Répartition à n Dimensions et Leursmarges; Publications de l’Institut de Statistique de l’Université de Paris: Paris, France, 1959. [Google Scholar]
- Alizadeh, H.; Mousavi, S.J.; Ponnambalam, K. Copula-Based Chance-Constrained Hydro-Economic Optimization Model for Optimal Design of Reservoir-Irrigation District Systems under Multiple Interdependent Sources of Uncertainty. Water Resour. Res. 2018, 54, 5763–5784. [Google Scholar] [CrossRef]
- Frees, E.W.; Valdez, E.A. Hierarchical Insurance Claims Modeling. J. Am. Stat. Assoc. 2008, 103, 1457–1469. [Google Scholar] [CrossRef]
- Genest, C.; Favre, A.C.; Beliveau, J.; Jacques, C. Metaelliptical copulas and their use in frequency analysis of multivariate hydrological data. Water Resour. Res. 2007, 43, W09401. [Google Scholar] [CrossRef]
- Salvadori, G.; De Michele, C. Frequency analysis via copulas: Theoretical aspects and applications to hydrological events. Water Resour. Res. 2004, 40, W12511. [Google Scholar] [CrossRef]
- Sun, C.Y. Bivariate Extreme Value Modeling of Wildland Fire Area and Duration. Forest Sci. 2013, 59, 649–660. [Google Scholar] [CrossRef]
- Cote, M.P.; Genest, C.; Omelka, M. Rank-based inference tools for copula regression, with property and casualty insurance applications. Insur. Math. Econ. 2019, 89, 1–15. [Google Scholar] [CrossRef]
- Masarotto, G.; Varin, C. Gaussian Copula Regression in R. J. Stat. Softw. 2017, 77, 1–26. [Google Scholar] [CrossRef]
- Parsa, A.R.; Klugman, S.A. Copula Regression. Var. Adv. Sci. Risk 2008, 5, 45–54. [Google Scholar]
- Murray, F.W. On the Computation of Saturation Vapor Pressure; Rand Corp.: Santa Monica, CA, USA, 1966. [Google Scholar]
- Nelson, R.B. Extremes of nonexchangeability. Stat. Pap. 2007, 48, 329–336. [Google Scholar] [CrossRef]
- Aas, K.; Czado, C.; Frigessi, A.; Bakken, H. Pair-copula constructions of multiple dependence. Insur. Math. Econ. 2009, 44, 182–198. [Google Scholar] [CrossRef]
- Bedford, T.; Cooke, R.M. Vines—A new graphical model for dependent random variables. Ann. Stat. 2002, 30, 1031–1068. [Google Scholar] [CrossRef]
- Mangiafico, S. Summary and Analysis of Extension Program Evaluation in R, Version 1. 2023. Available online: https://rcompanion.org/documents/RHandbookProgramEvaluation.pdf (accessed on 10 April 2024).
- Renard, B.; Lang, M. Use of a Gaussian copula for multivariate extreme value analysis: Some case studies in hydrology. Adv. Water Resour. 2007, 30, 897–912. [Google Scholar] [CrossRef]
- Vose, J.M.; Swank, W.T.; Harvey, G.J.; Clinton, B.D.; Sobek, C. Leaf Water Relations and Sapflow in Eastern Cottonwood (Populus deltoides Bartr.) Trees Planted for Phytoremediation of a Groundwater Pollutant. Int. J. Phytoremediat. 2000, 2, 53–73. [Google Scholar]
- Samuelson, L.J.; Stokes, T.A.; Coleman, M.D. Influence of irrigation and fertilization on transpiration and hydraulic properties of Populus deltoides. Tree Physiol. 2007, 27, 765–774. [Google Scholar] [CrossRef]
Copula Function | τ | p-Value |
---|---|---|
BB8 | 0.33 | 0.01 |
Clayton Copula | 0.21 | 0.01 |
Frank Copula | 0.16 | 0.01 |
Gumbel Copula | 0.04 | 0.01 |
Normal Copula | 0.59 | 0.01 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Ouyang, Y.; Sun, C. A Copula Approach for Predicting Tree Sap Flow Based on Vapor Pressure Deficit. Forests 2024, 15, 695. https://doi.org/10.3390/f15040695
Ouyang Y, Sun C. A Copula Approach for Predicting Tree Sap Flow Based on Vapor Pressure Deficit. Forests. 2024; 15(4):695. https://doi.org/10.3390/f15040695
Chicago/Turabian StyleOuyang, Ying, and Changyou Sun. 2024. "A Copula Approach for Predicting Tree Sap Flow Based on Vapor Pressure Deficit" Forests 15, no. 4: 695. https://doi.org/10.3390/f15040695
APA StyleOuyang, Y., & Sun, C. (2024). A Copula Approach for Predicting Tree Sap Flow Based on Vapor Pressure Deficit. Forests, 15(4), 695. https://doi.org/10.3390/f15040695