Tides on Other Earths: Implications for Exoplanet and Palaeo‐Tidal Simulations

A key controller of a planet's rotational evolution, and hence habitability, is tidal dissipation, which on Earth is dominated by the ocean tides. Because exoplanet or deep‐time Earth topographies are unknown, a statistical ensemble is used to constrain possible tidal dissipation rates on an Earth‐like planet. A dedicated tidal model is used together with 120 random continental configurations to simulate Earth's semidiurnal lunar tide. The results show a possible ocean tidal dissipation range spanning 3 orders of magnitude, between 2.3 GWto 1.9 TW (1 TW=1012 W). When model resolution is considered, this compares well with theoretical limits derived for the energetics of Earth's present‐day deep ocean. Consequently, continents exert a fundamental control on tidal dissipation rates and we suggest that plate tectonics on a planet will induce a time‐varying dissipation analogous to Earth's. This will alter rotational periods over millions of years and further complicate the role of tides for planetary evolution.


Introduction
The ocean tides drive and influence a range of geophysical processes. The tidal energy dissipated in shallow waters controls shelf sea stratification (Simpson et al., 1990) and sustains a vertical nutrient flux vital for marine primary productivity in shallow and deep environments (Hickman et al., 2012;Tuerena et al., 2019;Williams et al., 2018). In the deep ocean, tidally driven mixing is integral to the global overturning circulation (Kuhlbrodt et al., 2007;Srokosz et al., 2012;Wunsch & Ferrari, 2004), whereas the tidal range during particular geologic periods may exert an influence on the evolution of complex life (Balbus, 2014). The tides are also a first-order control on orbital evolution for the Earth-Moon system (Bills & Ray, 1999;Munk, 1968). This is illustrated through the current, anomalously high, (3.8 cm yr −1 ) recession rate of the Moon which does not match its radiometric age of 4.5 Gyr (Barnes, 2017;Dickey et al., 1994). The erroneous assumption of a constant modern globally integrated dissipation rate suggests an Earth-Moon system age of about 1.5 Gyr (Darwin, 1899;MacDonald, 1964;Munk, 1968). The missing factor may be the time-varying effect of the continental configuration on the tide, a condition shown to have driven long-term dissipation variability on Earth Kagan, 1997). Consequently, the implications of an ocean tide controlled by a planet's continental configuration is an important parameter when modeling the orbital evolution of an ocean-bearing planet.

