To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
The Arctic is undergoing increased warming compared to the global mean, with major implications for the mass balance of glaciers. Direct observations of mass balance in the Russian Arctic are sparse and remotely sensed volume changes do not provide information about climatic drivers. Here, we present simulations of the climatic mass balance and meltwater runoff from glaciers in Franz Josef Land and Novaya Zemlya from 1991 to 2022. Based on simulations of glacier climatic mass balance over the period 1991–2022, we present a first detailed view of mass balance evolution in Franz Josef Land and Novaya Zemlya. The simulations are conducted at a 2.5 km resolution using the CryoGrid model forced by the Copernicus Arctic Regional ReAnalysis (CARRA) product. Over the 30 year simulation period, the climatic mass balance of both Franz Josef Land (0.21 m w.e. a−1) and Novaya Zemlya (0.07 m w.e. a−1) is positive on average without a significant trend in annual climatic mass balance. There is still a tendency towards more frequent high-melt years after 2010 and the associated glacier runoff has intensified with record melt years occurring during the model period.
In 2015, a continuous 15.4 m snow/firn core was recovered from central South Georgia Island at ∼850 m a.s.l. All firn core samples were analyzed for major (Al, Ca, Mg, Na, K, Ti and Fe) and trace element concentrations (Sr, Cd, Cs, Ba, La, Ce, Pr, Pb, Bi, U, As, Li, S, V, Cr, Mn, Co, Cu and Zn) and stable water isotopes. The chemical and isotopic signal is well preserved in the top 6.2 m of the core. Below this depth, down to the bottom of the core, signal dampening is observed in the majority of the elemental species making it difficult to distinguish a seasonal signal. Thirteen elements (As, Bi, Ca, Cd, Cu, K, Li, Mg, Na, Pb, S, Sr and Zn) have crustal enrichment factor values higher than 10 suggesting sources in addition to those found naturally in the crust. While this study shows that 850 m a.s.l. is not high enough to preserve a record including recent years, higher-elevation (>1250 m a.s.l.) glaciers may be likely candidates for ice core drilling to recover better-preserved, continuous, recent to past glaciochemical records.
In this paper, we study the disturbance velocity and density fields induced by a sphere translating vertically in a viscous density-stratified ambient. Specifically, we consider the limit of a vanishingly small Reynolds number $(Re = \rho U a/\mu \ll 1)$, a small but finite viscous Richardson number $(Ri_v = \gamma a^3g/\mu U\ll 1)$ and large Péclet number $(Pe = Ua/D\gg 1)$. Here, $a$ is the sphere's radius, $U$ its translational velocity, $\rho$ an appropriate reference density within the framework of the Boussinesq approximation, $\mu$ the ambient viscosity, $\gamma$ the absolute value of the background density gradient, g is acceleration due to gravity and $D$ the diffusivity of the stratifying agent. For the scenario where buoyancy forces first become comparable to viscous forces at large distances, corresponding to the Stokes-stratification regime defined by $Re \ll Ri_v^{1/3} \ll 1$ for $Pe \gg 1$, important flow features have been identified by Varanasi & Subramanian (J. Fluid Mech., vol. 949, 2022, A29) – these include a vertically oriented reverse jet, and a horizontal axisymmetric wake, on scales larger than the primary (stratification) screening length of ${O}(aRi_v^{-1/3})$. Here, we study the reverse-jet region in more detail, and show that it is only the central portion of a columnar structure with multiple annular cells concentric about the rear stagnation streamline. In the absence of diffusion, corresponding to $Pe = \infty$$( \beta _\infty = Ri_v^{1/3}Pe^{-1} = 0)$, this columnar structure extends to downstream infinity with the number of annular cells diverging in this limit. We provide expressions for the boundary of the structure, and the number of cells within, as a function of the downstream distance. For small but finite $\beta _\infty$, two length scales emerge in addition to the primary screening length – a secondary screening length of ${O}(aRi_v^{-1/2}Pe^{1/2})$ where diffusion starts to smear out density variations across cells, leading to exponentially decaying density and velocity fields; and a tertiary screening length, $l_t \sim {O}(aRi_v^{-1/2}Pe^{1/2}[\zeta + \frac {13}{4}\ln {\zeta } + ({13^2}/{4^2})({\ln \zeta }/{\zeta })])$ with $\zeta = \frac {1}{2}\ln ({\sqrt {{\rm \pi} }Ri_v^{-1}Pe^3}/{2160})$, beyond which the columnar structure ceases to exist. The latter causes a transition from a vertical to a predominantly horizontal flow, with the downstream disturbance fields reverting from an exponential to an eventual algebraic decay, analogous to that prevalent at large distances upstream.
As snowlines retreat, the bare ice of Central Asian glaciers is increasingly exposed to short-wave radiation and high temperatures. The importance of bare-ice albedo for glacier melt rates is thus rising. Little is known about the variability of bare-ice albedo, its drivers or its implications for glacier melt. We address this gap by presenting the sub-seasonal and interannual variability of bare-ice albedo of Abramov Glacier in Kyrgyzstan between 1999 and 2022. We derived albedo products from Landsat surface reflectance data, investigated the relationship between air temperature and bare-ice albedo variability and explored the implications of this variability for glacier melt. Our results indicate that bare-ice albedo undergoes a sub-seasonal cycle controlled by air temperature and elevation-dependent refreezing events. Bare-ice albedo decreased over the tongue in July and August between 1999 and 2017, while, in 2018, a lateral displacement of the ice resulted in a shift in the patterns of bare-ice albedo. We found significant correlations between bare-ice albedo variability and both temperature and glacier melt at various timescales. Rising temperatures are thus expected to lead to darker bare ice and amplified feedback melt cycles. Integrating albedo variability into glaciological models is thus crucial for accurate predictions of accelerated glacier response to intensifying climate change.
We carry out a linear stability analysis of the flow of a thin layer of Newtonian fluid with a deformable free surface bounded at the bottom by a horizontal wall subjected to quasi-periodic oscillation in its own plane. Or's model (J. Fluid Mech., vol. 335, 1997, pp. 213–232), using a periodic oscillation, is extended to the configuration where oscillation has two incommensurate frequencies, $\omega _1$ and $\omega _2$, with an irrational ratio $\omega ={\omega _2}/{\omega _1}$. Using the long-wave expansion, we derive the asymptotic function involved in the long-wave instability criterion while taking into account the frequency ratio. It turns out that the maximum of this asymptotic function, as well as the frequency parameter at which long-wave instabilities occur, depend strongly on the frequency ratio. For arbitrary wavenumbers, the equations governing the problem under consideration are solved in space using Chebyshev's spectral collocation method, while the temporal resolution is performed using Floquet theory, knowing that an irrational number can be approximated by a rational number. For a large frequency ratio and for a velocity amplitude ratio equal to unity, we obtain, as in Or's work (J. Fluid Mech., vol. 335, 1997, pp. 213–232) considering the same frequency parameter interval, an alternation between the U shape and oblique shape referring respectively to instabilities of long wavelength and finite wavelength appearing in the diagram representing Reynolds number as a function of frequency parameter. By decreasing the frequency ratio towards $1/\sqrt {37}$, the three initial U-shaped and three oblique instabilities merge into a single U-shaped and a single oblique instability. This merging phenomenon also occurs when the ratio of the amplitudes of the superimposed velocities, linked to the introduction of the second frequency, increases from small values to unity. For a fixed frequency parameter, the effect of frequency ratio and velocity amplitude ratio on the marginal stability curves in terms of Reynolds number versus wavenumber is also investigated, focusing on the appearance of long wavelength instability and finite wavelength instability.
When one fluid is injected into a confined geometry such as a porous medium filled with another immiscible fluid, even at an extremely low injection speed, rapid filling of several pore spaces accompanied by retraction of multiple fluid–fluid interfaces can be observed. Such processes with fast liquid redistribution within the solid structure, called Haines jumps, are ubiquitous in many multiphase flow systems, which can impact fluid trapping, energy dissipation and hysteretic saturation in various engineering applications. Inspired by this mechanism, here, we propose a dual-channel structure to realise controlled Haines jumps during fluid displacement processes. Via theoretical analysis and numerical simulations, we show that the dynamics of fluid interfaces during Haines jumps can be quantitatively correlated with the driving capillary pressure and dissipating viscous stress, which enables simultaneous determination of the fluid viscosity and interfacial tension in the dual-channel multiphase system.
We investigate nonlinear energy transfer for channel flows at friction Reynolds numbers $Re_{\tau }=180$ and $590$. The key feature of the analysis is that we quantify the energy transferred from a source mode to a recipient mode, with each mode characterised by a streamwise wavenumber and a spanwise wavenumber. This is achieved through an explicit examination of the triadic interactions of the nonlinear energy transfer term in the spectral turbulent kinetic energy equation. First, we quantify the nonlinear energy transfer gain and loss for individual Fourier modes. The gain and loss cannot be obtained without expanding the nonlinear triadic interactions. Second, we quantify the nonlinear energy transfer budgets for three types of modes. Each type of mode is characterised by a specific region in streamwise–spanwise wavenumber space. We find that a transverse cascade from streamwise-elongated modes to spanwise-elongated modes exists for all three types of modes. Third, we quantify the forward and inverse cascades between resolved scales and subgrid scales in the spirit of large-eddy simulations. For the cutoff wavelength range that we consider, the forward and inverse cascades between the resolved scales and subgrid scales result in a net forward cascade from the resolved scales to the subgrid scales. The shape of the net forward cascade curve with respect to the cutoff wavelength resembles the net forward cascade predicted by the Smagorinsky eddy viscosity.
The two main large-scale features of Arctic sea-ice drift are the Beaufort Gyre and the Transpolar Drift Stream. They exhibit strong intraseasonal and interannual variability. Winter 2016/17 showed increased cyclone activity, leading to the collapse of the Beaufort Sea high and the reversal of the Beaufort Gyre. Winter 2020/21 displayed decreased cyclone activity and intense anticyclonic ice transport in the Beaufort Gyre. Here we show that the European Centre for Medium-Range Weather Forecasts’s (ECMWF) extended-range (46 days) retrospective forecasts were able to predict the ice motion during these cases. The initial contrasts in sea level pressure, surface winds and ice drift were well captured, and their temporal evolution—including the reversal of the usual drift direction—well reproduced by the forecasts initialized about a week before the event. Sea-ice thickness in the forecast exhibited initial errors even greater than 1 m that persisted throughout the forecast and negatively affected the ice speed forecast. Despite these shortcomings, the dynamic forecast outperformed the persistence and climatology forecast and represented the observed relation between surface winds and ice drift well. The benefit of dynamic forecasts is especially clear in cases that differ from climatology, like the one we focus on.
Data assimilation is a core component of numerical weather prediction systems. The large quantity of data processed during assimilation requires the computation to be distributed across increasingly many compute nodes; yet, existing approaches suffer from synchronization overhead in this setting. In this article, we exploit the formulation of data assimilation as a Bayesian inference problem and apply a message-passing algorithm to solve the spatial inference problem. Since message passing is inherently based on local computations, this approach lends itself to parallel and distributed computation. In combination with a GPU-accelerated implementation, we can scale the algorithm to very large grid sizes while retaining good accuracy and compute and memory requirements.
Azooxanthellate corals of Rhizangiidae, with their distinctive morphological and ecological features, are widely distributed across global oceans but remain under-studied due to identification challenges. Comprehensive underwater surveys across diverse marine habitats such as rocky reefs and submerged shipwrecks were undertaken, and findings highlighted the exclusive presence of Culicia stellata in natural rocks, off the coast of Kaup for the first time in Indian EEZ. The study also records the presence of Cladangia exusta, nearly six decades after the first record from off Cochin, Laccadive Sea. This research presents a detailed study of the taxonomic description, distribution, and ecological preferences of Cladangia exusta and Culicia stellata in the Laccadive Sea. This study underscores the importance of accurate species identification for effective conservation strategies and enriching biodiversity records.
Morphological and molecular methods were used to describe a new species of Trapania Pruvot-Fol, 1931 from shallow water kelp forests on the north-central coast of Peru. The new species, Trapania huarmeyana sp. nov., is distinguished from other species along the Eastern Pacific by external morphological characters such as its translucent white body with brown stripes and small spots on the dorsum, blotches on the base of the extra-branchial processes, extra rhinophoral processes and gill branches. Internally, T. huarmeyana sp. nov. is distinguishable by several morphological characteristics of the radula, jaws and genital organs. Phylogenetic trees recovered using Bayesian Inference and Maximum Likelihood analysis of DNA sequences support its distinct status and clarify its relationship to other species from the Eastern Pacific. This new species constitutes the first record of Trapania from the Humboldt Current Ecosystem, contributing to our understanding of the distribution of the genus in the South-eastern Pacific.
Although nudibranchs are common and attractive animals, our understanding of these marine gastropods in Vietnam remains limited. Prior research has suggested that combining morphological examination with molecular analysis results in more accurate identification of nudibranchs. However, previous studies in Vietnam have typically relied solely on morphological methods for nudibranch identification. In this study, the nudibranch species Halgerda batangas was recorded in Vietnam for the first time based on both morphological and molecular approaches. Halgerda batangas was characterized by a network pattern consisting of orange lines, relatively low dorsal tubercles with red-orange caps and white basal rings, and an orange line along the foot margin. Molecular analysis corroborated the morphological findings. These results suggest that integrating morphological and molecular methods is an effective approach for identifying nudibranchs.
The small-scale velocity gradient is connected to fundamental properties of turbulence at the large scales. By neglecting the viscous and non-local pressure Hessian terms, we derive a restricted Euler model for the turbulent flow along an undeformed free surface and discuss the associated stable/unstable manifolds. The model is compared with the data collected by high-resolution imaging on the free surface of a turbulent water tank with negligible surface waves. The joint probability density function (p.d.f.) of the velocity gradient invariants exhibits a distinct pattern from the one in the bulk. The restricted Euler model captures the enhanced probability along the unstable branch of the manifold and the asymmetry of the joint p.d.f. Significant deviations between the experiments and the prediction are evident, however, in particular concerning the compressibility of the surface flow. These results highlight the enhanced intermittency of the velocity gradient and the influence of the free surface on the energy cascade.
The yellownose skate (Dipturus chilensis) and roughskin skate (Dipturus trachyderma) are the only two elasmobranch species targeted by commercial fishing operations in Chile. Despite their importance, much of their biology and ecology remain poorly understood. This research aimed to evaluate the feasibility of tagging these species. In 2021, a pilot study was conducted at two locations, utilizing Petersen discs, acoustic transmitters, and pop-up satellite transmitters on both species. The results revealed a 6% recovery rate from the 50 skates tagged with Petersen discs, while 29.4% of those tagged with acoustic transmitters were successfully detected. Additionally, data from all ten satellite transmitters were successfully transmitted and recovered. The results revealed a maximum horizontal movement of 35.9 km, with the duration of liberty ranging from 8 to 275 days. Stocks of both species are currently depleted, and fishery management relies on closures and total allowable catches, where fishing effort is concentrated in short spatial and temporal windows. These particularities present significant challenges for implementing a national tagging programme, especially in terms of tag recovery. The main conclusion of this research is that the implementation of a tagging programme for both species is feasible. Satellite tagging provides the best results, but its higher implementation cost and limitations in use for relatively small skates could be mitigated by combining it with Petersen discs. Establishing a long-term tagging programme is essential for enhancing the understanding of distribution and migration patterns, which is crucial for enhancing conservation and management efforts for these skates in Chile.
While the giant anemone, Relicanthus daphneae, has been described as a characteristic inhabitant of the East Pacific Ocean since 1991, there are relatively few published occurrences worldwide. Here, we present the discovery and molecular verification of R. daphneae along the southern Central Indian Ridge, at the Rodriguez Triple Junction, and along the northern Southeast Indian Ridge within the BGR contract area for the exploration of marine massive sulphide deposits in the Indian Ocean. Individuals were solitary and attached exclusively to basalt hard substrates on the periphery of hydrothermal vent fields, at distances from active vents between 66 and 710 m. We report megafauna observed in close proximity to R. daphneae and, in one case, polychaetes on its tentacles and oral disc. For the first time, the giant anemone was observed capturing prey, a shrimp of the species Rimicaris kairei. Beyond this remark on the diet of these anemones, we also report other behavioural aspects for this species.
The concept of vortex lock-in for a single circular cylinder in an oscillating flow, induced through acoustic forcing, is revisited. Multiple cylinder diameters are investigated over a Reynolds number range between 500 and 7200. The lock-in behaviour is investigated quantitatively through hot-wire anemometry and planar particle image velocimetry measurements. The results corroborate previous findings describing the frequency range over which vortex lock-in occurs. It is found that the cylinder location in a standing wave (pressure node or velocity node) had a significant influence on the lock-in behaviour. A novel scaling which captures the onset of vortex lock-in is proposed which demonstrates that the Strouhal number is important in predicting the amplitude of the velocity fluctuations required to induce lock-in. Velocity fields also reveal the existence of bimodal vortex shedding during lock-in. This is confirmed using snapshot proper orthogonal decomposition which demonstrates that symmetric and alternate shedding modes are simultaneously present during lock-in and that symmetric shedding is inherent to the near wake region only. Reduced-order reconstruction of the instantaneous velocity fields confirmed that features associated with the forcing frequency control the shear layer roll-up up to $x/d=2.1$ while the influence of the asymmetric mode is simply to skew the trajectory of the vortex pair. Since vortex roll-up and the cylinder wake ends at $x/d=2.1$, the emergence of spectral content at $0.5f_e$ is attributed to a wavelength doubling measured between the vortical structures in the flow field.
Wind speed at the sea surface is a key quantity for a variety of scientific applications and human activities. For its importance, many observation techniques exist, ranging from in situ to satellite observations. However, none of such techniques can capture the spatiotemporal variability of the phenomenon at the same time. Reanalysis products, obtained from data assimilation methods, represent the state-of-the-art for sea-surface wind speed monitoring but may be biased by model errors and their spatial resolution is not competitive with satellite products. In this work, we propose a scheme based on both data assimilation and deep learning concepts to process spatiotemporally heterogeneous input sources to reconstruct high-resolution time series of spatial wind speed fields. This method allows to us make the most of the complementary information conveyed by the different sea-surface information typically available in operational settings. We use synthetic wind speed data to emulate satellite images, in situ time series and reanalyzed wind fields. Starting from these pseudo-observations, we run extensive numerical simulations to assess the impact of each input source on the model reconstruction performance. We show that our proposed framework outperforms a deep learning–based inversion scheme and can successfully exploit the spatiotemporal complementary information of the different input sources. We also show that the model can learn the possible bias in reanalysis products and attenuate it in the output reconstructions.
An efficient compression scheme for modal flow analysis is proposed and validated on data sequences of compressible flow through a linear turbomachinery blade row. The key feature of the compression scheme is a minimal, user-defined distortion of the mutual distance of any snapshot pair in phase space. Through this imposed feature, the model reduction process preserves the temporal dynamics contained in the data sequence, while still decreasing the spatial complexity. The mathematical foundation of the scheme is the fast Johnson–Lindenstrauss transformation (FJLT) which uses randomized projections and a tree-based spectral transform to accomplish the embedding of a high-dimensional data sequence into a lower-dimensional latent space. The compression scheme is coupled to a proper orthogonal decomposition and dynamic mode decomposition analysis of flow through a linear blade row. The application to a complex flow-field sequence demonstrates the efficacy of the scheme, where compression rates of two orders of magnitude are achieved, while incurring very small relative errors in the dominant temporal dynamics. This FJLT technique should be attractive to a wide range of modal analyses of large-scale and multi-physics fluid motion.
In recent years, the generalised quasilinear (GQL) approximation has been developed and its efficacy tested against purely quasilinear (QL) approximations. GQL systematically interpolates between QL and fully nonlinear dynamics by employing a generalised Reynolds decomposition. Here, we examine an exact statistical closure for the GQL equations on the doubly periodic $\beta$-plane. Closure is achieved at second order using a generalised cumulant approach which we term GCE2. GCE2 is shown to yield improved performance over statistical representations of purely QL dynamics (CE2) and thus enables direct statistical simulation of complex mean flows that do not entirely fall within the remit of pure QL theory. Despite the existence of an exact closure, GCE2 like CE2 admits the possibility of a rank instability that leads to differences with statistics obtained from GQL. Recognition of this instability is a necessary step before further progress can be made with the GCE2 statistical closure.