Tidal regime shift in Lingdingyang Bay,the Pearl River Delta: An identification and assessment of driving factors

Both natural changes (e.g., tidal forcing from the ocean and global sea level rise) and human-induced changes (e.g., dredging for navigation, sand excavation, and land reclamation) exert considerable influences on the long-term evolution of tidal regimes in estuaries. Evaluating the impacts of these factors on tidal-regime shifts is particularly important for the protection and management of estuarine environments. In this study, an analytical approach is developed to investigate the impacts of estuarine morphological alterations (mean water depth and width convergence length) on tidal hydrodynamics in Lingdingyang Bay, Southeast China. Based on the observed tidal levels from two tidal gauging stations along the channel, tidal wave celerity and tidal damping/amplification rate of different tidal constituents are computed using tidal amplitude and phase of tidal constituents extracted from a standard harmonic analysis. We show that the minimum mean water depth for the whole estuary occurred in 2006, whereas a shift in tidal wave celerity for the M₂ tide component occurred in 2009. As such, the study period (1990–2016) could be separated into pre-human (1990–2009) and post-human (2010–2016) phases. Our results show that the damping/amplification rate and celerity of the M₂ tide have increased by 31% (from 7 to 9.2 m⁻¹) and 28% (from 7 to 9 m·s⁻¹) respectively, as a consequence of the substantial impacts of human interventions. The proposed analytical method is subsequently applied to analyse the historical development of tidal hydrodynamics and regime shifts induced by human interventions, thus linking the evolution of estuarine morphology to the dominant tidal hydrodynamics along the channel. The observed tidal regime shift is primarily caused by channel deepening, which substantially enlarged the estuary and reduced effective bottom friction resulting in faster celerity and stronger wave amplification. Our proposed method for quantifying the impacts of human interventions on tidal regime shifts can inform evidence-based guidelines for evaluating hydraulic responses to future engineering activities.     

Modelling the global ocean tides: modern insights from FES2004

During the 1990s, a large number of new tidal atlases were developed, primarily to provide accurate tidal corrections for satellite altimetry applications. During this decade, the French tidal group (FTG), led by C. Le Provost, produced a series of finite element solutions (FES) tidal atlases, among which FES2004 is the latest release, computed from the tidal hydrodynamic equations and data assimilation. The aim of this paper is to review the state of the art of tidal modelling and the progress achieved during this past decade. The first sections summarise the general FTG approach to modelling the global tides. In the following sections, we introduce the FES2004 tidal atlas and validate the model against in situ and satellite data. We demonstrate the higher accuracy of the FES2004 release compared to earlier FES tidal atlases, and we recommend its use in tidal applications. The final section focuses on the new dissipation term added to the equations, which aims to account for the conversion of barotropic energy into internal tidal energy. There is a huge improvement in the hydrodynamic tidal solution and energy budget obtained when this term is taken into account.     

Effect of channel deepening on tidal flow and sediment transport: part I—sandy channels

Natural tidal channels often need deepening for navigation purposes (to facilitate larger vessels). Deepening often leads to tidal amplification, salinity intrusion, and increasing sand and mud import. These effects can be modelled and studied by using detailed 3D models. Reliable simplified models for a first quick evaluation are however lacking. This paper presents a simplified model for sand transport in prismatic and converging tidal channels. The simplified model is a local model neglecting horizontal sand transport gradients. The latter can be included by coupling (as post-processing) the simplified model to a 2DH or 3D flow model. Basic sand transport processes in stratified tidal flow are studied based on the typical example of the tidal Rotterdam Waterway in The Netherlands. The objective is to gain quantitative understanding of the effects of channel deepening on tidal penetration, salinity intrusion, tidal asymmetry, residual density-driven flow, and the net tide-integrated sand transport. We firstly study the most relevant tidal parameters at the mouth and along the channel with simple linear tidal models and numerical 2DH and 3D tidal models. We then present a simplified model describing the transport of sand (TSAND) in tidal channels. The TSAND model can be used to compute the variation of the depth-integrated suspended sand transport and total sand transport (incl. bed-load transport) over the tidal cycle. The model can either be used in stand-alone mode or with computed near-bed velocities from a 3D hydrodynamic model as input data.     

Effect of a secondary channel on the non-linear tidal dynamics in a semi-enclosed channel: a simple model