10.1029/2019GL085746
Key Points: • The possible tidal dissipation rates on Earth can span at least 3 orders of magnitude, with bathymetry this range narrows • Geologic features currently unique to Earth are crucial factors in the tidal energy budget • Time varying dissipation due to surface evolution such as tectonics should be considered in models of evolving planetary rotation To date, there are over 4,000 confirmed exoplanets, of which more than 630 are in multiplanet systems (see http://exoplanet.eu/catalog/ for the latest estimate). The masses of these planets range from gas giants of the order of ∼1,000 M ⊕ (Earth masses), to less massive objects that are often characterized as terrestrial planets, thought to be similar in composition and scale to the terrestrial planets in our solar system (exoplanets.eu). Concurrently, attempts have been made to estimate the "habitability" of these planets, a condition based on the probability of a given planet being able to support liquid water at its surface (Kasting et al., 1993)-one of many definitions for habitability. This is a function of host star irradiance, planet-star separation and the chemical composition of the planet and its climate system (Seager, 2013). These conditions determine the extent of the Habitable Zone, the theoretical shell around a given star at which liquid water could exist, given the appropriate climatic conditions (Kopparapu et al., 2013).
As the number of confirmed planets has grown, planets with potential oceans emerge (e.g., Tsiaras et al., 2019), and the number of planets that appear to be within their habitable zone is expected to increase . Estimates put the number of planets in the Milky Way as high as ∼10 10 (Dressing & Charbonneau, 2015). Since smaller, potentially terrestrial, planets far outnumber the larger gaseous planets (Cassan et al., 2012;Fressin et al., 2013;Howard, 2013), habitable planets may be common.
Little is known about these planets, which has led to the widespread use of numerical models of potential climatic conditions to determine a probability of a given body-hosting water (Seager, 2013), but there are large uncertainties in this approach. For example, varying host star irradiances and masses, and planetary atmospheric compositions, masses, rotation rates, eccentricities, and obliquities will drastically alter the conditions at the surface (Colose et al., 2019;Kasting et al., 1993;Way & Georgakarakos, 2017;Way et al., 2018Way et al., , 2016Williams & Pollard, 2002, 2003Yang et al., 2014). Few of these simulations, however, consider the effect of a water ocean tide, despite tidal friction being a key controller on orbital evolution, and hence habitability (Bills & Ray, 1999;Egbert et al., 2004;Green et al., 2019b;Lingam & Loeb, 2018).
If a planet's spin frequency is greater than its moon's orbital frequency, tidal friction will increase the semimajor axis length and rotational period of the planet, leading to a transfer of angular momentum due to the torques exerted on the body (Darwin, 1899). This process leads to an increase of the angular momentum of the tide raiser and a decrease in the planet's rotation rate, so the rotational evolution depends on the amount of tidal energy dissipation. This is in turn a function of the height of the tidal bulges and the rheology of the planet's interior (Renaud & Henning, 2017), including the fraction of water on it. On Earth today, the solid body dissipation for the dominating M 2 lunar tide is 0.08 TW (Ray et al., 1996), whereas the associated total ocean tidal dissipation rate is 2.4 TW (Egbert & Ray, 2001). Obviously, a liquid ocean can provide a more energetic tidal response than a solid body alone because it is easier to excite a tide in the ocean than in the solid Earth.
There are links between a planet's potential habitability and its rotation rate, for example, full spin-orbit synchronization at lower rates and the reduction of meridional atmospheric convection at higher rotation rates due to a stronger Coriolis force (Yang et al., 2014). To properly constrain a planet's climatology, and therefore habitability, a range of variables, including the rotation rate, topography and land/ocean mask, must be known (Colose et al., 2019;Way et al., 2016;Way & Del Genio, 2020;Yang et al., 2014). So a planet's past, present, and future total tidal dissipation rate must be quantified, to improve estimates of those other dependent properties.
Here, we constrain the potential range of tidal dissipation on an Earth-like planet using a dedicated numerical tidal model and a large sample of random continental configurations. This will allow us to produce bounds on the potential dissipation rates that can then be used for rotational evolution simulations and to provide error bounds on simulations of deep-time tides when Earth's surface looked very different from today's.

Methods
We use OTIS, the Oregon State University Tidal Inversion Software, a well established numerical tidal model that has been used to simulate deep-time past, present, and future tides on Earth (Egbert et al., 2004;Green et al., 2017Green et al., , 2018Wilmes et al., 2017) and on ancient Venus (Green et al., 2019b). The model has been bench-marked against other global nonassimilating tidal models and demonstrated to reproduce Earth's 10.1029/2019GL085746 present-day tide with a high degree of accuracy (Stammer et al., 2014). It solves the linearized shallow-water equations, where U is the depth-integrated volume transport (the current velocity u times the depth H), g is the gravitational constant, f the Coriolis vector, is the tidal amplitude. SAL is the elevation due to self-attraction and loading (SAL; here set to 8% of the amplitude following Ray, 1997) and EQ the equilibrium tidal elevation. The model is forced by the astronomical tide generated by the Lunar and Solar gravitational potential. F = F B + F w is the total loss of energy to bed friction (F B ) and tidal conversion (F w , describing the generation of an internal tide). The scalar product F · u thus gives the tidal dissipation rate, D. Bed friction is given by The drag coefficient, C d = 0.003, represents mean seabed roughness and is based on an appropriate value for present-day Earth (see, e.g., Taylor, 1920), and u is the combined velocity vector of all the tidal constituents. Note that the model is insensitive to the chosen value of C d , and sensitivity simulations (not shown) with it varying by a factor 3 did not significantly change the results. However, this may not be the case for aqua-planets (no continents), where the dissipation may scale directly with C d due to a lack of topography. Futhermore, for Earth-like ocean with fine sediments we do not expect C d to vary beyond this parameter range.
Tidal conversion is important in simulations with topography and is given by F w = C|U|, where C is the conversion coefficient given by (Zaron & Egbert, 2006): Here, = 50 is a scaling factor, N b is the buoyancy frequency at the seabed,N is the average buoyancy frequency in the vertical, and is the frequency of the constituent under investigation. Given that the surface properties of exoplanets are unknown, the buoyancy frequency is computed from a statistical fit based on observations from present-day Earth, or N(x, ) = 0.00524 exp(−z∕1,300), where z is the vertical coordinate, and the constants 0.00524 and 1,300 have units of s −1 and m, respectively. Tidal conversion will differ for other fluids (i.e., liquid Methane), but we assume a water ocean here.
All simulations were performed for three dominating tidal constituents, but we focus our discussion on those representing the principal semidiurnal lunar (M 2 ) and solar (S 2 ) tides. Their respective constituent periods for Earth today are 12.42, 12, and 23.93 hr. Note that OTIS handles tides only; there is no other forcing included in the model and although inertial oscillations are present we only discuss the ocean tides here.

Continental Configurations and Sensitivity Simulations
The model grids were generated using Planet Generator (https://topps.diku.dk/torbenm/maps.msp), a fractal map generator capable of simulating complex planet surface features, which is feasible because the shape and distribution of continents are approximately fractal-like in nature (Goodchild & Mark, 1987;Mandelbrot, 1982). The generator produces gray-scale images, where each pixel value is between 1 (white) and 0 (black). Here, the images were 360 × 180 pixels in size, and a flat-bottomed "bathtub" ocean grid, at a horizontal resolution of 1 • × 1 • , was built by allocating a depth of 4,500 m to all gray-scale pixels with a value below 0.2. The aim here is to evaluate the effect of the position and complexity of the continents on the tide, and a further series of sensitivity simulations on the boundary values was performed, with the coastline at gray-scale values of 0.1 and 0.4, as well as simulations with bottom topography (see below).
In total, 120 different continental configurations were evaluated, including an equilibrium tide simulation without continents (see table in the supporting information for a summary of all the simulations). As each configuration is unique, all configurations are independent of one another, so their general geometry and relative positions are heterogeneous. This represents (to some degree) the variation that would be expected between planets and between geological epochs on the same planet. To show the total effect of continent 10.1029/2019GL085746 size, ocean basin size, and coastal complexity, we use the nondimensional value R = L tot ∕ √ A ocean . L tot is the total coastline length (km) and A ocean (km 2 ) is the total ocean area. L tot was found by identifying the coordinates of the coastal boundary between individual continents and the ocean. The distance between neighboring coordinates were then calculated using a Great-circle technique: The coastal length for each continent was found, and the total global distance was the sum of these values. This is justified as there is little point considering distances smaller than the smallest horizontal resolution in the grid, although Earth's real tides are influenced by small, nonrandom topographic features.
We also completed two sets of simulations with ocean bathymetry for eight maps selected from the population (see below). The first has a set global average depth (SGAD), whereas the second set had a set maximum depth (SMD). Both were computed by multiplying the gray-scale values with a constant. It was equal to 4,500 for the SMD, so the deepest point will be 4,500 m and the shallowest 900 m. The SGAD constant was individually chosen so the mean depth of the ocean was 4,500 m. This gave a maximum depth between 9,500 and 42,000 m depending on the configuration. Depths of tens of kilometers, for a planet of Earth's radius, is approaching the limit of the shallow-water approximation employed in OTIS. The simulations are summarized in the supporting information Table S1; note that the extreme cases are not discussed further.
These eight maps, with bathymetry, were also used to test the sensitivity to stratification through simulations with enhanced ( = 500 in Equation 4) or reduced stratification (using = 5). The eight maps, but with bathtub bathymetry, were also subject to tests of the effect of the planet's rotation rate by adjusting Earth's day length to 3 hr and 8 days, respectively, so that the new M 2 periods are 1.51 and 124.20 hr. Only Earth's rotation period is changed, no adjustment is made to the Lunar period of 27.3 days. Note that a rotation period of 8 days still puts the planet within the fast-rotator regime of its climate dynamics (Way et al., 2018).

Present-Day Earth Sensitivity Simulations
Three simulations were run for present-day Earth. The control used present-day ocean bathymetry, the first of two sensitivity simulations had a bathtub ocean with the depth set to 4500 m everywhere (see Figures S1a and S1b), whereas the second had Earth's shelf seas removed by setting any ocean initially shallower than 1,000 to 1,800 m to highlight the importance of shelf dissipation. In the control, at the low resolution of 1 • × 1 • , OTIS overestimates the M 2 dissipation rate by a factor 2.5 compared to observed dissipation rates (Egbert & Ray, 2001; see our Figure S1), whereas S 2 is overestimated by a factor 3. The bathtub configuration is far less dissipative (Figures S1b and S1c) and underestimates the M 2 and S 2 dissipation rates by factors 25 and 200, respectively. Removing shelf seas gives 0.77 TW of M 2 dissipation, compared to the 6.3 TW of the control. These results make sense dynamically. In a coarse-resolution simulation not enough energy dissipates in the deep ocean due to underrepresented topography. Because the tidal energy is not lost in the deep ocean, the shelf sea currents are overestimated and, with frictional losses being proportional to the cube of the speed, the dissipation is overestimated. The same happens in the bathtub runs, but there is then no shelf to dissipate the energy on, and instead the global rate drops to deep abyssal values. The stratification sensitivity tests can be found in the supporting information.

Dissipation
As anticipated, the addition of continents of increasing complexity generates a more energetic tide, highlighted for the eight examples in Figure 1 and summarized in Figure 2; these are labeled A-H in the following and chosen because they cover the parameter space in Figure 2. The globally integrated dissipation in Figure  2a scales with the magnitude of continental configuration complexity, quantified here as a ratio, R, between the coastline length and the ocean basin area. As R increases, the M 2 tidal dissipation rate follows and soon reaches values 2 orders of magnitude larger than the rate in the equilibrium tide (i.e., the 2.3 GW at R = 0). The dissipation range for a certain value of R can span nearly an order of magnitude, but 36% of the configurations occur within the 0.1-1 TW range. The S 2 response is generally weak (not shown), with an amplitude of about 45% of M 2 , although some configurations produce an S 2 dissipation that is larger than M 2 because of S 2 resonances, for example, configuration E.
Both M 2 and S 2 reach dissipation maxima at continental distributions which are qualitatively among the most fragmentary, with several small-scale topographic features introducing local flow acceleration. For M 2 , the largest dissipation rate is 1.9 TW (Figure 1e), while S 2 peaks at 0.67 TW-values 19 and 260 times  The point at R = 0 is the control water world, the Earth controls are shown as a circle (bathtub) and triangle (with bathymetry), and the star is Earth without shelf seas. The color scale of continent area has been normalized against Earth's present continental area of 5.1 × 10 7 km 2 . Additional sensitivity simulations can be found in Table S1. Surfaces with some large number of small disconnected basins were not tested. (b) The same configurations, with a form of random bathymetry, Earth without shelf seas and the control are included for comparison. larger than in our bathtub Earth simulation and on par with the observed dissipation rates on Earth today (Egbert & Ray, 2001). As the coastline length continues to grow (i.e., R increases), the continent area must also increase. This limits the total ocean area and therefore the potential for large integrated tidal dissipation rates as there is a smaller ocean to dissipate energy in. In the extreme case, we have a planet with a large number of basins too small to host a tide. Thus, as R increases, we can expect dissipation to reduce, such as on a surface with many disconnected basins, smaller than usual in our ensemble. This form was not explicitly tested, but it would have quiescent tides under the conditions we use, potentially expanding this range to a greater number of orders of magnitude. We leave exploration of this limit to future work.
There is not a specific configuration, quantified by R, that is especially dissipative; two configurations can have similar R values, but be qualitatively dissimilar in the actual position of the landmasses (e.g., see configurations D and E in Figure 1). Given that global dissipation rates are dependent on oceanic area, a configuration that elevates dissipation closer to the equator should result in a larger global value than one with a continent of the same shape at higher latitudes. However, a number of the simulations possess a polar ocean,for example, H, or in some cases a single large hemispheric ocean with a large tide. The dissipation in these instances can have a large deviation from what could be expected, demonstrated in the increasing variability above very small R values.
On present-day Earth, deep ocean tidal dissipation is driven by conversion at bathymetric features (Egbert & Ray, 2001). The SMD configurations generally have an enhanced dissipation (see Figure 2b) influenced by the continental configuration. The two sets-bathtub and SMD-are significantly different (F 1,119 = 0.49, p < 0.05), with the bathtub dissipation values being poor predictors of the SMD values. The greatest increases are among the configurations that were the least energetic without bathymetry, while the bathtub maximum of 1.83 TW can be achieved with smaller and less complex continents. Because our results span 3 orders of magnitude, any upper limit on Earth's deep ocean dissipation remains highly uncertain; however, as our ensemble maxima are 1.83 and 1.92 TW, such a limit may lie close to these values at this resolution. While the variance is large (s 2 = 7.24 × 10 10 TW 2 ), adding bathymetry narrows the range of dissipation values compared to the runs without: 92% occur within 1 standard deviation (0.39 TW) of their mean (0.46 TW) with bathymetry, as opposed to 76% without bathymetry, although this still spans 2 orders of magnitude.

Sensitivity Simulations
In the SGAD set, the dissipation rises 3 orders of magnitude at the lowest R value (see Table S1). A consistent feature in these simulations is the convergence toward a common dissipation value as the stratification is strengthened by a factor of 10: The most energetic configurations are dampened, whereas those at the scatter's lower bound are enhanced. In the SMD set, shown in Figure 2b, the trend is similar; however, the dissipation shows less variability within the R scale, while the examples lose their dependence on R. The spatial distribution of the SMD runs are shown in Figure S4.
The planet's rotation rate exerts a fundamental control on the tide (Green et al., 2019b) by setting resonant periods. Increasing Earth's current rotation rate to 3 hr also means the number of semidiurnal periods in one present-day Earth day rises from 2 to 8, thereby potentially increasing the energy dissipation over a certain fixed period of time. This is indeed the case in the simulations herein, with dissipation enhanced in each example, setting an upper bound on the scatter. The large increases for maps G and F reflect the contribution of several smaller, resonant basins. With an 8-day day length, dissipation is reduced across all configurations. Crucially, the influence of continental configuration is still apparent, with the rotation adjustment conforming to the upper and lower bounds of the scatter seen in Figure 2a. Our three frequency sensitivity simulations show, unsurprisingly, that dissipation changes when the frequency is changed due to change in the resonant properties of the basins. However, there is not necessarily a continuous trend in the results, nor can we state if we have hit maximum global resonance in terms of rotation. For very fast or slow spin rates resonances become more localized and will not contribute to globally integrated dissipation rates to the same extent as at more moderate rates-see Green et al. (2019a) for a discussion.

An Influence on the Rotational Evolution
The slowing of the planets rotation rate is tied to the magnitude of tidal dissipation. In Table 1 we give an estimate of the variability introduced in the day length, LOD, as a function of just the M 2 frequency after 10 6 years. Eccentricity (e) and obliquity are 0, and we make the assumption that de/dt = 0. From conservation of energy and momentum, for a small moon orbiting a rapidly rotating planet such as the Earth-Moon system, the globally integrated dissipation, D, is related to the evolution of the Earth-Moon separation, a, as where M M is the Moon's mass (7.347 × 10 22 kg), is the Earth's angular rotation rate, a is the Earth-Moon separation, n = √ G(M E + M M )∕a 3 is the mean motion of the Moon (G = 6.674 × 10 −11 m 3 kg −1 s −2 is the gravitational constant and M E = 5.972 × 10 24 kg is Earth's mass). The angular deceleration, d /dt, is related to D as in which I is the rotating Earth's moment of inertia (8.04 × 10 37 kg m −2 ). An alternative method is given by MacDonald (1964). It is reasonable to assume that the values computed here and shown in Table 1 are underestimated by up to an order of magnitude because of the relatively coarse resolution used in our model. An elevated dissipation of course leads to an increased LOD and the linear dependence clearly shows the importance of particular geologic features as the values in Table 1 do not consider shelf seas. Alone, the real deep ocean dissipation (0.9 TW) gives a LOD of 5.1 s from Equation 6 over a million years, the total real integrated dissipation (3.6 TW) gives 24 s, or for just M 2 , 16 s. Our simulation without shelf seas slows by ∼5 s, compared to 42 s for the full bathymetry, illustrating the outsized effect of shallow water energy losses on the planet's rotational evolution.

Discussion
Continental configuration complexity matters in discussions of a planet's tide because, for any given Earth-like planet, increasing the number of basins and/or embayment type features will increase the probability of resonant amplification. Furthermore, increasing the complexity of coastlines generates further Geophysical Research Letters 10.1029/2019GL085746 energy loss due to acceleration around peninsula-type features. Numerical simulations have been used to test the effect of Earth's former continental configurations and bathymetry on the tide (Bjerrum et al., 2001;Green et al., 2017); however, to our knowledge there are no studies attempting to test a general influence of continental complexity on the ocean tide. This is difficult to do rigorously, given the problems in classifying complex non-Euclidean surfaces such as realistic continents, shapes that were once described as "monstrosities" (Mandelbrot, 1982).
Our results put a constraint on the rotational evolution of planets with oceans and are of use for exoplanet estimates and deep time Earth tidal simulations alike. They also highlight the fundamental effect topography has on ocean tides: Even with a few small continents or ocean bathymetry the dissipation can deviate by an order of magnitude from the water world estimate. There is also a strong link between rotation rate and dissipation (Green et al., 2019b), and exoplanet rotation is currently difficult to determine from observation.
With an equilibrium tide, Kepler-22b (initial rotation period of 1 day) would have a rotation rate of about 2 Earth days after 4.5 Gyr, while Proxima Centauri b (initial rotation period of 1 day) would experience a rotation rate of some 12 days (Barnes, 2017). We have shown that even with simple coastlines and a lack of bathymetry, the spin down could vary with an order of magnitude, and be order(s) of magnitude larger than in the equilibrium case.
As the continents aggregate and disperse during the 400-500 Ma-long, supercontinent cycle (Matthews et al., 2016), the tides become anomalously energetic during periods when the basins meet the conditions for resonance. A feature of this tectonically driven supertidal cycle (Green et al., 2018) is a much less energetic tide during supercontinental periods and a brief, more energetic state, for a few 10 Ma in between the subcontinents (Gotlib & Kagan, 1985;Green et al., 2017). The trend for more complex configurations to be more energetic in our results is analogous to this process and the tectonic setting thus exerts a first-order control on the tidal dissipation. For example, the basin configurations in Figures 2e and 2g are qualitatively similar, yet Figure 2e is substantially more dissipative due to its resonant state, whereas the large meridional continent in Figure 2g blocks the tide. Equally, the range of the dissipation scatter at any point along the complexity scale (R) is roughly an order of magnitude (see Figure 2a), which is comparable to the supertidal fluctuation. Higher R values tend to elevate dissipation, but not indefinitely. As continent size increases, a threshold must be reached where ocean area declines and suppresses dissipation. This is a key geologic control on the rotational evolution (see Table 1).
Earth's earliest continental crust is considered to have mostly formed between 3-2.5 Ga (Korenaga, 2013(Korenaga, , 2018, although this is debated. These early landmasses were probably small (Goodwin, 1996;Scotese, 2004): Estimates for the Late-Archean suggest that ∼12% of Earth's surface area was continental crust, with only 2-3% as emerged land (Flament et al., 2008). This is compared to 42.5% of the surface being continental crust today, of which 27.5% is emerged land area (Cawood et al., 2013;Flament et al., 2008;Schubert & Reymer, 1985). In the microcontinental configurations in Figures 1a-1c, a small semidiurnal tide would be expected, even with the shorter tidal period and day length at the time (Spalding & Fischer, 2019), but they still elevate the global dissipation by nearly an order of magnitude above the equilibrium tide ( Figure 2a). However, this is still 2 orders of magnitude smaller than in simulations with larger, more numerous landmasses that generate more complex basins. It is likely that while Earth had a tidal dissipation rate much greater than the equilibrium tide early on in its history, it may have become more energetic with the formation of larger continents. However this result is not a model for an Archean tide; as shown, shelf-sea area and deep-ocean bathymetry, both unknown for the period, would be key controls.
The amount of energy dissipated in our simulations suggests an upper limit near 1.9 TW for the deep-ocean dissipation, given the similarity between the simulations with and without bathymetry. If the overestimation due to our coarse resolution is considered, a value of ∼0.76 TW may be more realistic, supported by the dissipation value produced without shelf seas (Figure 2). These values do of course depend upon the topographic properties: Even with a randomized bathymetry the influence of continental configuration is clear but is not necessarily a good predictor of the tidal dissipation. Also, the smallest bathymetric features in our grids are of the order of ∼100 km horizontally, and are of a higher relief than the Earth grid. Grid A for example possesses a canyon many thousands of kilometers long. Given the nonlinear dependence on bathymetric slope in Equation 4), it is possible that the conversion fraction is overestimated, but this is an avenue for future work.

10.1029/2019GL085746
The broader implication for Earth-like exoplanets is that while an ocean may raise global tidal dissipation rates, a basin (or basins) with time-varying dimensions due to tectonics will, in a similar fashion to Earth, cause this value to fluctuate considerably over the planet's history. We suggest that it could span 3 orders of magnitude solely due to continental configuration. Planets with complex ocean bathymetry are likely to have larger tidal dissipation rates, and hence spin down much quicker, than more bathtub-like oceans. The apparent convergence generated by adding the SMD bathymetry is an argument for including even a randomized ocean floor in planetary tidal models.
The tidal modulation by Earth's tectonics gives these results a robust standard to compare against, but the complexity metric may also have implications for Earth's early history, for periods with limited records of continental position and size.