Resumen
One of the ways to quantify uncertainty of deterministic forecasts is to construct a joint distribution between the forecast variable and the observed variable; then, the uncertainty of the forecast can be represented by the conditional distribution of the observed given the forecast. The joint distribution of two continuous hydrometeorological variables can often be modeled by the bivariate meta-Gaussian distribution (BMGD). The BMGD can be obtained by transforming each of the two variables to a standard normal variable and the dependence between the transformed variables is provided by the Pearson correlation coefficient of these two variables. The BMGD modeling is exact provided that the transformed joint distribution is standard normal. In real-world applications, however, this normality assumption is hardly fulfilled. This is often the case for the modeling problem we consider in this paper: establish the joint distribution of a forecast variable and its corresponding observed variable for precipitation amounts accumulated over a duration of 24 h. In this case, the BMGD can only serve as an approximate model and the dependence parameter can be estimated in a variety of ways. In this paper, the effect of tuning this parameter is studied. Numerical simulations conducted suggest that, while the parameter tuning results in limited improvements in goodness-of-fit (GOF) for the BMGD as a bivariate distribution model, better results may be achieved by tuning the parameter for the one-dimensional conditional distribution of the observed given the forecast greater than a certain large value.