Non-linear tidal dynamics are investigated in a network that consists of a semi-enclosed main channel and a secondary channel at an arbitrary position. The water motion, governed by the one-dimensional shallow water equations, is forced by an incoming tidal wave. Solutions are obtained with the method of characteristics. The overall aim is to quantify and understand the spatial structure of different tidal harmonics (the principal tide and its non-linear overtides) and of tidal asymmetry for both the vertical and the horizontal tide in the main channel for different locations of the secondary channel. This is of practical interest in the context of possible construction of secondary channels to reduce tidal range in estuaries. Moreover, tidal asymmetry is an important factor in driving net sediment transport. Analysis of the different tidal harmonics shows that their characteristics are similar to those obtained with an earlier linear model. In particular, amplitudes of the harmonics are reduced landward of the secondary channel if the latter is positioned less than a quarter wavelength of the respective tidal wave away from the landward boundary. Thus, the distortions of the tide due to the presence of the secondary channel are generated locally and afterwards propagate through the network. Tidal asymmetry is quantified by examining tidal range, flood-to-ebb ratio and the duration of the falling tide and the duration between maximum flood and maximum ebb. A spatial non-uniform reduction in tidal range is observed that shows very localised increase and decrease depending on the position of the secondary channel. The changes in the velocity characteristics induce changes in net sediment transport. It turns out that the direction of the peak current, derived from the flood-to-ebb ratio, is not sensitive to the position of the secondary channel, whereas the duration between flood and ebb can change from more to less than half the tidal cycle. However, the changes in the velocity asymmetries are confined to a small region.     

Analytical study of the transverse distribution of along-channel and transverse residual flows in tidal estuaries

We present an analytical model to decompose complex along-channel and transverse residual flows into components induced by individual mechanisms. The model describes the transverse distribution of residual flows in tidally dominated estuaries. Scaling and perturbation techniques are used to obtain analytical solutions for residual flows over arbitrary across-channel bed profiles. The flows are induced by horizontal density gradients, tidal rectification processes, river discharge, wind, channel curvature and the earth's rotation. These rectification processes induce residual flows that are up-estuary to the right and down-estuary to the left of an estuarine channel (looking up-estuary in the northern hemisphere). The tidal rectification processes fundamentally change the transverse structure of along-channel residual flows in many tidal estuaries, as these processes cause the flows to be internally asymmetric about the mid-axis of the channel for relatively large tidal velocities, steep channels or narrow estuaries. In addition, velocity scales are derived from the analytical solutions to estimate the relative importance of the various residual flow mechanisms from estuarine parameters. A case study of a transect across the Upper Chesapeake Bay showed that important features of the residual flow observed in that transect are reproduced and explained by the analytical model. The velocity scales were able to identify the relevant residual flow mechanisms as well. The tidal rectification processes considered here result from advection of along-channel tidal momentum by Coriolis-induced transverse tidal currents.     

A simple approach to adjust tidal forcing in fjord models

To model currents in a fjord accurate tidal forcing is of extreme importance. Due to complex topography with narrow and shallow straits, the tides in the innermost parts of a fjord are both shifted in phase and altered in amplitude compared to the tides in the open water outside the fjord. Commonly, coastal tide information extracted from global or regional models is used on the boundary of the fjord model. Since tides vary over short distances in shallower waters close to the coast, the global and regional tidal forcings are usually too coarse to achieve sufficiently accurate tides in fjords. We present a straightforward method to remedy this problem by simply adjusting the tides to fit the observed tides at the entrance of the fjord. To evaluate the method, we present results from the Oslofjord, Norway. A model for the fjord is first run using raw tidal forcing on its open boundary. By comparing modelled and observed time series of water level at a tidal gauge station close to the open boundary of the model, a factor for the amplitude and a shift in phase are computed. The amplitude factor and the phase shift are then applied to produce adjusted tidal forcing at the open boundary. Next, we rerun the fjord model using the adjusted tidal forcing. The results from the two runs are then compared to independent observations inside the fjord in terms of amplitude and phases of the various tidal components, the total tidal water level, and the depth integrated tidal currents. The results show improvements in the modelled tides in both the outer, and more importantly, the inner parts of the fjord.     

A model for simulating the response of runoff from the mountainous watersheds of inland river basins in the arid area of northwest China to climatic changes

