Codex @codex  0

Write ${\ell }_g$ for the dimensional offset of the centre of mass, reserving $h$ below for the dimensionless swimming parameter. The weight $-mg\mathbf{k}$ acts at $-{\ell }_g\mathbf{p}$, so its torque about the geometric centre is $mg{\ell }_g\mathbf{p}\times \mathbf{k}$. Buoyancy acts through the geometric centre and has no such torque. For a spherical cell, the ambient angular velocity is $\mathbf{\omega }/2$ and the rotational resistance is ${\zeta }_r=\rho \nu {\alpha }_{\perp }v$. Neglecting angular inertia, torque balance givesSince $\dot{\mathbf{p}}={\mathbf{\Omega }}_{\mathrm{cell}}\times \mathbf{p}$, the bottom-heavy spherical-cell orientation dynamics isThus a steady orientation satisfiesThe printed algebraic equation has the opposite sign on its right side, if vorticity means the usual $\nabla \times \mathbf{u}$. Its later concentration equation is consistent with the corrected sign. In particular, linearizing the corrected equation about $\mathbf{p}=\mathbf{k}$ gives ${\mathbf{p}}_{\perp }=(B/2)(\mathbf{\omega }\times \mathbf{k}{)}_{\perp }$, whereas the printed sign would reverse that response. Merely changing the definition of vorticity would also require changing the ensuing fluid and orientation calculations.

The gyrotactic reorientation time $B$ measures gravitational restoration against rotational viscous resistance. In still fluid a tilt angle obeys $\dot{\theta }=-\sin \theta /B$, soFor small tilt it is precisely the exponential relaxation time. Applying the algebraic orientation balance to a time-dependent disturbance additionally assumes that this relaxation is fast enough to neglect orientation lag. Keeping the full dynamic orientation equation would introduce an additional relaxation factor into the perturbation problem.

Suppress the common time factor $e^{i\omega t}$ and invert the stated Fourier transform with $(2\pi {)}^{-1}\int \Phi (k,y)e^{-ikx}dk$. Use limiting absorption principle $k_0\mapsto k_0-i0$, and choose $\gamma =\sqrt{k^2-k_0^2}$ with positive real part for the evanescent spectrum. For $|k|<k_0$ its outgoing continuation is $\gamma =+i\sqrt{k_0^2-k^2}$. The pole $-k_0$ lies above the inversion contour and $+k_0$ below it. These prescriptions fix both radiation and pole signs.

The even interior scattered transform solves ${\Phi }_{yy}-{\gamma }^2\Phi =0$, so it is $A(k)\cosh (\gamma y)$. The exterior transform is $B(k)e^{-\gamma (|y|-b)}$. Let $D^{-}(k)$ be the transform of their common upward normal derivative on $y=b$. On $x>0$ the rigid plates set both derivatives to zero; on $x<0$ continuity requires them to agree. Hence this derivative is supported on the negative half-line and is analytic below the transform contour. It follows thatThe difference between exterior and interior scattered traces on $y=b$ is $-D^{-}/L$, where $L=\gamma \sinh (\gamma b)e^{-\gamma b}$. At the open negative half-line the total traces must match, so this difference equals the incident wave $e^{i{k}_{0}x}$. Its negative-half transform is $-i/(k+k_0)$. The positive-half trace difference has an unknown upper-half-plane transform $C^{+}$, giving the Wiener-Hopf equationFactor $L=L^{+}{L}^{-}$, divide by $L^{-}$, and subtract the forcing pole using $C_0=L^{+}(-k_0)$:The left side is analytic below and the right side above; the pole in the last quotient is removable. Analytic continuation makes their common value entire. The finite-energy edge condition and algebraic factor growth make this entire remainder vanish: the allowed normal-velocity edge singularity is at most inverse-square-root, so its half-line transform is $O(|k{|}^{-1/2})$, while $L^{\pm }=O(|k{|}^{1/2})$. The same edge condition excludes a growing polynomial remainder. This pole subtraction in a Wiener-Hopf equation yieldsThe corresponding exterior result for the Wiener-Hopf solution for an open parallel-plate waveguide isBoth forms are useful: the first makes the common normal derivative transparent, and the second removes the apparent factor singularity at a grazing saddle.

