Multivariate Drought Risk Analysis for the Weihe River: Comparison between Parametric and Nonparametric Copula Methods

: This study analyzed the multivariate drought risks for the Wei River basin by characterizing the interdependence between the standardized precipitation index (SPI) and the standardized precipitation evapotranspiration index (SPEI). Both parametric and nonparametric copulas were adopted to quantify the dependence between the SPI and SPEI. The results indicated that the Gaussian copula demonstrated the best ﬁ t in most cases, while the nonparametric copula method showed superiority over the parametric models at only one out of eighteen meteorological stations. The joint return periods (T OR , T AND , and T Kendall ) were computed through copula modeling, providing valuable insights into the co-occurrence of extreme drought events. For the SPI and SPEI with a 50-year return period, the T OR values range from 25.5 to 37.9 years, the T AND values ﬂ uctuate between 73.4 and 1233 years, and the T Kendall values range from 60.61 to 574.71 years, indicating a high correlation between the SPI and SPEI in the study area. The spatial analysis revealed varying pa tt erns across the basin with some regions more prone to experiencing simultaneous drought conditions characterized by both the SPI and SPEI. Furthermore, our results indicated that the SPEI exhibited more severity in drought characterization than the SPI due to its consideration of temperature e ﬀ ects. The disparities in the spatial features of the SPI and SPEI underscore the importance of incorporating multiple meteorological factors for a comprehensive drought risk analysis. This research contributes to a be tt er understanding of the drought pa tt erns and their joint risks in the Wei River basin, o ﬀ ering valuable information for drought preparedness and water resource management.


Introduction
Drought is a natural disaster that poses significant challenges to societies, economies, and ecosystems worldwide.In recent years, the increasing frequency and severity of drought events have raised concerns among researchers and policymakers [1,2].Monitoring and assessing drought risk are crucial for the development of effective mitigation and adaptation strategies.The standardized precipitation index (SPI) and standardized precipitation evapotranspiration index (SPEI) are two widely used drought indices in the fields of hydrology and meteorology [3,4].The SPI characterizes meteorological drought by measuring the deviation of precipitation from its long-term average, while the SPEI incorporates both precipitation and potential evapotranspiration, offering a comprehensive view of hydrological drought [4,5].These indices have been proven valuable for drought monitoring and early warning systems, enabling decision-makers to implement timely measures to alleviate the impacts of drought [6][7][8][9].However, assessing drought risk based solely on individual indices may not provide a complete understanding of the complex interactions and dependencies between meteorological and hydrological drought events.This limitation calls for a more integrated approach that considers the joint behavior of the SPI and SPEI to provide a more comprehensive characterization of drought risk [10].
Copula analysis has emerged as a powerful tool in the fields of hydrology and drought risk assessment, allowing researchers to model the joint behavior of multiple variables accurately.The application of copulas in hydrology dates back to the early 1990s, but it has gained significant traction in recent years due to its ability to capture the complex dependencies between hydrological processes and climatic variables [11,12].Copulabased methods have proven valuable in analyzing extreme events, providing a more comprehensive understanding of droughts and their interactions with other hydroclimatic variables.In the context of drought risk assessment, copulas have been applied to investigate the joint behavior of various drought indices, such as the SPI and SPEI, and other relevant hydroclimatic variables.Wang et al. [13] used a new copula-based standardized precipitation evapotranspiration streamflow index (SPESI) to characterize meteorological and hydrological drought in the Yellow River basin.Their findings highlighted the potential benefits of considering the joint probabilities of drought events, leading to more effective drought management strategies.Copula-based approaches have also been employed in the estimation of drought return periods and frequency analysis, providing valuable insights into the long-term behavior of drought events [14,15].Additionally, copula methods have been used to explore the dependence structure of other correlated variables, such as compound flood risks, climate downscaling, and so on [16][17][18].The flexibility of copulas allows researchers to select suitable copula families based on the data characteristics and research objectives.Archimedean copulas, such as the Clayton, Gumbel, and Frank copulas, are commonly used in drought risk analysis due to their simplicity and ability to model tail dependencies [19].In addition to Archimedean copulas, Gaussian copulas have also been widely applied in multivariate modeling of environmental variables and provide greater flexibility in capturing complex, non-monotonic relationships between variables [20,21].
Despite its many advantages, copula analysis in drought risk assessment faces certain challenges and limitations.The selection of an appropriate copula function and the determination of the optimal parameters can be sensitive to the dataset and require careful consideration [19].Recently, nonparametric copula methods have been proposed for coastal flood risk analysis to avoid misspecification issues associated with parametric copula models [22,23].These studies have demonstrated the applicability of nonparametric copula methods in capturing the interdependence among flood variables.However, the demonstration of nonparametric copula methods for drought risk analysis remains limited, and there is a lack of performance comparison between parametric and nonparametric copulas for multivariate drought risk inference.
Consequently, the objective of this study is to compare the performances of both parametric and nonparametric copulas on multivariate drought risk analysis for the Wei River basin.The drought events will be characterized by the SPI and SPEI values based on monthly rainfall and temperature observations.The annual minimum SPEI and the corresponding SPI would be further identified to reflect the most severe drought event in each year.The interdependence between the SPI and SPEI will be modeled using both parametric and nonparametric copulas, and their performances will be evaluated based on goodness-of-fit measures, such as the Akaike information criterion (AIC).Here, the AIC would be used to identify the most appropriate marginal and copula models since this method is able to evaluate the fitting effect of different distributions and penalize the model with more parameters.The best copula models will then be employed to reveal the joint risks of compound droughts, as reflected by the SPI and SPEI, for the Wei River basin.
This comparative analysis of copula methods for drought risk assessment will contribute to the existing literature by shedding light on the suitability of both parametric and nonparametric approaches for modeling drought dependencies.The findings of this study will provide valuable insights for drought management and decision-making in the Wei River basin and other regions facing similar hydroclimatic challenges.

