Resumen de Development of high-resolution l4 ocean wind products
Heat, moisture, gas, and momentum exchanges at the oceanic and atmospheric interface modulate, inter alia, the Earth's heat and carbon budgets, global circulation, and dynamical modes. Sea surface winds are fundamental to these exchanges and, as such, play a major role in the evolution and dynamics of the Earth’s climate. For ocean and atmospheric modeling purposes, and for their coupling, accurate sea-surface winds are therefore crucial to properly estimate these turbulent fluxes. Over the last decades, as numerical models became more sophisticated, the requirements for higher temporal and spatial resolution ocean forcing products grew. Sea surface winds from numerical weather prediction (NWP) models provide a convenient temporal and spatial coverage to force ocean models, and for that they are extensively used, e.g., the European Centre for Medium-range Weather Forecasts (ECMWF) latest reanalysis, ERA5, with ubiquitous hourly estimates of sea-surface wind available globally on a 30-km spatial grid. However, local systematic errors have been reported in global NWP fields using collocated scatterometer observations as reference. These rather persistent errors are associated with physical processes that are absent or misrepresented by the NWP models, e.g., strong current effects like the Western Boundary Current Systems (highly stationary), wind effects associated with the oceanic mesoscale (sea surface temperature gradients), coastal effects (land see breezes, katabatic winds), Planetary Boundary Layer parameterization errors, and large-scale circulation effects, such as those associated with moist convection areas. In contrast, the ocean surface vector wind or wind stress derived from scatterometers, although intrinsically limited by temporal and spatial sampling, exhibits considerable spatial detail and accuracy. The latter has an effective resolution of 25 km while that of NWP models is of ~150 km. Consequently, the biases between the two mostly represent the physical processes unresolved by NWP models. In this thesis, a high-resolution ocean surface wind forcing, the so-called ERA$^*$, that combines the strengths of both the scatterometer observations and of the atmospheric model wind fields is created using a scatterometer-based local NWP wind vector model bias correction. ERA$^*$ stress equivalent wind (U10S) is generated by means of a geolocated scatterometer-based correction applied separately to two different ECMWF reanalyses, the nowadays obsolete ERA-interim (ERAi) and the most recent ERA5. Several ERA$^*$ configurations using complementary scatterometer data accumulated over different temporal windows (TW) are generated and verified against independent wind sources (scatterometer and moored buoys), through statistical and spectral analysis of spatial structures. The newly developed method successfully corrects for local wind vector biases in the reanalysis output, particularly in open ocean regions, by introducing the oceanic mesoscales captured by the scatterometers into the ERAi/ERA5 NWP reanalyses. However, the effectiveness of the method is intrinsically dependent on regional scatterometer sampling, wind variability and local bias persistence. The optimal ERA$^*$ uses multiple complementary scatterometers and a 3-day TW. Bias patterns are the same for ERAi and ERA5 SC to the reanalyses, though the latter shows smaller bias amplitudes and hence smaller error variance reduction differences in verification (up to 11$\%$ globally). However, because of ERA5 being more accurate than ERAi, ERA$^*$ derived from ERA5 turns out to be the highest quality product. ERA$^*$ ocean forcing does not enhance the sensitivity in global circulation models to highly localized transient events, however it improves large-scale ocean simulations, where large-scale corrections are relevant. Besides ocean forcing studies, the developed methodology can be further applied to improve scatterometer wind data assimilation by accounting for the persistent model biases. In addition, since the biases can be associated with misrepresented processes and parmeterizations, empirical predictors of these biases can be developed for use in forecasting and to improve the dynamical closure and parameterizations in coupled ocean-atmosphere models.
