Wave Patterns of Gravity-Capillary Waves from Moving Localized Sources

We study wave patterns of gravity-capillary waves from moving localized sources within the classic setup of the problem of ship wakes. The focus is made on the co-existence of two wave systems with opposite signatures of group velocity relative to the localized source. It leads to the problem of choice of signs for phase functions of the gravity (“slow”) and capillary (“fast”) branches of the dispersion relation: the question generally ignored when constructing phase patterns of the solutions. We detail characteristic angles of the wake patterns: (i) angle of demarcation of gravity and capillary waves – “the phase Mach” cone, (ii) angle of the minimal group velocity of gravity-capillary waves – “the group Mach” cone, (iii, iv) angles of cusps of isophases that appear after a threshold current speed. The outer cusp cone is naturally associated with the classic cone of Kelvin for pure gravity waves. The inner one results from the effect of capillarity and tends to the “group Mach” pattern at high speeds of current. Amplitudes of the wave patterns are estimated within the recently proposed approach of reference functions for the problem of propagation of packets of linear dispersive waves. The effect of shape is discussed for elliptic reference sources.


Introduction
The problem of wavefield behind a moving source is well-known in various domains of modern physics and, in particular, in the dynamics of the atmosphere and ocean [1]. The brightest and deepest for understanding rich physics is the example of ship waves, first considered in the works of Lord Kelvin [2][3][4][5]. At the turn of the 20th century, Kelvin and his followers created the basis for many works developing the classical results (e.g., [6,7]).
Patterns of phases of linear waves can be constructed in a relatively simple and general way [8]. In what follows, we assume a motionless source and the flow velocity U = (U, 0, 0) directed along the x−axis in the positive direction. We restrict ourselves to the case of two-dimensional waves, which include surface water waves, Rossby waves, and internal gravity waves in a depth-limited ocean (atmosphere). Stationarity condition in the reference frame associated with the source ω(k) + kU = 0 (1) allows one to determine the transverse component of the wave vector k = (k, l) by specifying the longitudinal component k, which is usually associated with the size of the source. The phase ϕ = kx at an arbitrary point r = (x, y) is determined from the condition of stationarity (1), i.e., for the conditional extremum of the functionφ = kx − t (ω(k) + kU) The Lagrange multiplier t in (2) can be interpreted as the propagation time of a wave packet with a constant group velocity [9], since the conditional extremum of the phase ϕ is realized at The isophases ϕ = const provide a visual picture of the distribution of the wave troughs and crests in a stationary wake created by the source. Note that the straightforward implementation of this approach can lead to non-physical results when the wave pattern appears where it should not be. The reason for this may be the arbitrariness of the choice of the sign of the phase ϕ. Formally, ϕ + = kx and ϕ − = −kx are equal in rights. Usually, they refer to the principle of causality for choosing the sign of the phase and linking it with the positiveness of the time t.
The situation becomes somewhat more complicated when the stationarity condition (1) has two (or more) solutions. Lord Rayleigh [10] considered one of such cases for standing gravity-capillary waves on the surface of running water. For the flow speed higher than the minimum wave velocity C m = 2g/σ ≈ 23 cm · s −1 (phase and group speeds are the same) the Equation (1) has two roots corresponding to "slow" gravitational and "fast" capillary waves. The group velocities of these waves relative to a stationary source have opposite signatures: the energy of capillary waves runs upstream, while the energy of the gravitational waves lags behind. Rayleigh [10] accounted for this feature in the one-dimensional problem while discussed the solution for the two-dimensional case "by analogy".
The group velocity of waves is considered as an essential quantity of energy transfer since works of Hamilton [11]. It does not enter explicitly into the well-known expressions for the classic Kelvin's solution of pure gravity ship waves [3,4]. The Lighthill [8] interpretation of the Kelvin angle operates with a fixed ratio of 1/2, thus emphasizing the physical difference of the phase and group propagation. This fact alone results in the occurrence of the universal stationary wake patterns in the form of arcs of chevron.
The dispersion of pure capillary waves (ω(k) ∼ |k| 3/2 ) is positive, i.e., group velocity is larger than the phase one, and the corresponding wave patterns are not so peculiar as the Kelvin wake ones with their two systems of transverse and divergent waves. In a sense, the wave patterns look trivial. This visual simplicity disappears when passing to gravity-capillary waves where the above wave patterns represent asymptotics at very short and very long waves. The "intermediate" (in words of [12]) range of gravity-capillary waves turns out to be rich in non-trivial kinematic and dynamical effects. Two characteristic scales, minima of phase and group velocities, provide the dimensionless quantities of the Mach numbers. The "phase Mach number" U * = U/C m demarcates gravity and capillary branches of the ship wake. At relatively high U * 1.938 singularities of the phase patterns of the ship wake arise as cusps resembling the well-known gussets of Kelvin. A new family of cusps emerges as a feedback of the wave speed minimum, that is, as a result of the coupling of gravity and capillarity effects. The "group Mach number" U/C g,min (C g,min is minimal group velocity of gravity-capillary waves) forms a cone where diverging capillary and gravity waves cannot enter. This Mach cone is observed in laboratory and field experiments, e.g., ([13]).
The reference solution approach [9,14] resolves the problem of assessment of wake amplitudes under quite mild physical assumptions including wave caustics. At loci of the group velocity minima this method predicts the effect of quasi-dispersion [14], i.e., slow decay of wave amplitude with distance from the source. The observability of this effect is questionable under typical scales of wake but still worth to be mentioned in the context of the previously used asymptotic approaches (e.g., [13,15]) and the effect of the "group Mach" [13].
In contrast to the "point-wise" estimates of wave amplitudes with conventional asymptotic approaches (e.g., the stationary phase method), the reference solution method describes an evolution of a finite perturbation of elliptic shape. Thus, the effect of the shape on the wake patterns can be discussed with explicit formulas in hand.
In Section 2 the mathematical and physical background for the problem of gravity-capillary waves on running water is presented. The corresponding solutions will be formulated in Section 3 for phase patterns and asymptotic dependencies of wave amplitudes. Discussion of Section 4 is devoted to general aspects of the problem of linear wave propagation and their implications in geophysical hydrodynamics.