Copula Method
Copula analysis is a powerful statistical technique used to model the dependence structure between multivariate random variables while preserving their individual marginal distributions.Consider d correlated random variables X1, X2, …, Xd with their marginals denoted as F1, …, Fd, their joint cumulative distribution function (CDF) F(.) can be constructed through a copula function as follows [24]: where C(.) is the copula function, which is unique if Fi (i = 1, 2, …, d) is continuous.The multivariate probability density function (PDF) f(.) would also be formulated as follows [25]: where ( ,  , … ,  ) is the copula density function, and ui = Fi(xi). ( ) is the marginal PDF for random variable Xi.

Parametric Copulas
The Archimedean copulas are a class of copulas that are commonly used due to their simplicity and tractability.A bivariate Archimedean copula is defined as follows [26]: where ϕ is the copula generator that is a convex decreasing function with ϕ(1) = 0 and ϕ −1 (.) = 0 when u2 ≥ ϕ(0); θ is the parameter hidden in the generating function [26].The most commonly used Archimedean copulas include the Clayton, Gumbel, and Frank copulas.
In addition to the Archimedean copulas, the Gaussian copula is also commonly used when dealing with continuous variables.The joint CDF of the Gaussian copula is given as follows: where Φ denotes the multivariate Gaussian CDF with correlation matrix Σ, and Φ −1 is the inverse standard Gaussian CDF.
Table 1 presents the formulations for the Gaussian, Frank, Gumbel, and Joe copulas as well as their basic properties.All the copulas are characterized by specific mathematical functions with a fixed number of parameters.These copulas assume a predefined functional form, which simplifies the modeling process and parameter estimation.In parametric copula modeling, the selection of the appropriate copula function and the estimation of its parameters play a crucial role in accurately capturing the dependence structure between variables.Even though there is a large variety of parametric copula models, they notoriously lack flexibility and bear the risk of misspecification [27].

Nonparametric Copulas
A specific class of nonparametric density estimators are kernel estimators, which have been applied for exploratory data analysis and widely used in many disciplines [27,28].Consider n observations (Ui1, Ui2), i = 1, …, n, from a bivariate copula C, the corresponding density function c(u1, u2) can be estimated through the kernel density estimator as follows [27]: where the notation The kernel function K is typically a symmetric, bounded probability density function on R 2 , and bn > 0 is the smoothing or bandwidth parameter.
The estimator in Equation ( 5) will result in a considerable amount of probability mass outside the unit square, leading to ̂ ( ,  ) not being a valid density function on [0, 1] 2 due to its integral not equaling one [27].Furthermore, the estimator will also suffer from severe bias at the boundaries [27].Some approaches have been developed to tackle the above challenges, including the mirror-reflection method, the beta kernel method, and the transformation method [29,30].The kernel density estimator allows us to estimate the bivariate copula density function nonparametrically, avoiding the need to assume a specific functional form and providing a flexible approach to capture the underlying dependence structure of the data.

