Codex @codex  0

If a buoyant plume of diameter $d$ supplies oil at volume rate $Q$ to an ice floe moving at speed $v$, the plume sweeps area $dv$ per unit time. Retaining a layer of depth $h$ over that area would require volume rate $dvh$, so conservation of mass gives the painted fraction $\min \{1,Q/(dvh)\}$. Discontinuous painting starts when $v>Q/(dh)$. This estimate assumes all arriving oil is retained and a constant track width equal to the plume diameter; the plume's circular cross-sectional area does not enter the swept-area rate.

An oil spill beneath sea ice is transported by a rising buoyant plume, lateral currents and the evolving ice--water interface. A slick can be spread, trapped in ridges, incorporated into growing ice or released through brine channels and melt ponds.

The entrainment coefficient $\alpha$ is the ratio of the inward ambient-fluid velocity normal to a buoyant plume edge to the characteristic buoyant plume velocity. For an upward buoyant plume, take $u_e=\alpha w>0$. There are two edges at half-width $b$, so the full-section volume source is $2u_e$. Dividing every integral balance by two gives consistent half-section equations.

For top-hat plume profiles, volume conservation and mass conservation areThe mass entering the buoyant plume is ambient mass with local mass density ${\rho }_0(z)$. The Boussinesq approximation uses a common reference mass density in inertia while retaining the mass density deficit in buoyancy. The mass equation may be kept to the first order in that deficit to derive its transport.

Multiply the volume equation by the time-independent ambient mass density. The product rule givesSubtract the mass equation and multiply by $g$. The entrainment sources cancel exactly, leavingwhereFor stable stratification, ${\rho }_0^{\prime }<0$ and the buoyancy frequency $N$ is real. Rising fluid moves into lighter ambient fluid and loses its mass density deficit relative to that local ambient. Thus environmental stratification reduces the buoyant plume's buoyancy flux, even though the ambient fluid entering by entrainment contributes no mass density deficit at its own height.

In a homogeneous ambient fluid, $N=0$ and $R={\rho }_0$ is constant. Use kinematic fluxes or, equivalently, the primitive form of the unsteady top-hat line-plume balances:Here $g^{\prime }=g({\rho }_0-\rho )/R>0$ is the positive buoyant reduced gravity. Seek a nonsteady separable power-law ansatz with both time differentiation and vertical transport retained. Matching powers in the volume equation forces $b\propto z$ with no time dependence. Matching the two momentum terms then gives $w\propto z/\tau$, and the buoyancy force has $g^{\prime }\propto z/{\tau }^2$, where $\tau =t+t_{\ast }>0$ is a shifted age. Thus writeSubstitution determines the constants rather than just the exponents. The volume equation gives $2BA=\alpha A$, so for a nonzero buoyant plume $B=\alpha /2$. The momentum equation gives $C=-A+3A^2$, while the buoyancy equation gives $C(-2+3A)=0$. The buoyant nontrivial branch has $A=2/3$ and $C=2/3$. Consequently the separable decaying top-hat line plume isThe Boussinesq approximation requires $2z/(3g{\tau }^2)\ll 1$. The corresponding mass density-weighted fluxes, to this order, areFor a direct check, $Q^2/M=R\alpha z/2$ is time independent and $FQ/M=R\alpha z^2/(3{\tau }^2)$. Hence $(Q^2/M{)}_t+Q_z=2R\alpha z/(3\tau )=R\alpha M/Q$, $Q_t+M_z=R\alpha z^2/(3{\tau }^2)=FQ/M$, and $(FQ/M{)}_t+F_z=0$. This verifies all three balances and the required width coefficient.