Ray Equations and Isophase Patterns
This part follows the well-known results of ship wave kinematics used by Thomson [4] himself and his followers (e.g., [1,8]). We consider a point-wise source of perturbation of a stationary current with constant speed U directed along the axis x in positive direction. The source generates gravity-capillary waves which frequencies and wavenumbers obey the well-known dispersion relation [16] with the Doppler shift taken into account Gravity acceleration g = 981 cm · s −2 and kinematic water surface tension (δ = γ/ρ ≈ 74 cm 3 · s −2 γ ≈ 74 dyn · cm −1 , ρ = 1 g · cm −3 is water density) provide spatial and temporal scales of the problem L = δ 1/2 g −1/2 = 0.275 cm; T = δ 1/4 g −3/4 = 0.0167 s (5) associated with the wavelength λ m = 2πL = 1.73 cm and wave period T m = √ 2π/T = 0.074 s of minimal wave velocity C m = √ 2(gδ) 1/4 = √ 2L/T = 23.2 cm · s −1 [17]. The dimensionless wavenumber κ and frequency σ are defined as The corresponding dispersion relation (4) becomes with dimensionless velocity (phase Mach number) and angle θ counted clockwise from negative direction of axis x. The condition of stationarity (1) can be parameterized by θ as follows Dimensionless coordinates as functions of dimensionless time τ ξ = x/L; η = y/L; τ = t/T of the stationary phase can be calculated from the ray Equations (3) as Define the phase function ϕ ± = ±(kx + ly) = ± κ x ξ + κ y η and write it accounting for the condition of stationarity (1) Trajectories (11) and (12) can be written as functions of phase ϕ ± Expressing wavenumber κ in terms of angle θ one gets from (9) two roots Each root requires special choice of signature of the phase function ϕ (13). This step is usually associated with the principle of casuality. An equivalent formulation can be proposed based on the group velocity signatures corresponding to the roots κ 1,2 . Note that the critical value κ m = 1 determines minimal value of phase velocity, i.e., merging of the roots (16) along direction θ in wavevector space. For the gravity branch (κ 1 < 1) the group velocity is directed downstream that dictates the choice of sign '−' in ϕ (13). Alternatively, the sign '+' appears for the capillary branch (κ 2 > 1).
The line |U * cos θ| = 1 demarcates two families of isophases in wavevector space. In the coordinate space (ξ, η) the merging of gravity and capillary branches occurs at the asymptotes of isophases ϕ 1,2 = const 1,2 [1] tan 2.1.2. Minimal Group Velocity and Cusps of Isophases ϕ ± = const One more characteristic scale of the wake patterns can be related to the minimum of group velocity of the gravity-capillary waves. For water waves it gives C g,min = min(C g ) = 1 2 for dimensionless wavenumber κ min = 2/ √ 3 − 1. The corresponding phase velocity and, hence, C m ≈ 28.1 cm · s −1 is just 20% higher than the minimal C m . The minimum of the group velocity is closely related to the condition of singularity of parametric curves of isophases (14) and (15) dξ It leads to the polynomial equation in our dimensionless variables (cf. [15] Equation (20)) or ( [12] Equation (16)) that can be re-written in quite transparent form as a dependency of dimensionless speed U * on wavenumber in the cusp Dependencies Z 6 and U * (cusps) are shown in Figure 1a,b. The cusps appear when U * exceeds a threshold value U * ≈ 1.938 at the corresponding dimensionless wavenumber κ thr ≈ 0.276 (water wave wavelength λ cusp ≈ 6.27 cm) that is a root of the cubic equation For a given U * , two solutions for κ cusp exist. The longer wave solution correspond to the cusp of the classic Kelvin pattern for deep water gravity waves where two system of transversal and diverging waves merge. The shorter wave solution is limited by the condition of vanishing the denominator in (21). It is easy to show that this is exactly the condition of the group velocity minimum ∂ 2 σ/∂κ 2 Below we refer to this family of cusps as capillary one as far as it results from the capillarity effect. Plugging (21) into ray Equations (14) and (15) one gets directions of wavevectors at the cusp points (κ c is the wavenumber at the cusp point) and the direction of the cusp line in the coordinate space Dependencies of angles θ cusp , β cusp are shown in Figure 1c. The upper branch of the blue curve and lower one of the red correspond to the well-known cusps of gravity waves that conform the classic Kelvin pattern. In the limit U * → ∞ the cusps are located at the straight lines with the Kelvin angle β Kelvin = 19.47 • . The wavevectors at the cusps are tending to the universal angle θ Kelvin The remaining branches in Figure 1c are tending to "trivial" values β cusp → 0 and θ cusp → 90 • . The dependency of the inner (capillary) cusp cone angle β cusp on the current speed is an intriguing issue. In the limit of high speeds U * → ∞ it is tending to the inverse Mach number constructed of the current speed and the minimal group velocity as shown in Figure 1d tan As shown below, the direction of minimal group velocity C g,min (18) is one of quasi-dispersion (51) with peculiar dependence of wave amplitudes on distance [9]. Finally, we have three characteristic angles of the gravity-capillary wake patterns: 1.
Angle β m (17) demarcates families of gravity and capillary waves at low speeds U * 1.938. After the appearance of cusps of isophases, these two families can overlap as shown below. This angle can be naturally related to the Mach number constructed of the minimal phase velocity C m ; 2.
Cusps of isophases occur at relatively high speeds U * 1.938. In addition to the classic Kelvin pattern where diverging and transversal systems merge a new family of cusps, fan waves, in the words of [15], conform a divergent-fan pattern. Loci of the cusp is a cone with one of the angle β cusp (Equation (24), when Equation (20) has two solutions), while the fan wave fronts aligned along the demarcation line (angle β m ). Asymptotically, this cone tends to the Mach one constructed of the minimal group velocity of the gravity-capillary waves C g,min . It should be stressed that this cone is due to the capillarity effect in the domain of gravity waves; 3.
At the limit U * → ∞ the cone where transversal and divergent gravity waves are merging is tending to its classical Kelvin value β K = arcsin(1/3) ≈ 19.47 • . Thus, there is a continuous transition of this pattern between pure gravity and gravity-capillary waves.

