No Arabic abstract
We study the impact of three stochastic parametrizations in the ocean component of a coupled model, on forecast reliability over seasonal timescales. The relative impacts of these schemes upon the ocean mean state and ensemble spread are analyzed. The oceanic variability induced by the atmospheric forcing of the coupled system is, in most regions, the major source of ensemble spread. The largest impact on spread and bias came from the Stochastically Perturbed Parametrization Tendency (SPPT) scheme - which has proven particularly effective in the atmosphere. The key regions affected are eddy-active regions, namely the western boundary currents and the Southern Ocean. However, unlike its impact in the atmosphere, SPPT in the ocean did not result in a significant decrease in forecast error. Whilst there are good grounds for implementing stochastic schemes in ocean models, our results suggest that they will have to be more sophisticated. Some suggestions for next-generation stochastic schemes are made.
Many major oceanographic internal wave observational programs of the last 4 decades are reanalyzed in order to characterize variability of the deep ocean internal wavefield. The observations are discussed in the context of the universal spectral model proposed by Garrett and Munk. The Garrett and Munk model is a good description of wintertime conditions at Site-D on the continental rise north of the Gulf Stream. Elsewhere and at other times, significant deviations in terms of amplitude, separability of the 2-D vertical wavenumber - frequency spectrum, and departure from the models functional form are noted. Subtle geographic patterns are apparent in deviations from the high frequency and high vertical wavenumber power laws of the Garrett and Munk spectrum. Moreover, such deviations tend to co-vary: whiter frequency spectra are partnered with redder vertical wavenumber spectra. Attempts are made to interpret the variability in terms of the interplay between generation, propagation and nonlinearity using a statistical radiative balance equation. This process frames major questions for future research with the insight that such integrative studies could constrain both observationally and theoretically based interpretations.
A~machine learning framework is developed to estimate ocean-wave conditions. By supervised training of machine learning models on many thousands of iterations of a physics-based wave model, accurate representations of significant wave heights and period can be used to predict ocean conditions. A model of Monterey Bay was used as the example test site; it was forced by measured wave conditions, ocean-current nowcasts, and reported winds. These input data along with model outputs of spatially variable wave heights and characteristic period were aggregated into supervised learning training and test data sets, which were supplied to machine learning models. These machine learning models replicated wave heights with a root-mean-squared error of 9cm and correctly identify over 90% of the characteristic periods for the test-data sets. Impressively, transforming model inputs to outputs through matrix operations requires only a fraction (<1/1,000) of the computation time compared to forecasting with the physics-based model.
Atmosphere and ocean are coupled via air-sea interactions. The atmospheric conditions fuel the ocean circulation and its variability, but the extent to which ocean processes can affect the atmosphere at decadal time scales remains unclear. In particular, such low-frequency variability is difficult to extract from the short observational record, meaning that climate models are the primary tools deployed to resolve this question. Here, we assess how the oceans intrinsic variability leads to patterns of upper-ocean heat content that vary at decadal time scales. These patterns have the potential to feed back on the atmosphere and thereby affect climate modes of variability, such as El Ni~no or the Interdecadal Pacific Oscillation. We use the output from a global ocean-sea ice circulation model at three different horizontal resolutions, each driven by the same atmospheric reanalysis. To disentangle the variability of the oceans direct response to atmospheric forcing from the variability due to intrinsic ocean dynamics, we compare model runs driven with inter-annually varying forcing (1958-2018) and model runs driven with repeat-year forcing. Models with coarse resolution that rely on eddy parameterizations, show (i) significantly reduced variance of the upper-ocean heat content at decadal time scales and (ii) differences in the spatial patterns of low-frequency variability compared with higher resolution models. Climate projections are typically done with general circulation models with coarse-resolution ocean components. Therefore, these biases affect our ability to predict decadal climate modes of variability and, in turn, hinder climate projections. Our results suggest that for improving climate projections, the community should move towards coupled climate models with higher oceanic resolution.
The Atlantic Meridional Overturning Circulation (AMOC) transports substantial amounts of heat into the North Atlantic sector, and hence is of very high importance in regional climate projections. The AMOC has been observed to show multi-stability across a range of models of different complexity. The simplest models find a bifurcation associated with the AMOC `on state losing stability that is a saddle node. Here we study a physically derived global oceanic model of Wood {em et al} with five boxes, that is calibrated to runs of the FAMOUS coupled atmosphere-ocean general circulation model. We find the loss of stability of the `on state is due to a subcritical Hopf for parameters from both pre-industrial and doubled CO${}_2$ atmospheres. This loss of stability via subcritical Hopf bifurcation has important consequences for the behaviour of the basin of attraction close to bifurcation. We consider various time-dependent profiles of freshwater forcing to the system, and find that rate-induced thresholds for tipping can appear, even for perturbations that do not cross the bifurcation. Understanding how such state transitions occur is important in determining allowable safe climate change mitigation pathways to avoid collapse of the AMOC.
A primary goal of the National Oceanic and Atmospheric Administration (NOAA) Warn-on-Forecast (WoF) project is to provide rapidly updating probabilistic guidance to human forecasters for short-term (e.g., 0-3 h) severe weather forecasts. Maximizing the usefulness of probabilistic severe weather guidance from an ensemble of convection-allowing model forecasts requires calibration. In this study, we compare the skill of a simple method using updraft helicity against a series of machine learning (ML) algorithms for calibrating WoFS severe weather guidance. ML models are often used to calibrate severe weather guidance since they leverage multiple variables and discover useful patterns in complex datasets. indent Our dataset includes WoF System (WoFS) ensemble forecasts available every 5 minutes out to 150 min of lead time from the 2017-2019 NOAA Hazardous Weather Testbed Spring Forecasting Experiments (81 dates). Using a novel ensemble storm track identification method, we extracted three sets of predictors from the WoFS forecasts: intra-storm state variables, near-storm environment variables, and morphological attributes of the ensemble storm tracks. We then trained random forests, gradient-boosted trees, and logistic regression algorithms to predict which WoFS 30-min ensemble storm tracks will correspond to a tornado, severe hail, and/or severe wind report. For the simple method, we extracted the ensemble probability of 2-5 km updraft helicity (UH) exceeding a threshold (tuned per severe weather hazard) from each ensemble storm track. The three ML algorithms discriminated well for all three hazards and produced more reliable probabilities than the UH-based predictions. Overall, the results suggest that ML-based calibrations of dynamical ensemble output can improve short term, storm-scale severe weather probabilistic guidance