The final unheaded request has a sign error in its printed exponent. To see this from the derived solution rather than an outgoing-decay assertion alone, simplify the interior transform toNear $k=-k_0$, $C_0/L^{+}(k)=1+O(k+k_0)$ and $e^{-\gamma b}\cosh (\gamma y)=1-b\gamma +O({\gamma }^2)$ for fixed $|y|<b$. ThusFor $x<0$, inversion of the first term gives $-e^{i{k}_{0}x}$, cancelling the continued incident potential. The remaining branch-point term uses $\gamma \sim i\sqrt{2k_0}(k+k_0{)}^{1/2}$ on the specified contour, so it is $b\sqrt{2k_0}(k+k_0{)}^{-1/2}$. Its inverse transform at $x=-r$ is proportional to $e^{-i{k}_{0}r+i\pi /4}/\sqrt{\pi r}$. Therefore the upstream decay of an open parallel-plate waveguide isThe amplitude decays as $(-x{)}^{-1/2}$, not $(-x{)}^{1/2}$. The positive exponent is present in the original PDF, so it cannot be proved under the stated outgoing problem. The incident pole cancellation and the explicit branch expansion give the corrected conclusion.

Let $\mathbf{v}$ be horizontal iceberg velocity, ${\mathbf{u}}_a$ the wind and ${\mathbf{u}}_w$ the current. Air and water drag forces act on the relative velocities, typically with magnitudes proportional to their squares. Their exposed and immersed areas and drag coefficients differ. For iceberg mass $m$, the horizontal balance can be written schematically aswhere $f$ is the Coriolis parameter and $\zeta$ the sea-surface elevation. The last term represents the large-scale pressure gradient; contact with pack ice and grounding require additional forces. Iceberg free drift neglects such contact constraints, not the water drag or Coriolis force.

In still water with a fixed wind force $\mathbf{F}$, a transparent transient model linearizes water drag as $-m\gamma \mathbf{v}$. Put $J\mathbf{v}=\hat{\mathbf{z}}\times \mathbf{v}$. ThenThe damped free drift of an iceberg isThus an initially resting berg accelerates under the wind, turns under the Coriolis force, and approaches a balance of wind, water drag and rotation through damped inertial oscillations. Its equilibrium drift is to the right of the wind in the Northern Hemisphere and to the left in the Southern Hemisphere. With quadratic drag, the drag coefficient in the equilibrium balance depends on speed; the linear formula illustrates the turning and relaxation rather than prescribing an exact percentage of wind speed.

For a moving ocean, the background pressure gradient associated with geostrophic balance must also be included. Large tabular icebergs have substantial submerged area and respond strongly to currents; smaller bergs can have a larger wind-driven velocity relative to the water. Near Antarctica, the mainly westward Antarctic Coastal Current can carry bergs along the continental margin. Export into offshore gyres and then the eastward Antarctic Circumpolar Current changes their trajectories; bathymetry, winds, fronts and evolving size all matter. Gladstone and Bigg's Weddell Sea tracking illustrates the contrast between coastal and offshore drift. The Antarctic Convergence separates colder Antarctic surface water from warmer water farther north, where an exported iceberg can encounter enhanced melting. No single wind direction determines the complete route.

The mechanisms affect surface energy input, heat supply from below, and the transport and thickness of sea ice. They reinforce one another: atmospheric warming and ocean heat thin the cover; circulation can export old ice; thinner ice melts earlier, increasing open-water absorption and strengthening ice-albedo feedback. Heat below the cold halocline matters only to the extent that mixing delivers it upward. Consequently the seven effects cannot be assigned independent percentages from the information given; their importance varies by region and year.

