Temperature uncertainty models for land and sea surfaces can be developed based on statistical methods. In this paper, we developed a novel time-series temperature uncertainty model, which is the autoregressive moving average (ARMA) (1,1) model. The model was developed for an observed annual mean temperature anomaly X(t), which is a combination of a true (latent) global anomaly Y(t) for a year (t) and normal variable w(t). The uncertainty is taken as the variance of w(t), which was divided into land surface temperature (LST) uncertainty, sea surface temperature (SST) uncertainty, and the corresponding source of uncertainty. The ARMA model was analyzed and compared with autoregressive (AR) and autoregressive integrated moving average (ARIMA) for the data obtained from the NASA Goddard Institute for Space Studies Surface Temperature (GISTEMP) Analysis. The statistical analysis of the autocorrelation function (ACF), partial autocorrelation function (PACF), normal quantile–quantile (normal Q-Q) plot, density of the residuals, and variance of normal variable w(t) shows that ARMA (1,1) fits better than AR (1) and ARIMA (1, d, 1) for d = 1, 2.

Non-parametric estimation of copula density functions presents significant challenges. One issue is the unboundedness of certain copula density functions and their derivatives at the corners of the unit square. Another is the boundary bias inherent in kernel density estimation. This paper presents a kernel-based method for estimating bounded copula density functions, addressing boundary bias through the mirror reflection technique. Optimal smoothing parameters are derived via Asymptotic Mean Integrated Squared Error (AMISE) minimization and cross-validation, with theoretical guarantees of consistency and asymptotic normality. Two kernel smoothing strategies are proposed: the rule-of-thumb approach and least squares cross-validation (LSCV). Simulation studies highlight the efficacy of the rule-of-thumb method in bandwidth selection for copulas with unbounded marginal supports. The methodology is further validated through an application to the Wisconsin Breast Cancer Diagnostic Dataset (WBCDD), where LSCV is used for bandwidth selection.

Local likelihood methods are widely used to estimate calibration functions in conditional copula models. Recent work has established uniform stochastic equicontinuity and uniform convergence rates for local likelihood estimators of covariate-dependent copula parameters, yielding global consistency guarantees and supporting the stability of local optimization routines. This paper complements those results by deriving minimax lower bounds for uniform estimation over Hölder classes of calibration functions. Under mild regularity conditions on the copula family and the covariate design, we show that the minimax sup-norm risk over a compact covariate region is bounded below by the classical nonparametric rate for smooth functions on an s-dimensional domain. The proof combines a localized packing construction with a Fano–Le Cam testing argument, using second-order expansions of the conditional copula likelihood to control information distances. As a consequence, local polynomial likelihood estimators achieve the minimax rate up to the logarithmic factors inherent to uniform estimation, providing a sharp optimality justification for their use in conditional copula modeling.