The Reference Solutions Approach
Estimates of wave amplitudes in gravity-capillary wake follow the previous authors' papers on the so-called approach of reference solutions [9,14]. Here, we fix main points of the method.
The idea of the method has been exploited many times since Kelvin [2]. The initial condition for a wave variable f (below we assume f to be free surface elevation) is set in the form i.e., as a quasi-monochromatic wave packet with a carrier wavevector k 0 and "slow" modulation function F(x). Then, a special "reference function" F(x) is chosen for effective manipulations with linear solutions (e.g., [18], for one-dimensional case). The gaussian shape of F allows one for treating (26) as a superposition of linear harmonics with the gaussian distribution of amplitudes. The temporal evolution of this distribution is governed by the effect of dispersion of harmonics, that is, by properties of particular dispersion relationship (in our case, by Equation (4)). The distribution is assumed to be narrow in wavevector space and the dispersion effect is accounted up to the quadratic term in the formally small bandwidth parameter. The quadratic form appears as an argument of the corresponding exponent for an arbitrary time. It allows one for presenting wave packet evolution both in coordinate and wavevector spaces as one of the gaussian shape with parameters of the quadratic form being dependent on time.
The approximate solution can be written as follows [14] f (x, t) = 1 where "fast" x and "slow" ζ coordinates are used to discriminate "fast" effect of wave propagation in x-space and a "slow" dispersion of the wave packet of a finite bandwidth relatively to the carrier harmonics k 0 . The "fast" amplitude factor describes the overall changes of wave packet amplitude due to its finite bandwidth and can be seen as a multi-dimensional generalization of the corresponding amplitude factor in the method of stationary phase [14]. The second differentials appear in the denominator of (27) quite similarly to the stationary phase method. The function D in the argument of exponential function (27) provides an additional wave packet blurring. The "slow" coordinate describes deviations of harmonics relative to the carrier wave harmonic k 0 that is absent in the stationary phase method. The quadratic form Λ(k 0 , ζ, t) of exponential function in (27) preserves gaussian shape of the wave train modulation. Its isolines are ellipses, the parameters of which (orientation and eccentricity) vary with time.
One can restrict our consideration by the approximation of function f in (27) f In terms of dimensionless variables, one has for parameters of dispersion Note that dimensionless bandwidth components in (36)-(38) are scaled by the the capillarity length L (5) This is in contrast with the problem of pure gravity ship waves where the scaling is arbitrary and, generally, refers to the carrier wavelength.
The determinant D can be presented in the form where and With the stationarity condition (9) wave amplitude (determinant D) can be presented as function of angle θ (or angle β = arctan (y/x)) and dimensionless U * , the "phase Mach" number. The effect of such a transformation is different for gravity and capillary family of the wake. This justifies our appeal to the limiting cases of the dispersion relation (4), i.e., to the cases of pure gravity and pure capillary waves.

