Take horizontal motion towards the second projectile as positive. In projectile motion, both bodies have the same downward gravitational acceleration, so it cancels from their relative motion. Their coordinates are
For a collision at , equality of the horizontal coordinates requires . Equality of the heights then requires . Hence
The upward solution requires . A collision strictly above the ground further requires , or . The horizontal target distance does not independently affect the interception geometry. If the first shot is also required to hit that ground target, its parameters satisfy .
The specified rocket fuel fraction gives , so . Imposing yields
Then fixes
Numerically,
The associated constant mass loss rate is . The initial upward thrust acceleration already exceeds and increases as mass falls, so the downward speed decreases monotonically to zero. This calculation assumes constant lunar gravitational acceleration, constant relative rocket exhaust speed, and negligible drag.
Take downward distance and downward speed , with constant gravitational acceleration and quadratic drag coefficient . The printed vector square has no standard vector meaning; the directionally consistent drag force is . For this downward fall, Newton's second law becomes
The terminal velocity is obtained by setting the net force to zero:
From rest, the speed increases monotonically towards this limiting value. In fact ; thus it is a supremum approached at infinite time, and a tower of finite height need not permit it to be reached.
For the required distance, use the chain rule . Writing gives a linear ordinary differential equation:
Consequently,
Here has units of inverse length, so the result has the required units of distance and is independent of the falling mass.
Powered ascent. Use upward signed velocity, constant gravitational acceleration , constant relative rocket exhaust speed , vertical flight from rest, no drag force, and no impulse at stage separation. The constant- approximation is an explicit modelling assumption; a varying Newtonian gravitational field would require integrating a different ascent equation. Put and . During ascent , and momentum balance gives the rocket equation with gravity:
The launch rocket thrust exceeds , so immediate lift-off is consistent. At ascent burnout the entire attached vehicle still has mass , not ; the latter is the separated first stage. The burn duration and endpoint speed are
Coasting height. Integrating the ascent velocity from zero height gives
For example this follows by substituting in . At separation,
The separated stage then undergoes free fall in a constant gravitational field with initial upward speed . Conservation of mechanical energy during this coast gives
Its smaller mass after separation does not change that coasting height.
Reserve burn. On descending to , conservation of mechanical energy gives downward speed
Reset time to at reserve ignition. Now , and . Downward ejection supplies upward rocket thrust. The upward signed velocity is
Thus the formal free-flight burnout value is
corresponding to a formal downward speed of . Keeping throughout instead gives downward; premature rounding is especially damaging because three large velocity terms almost cancel.
The free-flight endpoint is below ground. Integrating the same reserve-burn velocity gives
At the stated endpoint this gives for . The surface is crossed at approximately , before reserve burnout. The stage has already hit the surface, so a physical post-impact burnout speed cannot be inferred without a ground-contact model. The boxed velocity is the intended ideal continuation, not an achievable free-flight landing result under these exact assumptions. In particular, the given numbers do not describe an exact soft landing. The constant-rate retro-rocket landing must satisfy both its velocity and height conditions.
Figure 1.
Constant-gravity reserve descent: the formal trajectory crosses ground before its nominal burnout
. Details.
For one wire, Ampère's circuital law on a concentric circle of radius gives
The other wire experiences the Lorentz force on its electric current, . Its direction is perpendicular to this magnetic field, so the force per unit length is
The right-hand rule gives attraction for parallel electric currents in the same direction and repulsion for opposite directions.
For the upper wire, the upward magnetic force per unit length balances its weight per unit length, , where is the gravitational acceleration. Therefore
This assumes the lower wire is held fixed and both wires remain straight. Within the vertical-motion model, increasing the separation weakens the upward magnetic force, so this equilibrium is vertically restoring.
A tidal force measures the difference of gravitational acceleration across an extended object. A freely falling reference frame removes a spatially uniform acceleration, but not its gradient. For a distant point mass at distance and direction from the reference centre, a small internal displacement gives
The tidal tensor has eigenvalues radially and in each transverse direction: radial stretching and transverse squeezing. Its trace vanishes in the source-free region, consistent with the gravitational Gauss law. An associated relative potential is .
The scale is across size , not the centre's absolute gravity . This distinction explains how a more distant massive object can exert a large force but a modest tide. Opposite-side stretching supports the elementary two-bulge picture; real ocean tides also involve rotation, basin geometry and fluid dynamics, so the static gradient alone does not predict coastal heights or timing. The same gradient stresses satellites and can disrupt a self-gravitating body when it competes with the body's own gravitational restoring acceleration.
Take upward as positive and interpret the stated as the total mass immediately before launch, including the harpoon. The remaining probe mass is . Conservation of momentum during the short firing interval gives
The harpoon speed used here is relative to the comet, as printed. From , the near-surface gravitational acceleration is
For the near-surface recoil height with constant downward thrust, the upward motion has downward acceleration . Therefore
Since , treating gravity as constant is consistent. This uses only the stated gravity and thrust after firing; a taut tether could reduce this free-flight height, but no tether length or elasticity is specified. If instead meant the probe mass after firing, the pre-firing mass would be and its momentum would leave the probe still moving downward, with no upward excursion. This illustrates nearly cancelling recoil momenta in payload ejection. The mass convention matters because the launch momenta nearly cancel.
Archimedes' principle states that the resultant hydrostatic pressure force on a body equals the weight of the fluid displaced, directed opposite to gravity. Let be the displaced fluid region, with outward normal from the body into the fluid. Since the fluid pressure force is inward,
Hydrostatic balance is , with downward gravitational acceleration . Therefore
for uniform fluid density and gravity. Uniform atmospheric pressure gives no net force on a closed surface. For a floating body, the same displaced-volume construction includes the free surface; any significant air buoyancy can be included as a second displaced fluid.
A freely floating body of mass satisfies , so its displaced volume is . A homogeneous body of density floats with submerged fraction ; if , hydrostatic buoyancy alone cannot support it. The buoyant force acts through the centre of buoyancy, the centroid of the displaced volume for constant density, while weight acts through the centre of mass. Their relative positions determine rotational equilibrium and stability. For a fully submerged body, a centre of mass below the centre of buoyancy gives a restoring torque; a floating body's changing displaced shape additionally requires the metacentric height criterion. Surface tension, externally driven flow and contact forces can alter the simple force balance.
The differential gravitational acceleration responsible for a tidal force is the difference between accelerations at separated points of a body. Relative to a freely falling reference point , a small offset experiences
The tidal tensor has units ; the vector called in the source is the resulting differential acceleration, in . An actual differential force is this vector multiplied by the test mass. Measurable effects include relative motion of separated freely falling masses, radial stretching, transverse squeezing, deformation and, when self-gravity or strength is insufficient, disruption. A spatially uniform acceleration alone produces no tidal deformation.
A spatially uniform external gravitational acceleration moves a body's centre of mass without stretching it. Tidal forces arise from the variation of that acceleration across the body. Let a distant mass lie at distance along unit vector , and let denote position relative to the body's centre. Subtract the centre's acceleration. To first order,
The tidal tensor is . It stretches along the line to the source with eigenvalue and compresses in the two perpendicular directions with eigenvalues . Its trace vanishes outside the source, consistently with Laplace's equation for the gravitational potential.
After removal of the constant and uniform-acceleration terms, the tidal potential is
An ideal equilibrium ocean on a sphere of radius follows an equipotential of the body's gravity plus this perturbation. For surface gravity , the leading displacement is
up to a constant fixing total water volume. The quadrupolar form produces bulges on both the near and far sides; the far-side bulge follows from differential acceleration, without postulating a repulsive gravitational force.
The Moon and Sun both raise terrestrial tides, with relative strength set mainly by , rather than by the total gravitational pull . Their contributions reinforce near new and full Moon to produce spring tides, and partly cancel near quadrature to produce neap tides. Earth's rotation relative to the tidal pattern gives a strong semidiurnal component in the simplest geometry. Real tide heights and timing depend on ocean depth, coastlines, basin resonances, Coriolis acceleration, friction and the declination of the forcing bodies. Tidal dissipation produces phase lags and torques, exchanging rotational and orbital angular momentum and gradually altering the spin and orbit. The equilibrium bulge is therefore a useful starting model, not an exact local tide forecast.
The uniform weight per unit length is , where is the mass density and the gravitational acceleration. The total weight is , so each support provides an upward force . The weight on the portion from to is acting at that portion's midpoint, a distance from the cut. In the same bending moment convention as part (c),
The units are . Its greatest magnitude occurs at , where .
Far from its ends, approximate the body by an infinite uniform cylinder. By cylindrical symmetry, the Newtonian gravitational field is radial. Applying Gauss's law for gravity to a coaxial cylindrical surface of radius and length outside the material gives
At the surface, , so
The gravitational acceleration points inward toward the axis. Its small magnitude reflects both the small size and the modest density. The infinite-cylinder approximation uses the specified distance from both ends; a spherical approximation would discard the relevant geometry. No rotational contribution is included, since no rotation rate is supplied.
Powered vertical flight. Take upward velocity as positive, and let be the downward exhaust speed relative to the rocket. During the rocket's mass changes by , and exhaust mass leaves with Earth-frame velocity . Compare the momentum of the rocket plus this newly expelled exhaust before and after the interval:
Discard the second-order product . The result is , or
The positive rocket thrust is . Integrating with and gives the rocket equation with gravitational loss:
Using for the rocket alone without the exhaust momentum flux would give the wrong equation for this open system.
Orbital speed under the stated gravity assumption. A circular orbit of radius requires centripetal acceleration . If the magnitude of gravitational acceleration is the stipulated constant even at the orbit,
This treats gravity as centrally directed for the orbital calculation, while retaining its constant magnitude. The subsequent flat-surface approximation treats its direction as locally vertical.
Burnout height. During the unpowered projectile motion, horizontal velocity is constant. At the top of the trajectory the vertical component vanishes, so the horizontal component must equal the orbital speed . At burnout, therefore,
These follow either by integrating or by the vertical work-energy theorem. Substituting the preceding constant- orbital result gives
The vertical component is initially positive, as required to reach the orbit after burnout.
Alternative inverse-square interpretation. A real Earth's exterior Newtonian gravitational field is not constant over a altitude change. If instead means the surface value , then Newton's law of universal gravitation gives and
Using this orbital speed but the paper's constant- flat projectile motion approximation for the coast would give
This alternative changes the global constant-gravity assumption; it is not the same model as the result. An exact curved Kepler orbit coast would also abandon the prescribed flat trajectory. The printed assumptions support the constant-magnitude answers above, while the inverse-square values explain the different answer obtained from the usual physical Earth-orbit formula.
The tunnel is taken thin enough that removal of its material does not significantly disturb spherical symmetry. For , the enclosed mass is
By the spherical shell theorem, the inward gravitational acceleration has magnitude
The released gravitational potential energy is the positive work done by gravity while the mass falls from radius to :
By conservation of mechanical energy,
Inside the cavity, the spherical shell theorem gives zero gravitational field, so the velocity stays constant. The original figure places B and C at opposite ends of a cavity diameter, a distance apart. Therefore
The quantity has units , so this is a time. There are no numerical values for or in the paper. The transit time is independent of and the falling mass . The cavity's gravitational potential is constant, not necessarily zero; only its gradient vanishes. The calculation is a hollow-sphere gravitational transit.
For a prograde equatorial synchronous circular orbit, equate gravitational acceleration and centripetal acceleration. With planetary radius , , and rotation period ,
The tabulated is a gravitational acceleration; its printed unit must be read as , rather than . Using gives
The orbital altitude is measured from the surface, rather than the centre.
At the visible horizon, the line of sight is tangent to the spherical surface. The limiting central angle satisfies . The visible spherical cap has area , giving
About of the whole surface is geometrically visible, with . This is a surface-area fraction, not an apparent-disc area ratio. Useful imaging near the horizon is further limited by foreshortening and the long atmospheric path.
For constant gravitational acceleration, hydrostatic balance gives the atmospheric column mass , with negligible pressure at the top. For the same surface temperature , the variable latent heat gives a smaller integrated , and hence
Multiplying by the same planetary surface area gives the same comparison for total mass. The models agree for a surface at or below the triple point, where they use the same sublimation branch. The comparison assumes the same surface temperature, rather than the same surface pressure.
Choose a common horizontal compensation depth beneath the deepest cooling-induced density anomaly, within mantle that can relax pressure differences over the timescale considered. In a normal oceanic plate model, this means below the old plate's thermal lithosphere, not at the base of the approximately 7 km thick oceanic crust.
The reason is that cooling changes density through both crust and lithospheric mantle. A datum at the Mohorovičić discontinuity would omit most of the changing column mass. A different datum for each age would also compare unlike columns. Once the common datum is sufficiently deep, extending it farther into laterally uniform mantle adds the same mass per unit area to each column and does not change the predicted relative subsidence.
There is no unique number fixed by the question. For an illustrative oceanic plate cooling model with a thermal thickness about 100–125 km, a common datum beneath that layer, of order 150–200 km below sea level, is a defensible choice provided the model has no substantial lateral density anomaly there. These are model choices, not measured universal boundaries. In a half-space cooling model, use the oldest age to choose a practical cutoff: makes the remaining temperature anomaly less than about 0.5% of its surface contrast. For and million years, this is about 225 km. A half-space model has no mathematically sharp finite thermal base.
Finally, equal column mass is the local, approximately constant-gravitational acceleration representation of equal pressure at the datum. A weak asthenosphere permits pressure relaxation through slow solid-state material creep; it need not be liquid. Short-wavelength loads, lithospheric flexure or dynamic topography can violate local column-by-column isostasy, requiring a more complete mechanical model.
Specify what “onset” means. I would distinguish the first local, perhaps short-lived subduction event from persistent, self-sustaining recycling in a connected global network of tectonic plates. Modern plate tectonics combines relatively rigid moving lithosphere, localized boundary deformation, seafloor spreading and sustained convergence. A sinking drip or a single imposed descending sheet is not sufficient evidence of that global regime. Different processes may become established at different times, so a transition interval may be a better hypothesis than an instantaneous worldwide switch.
Construct competing physical models. Solve conservation of mass, momentum and energy for evolving mantle convection and lithosphere, including temperature- and stress-dependent rheology, brittle yielding, plastic deformation, melting, compositional buoyancy and relevant phase changes. Use geodynamic modelling to compare a stagnant lid, episodic overturn, drip-dominated recycling and a sustained mobile lid. In a simplified incompressible creeping-flow framework,
Here is velocity, pressure, effective viscosity, the strain-rate tensor, density, gravitational acceleration, specific heat capacity, thermal conductivity, radiogenic power per unit mass, and dissipative heating per unit volume. These equations are a framework: compressibility, latent heat and chemical transport may require additional terms. A simple yield stress law represents cohesion , friction coefficient and pore pressure , illustrating why hydration and interface weakening can change the outcome.
Vary the early-planet controls rather than assuming the answer. Explore plausible initial mantle potential temperature, radiogenic heating, crust thickness and composition, water inventory, grain-size evolution, surface temperature, and boundary strength. Hotter mantle lowers viscosity but can produce thicker buoyant basaltic crust; it does not automatically make sustained subduction easier. A cold strong plate can transmit slab pull yet resist bending and initiation. Water weakening of rocks, melt-induced weakening, inherited faults and mechanically weak sediments can change this balance. Test proposed triggers such as mantle plume impingement or gravitational instability without presupposing an established global plate-force system.
Useful diagnostics include the Rayleigh number,
with thermal expansion coefficient , layer thickness , thermal diffusivity and temperature contrast , together with surface mobility, the localization of deformation, plate longevity, melt production and material recycling. A large Rayleigh number establishes vigorous convection under suitable conditions, not necessarily a mobile lid. A convincing model should sustain one-sided subduction, connected boundary evolution and seafloor spreading after any artificial trigger is removed. Compare initiation thresholds with continuation thresholds: hysteresis can allow an established regime to persist under conditions in which it would not spontaneously start.
As an example of the type of experiment, three-dimensional models have tested mantle plume-induced subduction initiation, finding sensitivity to plate buoyancy, magmatic weakening and hydration of the interface. The scientific use of such results here is to identify a feasible mechanism and its predicted signatures, not to adopt their success as a direct measurement of an ancient onset date.[ref]
Connect a model regime to a date through independent evidence. Couple the experiments to candidate planetary cooling histories, then predict observable pressure-temperature-time paths, magmatic compositions, deformation patterns and recycling rates for each regime. If independently dated rocks constrain a thermal or compositional state, ask which regimes can reproduce that state and the observed geological associations. Repeat over uncertain initial conditions and constitutive parameters, check numerical resolution, and compare two- and three-dimensional results where geometry matters. Models that prescribe a descending slab, fix plate velocity, or insert a weak boundary may be excellent forward models of a working subduction zone but do not demonstrate spontaneous initiation.
The result is a conditional statement: persistent plate tectonics is physically feasible under a range of conditions, and dated geological observations may favor that range at a particular time. Uncertain cooling histories, several feasible regimes and uncertain rheology prevent numerical feasibility alone from dating onset. The strongest constraint comes when the same parameter range explains multiple independently dated observations and alternative regimes fail those tests.
tidal forces measure differences in gravitational acceleration across an extended region. Neighboring freely falling particles each have zero four-acceleration, but their separation need not accelerate uniformly. Let be the tangent to a reference timelike geodesic and its infinitesimal separation from a neighboring one, chosen from an affinely parametrized geodesic congruence. Using the Riemann curvature tensor convention specified in Solution 3,
Changing the overall curvature convention changes this displayed sign. The covariant derivatives compare vectors by parallel transport, avoiding a coordinate-dependent subtraction of accelerations.
In the Newtonian limit, expand the force field at two nearby points:
For , a local orthonormal frame aligned with the radius gives
There is radial stretching and transverse compression; the trace vanishes outside the mass. Two freely falling test masses and a measurement of their changing separation implement gravity gradiometry. Their common acceleration drops out, while a relative acceleration proportional to separation remains. A local inertial frame removes Christoffel symbols at a point, but cannot remove nonzero curvature over a finite region.
A tidal force is a difference of gravitational acceleration across an extended region. In Newtonian gravity, . Two neighboring freely falling particles separated by have, to first order,
For a point mass, . In a local orthonormal radial/transverse frame,
A small freely falling cloud stretches radially and compresses transversely. Outside the source the three eigenvalues sum to zero, consistently with ; within matter the trace is . These are accelerations of separations, not the much larger common acceleration of the cloud's centre.
Figure 1.
Original local tidal map: a freely falling cloud stretches along the radial direction and contracts in both transverse directions
. Details.
In general relativity, freely falling test particles follow timelike geodesics. A freely falling frame can eliminate the connection at one event, but not the Riemann curvature tensor throughout a finite laboratory. Let be the geodesic tangent and a connecting vector for neighboring geodesics with . In the explicitly chosen curvature convention
the geodesic deviation equation is
It follows by commuting , using and zero torsion. In a parallel-propagated orthonormal frame, its spatial components directly describe measured relative acceleration. The weak-field, slow-motion limit recovers the Newtonian Hessian above. Reversing the curvature sign convention reverses the displayed curvature side, without changing the physical stretching. The equivalence principle removes locally uniform gravity, while geodesic deviation measures the residual tidal field. The connection between parallel transport, curvature and neighboring geodesics is developed systematically in the standard geometric construction.[ref]