The subaerial tidal sand area in the northern Jiangsu Province (Subei), stretching from Dongtai towards east with a fan shape, is an early developing stage of radial sand ridges distributed in the South Yellow Sea. Since 5000–6000 a BP, after the Holocene transgression maximum in the northern Jiangsu Province, subaqueous tidal sand bodies were exposed and changed into land gradually. The environmental magnetism analysis shows that subaerial tidal sand strata are formed by the convergent-divergent paleo-tidal current field. The sediment source of tidal sand strata came early from the Changjiang River and late from the Yellow River. Sea floor erosion by tidal currents also served as an important sand source. Drilling cores and ground-penetrating profile show that there exists no probability of sand supplying directly by a large river through the apical area of tidal sand ridges either on land or in the sea. Fluvial deposits supplied the tidal sand bodies by alongshore transportation, which corresponds to the conclusions obtained by the analyses of provenance and paleocurrent field. Project supported by the National Natural Science Foundation of China (Grant Nos. 43236120 and 49676288).     

中国大陆精密重力潮汐改正模型

利用理论和实验重力固体潮模型,充分考虑全球海潮和中国近海潮汐的负荷效应,建立了中国大陆的精密重力潮汐改正模型.结果表明,采用不同的固体潮模型会对重力潮汐结果产生相对变化幅度小于0.06％的差异;在沿海地区海潮负荷的影响约为整个潮汐的4％,而中部地区约为1％,其中中国近海潮汐模型的影响约占整个海潮负荷的10%,内插或外推潮波的负荷约占海潮负荷的3％.通过比较实测的重力数据表明,本文给出的重力潮汐改正模型的精度远远优于0.5×10^-8 m·s^-2,说明了本文构建的模型的实用性,可为中国大陆高精度重力测量提供有效参考和精密的改正模型.     

Local and remote response of the North Sea dynamics to morphodynamic changes in the Wadden Sea

The response of the tidal system in the southern North Sea to morphodynamic changes was investigated in a modelling study using fine resolution bathymetric observations available for 1982–2011. The Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM) was set up for the different sets of bathymetries. One set of bathymetry was compiled from a large number of bathymetric measurements over many years, while the other two reflected bathymetry state in the area of Wadden Sea during 2000 and 2011, respectively. The temporal and spatial evolution of bathymetry was dominated by migration of tidal channels. The M4 tide showed larger sensitivity to bathymetric change in the Wadden Sea than the M2 tide, whereas the structure of the latter remained rather robust. The largest change of the tidal wave due to the differences in bathymetries was located off the North Frisian Wadden Sea. Traces of changes were also found far away from the regions of their origin because the tidal waves in the North Sea propagate the local disturbances basin-wide. This illustrated an efficient physical mechanism of teleconnectivity, i.e. effecting the local responses to the larger-scale or remote change of ocean bottom caused by erosion and deposition. The tidal distortion resulting from the relatively small bathymetric changes was substantial, particularly in the coastal zone. This is a manifestation of the nonlinear tidal transformation in shallow oceans and is crucial for the sediment transport and the morphodynamic feedback, because of the altered tidal asymmetry.     

Determination of the freshwater budget of tidal flats from measurements near a tidal inlet

The freshwater budget of a tidal flat area is evaluated from long-term hydrographic time series from an observation pole positioned in a tidal channel in the Hörnum Basin (Germany). For each tidal cycle, the freshwater budget is calculated from the total imported and exported water volumes and the corresponding mean densities. The variability of the budget on a tidal scale is characterised by a period of twice the tidal period, exhibiting a minimum when the tidal flats are dry around daylight hours during the foregoing low tide, and a maximum when low tide occurs at night; enhanced evaporation on the flats at daylight hours is identified as the driving process. On the average over one year, while winter observations are missing, the freshwater budget is negative for the years 2002–2005 and positive only for 2006. The interannual mean is negative and amounts to a freshwater loss of about 2 mm day⁻¹, although the large-scale climate in this region is humid. The results demonstrate that the bulk parametrisations for the latent and sensible heat flux between the ocean and the atmosphere must not be applied for the tidelands.     

A conceptual model of depositional,rather than erosional,tidal channel development in the rapidly prograding Skagit River Delta (Washington,USA)