Limiting Cases of Wave Patterns. The Effect of Wave Dispersion on Wave Patterns
The classic case of pure gravity waves has been studied in detail since Lord Kelvin's works [2][3][4][5], while the alternative case of pure capillarity did not attract too much attention so far. This seems to be a serious gap since the change of signature of the wave dispersion leads to fundamental differences in the geometry of the patterns and the distribution of wave amplitudes. Here, we outline very few features of the issue leaving the extended discussion for future publications.
The limiting cases can be easily derived from general dimensionless relationships (9)-(16). Scaling parameters for pure gravity and capillary waves can be chosen arbitrarily for the power-like dispersion relations where specific temporal or spatial scales are absent. Thus, one can leave C m and, correspondingly, U * (9)- (16) or propose an alternative scaling for a particular problem to be solved. Both solutions are equally acceptable for the analysis of wave kinematics.
Stationarity condition (9) provides a pair of roots κ 1,2 obeying κ 1 κ 2 = 1. Thus, the formally small or large roots at U * → ∞ describe alternative limiting cases of pure gravity (κ 1 → (2U * cos θ) −2 → 0) and capillary waves (κ 2 → (2U * cos θ) 2 → ∞). Plugging these asymptotic expressions into trajectory Equations (14) and (15), one can easily get the well-known parametric equations for isophases for pure gravity waves (e.g., [9], Equation (26)) and their pure capillary counterparts within the generic scaling of this paper (5)-(8) One should emphasize a sort of similarity of pure gravity and capillary patterns: both the well-known Kelvin gusset of gravity waves and the parabolic pattern of capillary waves in Figure 2 are independent of the dimensionless current speed U * . The speed magnitude just re-scales these patterns in both directions but does not affect their geometry. This trivial conclusion results from the above mentioned fact of disappearance of characteristic scale (e.g., C m ) when passing to limiting cases of long (gravity) or short (capillary) waves.
The qualitative dissimilarity of the patterns in Figure 2 results, evidently, from different signatures of wave dispersion: negative for gravity and positive for capillary waves. For power-law dispersion in unmoving media ν = ακ n cf. (7) with a constant α and positive n the ratio of group and phase velocity is constant. At r < 1 (negative dispersion, case of gravity waves when energy propagates slower than phase) the Lighthill [8] illustration for the Kelvin cone becomes valid. For an arbitrary r < 1, the downstream cavity in the form of a cone with the corresponding opening angle is filled by waves.
In the alternative case of positive dispersion (capillary waves and waves with dispersion Equation (48) and n > 1) when energy propagates upstream the waves do not enter a downstream cavity of parabolic shape (Figure 2b). Dissimilarity in behavior of wave amplitudes requires a special discussion. Here, we just note that an additional physical scale, dimension of the source, in expressions (40)-(43) introduces an additional dimensionless parameter. The Froude number determined by the size of the source R x can be defined in the well-known way (e.g., [8]) for pure gravity waves Pure capillary waves require a special definition, say, as follows Irrespectively of the definitions (49)-(50), the amplitude dependencies on the current speed (40)-(42) are qualitatively different in the limiting cases: wake amplitude is decaying with speed for pure gravity while it is growing for the capillary counterparts. This delicate issue requires special study. The main focus of the paper is an intermediate range of gravity-capillary waves where both effects of gravity and capillarity compete on equal terms.