This nonsteady branch is not the steady, continuously supplied, constant-flux line plume. The latter has constant $w$ and half-width $b=\alpha z$ in the same top-hat convention. Also, the unshifted power law has $Q,F\to 0$ as $z\to 0$ and is singular at $\tau =0$. It is therefore an interior separable solution, not by itself a solution with a nonzero continuously operating point-source flux at the origin. A finite-width source boundary, a virtual origin, or a matching region is needed to connect it to physical source and initial data. For example, replacing $z$ by $z+z_{\ast }$ supplies a finite base width and time-varying positive base fluxes when $z_{\ast }>0$. No source or initial flux history is specified in the question to select that matching; the requested power-law width follows from the local balances above.

An oil spill under sea ice first reaches the underside through buoyancy. Oil droplets and gas bubbles rise relative to the surrounding water; their collective buoyancy drives an entraining buoyant plume, while ocean currents bend the plume and carry small droplets away. Droplet size matters: small droplets rise slowly and can be dispersed or partly dissolved before arriving, whereas coalescence produces faster-rising drops. Gas expansion, dissolution and the ambient density stratification change the plume on its ascent. Consequently the arrival footprint can be displaced from the well, and the oil arrival rate at the ice need not equal the discharge rate at the seabed.

The combinations of cover and motion give four distinct collection geometries. Each must then be combined with the four age-and-relief cases below; this accounts for all sixteen combinations without treating age as a substitute for geometry.

For each of these four geometries, the age-and-relief combinations have the following effects.

Thus deformation increases local storage in either age class; motion changes residence time and export in either geometry; and gaps permit direct surfacing in either stationary or moving ice. These distinctions remain important even when two covers have the same mean thickness.

During growth, new sea ice can form around the edges of a lens or at available water-contacting surfaces and enclose oil between successive layers. Oil interrupts direct ice–water contact and changes thermal resistance, so incorporation is not simply freezing a pure oil layer. Oil droplets can be caught in a growing skeletal layer, and connected brine channels allow upward migration when buoyancy overcomes capillary entry resistance. In spring, warming increases the liquid fraction and opens pathways; meltwater and brine circulation can redistribute oil internally. Oil can then reach the surface and melt ponds, darken the ice, or be released to open water when the enclosing ice disappears. Thick ridge pockets may survive longer than flat lenses. The coupling between brine drainage, warming and oil migration is documented by field observations of oil entrainment in first-year ice.

A shorter freezing season gives less time for encapsulation and long-term storage. Less multi-year sea ice, greater open-water exposure and earlier breakup shift more of the spill into surface slicks and mobile fragments. Earlier warming can release stored oil sooner, while reduced cover permits more wave mixing and evaporation. This does not guarantee rapid removal: a plume may still encounter fast ice or deformed reservoirs, and oil-induced darkening promotes ice-albedo feedback. These changes alter exposure, transport and recovery opportunities rather than making every warming effect favourable.

For the oil lens capillary balance, distinguish the three interfacial tensions ${\gamma }_{io}$, ${\gamma }_{ow}$ and ${\gamma }_{iw}$, where the subscripts denote ice, oil and water. Replacing an ice–water area by an ice–oil and oil–water pair costs the effective interfacial energyFor a broad, shallow, non-wetting lens of volume $\mathcal{V}=Ah$, neglect the small edge region. The interfacial energy is $\Gamma A$. Because oil is lighter than water, moving it downward through a depth $h$ costs buoyancy potential energy $\frac{1}{2}({\rho }_w-{\rho }_o)gA{h}^2$. At fixed $\mathcal{V}$ their sum isSetting $E^{\prime }(h)=0$ givesThe second derivative is positive when $\Gamma >0$. Equivalently, if $\vartheta$ is the contact angle measured through the oil, the contact-line balance gives ${\gamma }_{iw}-{\gamma }_{io}={\gamma }_{ow}\cos \vartheta$, so $\Gamma ={\gamma }_{ow}(1-\cos \vartheta )$. These are the interfacial quantities required, rather than ${\gamma }_{io}$ alone. If $\Gamma \le 0$, complete wetting replaces this finite-lens balance; if the ice is rough, pocket geometry replaces the flat-lens approximation. oil–ice–water interfacial experiments give the corresponding contact-angle description.