The origin and growth of blind tidal channels is generally considered to be an erosional process. This paper describes a contrasting depositional model for blind tidal channel origin and development in the Skagit River delta, Washington, USA. Chronological sequences of historical maps and photos spanning the last century show that as sediments accumulated at the river mouth, vegetation colonization created marsh islands that splintered the river into distributaries. The marsh islands coalesced when intervening distributary channels gradually narrowed and finally closed at the upstream end to form a blind tidal channel, or at mid‐length to form two blind tidal channels. Channel closure was probably often mediated through gradient reduction associated with marsh progradation and channel lengthening, coupled with large woody debris blockages. Blind tidal channel evolution from distributaries was common in the Skagit marshes from 1889 to the present, and it can account for the origin of very small modern blind tidal channels. The smallest observed distributary‐derived modern blind tidal channels have mean widths of 0·3 m, at the resolution limit of the modern orthophotographs. While channel initiation and persistence are similar processes in erosional systems, they are different processes in this depositional model. Once a channel is obstructed and isolated from distributary flow, only tidal flow remains and channel persistence becomes a function of tidal prism and tidal or wind/wave erosion. In rapidly prograding systems like the Skagit, blind tidal channel networks are probably inherited from the antecedent distributary network. Examination of large‐scale channel network geometry of such systems should therefore consider distributaries and blind tidal channels part of a common channel network and not entirely distinct elements of the system. Finally, managers of tidal habitat restoration projects generally assume an erosional model of tidal channel development. However, under circumstances conducive to progradation, depositional channel development may prevail instead. Copyright © 2006 John Wiley & Sons, Ltd.     

The effect of tide on the hydrology and morphology of a freshwater river

How does river hydrology and morphology change due to tidal influence? We contend that this is a question of particular consequence to many earth surface disciplines, but one that has not been adequately addressed. Previous studies have relied on gradients in channel morphology and stratigraphy to infer energy regime of channels. However, in tidal rivers geomorphology influences the energy regime while the energy regime influences morphology; thus, geomorphic and stratigraphic patterns do not fully resolve the mechanisms which lead to change. We addressed this problem by comparing measurements of hydraulic energy and channel morphology along a tidal gradient to predictions of these characteristics in the absence of tides, and attributed the differences to tidal processes. Measurements of discharge, channel area, and energy dissipation (in kJ day^–1) were made over a 24·8 hour period at four sites spanning the non‐tidal to tidal freshwater Newport River, NC. We then predicted those characteristics under non‐tidal conditions using hydraulic geometry relationships and literature values from coastal plain rivers. Discharge was enhanced more than 10‐fold by tide, and this tidal effect increased from upstream to downstream along the tidal gradient. Cross‐sectional area increased three‐fold due to tide. Energy dissipation measured in the upper tidal river was four‐fold lower than predicted to occur in the absence of tide because tides decreased average velocity and discharge. Energy dissipation measured downstream was similar to that predicted to occur without tides, although there was large uncertainty in predicted values downstream. While this limited dataset does not permit us to make broad generalizations for definitive models, it does provide a proof‐of‐concept for a new approach to addressing a critical problem at the interface of fluvial and coastal morphology. Copyright © 2013 John Wiley & Sons, Ltd.     

Morphology and dynamics of the intertidal floodplain along the Amazon tidal river

Depositional environments along the tidal river downstream of Óbidos have been proposed as important sinks for up to one third of the sediment discharge from the Amazon River. However, the morphology of the intertidal floodplain and the dynamics of sediment exchange along this reach have yet to be described. River-bank surveys in five regions along the Amazon tidal river reveal a distinct transition in bank morphology between the upper, central and lower reaches of the tidal river. The upper tidal-river floodplain is defined by prominent natural levees that control the transfer of water and sediment between the mainstem Amazon River and its floodplain. Greater tidal influence in the central tidal river suppresses levee development, and tidal currents increase sediment transport into the distal parts of the floodplain. In the lower tidal river, the floodplain morphology closely resembles marine intertidal environments (e.g. mud flats, salt marshes), with dendritic tidal channels incising elevated vegetated flats. Theory, morphology and geochronology suggest that the dynamics of sediment delivery to the intertidal floodplain of the Amazon tidal river vary along its length due to the relative influence and coupling of fluvial and tidal dynamics. © 2018 John Wiley & Sons, Ltd.     

Succession of macrobenthic fauna and nitrogen budget at two artificial tidal flats in Osaka Bay, Japan