Primary and Secondary Return Periods
If appropriate copula functions are selected to capture the joint probabilistic characteristics among correlated extreme variables, it becomes possible to derive conditional, primary, and secondary return periods.Specifically, joint (primary) return periods, referred to as "OR" and "AND" return periods, can be formulated as follows [26,[31][32][33]: where µ is the mean inter-arrival time of the two consecutive extreme events.Additionally, the secondary return period, referred to as the Kendall's return period, is defined as follows [26,[31][32][33]: where KC is Kendall's distribution associated with the theoretical Copula function C(.).

Overview of Wei River Basin
The Wei River basin, as shown in Figure 1, is one of the major river basins in China, encompassing an extensive area in the central and northwest regions.It is a crucial water resource for both agricultural and industrial activities in the region.The basin is characterized by a semi-arid to arid climate, experiencing irregular precipitation patterns with considerable variations in both temporal and spatial distributions.This region is prone to drought events due to its climate characteristics and limited water resources.Droughts in the region can have severe impacts on agriculture, water supply, and the ecosystem, posing significant challenges for water resource management and drought risk assessment [34].Over the past few decades, the Wei River basin has experienced varying degrees of drought severity, ranging from mild to severe drought events [35].These drought occurrences have highlighted the need for effective drought risk assessment and management strategies to mitigate the socio-economic and environmental impacts of drought in the region.

Data Collection and Drought Identification
The historical rainfall and temperature data between 1971 and 2016 were collected from the China Meteorological Data Service Center (http://data.cma.cn/en,accessed on 30 September 2022).The monthly rainfall and temperature data from a total of 18 meteorological stations located around the Wei River basin were utilized in this study for drought risk inference.The geographical positions of these stations are indicated in Figure 1 with the detailed locations and elevations presented in Table 2.Both the SPI and SPEI were adopted to characterize the historical drought events.For the SPI, the gamma probability distribution function was used to calculate the probability distribution of the precipitation in a specific month as follows [3,36]: / ,  > 0 (10a) where β > 0 and γ > 0 are scale and shape parameters, respectively, calculated as follows: where xi is the monthly precipitation, n refers to the sample size, and ̅ is the average of the precipitation samples.The probability of the random variable x less than x0 can be derived as follows: The gamma probability distribution in Equation ( 11) is normalized as follows: The gamma function does not include the case of x = 0, and the probability distribution with the x = 0 cases would be modified as follows: where q is the probability for zero precipitation, and F(x) is the probability from the gamma distribution.The SPI can finally be derived as follows [37]: where the constants are c0 = 2.515517, c1 = 0.802853, c2 = 0.010328, d1 = 1.432788, d2 = 0.189269, and d3 = 0.001308.The Gamma distributions for the monthly precipitation were tested with the Anderson-Darling test with their p-values provided in Table S1.The results suggest that the Gamma distributions can pass the statistical test for most months at each station except that some rejections occurred in December at some stations.However, the Gamma distributions would still be applicable in this study since only the annual minimum SPI values were analyzed, which seldom occurred in December.
The SPI has its limitations since it only considers the precipitation in its calculation process and cannot fully reflect the impact of climate warming on drought [4].Consequently, the SPEI has been developed based on the monthly water balance to reflect the impact of surface evaporation changes in order to compensate for the lack of the SPI [4,36].
In order to obtain the monthly water balance, the potential evapotranspiration (PET) was first calculated with the Hargreaves model as follows [38,39]: where Tmean is the average air temperature (°C); Tmax and Tmin are the maximum and minimum air temperatures (°C), respectively; and Ra is the daily net radiation on the land surface (MJ m −2 d −1 ).The water balance can then be obtained as follows: Where Di, Pi, and PETi respectively denote the water balance, monthly precipitation, and monthly potential evapotranspiration.In this study, the GEV distribution was employed to normalize the water balance series (i.e., Di) with the density function expressed as follows [38]: where z(x) = (x − µ)/σ; ξ, σ, and µ are the shape, scale, and location parameters, respectively.Based on the cumulative distribution function F(x) of GEV, the SPEI can be easily standardized as follows [38]: where W = −2ln () for P ≤ 0.5, and P = 1 − F(x).If P > 0.5, then P is replaced with 1 -P, and the sign of the resultant SPEI is reversed.The constants are c0 = 2.515517, c1 = 0.802853, c2 = 0.010328, d1 = 1.432788, d2 =0.189269, and d3 =0.001308.The GEV distributions were tested with the Anderson-Darling test with their p-values provided in Table S2 [40].The results suggest that the GEV distributions can pass the statistical test for all months at each station, indicating its applicability to derive the SPEI values.

Probability Estimation of Individual Drought Index
Figure 2 presents the general procedures to derive the SPI and SPEI values and then predict the multivariate drought risks.In this study, the one-month SPI and SPEI values were derived respectively from the Gamma and GEV distributions.The annual minimum SPEI values as well as the corresponding SPI values were identified for the period of 1971-2016.Table 3 presents some basic statistics of the SPI and SPEI values as well as their Kendall correlations at the different stations.The results suggest that the droughts respectively characterized by the SPI and SPEI are highly correlated, which raised the requirement to reveal the interdependence between these two kinds of drought.Figure 3 presents the temporal variations in the SPI and SPEI at stations 52986 and 52996, which also indicate the correlation between these two drought indicators.One of the major advantages of copula methods for multivariate drought risk analysis is that they allow for the quantification of marginal distributions of individual variables and their dependence structures in separate processes.In this study, multiple distribution methods were employed to quantify the random features of the annual minimum SPEI and SPI values.The goal was to identify the most appropriate distribution models for these two indices at different stations.
In detail, five distributions, namely the Gumbel, generalized extreme value (GEV), Gamma, lognormal, and Weibull distributions, were employed to quantify the distribution features of the SPI and SPEI.The performance of each distribution model was further evaluated using both the Kolmogorov-Smirnov (KS) test and the Akaike information criterion (AIC) to identify the most appropriate distributions for the SPI and SPEI.The KS test was used to assess how well each distribution fits the empirical data, while the AIC provided a measure to balance goodness-of-fit and model complexity, aiding in the selection of the best-fitting distributions for the drought indices.
Table 4 presents the marginal distribution selections for the SPI and SPEI at various stations within the Wei River basin.For the SPI, the generalized extreme value (GEV) distribution was predominantly identified as the most suitable model, which presented the lowest AIC values among the five candidate distributions for most stations.Also, the pvalues from the KS test are all larger than 0.05, showing the statistical applicability of the selected GEV models.Additionally, the Weibull distribution was chosen for three stations due to the lowest AIC values and high p-values at these stations.Similarly, for the SPEI, the GEV distribution was the primary choice, being deemed most appropriate for 13 out of the 18 stations.The Weibull distribution was the second option, considered suitable for three stations, while the choices for the remaining two stations were the Gamma and Gumbel distributions, respectively.Overall, the findings demonstrate that the GEV distribution serves as a robust choice for characterizing the marginal distributions of both the SPI and SPEI across the Wei River basin.However, it is evident that certain stations exhibited slight variations in their optimal distribution models, highlighting the spatial variability in drought characteristics within the basin.

Quantification of Interdependence between SPI and SPEI through Both Parametric and Nonparametric Copulas
As presented in Table 3, the SPI and SPEI are highly correlated with their Kendall's correlation larger than 0.3 for most stations and highest correlation approaching 0.7.In order to quantify the interdependence between the SPI and SPEI, both parametric copula models, as listed in Table 1, and nonparametric copula models were utilized.Similar to the selection process of marginal distributions, we employed the Kolmogorov-Smirnov (KS) test and the Akaike information criterion (AIC) to identify the most appropriate copula model at each station.For the nonparametric copula models, we calculated the effective number of parameters using the R package "kdecopula" [27], which helped us compute the AIC values for these models.Additionally, we used the root-mean-square error (RMSE) as an additional metric to evaluate the performance of the different copula models.
Table 5 presents the performance of both the best parametric copula and nonparametric copula models at different stations in the Wei River basin for multivariate drought risk analysis.The parametric copulas include the Gaussian, Gumbel, Frank, and Joe copulas, as shown in Table 1, in which the best parametric copulas with the lowest AIC values are presented in Table 5.In comparison, the nonparametric copula is expressed in Equation (5), and its performances at all stations are presented in Table 5.
The performance of the parametric copula models varied across the stations with different copula functions showing better fits in different cases.In general, the Gaussian copula was primarily selected, being deemed most appropriate for six out of the eighteen stations, followed by the Gumbel, Joe, and Frank copulas.
Regarding the nonparametric copula methods, they demonstrated competitive performance compared to the parametric models.However, it is important to note that the nonparametric model does not necessarily outperform the parametric copulas in quantifying the interdependence between the SPI and SPEI in the Wei River basin.Specifically, the nonparametric copula generated lower RMSE values at eight out of eighteen stations.Nevertheless, due to the inclusion of more effective parameters in the nonparametric copula, this model showed the least AIC value only at Station 53817.In other words, the nonparametric copula exhibited the best performance at only one station in this area.Table 6 presents the SPI and SPEI values with a 50-year return period (RP) and their corresponding joint return periods across the Wei River basin.Here, the 50-year return period is considered since (i) this RP would generally show severe drought events, and (ii) this RP has been analyzed in most drought research studies.Drought events with other RPs can be similarly generated with the proposed modelling method.
The analysis of the results reveals the severity of drought conditions in the region based on the SPI and SPEI values with a 50-year RP at various meteorological stations.The SPI values range from −2.14 to −2.64, while the SPEI values vary from −2.48 to −3.00, indicating the presence of significant and prolonged drought conditions across the basin.
Furthermore, the joint return periods, including T OR , T AND , and T Kendall , obtained through copula modeling provide valuable insights into the co-occurrence of extreme drought events.T OR represents the time of occurrence for droughts when either the SPI or SPEI falls below their respective 50-year RP thresholds.As shown in Equation ( 6), u1 = u2 = 0.98 since both the SPI and SPEI have an RP of 50 years.T OR can then be derived from the obtained copula model based on Equation (6).For instance, at station 52986, the Gumbel copula was selected; thus, we can have T OR = 1/[1 − exp {−[(−ln( )) + (−ln( ) ] / }] = 33.18years, where u1 = u2 = 0.98 and θ = 1.669, obtained in Section 4.2.In summary, the T OR values range from 25.5 to 37.9 years, indicating the timing of individual drought events below the SPI or SPEI with a 50-year RP.On the other hand, the T AND values fluctuate between 73.4 and 1233 years, representing the time of occurrence for droughts when both the SPI and SPEI are simultaneously below their 50-year RP thresholds.
Additionally, the T Kendall values, ranging from 60.61 to 574.71, indicate the likelihood of compound drought occurrences where both the SPI and SPEI experience extreme droughts simultaneously.These joint return periods emphasize the importance of considering the interdependence between the SPI and SPEI for a comprehensive understanding of drought risks in the Wei River basin. Figure 4 illustrates the spatial variations of the SPI and SPEI values with a 50-year RP in the Wei River basin, which were interpolated with the Kriging method based on station results.The results indicate that the SPEI appears to show more severity than the SPI as it takes into account the influence of temperature on drought conditions in addition to precipitation.Moreover, the SPI and SPEI exhibit different spatial features across the basin.In particular, the northwest region experiences the most severe droughts according to the SPEI, while the SPI identifies relatively milder drought severity in this area.Conversely, the central-south area exhibits the most severe droughts according to the SPI with drought severity characterized by the SPEI ranging between −2.6 and −2.7, approaching the least drought severity regions (i.e., central and southeast) as identified by the SPEI.
Figure 5 displays the spatial variations in the joint return periods in AND, OR, and Kendall for the SPI and SPEI with a 50-year RP in the Wei River basin, which are also interpolated with the Kriging method.Due to the distinct dependence patterns between the SPI and SPEI across the basin, the joint return periods in AND, OR, and Kendall exhibit different spatial variation features.As observed in Figure 5, the western and central-northern regions are more prone to experiencing simultaneous severe droughts characterized by both the SPI and SPEI with the T AND return periods less than 400 years.In contrast, the central-southeastern part has relatively fewer chances of experiencing simultaneous SPIand SPEI-based droughts with the T AND return period potentially exceeding 800 years.The spatial variations in TOR present a different feature compared to T AND .It indicates that the western, northeastern, and southeastern regions are less likely to encounter a 50-year drought represented solely by either the SPI or SPEI, while the central-southeastern and northwestern regions are more likely to experience a 50-year SPI or SPEI drought.
Regarding T Kendall , as shown in Figure 5c, its spatial variations are similar to those of T AND but with relatively shorter return periods.This suggests that the basin may experience compound drought events with shorter recurrence intervals, highlighting the possibility of concurrent extreme drought occurrences based on the interdependence between the SPI and SPEI.
These spatial variation features of the joint return periods emphasize the complex and region-specific nature of drought risks in the Wei River basin.Understanding these variations is crucial for developing targeted drought management strategies and implementing adaptive measures in different regions to enhance the basin's resilience to drought events.Policymakers and water resource managers can utilize this information to prioritize drought mitigation efforts and allocate resources effectively, considering both individual and joint drought occurrences in the region.

Conclusions
In this study, we conducted a comprehensive multivariate drought risk analysis in the Wei River basin using copula methods, specifically comparing the performances of parametric and nonparametric copulas.Our findings shed light on the interdependence between the standardized precipitation index (SPI) and the standardized precipitation evapotranspiration index (SPEI) and their implications for drought risk assessment in the region.
Through the selection of appropriate copula functions and marginal distributions, we quantified the joint probabilistic characteristics of the SPI and SPEI, enabling us to identify the most suitable copula models for each station.Our results demonstrated that the generalized extreme value (GEV) distribution was predominantly identified as the most appropriate marginal model for both the SPI and SPEI.The Gaussian copula emerged as the primary selection among the parametric copula models, while the nonparametric copulas only showed better performance at one out of eighteen meteorological stations in capturing the interdependence between the variables.
The joint return periods, including T OR , T AND , and Kendall's return period (T Kendall ), provided valuable insights into the co-occurrence of extreme drought events.We observed varying spatial patterns in the basin with certain regions more prone to experiencing concurrent drought conditions characterized by both the SPI and SPEI.Our spatial analysis also revealed that the SPEI exhibited more severity in drought characterization than the SPI, highlighting the significance of considering both precipitation and temperature factors in drought assessments.The disparities in the spatial features of the SPI and SPEI underscore the need for a comprehensive approach that incorporates multiple meteorological variables to enhance drought risk analysis accuracy.Overall, our study contributes to a better understanding of the drought patterns and their joint risks in the Wei River basin.The copula-based approach demonstrated its effectiveness in quantifying the interdependence between the SPI and SPEI, providing valuable information for water resource management and drought resilience planning in the region.The insights gained from this research can serve as a basis for informed decision-making and the development of targeted drought mitigation and adaptation strategies.Policymakers and water resource managers can utilize this knowledge to implement region-specific measures and policies to combat the increasing drought risks in the basin effectively.
In conclusion, our study showcases the value of copula methods, including both parametric and nonparametric copulas, in multivariate drought risk analysis and emphasizes the importance of considering joint drought occurrences to enhance drought preparedness and water management strategies in the Wei River basin.Moreover, parametric copulas would be sufficient to analyze multivariate drought risks in the Wei River basin for most cases.Additionally, the comparison of parametric and nonparametric copulas for quantifying the interdependence between the 3-month SPI and SPEI is presented in Table S3.The results also demonstrated the applicability of the parametric copulas on investigating long-term multivariate drought risks.As climate change continues to influence hydrological patterns, non-stationarity may exist in both the SPI and SPEI series.Consequently, more studies are required to explore whether parametric copulas would still have better performances than nonparametric models under nonstationary conditions.

Figure 1 .
Figure 1.The location of the Wei River basin.

Figure 2 .
Figure 2. The procedures for multivariate drought risk analysis.

Figure 3 .
Figure 3. Temporal variations in the SPI and SPEI at stations 52986 and 52996.

Figure 4 .
Figure 4. Spatial variation in the SPI and SPEI with a return period of 50 years.

Figure 5 .
Figure 5. Spatial variation in joint return periods for the SPI and SPEI with a 50-year RP.

Table 1 .
Basic properties of parametric copulas.

Table 3 .
Basic statistics of the SPI and SPEI values and their correlations at different stations.

Table 4 .
Selection of marginal distribution at different stations.

Table 6 .
SPI and SPEI values with a 50-year RP and their corresponding joint return periods.