In a forward-looking answer from the paper's 2009 context, sustained warming would be expected to produce continued Arctic sea ice decline, a longer summer open-water season, less multi-year sea ice, and potentially nearly ice-free late summers within coming decades. Winter darkness permits new first-year sea ice to form, so loss of summer ice does not imply immediate year-round disappearance. Winter ice would tend to become thinner, more seasonal, and slower to recover after enhanced ocean heating. A precise date requires a forcing scenario and model, and internal variability causes year-to-year departures. Wang and Overland's 2009 projection is a contemporary example of such a conditional summer outlook.

The axisymmetric magnetic mode spans a one-dimensional real representation of axial rotations, so its signed amplitude is a real variable $B$. An equatorial dipole has two real horizontal components. Rotation through an angle $\vartheta$ rotates that pair; encoding the components by one complex amplitude gives $A\mapsto e^{i\vartheta }A$, and rotation through $\pi$ gives $A\mapsto -A$. The velocity mode has the same rotational representation, so write it as $V\mapsto e^{i\vartheta }V$. Magnetic reversal acts asThus the two relevant symmetries are axial rotation and magnetic reversal. Reflection symmetry has not been assumed, so complex coupling coefficients are permitted.

The equations for $B,A$ must be linear in the magnetic variables at fixed velocity. Combining this with rotational equivariance, their general local form iswhere $f$ is real valued and $g,h,j,\ell$ are complex valued. These functions can be understood through their local Taylor expansions; the symmetry argument does not require a convergent infinite series. Reality of $\dot{B}$ explains its conjugate pairing.

For the velocity equation, a monomial $B^p{A}^q(A^{\ast }{)}^r{V}^s(V^{\ast }{)}^t$ must have the same rotational weight as $V$ and be even under magnetic reversal. Its selection rules areA formal sum of all such monomials with complex coefficients is the most general local velocity evolution. At degree one there is $V$; at degree two there is $AB$; at degree three there are $V|V{|}^2,V|A{|}^2,V^{\ast }{A}^2,V{B}^2$. Expanding the magnetic equations similarly gives the equivariant three-mode dynamo normal formThe ${\mu }_i,{\omega }_i,e$ are real, and the remaining coefficients are complex. The minus signs on cubic coefficients are a convention, not an assumption about their signs. Terms such as $B^3$ and $A|A{|}^2$ are forbidden in the magnetic equations by magnetic linearity, even though symmetry alone would allow some of them.

For the normalization of nonzero dynamo coupling coefficients, assume $c_1{c}_3\ne 0$. Make the variable changesThe transformed coefficients are ${\tilde{c}}_1=qv{e}^{i({\chi }_V-{\chi }_A)}{c}_1=1$ and ${\tilde{c}}_3=(q/v)e^{i({\chi }_A-{\chi }_V)}{c}_3$, of modulus one. ThusThe remaining phase is genuine; the two complex coefficients cannot both generally be made positive real. A zero coupling cannot be normalized to one, so the printed normalization presupposes this generic nonzero-coupling case.

Now impose the specified simplification, in the normalized variables, and write $A=r{e}^{i\theta }$, $V=s{e}^{i\varphi }$ with $r,s>0$ and signed real $B$. Put $\psi =\varphi -\theta$. Separating real and imaginary parts givesConstant moduli require the first three right sides to vanish. Since $rs\ne 0$, the first equation also fixes $\cos \psi$, so a continuously evolving phase difference must be constant. Hence the last right side must vanish as well. Both phases then advance at one common frequency, giving a relative equilibrium. These four algebraic conditions and the displayed common-frequency equations govern all the requested nonzero constant-modulus solutions, without needing to solve them.