Nitrogen budgets and predominant benthic organisms were examined at a matured artificial tidal flat of Osaka Nanko bird sanctuary and an artificial young tidal flat of Hannan Second District, Osaka Bay, Japan. At the Osaka Nanko bird sanctuary, indexes of the quality of the bottom sediment increased over the time course and no abiotic area was observed in the macrobenthic fauna. Dissolved inorganic nitrogen imported into the tidal flat changed its form and was exported in the form of dissolved organic nitrogen. As a result, total nitrogen was trapped at the tidal flat at an average rate of 138 mgN m(-2)day(-1). At the young tidal flat of Hannan Second District net exchange of nitrogen varied markedly among 2000, 2001, 2002 and 2003. In September 2000 when only 5 months had passed after construction, the tidal flat functioned as a site of source for nitrogen. However, it changed as a sink for nitrogen in 2001-2002 in relation to the growth of seaweeds (Ulva spp. and Gracilaria vermiculophylla) and clams (Tapes philippinarum). In 2003 nitrogen was again released from the artificial tidal flat due to the decrease of biomass of the dominant organisms. Comparing the nitrogen budget between the two sites on a per unit time and area basis, net exchanges of nitrogen by tidal exchange were apparently low and relatively constant at the tidal flat of Osaka Nanko bird sanctuary. This suggests that the aged artificial tidal flat of Osaka Nanko bird sanctuary has a stable ecosystem compared to the young artificial tidal flat of Hannan Second District.     

Sensitivity of tidal sand wavelength to environmental parameters: A combined data analysis and modelling approach

An integrated field data-modelling approach is employed to investigate relationships between the wavelength of tidal sand waves and four environmental parameters: tidal current amplitude, water depth, tidal ellipticity and median grain size. From echo sounder data at 23 locations on the Dutch continental shelf, the average wavelengths of observed sand waves are determined and compared with the wavelengths obtained with a process-based model. The latter describes the initial formation of these bedforms due to feedbacks between the tidal current and the erodible bed and uses environmental parameters for the 23 locations as input. Good agreement between observed and modelled wavelengths is found if the bottom stress experienced by tidal currents is adequately quantified. Model results show that the wavelength of sand waves increases with increasing water depth, tidal ellipticity and grain size (coarse sand), whilst it decreases with increasing tidal current amplitude and grain size (fine sand). Due to the limited number of stations and the fact that all four parameters change from location to location, the modelled relationships are only partly supported by the field observations.     

A numerical study of horizontal dispersion in a macro tidal basin

Tidal circulation in Cobscook Bay, a macro tidal basin, is simulated using the three-dimensional, nonlinear, finite element ocean model, QUODDY_dry. Numerical particles are released from various transects in the bay at different tidal phases and tracked for several tidal cycles. Initially, nearby particles in the main tidal channel experience a great deal of spreading and straining, and after a few tidal cycles, they are separated in different parts of the bay. The fundamental mechanism for particle dispersion is the chaotic advection that arises from long tidal excursions passing through many residual eddies. A loosely correlated, inverse relationship between the two dimensionless parameters, ν (the ratio of the residual current to the tidal current) and λ (the ratio of the tidal excursion to the main topographic scale), can be constructed for large values of ν. Several Lagrangian statistical measures are used to quantify and distinguish dispersion regimes in different parts of Cobscook Bay. It is found that the effective Lagrangian dispersion coefficient can be estimated using the product of the magnitude of residual currents and the tidal excursion.     

Tidal inlet velocity asymmetry in diurnal regimes

Tidal velocity asymmetry at inlets influences sediment transport pathways and the morphological evolution of estuaries/lagoons connected to these inlets. Generation of overtides is generally seen as the main cause of tidal velocity asymmetry. Whilst majority of studies examining tidal velocity asymmetry have concentrated on inlets located in semi-diurnal tidal regimes, here, attention is focused on the processes responsible for causing tidal velocity asymmetry at inlets located in diurnal tidal regimes. Using field data collected from three West Australian inlets, it is shown that tidal velocity asymmetry in this type of system is caused by the oceanic tidal conditions. It is also shown that in these systems, the occurrence of flood/ebb dominance can be determined using oceanic tidal elevations, which are more readily available than inlet current data. In contrast to semi-diurnal systems the flood/ebb dominance in diurnal systems varies throughout the year depending on the phase angle relationship between the significant oceanic tidal constituents. The net sediment transport in to/out of these systems, which determines the morphological evolution of the systems, is shown to be governed more by the degree of tidal velocity asymmetry rather than the number of occurrences or duration of flood/ebb-dominant periods.     

Tidally generated island wakes and surface water cooling over Izu Ridge