Anomalous Dispersion of Gravity-Capillary Ship Wake
Explicit Formulas (40)-(43) turn us back to the effect of minimum of group velocity C g,m . In the general case, A = 0 determines the amplitude decay f ∼ ϕ −1 or τ −1 that is consistent with the result for pure gravity waves (e.g., [9]). The capillarity changes this result qualitatively when A vanish at that is exactly at the wavenumber of the minimum of group velocity (cf. Equation (18)). It is interesting that the combination of the second derivatives responsible for dispersion of two-dimensional wave packet vanishes as well, i.e., (see Equation (36)) The case of quasi-dispersion has been first described in Gnevyshev and Badulin [14] for Rossby waves in the ocean with homogeneous vertical stratification (buoyancy frequency is constant). The vanishing of the second derivative Ω kk is just a special case when the omni-directional method of stationary phase predicts a slower dispersion of wave packets. Using a more general approach of the reference functions, one has an anomalously slow dispersion τ −1/2 in the two-dimensional problem.
It is easy to see that the anomalous (quasi-) dispersion occurs for the direction β group specified in the previous section (25) and for relatively high speeds U * > 2 √ 3 − 3 −1/4 ≈ 1.212. One should note that amplitudes are vanishing at the demarcation asymptote at β m where κ = 1 (17). At the line of quasi-dispersion β group (25) and (52) the amplitudes are not zero but asymptotically decaying slower than in close directions.

Results
In this section, we present illustrations of the above theoretical analysis for different parameters. The dimensionless speed U * , the Mach number related to the minimum of phase speed, determines characteristic wave patterns. Dimensionless coordinates ξ, η make distributions of wave phases and amplitudes universal, i.e., independent of the nature of fluids. The dimensionless variable wavenumber κ = L|k| can be treated as the inverse Bond number with wavelength λ as a characteristic scale. In our discussion of results, we mostly focus on an intermediate range of waves where both gravity and capillarity effects are equally important. For illustrations of wave phase patterns, we choose characteristic values of the current speed: 1.
The phase speed of the first Wilton resonance [20]; 2.
The phase speed of the minimum group velocity; 3.
The threshold value U cusps = 1.938 of the isophase cusps; 4. U * = 4 as a case when families of gravity and capillary branches essentially overlap.
When discussing wave amplitude distributions, we consider just three cases of the above list, avoiding a special case of wave caustics and superposition of different wave systems (gravity and capillary).

The First Wilton Resonance
A stationary wave pattern appears when the current speed exceeds the minimal phase speed of waves, i.e., for U * = U/C m > 1. One can start illustrations with the case of the first Wilton resonance [20]. Wavelengths of the Wilton harmonics for water waves λ Wilton1 = 2.44 cm and λ Wilton2 = 1.22 cm correspond to rather small excess of the current speed over the minimal speed C m ≈ 23 cm · s −1 , i.e., U * = 3/2 √ 2 ≈ 1.03. Figure 3a shows two systems of waves. Below, we refer to waves running upstream and shown by red isolines as a capillary family. Their downstream counterparts of gravity family are given in blue. The two families are independent at this current speed, that is, their domains are not overlapping. The demarcation line of the families, the asymptote (17), is shown in magenta. The classic Kelvin cone is given by a black solid line for reference.  The cone of "group Mach" U/C g,min is shown in cyan. It is wider than the Kelvin cone (black lines).

