This study presents a way to parametrically model and simulate multivariate distributions under potential non-normality, autocorrelation and heteroscedasticity and illustrates its application to agricultural risk analysis. Specifically, the joint probability distribution (pdf) for West Texas irrigated cotton, corn, sorghum, and wheat production and prices is estimated and applied to evaluate the changes in the risk and returns of agricultural production in the region resulting from observed and predicted price and production trends. The estimated pdf allows for time trends on the mean and the variance and varying degrees of autocorrelation and non-normality (kurtosis and right- or left-skewness) in each of the price and production variables. It also allows for any possible price-price, production-production, or price-production correlation.