The Kuroshio current is well known for generating cold wakes behind islands over Izu Ridge in Northwestern Pacific. Observational data from the geostationary Himawari-8 satellite for 2015–2017 revealed the occurrence of cold waters during the period when the Kuroshio current flows away from the islands. With a focus on tidal currents, this study presents an investigation of dynamical processes responsible for the formation of areas with low sea surface temperature (SST) through the adoption of a high-resolution numerical ocean model for an event that happened in July 2017. Areas with cold water emerged only when tidal currents are included in the numerical model. The model results indicate the cold surface waters are formed in the vicinity of the islands because of upwelling and vertical mixing. Qualitative features of the cold water formation for each island are found to depend on its size, topography, and ambient currents. Near Kozu Island, the tidal excursion is large enough to cause eddy shedding. These shed eddies are stirred by tidal currents to extend the surface cooling effect to wider areas. Near Hachijo Island, a persistent wake is formed by the ambient northward current. Inclusion of tidal currents destabilizes the wake, and consequently leads to the formation of a low SST area, although no clear eddy shedding is detected. The flow patterns around the islands are classified using an additional non-dimensional parameter, defined as the ratio between tidal excursion and island diameter.

     

Lunar seismicity and tectonics

Seismic signals from 300–700 deep moonquakes and about four shallow moonquakes are detected by the long-period seismometers of two or more of the Apollo seismic stations annually. Deep-moonquake activity detected by the Apollo seismic network displays tidal periodicities of 0.5 and 1 month, 206 d and 6 a. Repetitive moonquakes from 60 hypocenters produce seismograms characteristic of each. At each hypocenter, moonquakes occur only within an active period of a few days during a characteristic phase of the monthly lunar tidal cycle. An episode of activity may contain up to four quakes from one hypocenter. Nearly equal numbers of hypocenters are active at opposite phases of the monthly cycle, accounting for the 0.5-month periodicity. The 0.5- and 1-month activity peaks occur near times of extreme latitudinal and longitudinal librations and earth-moon separation (EMS). The 206-d and 6-a periodicities in moonquake occurrence and energy release characteristics are associated with the phase variations between the librations and EMS. Because of the exact relationship between tidal phases and the occurrence of deep moonquakes from a particular hypocenter, it is possible to predict not only the occurrence times from month to month, often to within several hours, but also the magnitudes of the moonquakes from that hypocenter. The predicted occurrence of large A₁ moonquakes in 1975, following a 3-a hiatus, confirms the correlation between A₁-moonquake activity and the 6-a lunar tidal cycle and implies a similar resurgence for all of the deep moonquakes. Because no matching shallow moonquake signals have been identified to date, tidal periodicities cannot be identified for the individual sources. However, shallow moonquakes generally occur near the times of extreme librations and EMS and often near the same tidal phase as the closest deep moonquake epicenters. With several possible exceptations, the deep-moonquake foci located to date occur in three narrow belts on the nearside of the moon. The belts are 100–300 km wide, 1,000–2,500 km long and 800–1,000 km deep and define a global fracture system that intersects in central Oceanus Procellarum. A fourth active, although poorly defined, zone is indicated. The locations of 17 shallow-moonquake foci, although not as accurate as the deep foci, show fair agreement with the deep-moonquake belts. Focal depths calculated for the shallow moonquakes range from 0–200 km. Deep-moonquake magnitudes range from 0.5 to 1.3 on the Richter scale with a total energy release estimated to be about 10¹¹ erg annually. The largest shallow moonquakes have magnitudes of 4–5 and release about 10¹⁵–10¹⁸ erg each. Tidal deformation of a rigid lunar lithosphere overlying a reduced-rigidity asthenosphere leads to stress and strain concentrations near the base of the lithosphere at the level of the deep moonquakes. Although tidal strain energy can account for the deep moonquakes in this model, it cannot account for the shallow moonquakes. The tidal stresses within the lunar lithosphere range from about 0.1 to 1 bar and are insufficient to generate moonquakes in unfractured rock, suggesting that lunar tides act as a triggering mechanism. The largest deep moonquakes of each belt usually occur near the same characteristic tidal phases corresponding to near minimum or maximum tidal stress, increasing tidal stress, and alignments of tidal shear stresses that correspond to thrust faulting along planes parallel to the moonquake belts and dipping 30–40°. With few exceptions, the shallow moonquakes occur at times of near minimum tidal stress conditions and increasing tidal stress that also suggest thrust faulting. The secular accumulation of strain energy required for the shallow moonquakes and implied by the uniform polarities of the deep moonquake signals probably results from weak convection. A convective mechanism would explain the close association between moonquake locations and the distribution of filled mare basins and thin lunar crust, the earth-side topographic bulge, and the ancient lunar magnetic field. The low level of lunar seismic activity and the occurrence of thrust faulting both at shallow and great depths implies that the moon is presently cooling and contracting at a slow rate.