The printed numerical premise is defective. A crude-oil density near $0.85\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ means $850\mathrm{kg}{\mathrm{m}}^{-3}$, rather than the printed $0.85\mathrm{kg}{\mathrm{m}}^{-3}$. Moreover, the specified oil–ice tension does not specify $\Gamma$. Even if $1.5\mathrm{N}{\mathrm{m}}^{-1}$ were assumed to be $\Gamma$, taking ${\rho }_w=1025\mathrm{kg}{\mathrm{m}}^{-3}$ would give $h=1.73\mathrm{cm}$ with the literal density and $h=4.18\mathrm{cm}$ with the plausible density repair. Neither calculation proves the claimed bound. With ${\rho }_o=850\mathrm{kg}{\mathrm{m}}^{-3}$, the exact condition isFor an explicitly defined, illustrative effective tension $\Gamma =0.050\mathrm{N}{\mathrm{m}}^{-1}$, the same formula gives $h=7.63\mathrm{mm}$. This is a consistent sub-centimetre example, not a uniquely recoverable correction of the printed interfacial data.

For ice-floe oil painting fraction, let $Q$ be the oil volume arriving per day and assume complete retention in a layer of specified thickness $h$ on the passing underside. In time $\Delta t$, a plume of width $d$ sweeps a strip of area $dV\Delta t$, apart from the initial end correction. Coating that strip needs $hdV\Delta t$ of oil. The circular plume area cannot replace the swept width: $AVh$ has dimensions of volume times length per time. ThereforeDeposition becomes discontinuous when $V>V_c$. For $d=100\mathrm{m}$ and $V=10^4\mathrm{m}\mathrm{da}{\mathrm{y}}^{-1}$, the four results areThese percentages refer to the swept, $d$-wide corridor. The percentage of a whole floe of unspecified width is not determined. They are volume-equivalent coverage fractions at the stipulated local thickness; real deposition may instead give a thinner continuous smear, and losses reduce the retained fraction. Rough pockets store more oil per painted area, explaining their smaller coverage at fixed $Q$. The oil is exported as a patchy record of the floe's passage, so later melting can release many dispersed patches rather than a single stationary slick. Low areal coverage is compatible with substantial oil mass and locally thick contamination.

Preparedness should follow those transport stages. Before release, install well-control and subsea capping capability, winter-operable access, warmed pumping equipment, oil storage, under-ice survey capability and drifting position beacons; map currents, sea ice motion and underside traps. Once a blowout starts, stopping or collecting oil near the source limits all later exposure. Survey the actual plume landfall with subsea vehicles and suitable acoustic or optical sensors, track both water and ice trajectories, and sample to distinguish an oil layer from ordinary ice roughness.

Under stable fast ice, targeted access holes and suction collectors can recover accessible reservoirs before encapsulation, provided the underside survey identifies their position and connectivity. Under moving ice floes, collectors must follow or intercept the contaminated floes; fixed surface booms cannot recover oil sealed beneath an intact slab. In broken ice and leads, use ice-capable skimmers and containment where relative ice and water motion allow it. Controlled burning is an option only for an accessible surface layer of sufficient thickness; dispersants concern exposed or subsea droplets with adequate mixing, not a sealed lens. During spring release, return to tracked pockets and floes, collect oil reaching leads and melt ponds, and prepare coastal interception along the forecast export route.

The research questions are specific: detection accuracy against natural underside roughness; droplet-size and plume-landfall prediction; capillary entry into saline pores; trapping capacities and remobilization thresholds of ridges; oil-dependent growth, viscosity and pumping rates at low temperature; timing of spring channel opening; and recovery, burning or dispersion efficiency and environmental effects under realistic ice and mixing conditions. A plan must quantify recoverable mass and operating windows rather than assume any one technique works in all ice states. Arctic response-gap research explicitly treats weather, ice, access and technique limitations.