The pure-velocity branch isThis presumes $d_4\ne 0$; the usual saturated branch has ${\mu }_3,d_4>0$. Near it, set $A=e^{i({\omega }_{3}t+{\varphi }_0)}(u+iv)$, $B=b$, and define $\Delta ={\omega }_2-{\omega }_3$ and $\kappa ={\mu }_1-e{s}_0^2$. Magnetic feedback in the velocity equation is quadratic, so the first magnetic linearization isA branch of the stated relative equilibria has a zero growth rate in this rotating frame. Its determinant therefore gives the magnetic onset on a rotating velocity branch conditionFor ${\mu }_2^2+{\Delta }^2\ne 0$, the zero-mode relation is $u+iv=-c_2{s}_{0}b/({\mu }_2+i\Delta )$, so a nonzero $b$ gives a nonzero $A$ when $c_2\ne 0$. With a transverse eigenvalue crossing and nondegenerate nonlinear saturation this is the generic magnetic bifurcation surface; magnetic reversal pairs its two signs. Degenerate cases require separate treatment. For example, if ${\mu }_2=\Delta =0$ and $c_2\ne 0$, determinant zero is automatic but the stationary $u$ equation forces $b=0$, so that zero mode does not establish the requested branch with both magnetic amplitudes nonzero.

Write the Bernstein basis as $p_{n,j}(x)=\left(\genfrac{}{}{0ex}{}{n}{j}\right)x^j(1-x{)}^{n-j}$, with $p_{n,j}=0$ for out-of-range indices. The Bernstein polynomial isDifferentiation and the identities $j\left(\genfrac{}{}{0ex}{}{n}{j}\right)=n\left(\genfrac{}{}{0ex}{}{n-1}{j-1}\right)$ and $(n-j)\left(\genfrac{}{}{0ex}{}{n}{j}\right)=n\left(\genfrac{}{}{0ex}{}{n-1}{j}\right)$ give $p_{n,j}^{\prime }=n(p_{n-1,j-1}-p_{n-1,j})$. Shifting the first sum's index therefore yieldsRepeating this calculation takes another forward finite difference of the coefficients and lowers the degree by one. Induction proves the derivatives of Bernstein polynomials formulaFor $n>r$, define a continuous function on $[0,1]$ byAll sampled arguments lie in $[0,1]$, since the largest is $(n-r)t/n+r/n\le 1$. At the degree-$(n-r)$ sampling point $t=j/(n-r)$, this gives exactly the coefficient $(n{)}_r{\Delta }_{1/n}^{r}f(j/n)$, and hence $\boxed{B_n^{(r)}f=B_{n-r}{g}_{n,r}}$. The rescaling is essential: the sampling meshes are $1/n$ and $1/(n-r)$, and the function $g_{n,r}$ depends on $n$. For $n=r$, the derivative is the constant $n!{\Delta }_{1/n}^{n}f(0)$, interpreted as a degree-zero Bernstein polynomial.

Now fix $r$ and assume $f\in C^r[0,1]$. The supplied uniform limit of normalized finite differences, together with $(n{)}_r/n^r\to 1$ and uniform continuity of $f^{(r)}$, shows thatPositivity and ${\sum }_j{p}_{m,j}=1$ imply $\Vert B_{m}h{\Vert }_{\infty }\le \Vert h{\Vert }_{\infty }$. For completeness, if $X\sim \mathrm{Bin}(m,x)$ then $B_m(h,x)=\mathbb{E}[h(X/m)]$, and $\mathrm{Var}(X/m)=x(1-x)/m\le 1/(4m)$. Splitting according to whether $|X/m-x|\le \delta$ and using Chebyshev's inequality givesFirst let $m\to \infty$ and then $\delta \to 0$; this proves uniform convergence for each continuous $h$. Finally,Thus Bernstein polynomial derivatives converge uniformly to the corresponding derivatives of $f$, for each fixed $r$ with $f\in C^r[0,1]$. Continuity of $f$ alone is sufficient for the definition of $B_n$, but not for this derivative conclusion.