The Minimum of the Group Velocity
The case of minimal group velocity shows nothing new in the phase distributions. This is the threshold value of U * when the inner Mach cone emerges and the condition of anomalous wave dispersion (52) starts to be satisfied. The "group Mach" cone is still essentially larger than the Kelvin one.

Cusps of Isophases at U * (cusps) ≈ 1.938
As mentioned above, the case of the appearance of singularities of isophases is of special interest. All the characteristic cones (demarcation-magenta, "group Mach" and Kelvin's one are well separated.

Merging Wave Systems. Inner Mach's cone
The domains of gravity and capillary families continue to be separated by the demarcation lines up to U * ≈ 3 when the fan waves match the classic Kelvin cone. At slightly lower U * , the domains of capillary and gravity families start to overlap. While the capillary family isophases remain smooth, the gravity family curves form two types of cusps that correspond to two roots of (20) and (21). The longer wave cusp can be naturally related to the well-known caustics of the Kelvin wake of pure gravity waves. The shorter wave cusps (fan waves in words of [15]) are genetically related to the capillarity effect as it has been noted by Moissy and Rabaud [13]. Asymptotically, when U * → ∞ this Mach cone is shrinking as (24). The dimensionless wavenumbers at the cone are limited by maximal value κ min = 2/ √ 3 − 1 ≈ 0.393 when the denominator of (21) vanishes. This wavenumber corresponds to the minimum of the group velocity of gravity-capillary waves as noted above. Thus, the inner cone of the anomalous dispersion (25) and (52) always remains inside the cone of capillary cusps.

Wake Amplitudes
The reference solution method provides analytical solutions (40)-(42) that can be easily converted to graphics. Presenting these figures, the authors are trying to accompany them by comments on possible physical effects governing the amplitude distributions.
The key amplitude factor D depends on two variables: the wavenumber κ that is determined by wave direction θ at a given speed U * (16) and the distance along the ray expressed by time τ or phase ϕ. The first remark to be made is the presence of two terms A and B in the expression of D (40): the first term on the right-hand side is proportional to phase ϕ in power 4 while the second one is just squared phase. Thus, vanishing of A qualitatively changes asymptotical behavior along wave rays (with distance from the source), while vanishing the term (1 − A) can lead to an extremum of amplitudes both as functions of distance and angle, and rather peculiar distributions of amplitudes in the ship wake.
The second remark concerns the demarcation line κ = 1 where amplitudes of wake vanish as it is seen from (40)-(42): A, B are singular in this case. The observability of this effect is an interesting point for experimentalists.
In this paper, we present distributions of wave amplitudes for relatively low speeds when these amplitudes are governed by the only wave system at a given coordinate point that is for U * ≤ 1.938. At higher speeds after singularity occurrence, two (or even three) wave systems can be easily superimposed as it has been done for pure gravity waves in [9]. Here, we restrict ourselves by the simplest case. Note that in this case there is no similarity with the case of pure gravity waves where two different overlapping wave systems (divergent and transversal) always co-exist irrespectively to the current speed value.

Wakes of Circular Sources
For the circular shape of the source when ∆ x = ∆ y , wave amplitudes (terms A, B, D) do not depend on angle θ explicitly. While wavenumber κ depends on θ monotonically basic features of distributions can be investigated qualitatively. Figure 4 shows the amplitude factor D −1/2 (k 0 , t) in (35) as function of dimensionless coordinates (ξ, η), that is, a relative decrease of particular wave harmonics with the distance along a ray. Wave patterns are shown in two columns: the left panels are zoomed images of the right ones. Colorbars are scaled by the exponent n of 2 −n .  (e) (f)  (17) is shown in magenta, the "group Mach cone" 0.7681C m /U * -cyan, the Kelvin cone-black. The patterns qualitatively resemble experimentally observed ones by (Rousseaux and Kellay [17], cf. Figures 7 and 10).
At the low speed of the first Wilton resonance U * ≈ 1.03 (Figure 4a,b), the isolines of amplitudes in the wake resemble slightly distorted circles. The amplitude vanishes at the demarcation line κ = 1. It decays approximately 15 times on the scale of the right panel of Figure 4 about 5 m. Figure 4c,d likely reflects an effect of the group velocity minimum. The source produces a sort of bulldozing pattern with a mushroom-shaped front and a track along the axis behind. The amplitude decay is minimal near the current axis: only 6 times at a distance of about 30 m (1600 · L) along the axis versus 30 times at oblique directions.
At higher threshold speed U * = 1.983 (Figure 4e,f), the capillary waves propagating upstream show slower decay than their gravity counterparts. The high amplitudes are concentrated near the "group Mach" cone (cyan lines). Amplitudes of the gravity waves inside this cone appear almost one order lower than in the previous case at the same distances.

Wakes of Elliptic Sources
The reference solution method, in contrast to its asymptotic counterparts, accounts for the shape of the initial perturbation. These are bandwidth parameters ∆ x , ∆ y in formulas for wave amplitudes defined by (39) and their reciprocal values 1/∆ x , 1/∆ y that characterize the initial perturbation dimension in the coordinate space. For correct comparison, two pairs of the bandwidth parameters that satisfy an evident condition have been chosen, i.e., two elliptic distributions of equal square areas in the coordinate and wavenumber spaces. Figure 5 shows wave patterns for two sets of the shape parameters for two values of dimensionless speed U * .
Case ∆ x = 0.2, ∆ y = 0.05 (Figure 5a,b,e,f) presents a source for which its along-current size is 4 times shorter than its front size. This source demonstrates a pronounced bulldozing effect.
(minimum of group velocity), the mushroom pattern is larger than for the circular source (cf. Figure 4c,d) and the track is essentially narrower. For higher speed U * = 1.938, a similar mushroom-shaped pattern is seen near the Mach cone (cyan lines).
Narrow source with high length-beam ratio ∆ y /∆ x = 4 shows peculiar features in Figure 5c,d,g,h. For (minimum of group velocity), the wake appears larger than in other cases (circle and wide sources). For U * = 1.938, the wake perturbations are concentrating near the Mach angle and strongly decaying both to the short capillary waves and to the interior of the Kelvin angle of gravity waves. (g) (h) Figure 5. Isolines of wake amplitudes for the elliptic sources in two resolutions (left and right columns).
It would be interesting to verify the above estimates of the wave amplitude patterns experimentally.

Discussion and Conclusions
The gravity-capillary waves play a special role in ocean studies. One of the issues is the remote sensing applications. Satellite altimeters (Ka-and Ku-band) and other space-and airborne systems operate in the intermediate range of wavelength where the effects of gravity and capillarity compete on equal terms. This is a reason for focusing on the specific range of scales and/or relatively low dimensionless speed of U * . In this case, wave patterns differ dramatically from the well-known pictures of pure gravity waves with the classic Kelvin cone.
We started with the analysis of the gravity-capillary wave kinematics following the conventional approach developed in previous works (e.g., [1,[21][22][23]). The feature of our approach is in the scale choice. The wave Bond number κ = |k|L provides a useful reference κ = 1 for capillary (κ > 1) and gravity (κ < 1) families of the ship wake in wavenumber and coordinate spaces. The counterpart of the Froude number, dimensionless speed of current U * = U/C m , is exploited for specifying a variety of patterns of gravity-capillary waves.
Parameters of wave patterns are specified analytically, thus extending the results of recent works (e.g., [13,24]). In particular, analytical expressions for parameters of cusps of isophases are presented (e.g., (22)- (24)) and related to phase and group velocity extremes, and characteristic angles (see Figure 1).
The analytical expressions for the wake amplitudes are derived with the reference solutions approach for the first time. The "phase Mach" cone is shown to be invisible because of vanishing amplitudes along the demarcation line of gravity and capillary families. On the contrary, the "group Mach" cone is emphasized by the effect of quasi-dispersion of two-dimensional waves when the corresponding dispersion factor (51) vanishes. This is in line with the analytical and experimental results of Moissy and Rabaud [13]. The amplitude effect of the "group Mach" cone is demonstrated for a limited range of dimensionless speed U * ≤ 1.938 when cusps are absent. It is natural to assume that the presence of the merging waves at the cusps enhances this effect and emphasizes the "group Mach" effect at higher speeds.
The reference solution approach goes beyond the point-wise paradigm of its asymptotic counterparts (e.g., stationary phase approach), thus accounting for the effect of the source shape.
The corresponding estimates have been carried out for elliptic sources of different orientation. The bulldozing effect of the shapes on the wake has been demonstrated in Figure 5 fairly well.
Author Contributions: Conceptualization, formal analysis, writing-original draft preparation-V.G. Methodology, formal analysis, visualization, writing-review and editing-S.B. Both authors have read and agreed to the published version of the manuscript.