Search

Geoelectric fields at the Earth's surface caused by magnetic storms constitute a hazard to the operation of electric power grids and related infrastructure. The ability to estimate these geoelectric fields in close to real time and provide local predictions would better equip the industry to mitigate negative impacts on their operations. Here we report progress toward this goal: development of robust algorithms that convolve a magnetic storm time series with a frequency domain impedance for a realistic three-dimensional (3-D) Earth, to estimate the local, storm time geoelectric field. Both frequency domain and time domain approaches are presented and validated against storm time geoelectric field data measured in Japan. The methods are then compared in the context of a real-time application.

Tides are the most predictable of oceanographic phenomena, due both to the simplicity and predictability of the astronomical forcing and to the near linearity of the ocean's dynamical response. In the classical and simplest scenario, tidal prediction is based on harmonic analysis of past measurements at a fixed location. Limits to predictability arise because isolated astronomical spectral lines are broadened into "cusps" of incoherent energy, for example through interactions with non-tidal flows. Tidal prediction at locations without past measurements has historically been a major challenge, but, owing to near-global observations of modern satellite altimeters, the empirical harmonic approach now yields reasonably accurate predictions throughout most of the open ocean. Advances in numerical modeling and data assimilation allow these predictions to be refined (especially in shallow seas where observations remain insufficient to directly constrain tidal wave structure) and extended to include tidal currents. We review recent progress in the development of global and regional-scale tidal prediction capabilities, summarize accuracy of available charts, and briefly consider outstanding issues. Satellite altimetry has also helped unravel the global tidal energy budget and has clarified the role of internal tides as a sink of tidal energy. We summarize these results, and then turn to the challenging problem of predicting internal tides. To the extent that low modes remain coherent with the surface tide, elevations can be directly mapped. We review current efforts in this direction, which are already producing charts with predictive capability. However, internal tides (especially higher modes) interact much more strongly with lower-frequency ocean flows, so a significant fraction of this tidal signal is intermittent and incoherent. We close with a brief review of ongoing efforts to model global tides in combination with wind-forced ocean motions.

Sustained accurate measurements of earth rotation are one of the prime goals of Global Geodetic Observing System (GGOS). We here concentrate on the fortnightly (Mf) tidal component of earth-rotation data to obtain new results concerning anelasticity of the mantle at this period. The study comprises three parts: (1) a new determination of the Mf component of polar motion and length-of-day from a multi-decade time series of space-geodetic data; (2) the use of the polar-motion determination as one constraint in the development of a hydrodynamic ocean model of the Mf tide; and (3) the use of these results to place new constraints on mantle anelasticity. Our model of the Mf ocean tide assimilates more than fourteen years of altimeter data from the Topex/Poseidon and Jason-1 satellites. The polar motion data, plus tide-gauge data and independent altimeter data, give useful additional information, with only the polar motion putting constraints on tidal current velocities. The resulting ocean-tide model, plus the dominant elastic body tide, leaves a small residual in observed length-of-day caused by mantle anelasticity. The inferred effective tidal 0 of the anelastic body tide is 90 and is in line with a omega-alpha frequency dependence with alpha in the range 0.2--0.3.

Knowledge of global oceanic tides has markedly advanced over the last two decades, in no small part because of the near-global measurements provided by satellite altimeters, and especially the long and precise Topex/Poseidon time series e.g. [2]. Satellite altimetry in turn places very severe demands on the accuracy of tidal models. The reason is clear: tides are by far the largest contributor to the variance of sea-surface elevation, so any study of non-tidal ocean signals requires removal of this dominant tidal component. Efforts toward improving models for altimetric tide corrections have understandably focused on deep-water, open-ocean regions. These efforts have produced models thought to be generally accurate to about 2 cm rms. Corresponding tide predictions in shelf and near-coastal regions, however, are far less accurate. This paper discusses the status of our current abilities to provide near-global tidal predictions in shelf and near-coastal waters, highlights some of the difficulties that must be overcome, and attempts to divine a path toward some degree of progress. There are, of course, many groups worldwide who model tides over fairly localized shallow-water regions, and such work is extremely valuable for any altimeter study limited to those regions, but this paper considers the more global models necessary for the general user. There have indeed been efforts to patch local and global models together, but such work is difficult to maintain over many updates and can often encounter problems of proprietary or political nature. Such a path, however, might yet prove the most fruitful, and there are now new plans afoot to try again. As is well known, tides in shallow waters tend to be large, possibly nonlinear, and high wavenumber. The short spatial scales mean that current mapping capabilities with (multiple) nadir-oriented altimeters often yield inadequate coverage. This necessitates added reliance on numerical hydrodynamic models and data assimilation, which in turn necessitates very accurate bathymetry with high spatial resolution. Nonlinearity means that many additional compound tides and overtides must be accounted for in our predictions, which increases the degree of modeling effort and increases the amounts of data required to disentangle closely aliased tides.

Two papers on tidal dissipation were completed and published. The first of these extended our earlier work, which focused on the dominant M2 constituent, to include 7 additional constituents. In addition to confirming a total deep water dissipation total very close to 1 TW, this study demonstrated significant differences in the distribution of dissipation between diurnal and semi-diurnal constituents. The second paper involved an extensive modeling study of tides in the present day and the last glacial maximum. In this study we showed that accuracy of tidal solutions for the present day Ocean were significantly improved by including a parameterization of internal tide drag over rough topography in the deep ocean. It was also demonstrated that a complete self-consistent treatment of Ocean self attraction and tidal loading was required for accurate solutions.

New empirical estimates of the long-period fortnightly (Mf) tide obtained from TOPEX/Poseidon (T/P) altimeter data confirm significant basin-scale deviations from equilibrium. Elevations in the low-latitude Pacific have reduced amplitude and lag those in the Atlantic by 30 deg or more. These interbasin amplitude and phase variations are robust features that are reproduced by numerical solutions of the shallow-water equations, even for a constant-depth ocean with schematic interconnected rectangular basins. A simplified analytical model for cooscillating connected basins also reproduces the principal features observed in the empirical solutions. This simple model is largely kinematic. Zonally averaged elevations within a simple closed basin would be nearly in equilibrium with the gravitational potential, except for a constant offset required to conserve mass. With connected basins these offsets are mostly eliminated by interbasin mass flux. Because of rotation, this flux occurs mostly in a narrow boundary layer across the mouth and at the western edge of each basin, and geostrophic balance in this zone supports small residual offsets (and phase shifts) between basins. The simple model predicts that this effect should decrease roughly linearly with frequency, a result that is confirmed by numerical modeling and empirical T/P estimates of the monthly (Mm) tidal constituent. This model also explains some aspects of the anomalous nonisostatic response of the ocean to atmospheric pressure forcing at periods of around 5 days.

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Maintenance of the large-scale thermohaline circulation has long presented an interesting problem. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Maintenance of the large-scale thermohaline circulation has long presented a problem to oceanographers. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.