Self-and Inter-Crossover Points of Jasons’ Missions as New Essential Add-on of Satellite Altimetry in the Sub-Arctic Seas and the Southern Ocean

: For decades, satellite altimetry is providing reliable data on the sea level, surface currents, tides, planetary and wind waves and sea ice. We propose a potential enhancement of collecting special data subsets for a better temporal and spatial sampling of dynamics of the Sub-Polar seas. We show that peculiarities of the orbits of the TOPEX/Poseidon and Jasons’ missions generate a regular grid of crossovers which locations are ﬁxed in space and time of the orbit cycles. Two types of these special crossovers are speciﬁed: self-crossovers (SC) occurring as intersections of tracks of the same satellite and inter-crossovers (IC) appearing during the phases of interleaved orbits when tracks of two different satellites are crossing. The time delay between two consecutive measurements in the high-latitude sites is less than 2 h for SC and 1 h for IC, thus providing a ground for monitoring large-scale currents (including current direction) and wind waves. The sets of the special crossovers can be regarded as networks of virtual ocean buoys. The Norwegian Sea case study illustrates the value of the special crossovers for the Sub-Polar seas and the feasibility of building the networks of SC and IC in other areas of the World Ocean.


Introduction
Satellite altimetry has successfully been developing for ocean studies during the past three decades.It shows numerous applications for the sea level, ocean dynamics, coastal oceanography, planetary waves, ocean tides, wind and wave, ice cover, Earth's gravity field and climatology research [1][2][3][4][5][6][7][8][9].The altimetry "adequate sampling . . .have revolutionized oceanography" by providing data previously impossible to obtain, in the words of the outstanding oceanographer Walter Munk [10].The TOPEX/Poseidon, "the most successful ocean experiment of all time" [10], together with its follow-on Jason's missions is operating continuously at the same reference orbits for almost thirty years.
The crossovers, that is, points of intersection of ground tracks of altimeters, play an exclusive role in maintenance of functionality and performance of satellite altimetry.Calibration and validation of measurements, assessment of the measurement corrections, retracking procedure and reconstruction of spatial and temporal fields of the measured physical quantities inevitably refer to the crossovers [3].For the primary task of satellite altimetry, measurements of the ocean surface, the acceptable delay between two measurements is assumed to be less than 10 days [11].For wind-wave studies the corresponding spatial mismatch of measurements 50 km and time delay 1 h is recommended (e.g., http://globwave.ifremer.fr/products/globwave-satellite-data/satellite-crossovers,last access on 11 February 2021).The latter implies severe limitations on the number of crossovers and orientation of the altimeter orbits.Generally, such kind of crossovers is possible to find only for different altimetry missions that dictates new difficulties to overcome.We mention only two issues among many others:

•
The crossovers of different missions require their careful intercalibration accounting for features of hardware (different types of radars) and data processing methods (including different types of applied corrections).• Incommensurability of their orbit periods excludes repetitive monitoring of the particular sites because the corresponding crossovers appear to be unique, i.e., are irregularly distributed in space and time.
TOPEX/Poseidon-Jasons' missions provide special crossover datasets that are free of the above-mentioned difficulties and generate an adequate spatial and temporal sampling of the ocean processes at high latitudes [12].The distance between ground tracks of these satellites varies from more than 300 km at the Equator to just several kilometers in high latitudes.Time intervals between two consecutive measurements at the same point (selfcrossover-self-intersection of the altimeter track) can also be rather short, close to one orbit period (less than two hours).The orbit repetition period of about 10 days provides a series of such cross-over points with time lag of the same period.The locations of these points are fixed, which makes them an essential component of calibration/validation procedures of the altimetry.The orbit inclination of 66.04 • of TOPEX/Poseidon-Jasons' missions implies the coverage of the major part of the ice-free area of the World Ocean by measurements of these satellites.
Based on the above listed features, we present a special subset of data of High Latitude Crossovers (HLC) with minimal time intervals between consecutive measurements and identify two of their types: • self-crossovers (SC) when the crossover point is provided by measurements of the same satellite; and • inter-crossovers (IC) when measurements in a point are made by different altimeters, e.g., by a pair TOPEX/Poseidon/Jason-1, Jason-1/Jason-2 or Jason-2/Jason-3 during phases of interleaved orbits.
Because of these, individual (for SC) and synchronized (for IC) orbits of Jason's altimeters generate regular in space and periodic in time measurements in a number of special fixed sites.The SC and IC altimeter data with regular space and time sampling can be regarded as virtual stationary sea stations, analog of continuous buoy measurements of the sea level, wind speed and wave height.There are evident analogs in satellite altimetry of rivers when virtual water level gauges are established at locations where satellite ground tracks cross riverbed or river tributary.It makes the topic of the new crossovers quite promising, especially for the Sub-Arctic seas.
The above arguments constitute an idea and the essence of this study: special subsets of altimetry data with improved sampling can provide additional information on ocean dynamics and new opportunities for monitoring the seas.
In Section 2, we describe useful features of the CS and IC of the TOPEX/Poseidon-Jason's missions.
Section 3 describes geographical distributions of SC and IC in the World Ocean and possible gains of using the data of these special crossovers, especially in the Sub-Arctic seas.
Section 4 presents a case study where advantages and problems of the data subset to be solved are discussed.The paper is finalized by Discussion and Conclusions.

Self-Crossovers: Standard Jason'S Orbits
Figure 1a shows standard tracks of the Jason's missions in the North Atlantic.The ascending track N165 is given in red as a reference one.The self-crossover (SC) of the orbit occurs at descending track N168 (plus 3 to the ascending track number) near the limiting latitude φ t ≈ 66.05 • in a time lag a bit longer that one orbit period (≈112 min).In one more orbit period (plus 5 to the reference track N165), a similar SC occurs at lower latitude.The new crossovers with the fixed ascending track (N165 in our case) continues to appear with every new orbit, thus travelling to the South Pole.The passage of crossovers between poles takes approximately 1 day (about 12 orbit periods) for the fixed reference track (track N165 in our case of Figure 1).
The selected area of the North Atlantic in Figure 1a contains four crossovers of this type (descending tracks are shown in green).Time delay of SC is growing from north to south for the fixed ascending track and from south to north for a fixed descending track.Variation of the ocean conditions between two consecutive measurements in SC, evidently, depends on time delay and on characteristic scales of variability in the particular site.Time delay is minimal in high latitudes where variability of ocean environment can be rather strong.On the other hand, low variability in tropical crossovers can essentially reduce the effect of relatively large time lag of about half-day.In this way, quality of a particular crossover is not easy to assess in advance.This important issue is detailed in the next section.Here, one should emphasize that SCs occur regularly at the same points and time of the satellite cycle (i.e., for fixed pairs of track numbers), thus comprising an abundant data subset.

Inter-Crossovers: The Jason's Interleaved Orbits
At the beginning of every new mission (TOPEX/Poseidon-Jason-1; Jason-1-Jason-2; Jason-2-Jason-3) two satellites conform a special configuration of the so-called interleaved orbits when "old" satellite is positioned midway between standard ground tracks and two satellites are placed on opposite sides of the Earth.These phases have been realized during mid-September 2002-October 2005 for the tandem TOPEX/Poseidon-Jason-1, mid-February 2009-March 2012 for Jason-1/2 and mid October 2016-April 2017 for Jason-2/3.Such a configuration of two orbits provides a new family of crossovers we call intercrossovers (IC).For IC, the shortest time lag between consecutive measurements is about one-half of the orbit period (56 min).It occurs near the latitudes approximately 65.95 • in the Northern and Southern Hemispheres, i.e., a bit closer to the limiting latitudes of the orbits than SCs.As a result, similar ICs with longer time delay (approximately 1 1  2 , 2 1 2 . . .orbit periods) conform a configuration resembling the SC one shown in Figure 1b for the same reference track of Jason-3 N165 and descending tracks of Jason-2 N40, 42, 44, 46.Table 1 presents meridional dependencies of variability of parameters measured by altimeter for the ICs with the fixed ascending track N165 of Jason-3 and Jason-2 tracks N40-62.The corresponding ICs latitudes vary from the most northern to the most southern limit during the one-day period.The estimates are given for all the duration of the phase of interleaved orbits (October 2016-April 2017, totally, 19 cycles).Similar estimates can be made for SCs in Figure 1a.
The most northern IC has a minimal time lag a bit less than 1 h.The delay of the most southern one N62 in Table 1 is about one day (23 h 28 min).Variability of significant wave height H s and normalized radar cross-section σ 0 directly measured by altimeters grows in moderate latitudes with the time lag.The standard deviation (STD) of difference of measurements of Jason-2 and Jason-3 in IC with time lag about 9 h (Jason-2 track 48, Line 5) is six times higher than one of minimal time lag for wave height H s (Column 5), more than three times for normalized radar cross-section σ 0 (Column 7) and 3.7 times for altimetry-derived near-surface wind speed U alt (Column 8).The variability is even stronger for extreme values of wave height H s (Column 6) exceeding one order of magnitude.The variability weakens in tropical and subtropical ICs (Lines 6 and 7) but remains higher than in the Sub-Arctic IC with the minimal time delay.An additional aspect of this variability is its effect on numerous corrections of altimetry measurements.Some of these corrections are directly linked with the time delay, e.g., tidal corrections.Other corrections, such as Sea State Bias (SSB), are estimated parametrically [13][14][15] as functions of the directly measured values (e.g., H s and σ 0 ).One can easily show that all these corrections are minimal for the HLCs with minimal time lags that nudges us to draw special attention to these crossovers located near the Polar Circles.These special crossovers provide an additional information on anomalies of ocean currents and characteristics of wind-driven waves and near-surface winds (altimetric wind speed, e.g., [16]).The angle between the tracks in the SCs and ICs is quite large (about 30 • and 15 • correspondingly) that can be used for calculation of quasi-instantaneous values of: • anomalies of geostrophic current velocities and their directions [17][18][19]; and • wave heights, periods and directions of wind-driven waves [15,20,21].
It should be stressed that an instantaneous assessment of the geostrophic currents governed by balance of the Coriolis force and the sea surface height anomaly (SSHA) gradient is, generally, unattainable in the classic altimetry.Usually, estimates of noncollinear components of the geostrophic velocity vector are widely spaced over time.Interpolation of altimetry data in space and time is used to construct a generalized field of the ocean circulation.It leads to undersampling of relatively fast and small-scale ocean dynamics, especially in the Sub-Polar seas where the Rossby radius of deformation is small.The instantaneous along-track assessment of geostrophic currents with SCs and ICs near the limiting latitudes becomes a worthy solution of the problem.

Crossover Clusters: Temporal Sampling of SCs and ICs
As mentioned above, inter-mission crossovers are irregularly distributed in space and time due to different periods of satellites.The standard orbits of the TOPEX/Poseidon and Jasons' missions gives an essential advantage providing fixed periodically spaced sites where crossovers occur each cycle of approximately 10 days.Different types of these crossovers can be grouped into clusters, as shown in Figure 2, for the phase of interleaved orbits of Jason-2/3 We introduce codes for particular crossovers: "22" and "33" mean SCs of Jason-2 and Jason-3, correspondingly.The IC "23" means that Jason-2 passes the crossover before Jason-3, and "IC 32" is used for the opposite case.Each cluster comprises two SCs (circles and squares) and two ICs (triangles) centered relative to the ICs latitude.Table 2 presents key characteristics of the cluster of crossovers with minimal time lags we are mostly interested in.The cluster's SCs are located at approximately the same latitude about 65.44 • N(S).IC 32 is approximately 30 nautical miles (nmi) to the north (south) close to the latitude 65.94 • and its counterpart IC 23 is 3.6 nmi more apart to the north (south) at approximately the same longitude.Time delays for SCs 22 and 33 are evidently very close to each other because of equivalence of individual orbit periods of satellites (Column 4 of Table 2).The difference of about 9 min of time lags of ICs is explained by shift of Jason-2 and 3 orbits slightly different than one-half of the orbit period.Taking the track number of Jason-3 of IC 32 as a reference, one can get the track numbers of all elements of the cluster following rules of two last columns of Table 2 (track number more than 254 means a shift on one cycle).
Figure 2 shows location of three clusters (Figure 2a) together with timeline of crossovers in Figure 2b.Codes of crossovers correspond to colors of the clusters.One can see that the group of three clusters provide at least one measurement per day for distance of approximately 212 nmi along latitude.As shown below, many crossovers with minimal time lag are located in ocean straits.Thus, the issue of clustering of SCs and ICs can be useful for assessment water mass transport across these areas.
All the above mentioned predetermines our special attention to the high latitude SCs and ICs with minimal time lags, their geographical distribution and features of dynamics.

Geography of High Latitude Crossovers: Issues of Monitoring Sub-Polar Seas
In this section, we discuss the geographic location of SC and IC points with minimal time lags in both hemispheres, which fortunately cross a series of crucially important water bodies (Figure 3): 1.
The Norwegian Sea 2.
The Denmark Strait 3.
The Davis Strait 4.
Foxe Basin north of Hudson Bay in the Canadian Arctic Archipelago 5.
The White Sea 6.
The Southern Ocean Sea ice conditions in the White Sea, Foxe Basin, the Bering Strait and along Antarctica can significantly affect quantity and quality of the altimetry data.Below, we show all these areas in their connection with specific tasks for their studies and monitoring.

The Norwegian Sea
In the Norwegian Sea, the chain of 18 SCs and 18 ICs forms nine full clusters in the southern part of the sea between Iceland and Norway (Figure 3a).This is an important area where calculations of water mass and heat transfer from the Atlantic Ocean to the Arctic Ocean are made to understand climate change in the Arctic Ocean.Large-scale circulation in the Norwegian Sea and adjacent areas largely determines the climate in the North Atlantic Ocean [22][23][24][25][26][27][28].
The hydrological regime of the Norwegian Sea is determined by the flux of the Atlantic water, transferred by the Norwegian Current (a branch of the North Atlantic Current) through the Faeroe-Shetland Strait.In addition, a small part of the Atlantic water penetrates between Iceland and the Shetland Islands.Besides, freshened water of the Baltic Current comes from the North Sea along the coast of Norway.The total volume of the Atlantic water coming to the Norwegian Sea amounts to about 100,000 km 3 /year, and the water carried out to the Atlantic Ocean is about 40,000 km 3 /year.The maximum inflow of water in the Norwegian Sea is observed in December and January, the minimum being in April and May.The magnitude of interannual fluctuations of the Atlantic water inflow reaches 40,000 km 3 /year [29,30].The chain of new crossovers crosses the following currents which determines water exchange between the North Atlantic Ocean and the Norwegian Sea (from west to east): (i) the Coastal Iceland branch of the Irminger Current; (ii) the East Icelandic Current, which flows southward around eastern coast of Iceland; (iii) the Norwegian Current; and (iv) the Norwegian Coastal Current (to the south of 65 • N it is called the Baltic Current), which flows northward into the Norwegian Sea [29,30].
A possibility of along-track estimates of geostrophic current magnitudes and directions at the high latitude SCs and ICs from instantaneous SSHA provided by altimetry represents a strong advantage in establishing permanent monitoring of fields of currents and calculation of precise water exchange anomalies between sub-basins in this region.Nine clusters provide at least three independent estimates of the current field across the cross-section per day (see Section 2.3 and Figure 2).Thus, currents and water exchange can be monitored with essentially better time resolution that is usually made on the basis of gridded altimetry products for SSHA.
The region is also characterized by intense wind wave variability, especially in the conditions of severe weather which is a characteristic feature of the region.In winter season, wave heights can exceed 20 m, while the duration of "weather windows" (periods with waves less than 1.5 m) is never more than three days [31].Wave models reproduce the wave and wind fields in the Norwegian Sea with high accuracy, showing a good correlation with buoy measurements [32,33].Nevertheless, permanent wave monitoring at the HLCs can further support the analysis of the sea state under ongoing climate change.Being concentrated near the Polar Circle away from stationary buoys, these points can also be used for the additional calibration and validation of wave models.

The Denmark Strait
In the Denmark Strait, the chain of four full clusters (eight SCs and eight ICs) and five more SCs in the near-shore zone is located between Greenland and Iceland, just southward of the southern border of the Greenland Sea (Figure 3b).
The strait stretches 480 km long and 290 km wide, and it has a sill of 191 m depth at the Greenland-Iceland Rise.Being relatively narrow, it plays an exceptional role in the ocean circulation and water transport in the North Atlantic Ocean.In the Denmark Strait, the East Greenland Current flowing southward meets the northern branch of the Irminger Current (also known as the North Icelandic Irminger Current) flowing around Iceland northward and forms the East Greenland Frontal Zone [29,30].That is why most hydrographic surveys and theoretical and model studies in this region are focused on the water transport between the Greenland Sea (Iceland Sea) and the North Atlantic Ocean [34][35][36].Repeated transects across the Denmark Strait based on local SCs and ICs allows one to derive new high accuracy time series of sea level, current velocity, water transport, significant wave height (SWH) and wind speed measurements and compare them with standard oceanographic (CTD, ADCP) and wave data (visual observations, model simulations and buoy measurements).
The Greenland ice sheet essentially conditions the features of storm activity in this area, in particular the phenomena of bora-types in the western part of the section called a "Piteraq".This is a very strong katabatic wind which rushes down the east coast of Greenland usually in the autumn and winter.Wind speeds typically reach 50-80 m/s (180-288 km/h).On 6 February 1970, the Tasiilaq village located southward of the Denmark Strait was hit by the worst Piteraq ever documented in Greenland with estimated wind speed of 90 m/s (325 km/h) causing severe damage.Since the beginning of 1970, special Piteraq warnings are issued by the Danish Meteorological Institute [37].Measurements of wind speed at the HLCs in the Denmark Strait can contribute to the understanding the characteristics and impact of the Piteraq on water area and estimate the distance of its impact from the coast.

The Davis Strait
The chain of three full clusters (six CSs and six ICs) and three more nearshore SCs is located in the narrowest part of the Davis Strait, between the Baffin Island and western coast of Greenland, right in the middle of the strait which stretches from 70 • to 60 • N where it meets the Labrador Sea (Figure 3c).The width of the strait is between 300 and 1000 km, and its depth varies from 350 to 500 m at the sill in the narrowest part along the Arctic Circle.The Davis Strait separates the Baffin Bay and the Labrador Sea.This area is poorly understood due to a lack of in-situ measurements.
Two currents pass through the strait in different directions.The relatively warm West Greenland Current carries water along the coast of the West Greenland to the north, and the cold Labrador Current is along the eastern coast of the Baffin Land to the south, into the Labrador Sea and further into the Atlantic Ocean.The cold current carries numerous icebergs breaking away from the Greenland Ice Sheet in Melville Bay.Thus, near the eastern coast of the strait, the conditions for navigation are noticeably better; the navigation season in this part of the strait lasts from mid-summer to late autumn, and it was there that fishing and whaling vessels traditionally went north.At this time of the year, the main ports of Greenland are in operation: Paamiut, Nuuk and Sisimiut.In the winter months, more than 50% of the water area is covered with ice, while, in the warm season, more than 90% of the water area is ice-free [38].More than 90% of the water area of the Davis Strait is used for navigation: from 50 to 70 large-tonnage vessels are constantly in the waters of the strait, and the total number of vessels passing through the strait ranges in different years from 400 to 700.The main types of vessels are fishing, cargo and passenger, including coastal ferries off the Greenland coast, tourist liners and tankers [38,39].
Regular satellite altimetry tracks crossing the Davis Strait allow assessing anomalies of water exchange between the Baffin Bay and the Labrador Sea.During the ice-free period, there is the possibility to establish a permanent monitoring of water exchange through the Davis Strait and of the West Greenland and Labrador currents.These measurements can contribute significantly to the oceanography of this region under climate change.
The non-Jasons' tracks are almost parallel to the strait axis and, thus, their use for instantaneous assessment of water mass transfer across the strait is questionable.On the contrary, the Jasons' satellite tracks are close to zonal directions, i.e., they better capture geostrophic currents along the strait axis.The instantaneous currents can be successfully estimated from two intersecting tracks with a short time lag.A similar approach was tested in the Caspian Sea when calculating water exchange between sub-basins from satellite altimetry data [18,19].Jasons' HLCs with minimal time lags make such periodic in time along-track estimates at fixed sites more adequate to the problem of monitoring water exchange.
The Arctic warming shows new opportunities for intensification of shipping via the Northwest Passage, which starts from the Davis Strait [39].The present and future shipping activities (including those related to offshore oil and gas exploration) require better knowledge on wind and wave regime in this area which can be gained from the advantages of satellite altimetry including those related to measurements of wave characteristics at SCs and ICs.

Foxe Basin
The length of the Foxe Basin in the Canadian Arctic Archipelago is of 600 km with a width of 500 km.The chain of four full clusters (eight SCs and eight ICs), two ICs in the main water body and one SC in the eastern nearshore area covers the basin between the Melville Peninsula on the west and the Baffin Island on the east (Figure 3d), the Frosen Strait and the Roes Welcome Sound.The basin is shallow with depths of 100 m in the northern and central parts and up to 400 m in the southern part, which is connected with the Hudson Bay on the west and Hudson Strait on the east via Foxe Channel.Foxe Basin is rarely ice-free even in summer season, which complicates all types of oceanographic observations.Tidal currents and strong winds keep the ice pack in constant motion and contribute to the formation of numerous polynyas and shore leads, which are found throughout the region, especially in the southern part of the Foxe Basin.Regional warming allows establishing monitoring of wind-wave characteristics and water exchange between Foxe Basin and Hudson Bay during 2-3 months with the help of satellite altimetry and contributing significantly in the unknown oceanography of this area with a unique wildlife.

The White Sea
In the White Sea, two incomplete clusters (three SCs and four ICs) are located in the middle of the sea (Figure 3e).The crossover events occur on seven different days of each 10-day Jason's cycle, thus providing a reasonably good temporal coverage of the area.The White Sea has a very complicated shoreline and complex water dynamics conditioned by wind forcing, water exchange with the Barents Sea in its northeastern part, strong tidal motions and river inflows [30,[40][41][42][43].The cyclonic activity over the White Sea is quite intense forming specific wind-wave climate.In winter, the wind-wave regime is determined mostly by the influence of the Icelandic Low.In summer, it is formed under the influence of Arctic invasions of cold air.The storminess activity in the White Sea is relatively low due to its small size (about 250 km ×250 km) and winter ice cover.
In general, the White Sea is well studied region covered by visual observations, stationary buoys and weather stations measurements complemented by numerous national research programs, sea expeditions [41] and global and regional wind-wave models [44].Nevertheless, three SCs and four ICs in the White Sea can be used as an alternative source of valuable data on sea level variability, tides, currents, water exchange with the Barents Sea, wind and wave regime and ice cover [45,46].Satellite altimetry is already applied to the White Sea in order to study sea level variability, tidal, wind and ice regime [17,18,[47][48][49][50][51].

The Bering Strait
In the Pacific Ocean, three SCs and two ICs fall exactly into the Bering Strait (Figure 3f), which is of great value for the assessment of water exchange between the Pacific and Arctic oceans.The Bering Strait is 90 km wide and 60 m deep and has a cross-section of 3.4 km 2 .The water transport mainly directed into the Arctic Ocean varies between 0.5 and 2.5 Sv in relation to atmospheric condition (local winds and atmospheric pressure).The strait is usually completely frozen over from October to June [30,52].The first extensive survey of the Northern Bering Sea was probably carried on in the mid-1980s for five years in the scope of the NSF ISHTAR Program [30,53].The observations were concentrated in the summer months with the objective of determining the main physical, chemical and biological characteristics of the Northern Bering Sea in typical summer situations.The ISHTAR Program included, beside extensive field surveys, the application of a three-dimensional General Circulation Model derived and adapted to the Northern Bering Sea from the threedimensional model developed at the GeoHydrodynamics and Environmental Research Laboratory (GHER) of the University of Liege, Belgium [30,[54][55][56][57].The objectives of the model were to gain a general understanding of the Northern Bering Sea dynamics in the summer season, including modeling of the currents in the Bering Strait which were compared with measurements at buoys deployed in summer 1985 and 1986.Now, a regular monitoring of water exchange can be based on the evaluation of the vectors of geostrophic velocity for every Jason cycle (i.e., about 10 days) at five crossovers (Figure 3f).This will provide unique regular information on water exchange between the Bering Sea and Chukchi Sea with a step of 10 days for the past 30 years since the beginning of the TOPEX/Poseidon mission in September 1992.Wind and wave measurements are also of high value taking in mind intensification of shipping activities along the Northern Sea Route with a decline of the Arctic sea ice cover.

The Southern Ocean
In the Southern Hemisphere, crossovers form a belt ranging between 65.5 • and 66.0 • S near the Antarctic Circle (Figure 3g).The points do not cross the Polar Front (PF) and the Antarctic Circumpolar Current (ACC), but in some places they cross the Southern Boundary (SB) of the ACC and the Antarctic Slope Front (ASF) located over the upper continental slope of Antarctica (for example, close to the Drake Passage from the side of the Pacific Ocean) [30,[58][59][60][61][62].In the ice-free conditions, it would be possible to measure ocean currents, eddies, calculate energy of jets and eddies.Regular measurements of wave heights and periods can fill the gaps in the sea state data in the Southern Ocean caused by a lack of buoys and visual observations, and can also be used to validate wave models.The experiment is complicated by ice conditions, but the regular coverage of the area with crossover points (497 by about 16,000 km) will doubtless be useful for better understanding of the ongoing metocean processes in the coastal zone of Antarctica.

High Latitude Crossovers in the Norwegian Sea
In this section, we present an analysis of HLCs in the Norwegian Sea based on open access altimetry data of the portal AVISO (http://www.aviso.altimetry.fr/en/home.html,last access on 11 February 2021).Ocean data reanalysis of Copernicus Marine Service (https://resources.marine.copernicus.eu/,last access on 11 February 2021) was used for comparison of large-scale (geostrophic) currents and wind wave fields retrieved in the crossovers.

High Latitude Crossovers for Assessment Currents and Wind Waves
Assessment of large-scale current velocities is a key output of the satellite altimetry for ocean studies.Assuming a balance of the Coriolis force and pressure gradient due to the deviation of sea surface from equipotential surface of the Earth gravity field, one has 2ΩUn sin φ = −g∇ηs (1) where 2Ω = f is the Coriolis parameter and φ is the latitude of the site of measurements.The unit vector s is tangential to the velocity vector U in Equation ( 1) while its orthogonal counterpart n is tangential to the pressure gradient ∇η.Thus, U and ∇η are always orthogonal to each other within the model of geostrophic balance.
Figure 4a provides the well-known scheme of altimetry measurements of currents in a crossover point.Non-parallel tracks (red is ascending as in our previous figures and green is descending) allow the components of the geostrophic current, which are perpendicular to track directions (dashed lines in Figure 4), to be estimated.The linearity of Equation (1) eases estimate of both wind values and associated errors.Local differences in current velocity can be used to estimate changes in water transport.Anomalies of current velocity, derived from the Sea Level Anomaly, can be used for calculation of anomalies of water transport, or they can be transformed into absolute current velocity (and then to absolute water transport) by adding the Mean Dynamic Topography derived from the difference between the mean sea surface and the geoid or from ocean drifters and hydrographic profiles.
Analogously to the above scheme, the physical model of wave periods from altimetry data [20] is also based on estimates of spatial gradient.The appearance of gradient of H s is associated with two important features of the model.First, the fetch-limited scenario of wind-wave evolution is considered as the most adequate to the dynamics of wind-driven sea [63][64][65].This allows the variations of wave height to be related with spectral flux of energy, which is documented by instantaneous sea state parameters: wave height and wave period.Second, these variations of H s (∇H s ) are hypothesized to be aligned with the direction of wave group velocity (general direction of wave propagation).Thus, the balance equation can be reduced to the one for the spatial derivative of total energy E = H 2 s /16 and the energy flux due to nonlinear wave-wave interactions P [66] C peak g(s) Assuming self-similarity of wave evolution [67,68], one can express the energy flux P as function of instantaneous wave height and period.The resulting relationship among H s , wave period T p and its along-propagation derivative ∂H s /∂s can be used to construct a method of wave period estimates from altimetry data [20] (s is the unit vector and C peak g(s) is the group velocity of the spectral peak component along the direction of wave propagation).
The result of the model can be expressed in a quite compact relationship between wave steepness and the spatial gradient [20,21] where α ss ≈ 0.67 is dimensionless parameter of the theory of weak turbulence of winddriven waves [66,69].Wave steepness is defined in terms of significant wave height H s and spectral peak period One should emphasize that the model "forgets" the wave direction when operating with the gradient projection in the left-hand side of (2).Despite this difficulty, the modulo of "full gradient" can be reconstructed as a maximal possible sum of vectors with given projections onto altimeter tracks and the corresponding wave direction can be estimated with the "back-forth" ambiguity.This trouble of inherently nonlinear model (3) for wind waves is in heavy contrast with the scheme of linear conversion of scalar field of sea surface height into vector field of pressure gradient and, then, to the vector of geostrophic current.
Figure 4b shows gradients of H s for estimates of wave steepness and, then, wave period in the spirit of the similar scheme for geostrophic currents in Figure 4a.Two different configurations for alternative signatures of along-track derivatives are shown by solid and dash-dotted vectors, thus illustrating the problem of ambiguity of assessment of direction of wind waves.This ambiguity can be resolved with a priori information on wind and waves in the area under study.Figure 4 also shows a possible problem of HLCs when angle between tracks is rather small.Errors of measurements can affect significantly estimates of directions.In the case of wind wave measurements, high nonlinearity of the relationship (3) mitigate the problem when angle between track and wave direction is not too close to π/2.Generally, the single track measurements give rather good approximation [21].On the contrary, linearity of the geostrophic relationship (1) becomes a problem when estimating currents at low angles α in Figure 4a.
In the following section, we discuss features of reconstruction of slowly varying largescale currents and sea waves as it is documented by reanalysis of the Copernicus Marine Service (https://marine.copernicus.eu/,last access on 11 February 2021) and reconstructed from along-track estimates of the HLC.The Norwegian Sea is a good area for such a study.Here, HLCs form a regular array of nine full clusters and, thus, provide at least three daily measurements in the crossover points.

Geostrophic Currents in the Area between Iceland and Norway
The GLORYS12V1 product has been used for assessment of large-scale currents.The data are provided at a standard regular grid at 1/12 • (approximatively 8 km) resolution.This eddy resolving reanalysis is based on altimetry data since 1993.We choose the period 17-27 October 2016 covering the whole cycle of Jasons' altimeters.The absolute values of the current velocity are shown in Figure 5a,b for the first and last (27 October) days of the period under study.The figures reflect the complicated structure of jets in this area described in Section 3.1.At the same time, the fine structure of the currents that visually correlates with directions and locations of Jasons' tracks look somewhat suspicious.This distrust is confirmed by trivial spectral analysis.Discrete Fourier decomposition of a current velocity component G as function of longitude λ is defined as follows (i is the imaginary unit):  Here, Fourier amplitude A n is calculated for the corresponding spatial nth frequency of the discrete set of longitudes λ n (n = 0, 1, 2 . . .N, N = 128 in this study) Longitudes λ W , λ E are western and eastern boundaries of the transect in Figure 5. Figures 6b,d and 7b,d show spectral densities of the current components S n (Λ) = A 2 n as function of the spatial period 1/Λ.For the latitude φ = 66 • , slightly below the limiting one for the Jasons' missions (Figure 6b,d), the spectral peaks at the period ≈1.42 • matches the spacing of Jason's 254 tracks (360 • /254 ≈ 1.4173 • ).Only 30 nmi (≈55.56 km) to the north, the period of spectral peak becomes 50% longer (≈2 • ).This example clearly shows possible problems of the reanalysis based on the altimetry data near the latitude where an essential fraction of data (Jasons' data) disappears.
The problem of the current reanalysis is confirmed by comparison with along-track estimates of geostrophic currents for SCs and ICs for the period under study in Figure 8.Both meridional and zonal components show rather strong dispersion.At the same time, the along-track estimate of zonal component of geostrophic current shows remarkably low variations relative to a median value close to zero.On the contrary, the reanalysis predicts rather strong currents in the zonal direction.Both patterns, reanalysis and alongtrack estimates, do not contradict observations, but the dominance of meridional currents (north-south direction) looks more realistic , in the authors' opinion.
In fact, there is a problem of along-track estimates as long as tracks are close to parallels in high latitudes and, thus, capture meridional currents (cross-track variations of SSHA) better than the reanalysis that smoothes data more or less isotropically.

Wind Waves. Wave Period Estimates by Altimetry
The short time delay between consecutive measurements in HLC is an essential gain for wind wave studies.In some cases, it allows one to capture periods and directions of wave field that are, generally, unavailable from altimetry data using along-track records.
The reanalysis WAVERYS on a regular 1/5 • grid with a 3 h time step of the Copernicus Marine Service has been used for demonstration features of the HLC data subset for the same time interval 17-27 October 2016.The altimetry data are assimilated by MFWAM (Meteo-France WAve Model) taking into account oceanic currents from the GLORYS12 physical ocean reanalysis.
An essential result is good agreement of two consecutive measurements in the crossovers as illustrated by Figure 9a.Both ICs (triangles) with a delay of about 1 h and SCs (circles and squares) with 2 h delay show dispersion of about 30 cm that is close to the declared accuracy of wave height measurements by altimeters and ocean buoys [70].As expected, the wave heights measured by altimeters agree perfectly well with the reanalysis (not shown here).
Wave periods cannot be measured by altimeters directly but can be recovered within different approaches.The most used empirical parametrical models (e.g., [71][72][73]) postulate dependence of wave period on significant wave height H s and the normalized radar crosssection (NRCS) σ 0 .These models are able to provide the accuracy better than 1 s [73].At the same time, validity of these models can be questionable for high latitudes where buoy measurements are absent and, thus, the models are not properly substantiated.Additionally, these parameterizations are mission dependent: every mission requires its own set of empirical coefficients.
Alternative models of wave periods (e.g., [20,74]) rely upon physical links of wave parameters and, in this way, are less dependent on features of particular altimeter and an instantaneous sea state environment.Figure 9b shows quite good correspondence of the physical model (Equations ( 3) and ( 4)) and the WAVERYS reanalysis for the spectral moments (0, 2), wave period  High deviation of the model from the reanalysis quantities occurs at relatively short periods less than 4 s (wavelength less than approximately 25 m) where the model validity is questionable.In Figure 9b, both track records in the crossovers are used to retrieve the "full" period in accordance with the scheme of Figure 4b.This "full" period does not differ significantly from its approximate counterpart when the full gradient |∇H s | is replaced by the along-track derivative of one of altimeters.We present the ascending track calculations in Figure 9c with less than 15% of the period overestimate.
Results from the comparison between the physical model of Badulin [20] and its empirical counterpart of Gommenginger et al. [71] are reported in Figure 9d and show a remarkable similarity.The excess of the period of the model [71] is essentially higher than one fixed for a one-year period at ocean buoy sites (cf. Figure 2d in [21]).Such difference can be explained by characteristics of altimeters (Envisat in [21] and Jasons in this study) or, as an alternative, by the effect of high latitudes where the empirical model [71] has never been tested.
The apparent correlation of the two period estimates demonstrates, in our opinion, a good potential of the parametric models regardless their quantitative fail.The pronounced linear correlation of two models in Figure 9d signifies correspondence of the parametric representation to the inherent physical feature of self-similarity of evolution of winddriven waves.

Discussion
The HLCs with fixed location and minimal time lags between consecutive measurements can be seen as a network of virtual ocean buoys with enhanced temporal and spatial sampling.Building and exploitation of the network may provide a solid ground for international collaboration in studies of currents, water exchange, wind and wave characteristics in the Sub-Arctic seas and the Southern Ocean.
One can propose a list of particular tasks for the HLC network.The tackling of these tasks will shed light on important problems of ocean dynamics in Sub-Polar seas.
The first set of tasks mentioned above is associated with studies of the areas in Section 3 considering their special features.For greater oceanography purposes, this means getting reliable high-quality data on currents, winds and waves from the remote regions of Sub-Polar seas.
The permanent monitoring of large-scale currents with regular temporal sampling (see Section 4.2) becomes self-evident for the Norwegian Sea and Denmark Strait as vital gates connecting the North Atlantic and Arctic oceans and, thus, essentially determining the climate changes in the Arctic Ocean and the North Atlantic.
Measurements of currents, wind and waves in the Davis Strait and Foxe Basin will significantly enrich our knowledge of sea dynamics in these remote regions with severe weather and ice conditions.We see the HLC network as a decisive add-on and extension of rare sea campaigns and in situ observations in these areas.
The Arctic warming shows new opportunities for intensification of shipping via the Northwest Passage which starts from the Davis Strait.Satellite altimetry monitoring at SCs and ICs will bring valuable knowledge to the future studies, exploitation and nature conservation of these areas.
The White Sea is well supported by a set of meteostations and regular expeditions organized by the Russian Academy of Sciences, Ministry of Natural Resources and Environment of the Russian Federation and universities.This regional sea can be regarded as a testing ground for calibration and validation of satellite altimetry against in-situ measurements in the open sea and near-shore environment.
The Bering Strait plays a key role in water exchange between the Bering Sea and Arctic Ocean.Several SCs and ICs are located right in the Strait and, thus, are of vital importance primarily for regular water exchange measurements through the Strait.
In the Southern Ocean, we have about 400 HLCs which can be regarded as a huge virtual buoy network since the launch of TOPEX/Poseidon in September 1992.This worldaround zonal loop allows for tracing seasonal and interannual variability of currents, wind, waves and swells.It is difficult to imagine that such a unique experiment could ever be carried out by traditional in-situ tools.
The second block includes specific problems of altimetry methods.We do not mention these problems in the paper intentionally because of their technical complexity and numerous logistical issues.Here, we restrict the discussion to very few remarks.
Fixed locations of SCs and ICs and rather short time intervals between consecutive measurements provoke a question on control of altimetry measurements and building regional models of numerous altimetry corrections based, firstly, on data of HLCs.White Sea and Foxe Basin with their high tidal activity are the first candidates for developing such an approach.
Interpolation (and, quite often, extrapolation) of altimetry data can lead to suspicious results, as demonstrated in Section 4.2.The instantaneous estimates in HLCs can help to fix possible artifacts due to global data processing.
The final point of the discussion is the value of the TOPEX/Poseidon-Jason program for studies of high-latitude areas of the World Ocean as compared to other altimetry missions.HLCs generated by inclined orbits of the program cover vital sea areas.On the contrary, near-pole orbits of other altimeters (e.g., Envisat, SARAL/AltiKa, Sentinel, etc.) provide an abundant series of HLCs in ice-covered seas above 80 • latitudes and continental part of the Antarctica.Further support and development of the Jasons' missions, especially of their special phases of interleaved orbits, could be a good contribution to studies and exploitation of Sub-Polar seas.

Conclusions
In the present paper, we propose new essential add-ons of satellite altimetry related to the peculiarities of the orbits of the TOPEX/Poseidon and Jasons (Jason-1, Jason-2, Jason-3) satellite missions which have not been mentioned before in the scientific literature.The first one is related with the use of the self-crossover (SC) points that occur near the latitude of the Jason's orbits turn just below the Polar Circles where the satellite intersects its own orbit with one orbit period of 112 min.The second one is related to the inter-crossover points (IC) that occurs as intersection of tracks of interleaved orbits of two subsequent altimetry satellites with the time lag of about 1 h (approximately one half of the orbit period).
The features of self-and inter-crossovers relevant for improving the measurements of sea conditions are listed below: 1.
Measurements at these crossovers are carried out by the same (in the case of selfcrossovers) or similar (in the case of inter-crossovers) radar altimeters with the subsequent use of compatible processing and retracking procedures.

2.
Points of both types of crossovers are the same within the accuracy of the altimeter tracks.It permits organizing repetitive monitoring of currents, water exchange, wind and wave characteristics at these points with periods of full cycle for the SCs (about 10 days) and half-cycle (about five days) for ICs.

3.
The short time lag of high latitude crossovers (HLC) less than 2 h for SCs and about 1 h for ICs (see Table 2) reduces the effect of variability of sea environment (sea state) even for such short-scale phenomenon as wind-driven waves.

4.
Angles between tracks in the intersection points are large enough (about 30 • for SCs and 15 • for ICs) to get full vector of geostrophic current velocity and wave direction from along-track records.In the conventional satellite altimetry, estimates of current velocity are usually based on interpolated fields of sea surface height anomaly that are associated with essential smoothing of instantaneous fields in time and space.The latter can essentially corrupt the resulting current field.

5.
The number of these crossovers in the vital ice-free Sub-Polar areas of the World Ocean is large enough to conduct permanent monitoring of key oceanographic parameters for ocean dynamics, marine meteorology and climate research.6.
In many regions, the spatial distribution of these crossovers allows the assessment of anomalies of water exchange between the most important sub-basins in the Sub-Arctic seas.It is interesting to apply this methodology to the Southern Ocean where these crossovers do not cross the Antarctic Circumpolar Current (ACC), its branches and its fronts because the chain of these crossovers goes along the ACC and all the crossovers are located between the ACC and the coast of Antarctica.
The Norwegian Sea case study shows important advantages in the application of the new technology for ocean dynamics, sea state and climate studies which was unattainable for conventional satellite altimetry.Our preliminary results show a good correspondence of wave characteristics derived from HLC and reanalysis data.For the first time, zonal and meridional components of currents in the Norwegian Sea were estimated and found consistent with a classical scheme of currents on the transection between Iceland and Norway.The authors encourage the altimetry community to pay more attention to the specific subsets of altimetry data presented in this study as an important add-on to conventional approaches.
We have to note that the precise orbits of Jason satellites are fundamental for altimetry records.Jason satellites as well as other altimetry satellites are equipped with three main instruments for precise orbit determination: (i) the Global Navigation Satellite System (GNSS) or Global Positioning System (GPS) receivers; (ii) Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) receiver, and (iii) Satellite Laser Ranging (SLR) or Laser Retroreflector (LRR) retroreflectors.
The main progress in the Low Earth Orbit (LEO) modeling in terms of GNSS orbits is related to employing proper GNSS antenna calibrations and fixed phase GNSS ambiguities [75].The DORIS receivers now have more channels which allow improvements of the ground network in the DORIS-based Jason orbits [76].It was shown that the Terrestrial Reference Frame 2014 (TRF2014) SLR/DORIS complement impact on the ITRF2008 orbit is only 1-2 mm RMS radial difference during 1992-2009, and it increases after 2009, up to 5 mm RMS radial difference in 2016.Crossover data show the Deutsche Geodätisches Forschungsinstitut DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and those of the Jet Propulsion Laboratory JTRF2014 orbits for the Jason-1 period [76].
Arnold et al. [77] presented a parameter estimation approach to investigate systematic orbit errors and showed that SLR validation of the LEO satellites is not only able to detect radial but also along-track and cross-track offsets.Using the residual-based parameter estimation approach, corrections to the International Laser Ranging Service (ILRS) station coordinates, range biases and timing offsets were derived.As a result, RMS residuals of 5-10 mm were obtained over one-year data in 2016 using observations from a subset of high-performance stations and ambiguity-fixed orbits of four LEO missions [77].
Melachroinos et al. [78] investigated the effect of geocenter motion, which is a nonstationary process that includes secular trends, on Jason-2 orbits and the mean sea level.They found that, for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a −4.67 ± 3.40 mm error in the Z-component of the orbit frame, which creates the systematic error in the Mean Sea Level (MSL) of 1.06 ± 2.66 mm mainly due to the uneven distribution of the oceans between the Northern and Southern hemispheres [78].
The recent improvement in the precise orbit determination of Jason orbits includes the improvement of the orbit models, for instance, the improved macromodels such as box-wing models used to properly handle the non-gravitational forces emerging from the solar radiation pressure, albedo and atmospheric drag [79][80][81] The present study was based on Level 2 altimetry products where recent results of all the techniques, orbit correction and modeling are converted into physical parameters of Sea Surface Height (SSH), Sea Level Anomaly (SLA), Significant Wave Height (SWH), etc.This is why the discussion of the issues related to the modeling of the Jason orbits, including the LEO (Low Earth Orbit) modeling issues of Jason satellites, are beyond the scope of the paper.
Finally, we have to note that the more precise determination of the satellite orbit will not impact the new idea and results presented in the paper because we draw attention to new crossovers which have never been mentioned before whatever the accuracy of the Jason satellites is.The present spatial accuracy of the Jason satellites orbit of about ±1 km [70] is

Figure 1 .
Figure 1.Crossovers of Jasons' orbits in the North Atlantic.The ascending track N165 of Jason-3 is given in red, and green trajectories show its descending counterparts generating crossovers within the one-day delay: (a) self-crossovers (SCs) of Jason-3 standard orbits; and (b) inter-crossovers (ICs) of Jason-3 and Jason-2 during the phase of interleaved orbits October 2016-April 2017.

Figure 3 .
Figure 3. SCs and ICs of the phase of interleaved orbits of Jason-2/3: (a) the Norwegian Sea; (b) the Denmark Strait; (c) the Davis Strait; (d) the Foxe Basin; (e) the White Sea; (f) the Bering Strait; and (g) the Southern Ocean.

Figure 4 .
Figure 4. Setups of altimetry measurements in a crossover point.(a) Measurements of large-scale currents by satellite altimeter.Along-track differences of the sea surface height anomaly η 12 − η 11 and η 22 − η 21 are used to get transversal projection of the current velocity.(b) Measurements of wave steepness.Along-track differences of the significant wave height H s (H s12 − H s11 ) and (H s22 − H s21 ) are used to estimate along-track components of ∇H s with "back-forth" ambiguity.Two possible configurations are shown by solid and dash-dotted vectors.

Figure 5 .
Figure 5. Geostrophic currents in the Norwegian Sea from the GLORYS12V1 reanalysis.Absolute value of current velocity in m/s for: 17 October (a); and 27 October (b).White circles show crossover sites (for details, see Figure 3a).

Figures 6 and 7
Figures 6 and 7 present two cross-sections along latitudes φ = 66 • and φ = 66.5 • , i.e., only 30 nmi to the north derived from the GLORYS12V1 reanalysis.The zonal component U and the meridional one V demonstrate general feature of circulation in the area shown in Figure 5.The meridional component of the current V is directed to the North along the main eastern part of the transect between Iceland and Norway in full agreement with distribution of the Coastal Iceland branch of the Irminger Current and the East Icelandic Current, which flow southward around eastern coast of Iceland (see western side of the transect), and the Norwegian Current and the Norwegian Coastal Current, which flow northward into the Norwegian Sea.Several peaks directed northward represent known branches of the Norwegian Current and the Norwegian Coastal Current at coast of Norway.Discrete Fourier decomposition of a current velocity component G as function of longitude λ is defined as follows (i is the imaginary unit):

Figure 6 .
Figure 6.Transection along latitude 66 • N on 17 October 2016: (a,c) the current velocity components; and (b,d) spatial spectral densities of the current components vs. spatial period (degrees).

Figure 7 .
Figure 7. Same as in Figure 6 at transection along 66.5 • N.

Figure 8 .
Figure 8. Scatterplot of current velocity components of the reanalysis vs. along-track estimates of geostrophic currents at ISs and SCs: (a) zonal; and (b) meridional components.

Table 1 .
Variability of the altimetry measurements in ICs during the phase of interleaved orbits of Jason-2/3 (mid October 2016-April 2017).Ascending track N165 of Jason-3 is taken as a reference one, descending track numbers of Jason-2 are shown in the first column.Estimates for last two lines (tracks 60 and 62) are absent because of low data quality.

Table 2 .
Spacing of crossovers with minimal time lag.