The accuracy of state‐of‐the‐art global barotropic tide models is assessed using bottom pressure data, coastal tide gauges, satellite altimetry, various geodetic data on Antarctic ice shelves, and ...independent tracked satellite orbit perturbations. Tide models under review include empirical, purely hydrodynamic (“forward”), and assimilative dynamical, i.e., constrained by observations. Ten dominant tidal constituents in the diurnal, semidiurnal, and quarter‐diurnal bands are considered. Since the last major model comparison project in 1997, models have improved markedly, especially in shallow‐water regions and also in the deep ocean. The root‐sum‐square differences between tide observations and the best models for eight major constituents are approximately 0.9, 5.0, and 6.5 cm for pelagic, shelf, and coastal conditions, respectively. Large intermodel discrepancies occur in high latitudes, but testing in those regions is impeded by the paucity of high‐quality in situ tide records. Long‐wavelength components of models tested by analyzing satellite laser ranging measurements suggest that several models are comparably accurate for use in precise orbit determination, but analyses of GRACE intersatellite ranging data show that all models are still imperfect on basin and subbasin scales, especially near Antarctica. For the M2 constituent, errors in purely hydrodynamic models are now almost comparable to the 1980‐era Schwiderski empirical solution, indicating marked advancement in dynamical modeling. Assessing model accuracy using tidal currents remains problematic owing to uncertainties in in situ current meter estimates and the inability to isolate the barotropic mode. Velocity tests against both acoustic tomography and current meters do confirm that assimilative models perform better than purely hydrodynamic models.
Key Points
Tide model accuracy assessmentImproved accuraciesTidal current estimates
Global comparisons of barotropic and internal tides generated in an eddy‐resolving ocean circulation model are made with tidal estimates obtained from altimetric sea surface heights and an ...altimetry‐constrained tide model. As far as we know, our Hybrid Coordinate Ocean Model (HYCOM) simulations shown here and in an earlier paper are the only published high‐resolution global simulations to contain barotropic tides, internal tides, the general circulation, and mesoscale eddies concurrently. Comparing the model barotropic tide with a global data‐assimilative shallow water tide model shows that the global tidal elevation differences are approximately evenly split between discrepancies in tidal amplitude and phase. Both the model and observations show strong generation of internal tides at a limited number of “hot spot” regions with propagation of beams of energy for thousands of kilometers away from the sources. The model internal tidal amplitudes compare well with observations near these energetic tidal regions. Averaged over these regions, the model and observation internal tide amplitude estimates agree to approximately 15% for the four largest semidiurnal constituents and 23% for the four largest diurnal constituents. Away from the hot spots, the comparison between the model and altimetric amplitude is not as good due, in part, to two problems, errors in the model barotropic tides and overestimation of the altimetric tides in regions of strong mesoscale eddy activity. Examining the general energy distribution of the simulated internal tide is an important first step in the evaluation of internal tides in HYCOM.
Key Points
Model and observations show generation of internal tides in limited regions
Our global model is able to generate internal waves consistent with observations
Barotropic phase errors are a major source of errors in the model internal tide
•We tune several linear (modified) wave drag schemes in a tidal barotropic model.•We compare with TPXO elevation and dissipation rates, horizontally and with depth.•A scalarized tensor wave drag ...performs as well or better than the tensor equivalent.•Reducing the wave drag strength at supercritical slopes improves the solution.•The inclusion of artificial abyssal hill roughness yields a modest improvement.
A global tuning experiment for the semidiurnal tide is performed with a barotropic model. The model is forced with the M2 equilibrium tide and accounts for the self-attraction and loading (SAL) term. In addition to a quadratic drag, various linear internal wave drag terms adjusted by a scale factor of O(1) are applied. The drag terms include the original Nycander (2005) tensor scheme, the Nycander tensor scheme reduced at supercritical slopes, and their scalar sisters, a Nycander scalar scheme computed for additional abyssal hill roughness, and the Jayne and St. Laurent (2001) scalar scheme. The Nycander scheme does not have a tunable parameter, but to obtain the best tidal solutions, it is demonstrated that some tuning is unavoidable. It is shown that the scalar Nycander schemes yield slightly lower root-mean square (RMS) elevation errors vs. the data-assimilative TPXO tide model than the tensor schemes. Although the simulation with the optimally tuned original Nycander scalar yields dissipation rates close to TPXO, the RMS error is among the highest. The RMS error is lowered for the reduced schemes, which place relatively more dissipation in deeper water. The inclusion of abyssal hill roughness improves the regional agreement with TPXO dissipation rates, without changing the RMS errors. It is difficult to have each ocean basin optimally tuned with the application of a constant scale factor. The relatively high RMS error in the Atlantic Ocean is reduced with a spatially varying scale factor with a larger value in the Atlantic. Our best global mean RMS error of 4.4cm for areas deeper than 1000m and equatorward of 66° is among the lowest obtained in a forward barotropic tide model.
Abstract
Motivated by the potential of oceanic mesoscale eddies to drive intrinsic low-frequency variability, this paper examines geostrophic turbulence in the frequency–wavenumber domain. ...Frequency–wavenumber spectra, spectral fluxes, and spectral transfers are computed from an idealized two-layer quasigeostrophic (QG) turbulence model, a realistic high-resolution global ocean general circulation model, and gridded satellite altimeter products. In the idealized QG model, energy in low wavenumbers, arising from nonlinear interactions via the well-known inverse cascade, is associated with energy in low frequencies and vice versa, although not in a simple way. The range of frequencies that are highly energized and engaged in nonlinear transfer is much greater than the range of highly energized and engaged wavenumbers. Low-frequency, low-wavenumber energy is maintained primarily by nonlinearities in the QG model, with forcing and friction playing important but secondary roles. In the high-resolution ocean model, nonlinearities also generally drive kinetic energy to low frequencies as well as to low wavenumbers. Implications for the maintenance of low-frequency oceanic variability are discussed. The cascade of surface kinetic energy to low frequencies that predominates in idealized and realistic models is seen in some regions of the gridded altimeter product, but not in others. Exercises conducted with the general circulation model suggest that the spatial and temporal filtering inherent in the construction of gridded satellite altimeter maps may contribute to the discrepancies between the direction of the frequency cascade in models versus gridded altimeter maps seen in some regions. Of course, another potential reason for the discrepancy is missing physics in the models utilized here.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The dynamics of the flow field surrounding New Zealand are investigated using a series of global ocean models. The physical mechanisms governing the direction, magnitude, and location of the East ...Australian Current (EAC), the Tasman Front, the East Auckland Current (EAUC), and the East Cape Current (ECC) are studied using numerical simulations whose complexity is systematically increased.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Global maps of the mesoscale eddy available potential energy (EAPE) field at a depth of 500 m are created using potential density anomalies in a high‐resolution 1/12.5° global ocean model. Maps made ...from both a free‐running simulation and a data‐assimilative reanalysis of the HYbrid Coordinate Ocean Model (HYCOM) are compared with maps made by other researchers from density anomalies in Argo profiles. The HYCOM and Argo maps display similar features, especially in the dominance of western boundary currents. The reanalysis maps match the Argo maps more closely, demonstrating the added value of data assimilation. Global averages of the simulation, reanalysis, and Argo EAPE all agree to within about 10%. The model and Argo EAPE fields are compared to EAPE computed from temperature anomalies in a data set of “moored historical observations” (MHO) in conjunction with buoyancy frequencies computed from a global climatology. The MHO data set allows for an estimate of the EAPE in high‐frequency motions that is aliased into the Argo EAPE values. At MHO locations, 15–32% of the EAPE in the Argo estimates is due to aliased motions having periods of 10 days or less. Spatial averages of EAPE in HYCOM, Argo, and MHO data agree to within 50% at MHO locations, with both model estimates lying within error bars observations. Analysis of the EAPE field in an idealized model, in conjunction with published theory, suggests that much of the scatter seen in comparisons of different EAPE estimates is to be expected given the chaotic, unpredictable nature of mesoscale eddies.
Key Points
Global maps of the mesoscale eddy available potential energy are made from a HYCOM simulation and reanalysis
Modeled eddy available potential energy compares well to Argo observations globally, and to moored instruments locally
Model‐data comparisons of eddy available potential energy exhibit intrinsic scatter
The autocovariance of the semidiurnal internal tide (IT) is examined in a 32 d segment of a global run of the HYbrid Coordinate Ocean Model (HYCOM). This numerical simulation, with 41 vertical layers ...and 1/25∘ horizontal resolution, includes tidal and atmospheric forcing, allowing for the generation and propagation of ITs to take place within a realistic eddying general circulation. The HYCOM data are in turn compared with global observations of the IT around 1000 dbar, from Argo float park-phase data and mooring records. HYCOM is found to be globally biased low in terms of the IT variance and decay of the IT autocovariance over timescales shorter than 32 d. Except in the Southern Ocean, where limitations in the model cause the discrepancy with in situ measurements to grow poleward, the spatial correlation between the Argo and HYCOM tidal variance suggests that the generation of low-mode semidiurnal ITs is globally well captured by the model.
Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most ...critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally.
With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between ...spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.
Comb polyelectrolytes consisting of a polysiloxane backbone with tetraglyme and lithium-sulfonate terminated perfluoroether side chains have been synthesized. Both types of side chains were ...independently prepared with allyl functionality and attached to poly(methylhydrosiloxane) by hydrosilylation. Several methods are described to form lithium sulfonate from sulfonate esters. Spectroscopic and thermal characterizations of the resulting polyelectrolyte are described.