Soft matter physics of the ground beneath our feet.

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

Wavier jet streams driven by zonally asymmetric surface thermal forcing.

Recent studies have argued that global warming is responsible for a wavier jet stream, thereby driving midlatitude extreme flooding and drought. Polar amplification-the relative enhancement of high-latitude temperatures under global warming-is argued to be the principal climate state driving midlatitude extremes. Namely, the decreased meridional temperature gradient suppresses the mean zonal winds, leading to wavier midlatitude jets. However, although observations are consistent with such a linkage, a detailed dynamical mechanism is still debated. Here, we argue that the Northern Hemisphere land-sea thermal forcing contrast that underlies zonally asymmetric forcing drives a response in the planetary geostrophic motion, which provides balanced mean fields for synoptic eddies in midlatitudes and thus for wavier jet streams. We show that when the barotropic zonal mean wind U is smaller than a threshold, proportional to the ß-plane effect and dry static stability, the flow field exhibits a dramatic transition from a response confined near the surface to one reaching the upper atmosphere. As global warming enhances polar amplification, the midlatitude jet stream intensity is suppressed. The confluence of these effects leads to wavier jet streams.

Coupling functions in climate.

We examine how coupling functions in the theory of dynamical systems provide a quantitative window into climate dynamics. Previously, we have shown that a one-dimensional periodic non-autonomous stochastic dynamical system can simulate the monthly statistics of surface air temperature data. Here, we expand this approach to two-dimensional dynamical systems to include interactions between two sub-systems of the climate. The relevant coupling functions are constructed from the covariance of the data from the two sub-systems. We demonstrate the method on two tropical climate indices, the El-Niño-Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), to interpret the mutual interactions between these two air-sea interaction phenomena in the Pacific and Indian Oceans. The coupling function reveals that the ENSO mainly controls the seasonal variability of the IOD during its mature phase. This demonstrates the plausibility of constructing a network model for the seasonal variability of climate systems based on such coupling functions. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Hydrodynamic interactions and the diffusivity of spheroidal particles.

It is intuitive that the diffusivity of an isolated particle differs from those in a monodisperse suspension, in which hydrodynamic interactions between the particles are operative. Batchelor [J. Fluid Mech. 74, 1-29 (1976) and J. Fluid Mech. 131, 155-175 (1983)] calculated how hydrodynamic interactions influenced the diffusivity of a dilute suspension of spherical particles, and Russel et al. [Colloidal Dispersions (Cambridge University Press, 1991)] and Brady [J. Fluid Mech. 272, 109-134 (1994)] treated nondilute (higher particle volume fraction) suspensions. Although most particles lack perfect sphericity, little is known about the effects of hydrodynamic interactions on the diffusivity of spheroidal particles, which are the simplest shapes that can be used to model anisotropic particles. Here, we calculate the effects of hydrodynamic interactions on the translational and rotational diffusivities of spheroidal particles of arbitrary aspect ratio in dilute monodisperse suspensions. We find that the translational and rotational diffusivities of prolate spheroids are more sensitive to eccentricity than for oblate spheroids. The origin of the hydrodynamic anisotropy is that found in the stresslet field for the induced-dipole interaction. However, in the dilute limit, the effects of anisotropy are at the level of a few percent. These effects have influence on a vast range of settings, from partially frozen colloidal suspensions to the dynamics of cytoplasm.

Fluctuations in Arctic sea-ice extent: comparing observations and climate models.

The fluctuation statistics of the observed sea-ice extent during the satellite era are compared with model output from CMIP5 models using a multifractal time series method. The two robust features of the observations are that on annual to biannual time scales the ice extent exhibits white noise structure, and there is a decadal scale trend associated with the decay of the ice cover. It is shown that (i) there is a large inter-model variability in the time scales extracted from the models, (ii) none of the models exhibits the decadal time scales found in the satellite observations, (iii) five of the 21 models examined exhibit the observed white noise structure, and (iv) the multi-model ensemble mean exhibits neither the observed white noise structure nor the observed decadal trend. It is proposed that the observed fluctuation statistics produced by this method serve as an appropriate test bed for modelling studies.This article is part of the theme issue 'Modelling of sea-ice phenomena'.

Fluctuation spectra and force generation in nonequilibrium systems.

Many biological systems are appropriately viewed as passive inclusions immersed in an active bath: from proteins on active membranes to microscopic swimmers confined by boundaries. The nonequilibrium forces exerted by the active bath on the inclusions or boundaries often regulate function, and such forces may also be exploited in artificial active materials. Nonetheless, the general phenomenology of these active forces remains elusive. We show that the fluctuation spectrum of the active medium, the partitioning of energy as a function of wavenumber, controls the phenomenology of force generation. We find that, for a narrow, unimodal spectrum, the force exerted by a nonequilibrium system on two embedded walls depends on the width and the position of the peak in the fluctuation spectrum, and oscillates between repulsion and attraction as a function of wall separation. We examine two apparently disparate examples: the Maritime Casimir effect and recent simulations of active Brownian particles. A key implication of our work is that important nonequilibrium interactions are encoded within the fluctuation spectrum. In this sense, the noise becomes the signal.

A unified nonlinear stochastic time series analysis for climate science.

Earth's orbit and axial tilt imprint a strong seasonal cycle on climatological data. Climate variability is typically viewed in terms of fluctuations in the seasonal cycle induced by higher frequency processes. We can interpret this as a competition between the orbitally enforced monthly stability and the fluctuations/noise induced by weather. Here we introduce a new time-series method that determines these contributions from monthly-averaged data. We find that the spatio-temporal distribution of the monthly stability and the magnitude of the noise reveal key fingerprints of several important climate phenomena, including the evolution of the Arctic sea ice cover, the El Nio Southern Oscillation (ENSO), the Atlantic Nio and the Indian Dipole Mode. In analogy with the classical destabilising influence of the ice-albedo feedback on summertime sea ice, we find that during some time interval of the season a destabilising process operates in all of these climate phenomena. The interaction between the destabilisation and the accumulation of noise, which we term the memory effect, underlies phase locking to the seasonal cycle and the statistical nature of seasonal predictability.

Roughness as a Route to the Ultimate Regime of Thermal Convection.

We use highly resolved numerical simulations to study turbulent Rayleigh-Bénard convection in a cell with sinusoidally rough upper and lower surfaces in two dimensions for Pr=1 and Ra=[4×10^{6},3×10^{9}]. By varying the wavelength λ at a fixed amplitude, we find an optimal wavelength λ_{opt} for which the Nusselt-Rayleigh scaling relation is (Nu-1∝Ra^{0.483}), maximizing the heat flux. This is consistent with the upper bound of Goluskin and Doering [J. Fluid Mech. 804, 370 (2016)JFLSA70022-112010.1017/jfm.2016.528] who prove that Nu can grow no faster than O(Ra^{1/2}) as Ra→∞, and thus with the concept that roughness facilitates the attainment of the so-called ultimate regime. Our data nearly achieve the largest growth rate permitted by the bound. When λ≪λ_{opt} and λ≫λ_{opt}, the planar case is recovered, demonstrating how controlling the wall geometry manipulates the interaction between the boundary layers and the core flow. Finally, for each Ra, we choose the maximum Nu among all λ, thus optimizing over all λ, to find Nu_{opt}-1=0.01×Ra^{0.444}.

Surface tension and a self-consistent theory of soft composite solids with elastic inclusions.

The importance of surface tension effects is being recognized in the context of soft composite solids, where they are found to significantly affect the mechanical properties, such as the elastic response to an external stress. It has recently been discovered that Eshelby's inclusion theory breaks down when the inclusion size approaches the elastocapillary length L≡γ/E, where γ is the inclusion/host surface tension and E is the host Young's modulus. Extending our recent results for liquid inclusions, here we model the elastic behavior of a non-dilute distribution of isotropic elastic spherical inclusions in a soft isotropic elastic matrix, subject to a prescribed infinitesimal far-field loading. Within our framework, the composite stiffness is uniquely determined by the elastocapillary length L, the spherical inclusion radius R, and the stiffness contrast parameter C, which is the ratio of the inclusion to the matrix stiffness. We compare the results with those from the case of liquid inclusions, and we derive an analytical expression for elastic cloaking of the composite by the inclusions. Remarkably, we find that the composite stiffness is influenced significantly by surface tension even for inclusions two orders of magnitude more stiff than the host matrix. Finally, we show how to simultaneously determine the surface tension and the inclusion stiffness using two independent constraints provided by global and local measurements.

Confinement effects in premelting dynamics.

We examine the effects of confinement on the dynamics of premelted films driven by thermomolecular pressure gradients. Our approach is to modify a well-studied setting in which the thermomolecular pressure gradient is driven by a temperature gradient parallel to an interfacially premelted elastic wall. The modification treats the increase in viscosity associated with the thinning of films, studied in a wide variety of materials, using a power law and we examine the consequent evolution of the confining elastic wall. We treat (1) a range of interactions that are known to underlie interfacial premelting and (2) a constant temperature gradient wherein the thermomolecular pressure gradient is a constant. The difference between the cases with and without the proximity effect arises in the volume flux of premelted liquid. The proximity effect increases the viscosity as the film thickness decreases thereby requiring the thermomolecular pressure driven flux to be accommodated at higher temperatures where the premelted film thickness is the largest. Implications for experiment and observations of frost heave are discussed.

Surface tension and the Mori-Tanaka theory of non-dilute soft composite solids.

Eshelby's theory is the foundation of composite mechanics, allowing calculation of the effective elastic moduli of composites from a knowledge of their microstructure. However, it ignores interfacial stress and only applies to very dilute composites-i.e. where any inclusions are widely spaced apart. Here, within the framework of the Mori-Tanaka multiphase approximation scheme, we extend Eshelby's theory to treat a composite with interfacial stress in the non-dilute limit. In particular, we calculate the elastic moduli of composites comprised of a compliant, elastic solid hosting a non-dilute distribution of identical liquid droplets. The composite stiffness depends strongly on the ratio of the droplet size, R, to an elastocapillary lengthscale, L. Interfacial tension substantially impacts the effective elastic moduli of the composite when [Formula: see text]. When R<3L/2 (R=3L/2) liquid inclusions stiffen (cloak the far-field signature of) the solid.

Segregated Ice Growth in a Suspension of Colloidal Particles.

We study the freezing of a dispersion of colloidal silica particles in water, focusing on the formation of segregated ice in the form of ice lenses. Local temperature measurements in combination with video microscopy give insight into the rich variety of factors that control ice lens formation. We observe the initiation of the lenses, their growth morphology, and their final thickness and spacing over a range of conditions, in particular the effect of the particle packing and the cooling rate. We find that increasing the particle density drastically reduces the thickness of lenses but has little effect on the lens spacing. Therefore, the fraction of segregated ice formed reduces. The effect of the cooling rate, which is the product of the temperature gradient and the pulling speed across the temperature gradient, depends on which parameter is varied. A larger temperature gradient causes ice lenses to be initiated more frequently, while a lower pulling speed allows for more time for ice lenses to grow: both increase the fraction of segregated ice. Surprisingly, we find that the growth rate of a lens does not depend on its undercooling. Finally, we have indications of pore ice in front of the warmest ice lens, which has important consequences for the interpretation of the measured trends. Our findings are relevant for ice segregation occurring in a wide range of situations, ranging from model lab experiments and theories to geological and industrial processes, like frost heave and frozen food production.

Interfacial tension and a three-phase generalized self-consistent theory of non-dilute soft composite solids.

In the dilute limit Eshelby's inclusion theory captures the behavior of a wide range of systems and properties. However, because Eshelby's approach neglects interfacial stress, it breaks down in soft materials as the inclusion size approaches the elastocapillarity length L≡γ/E. Here, we use a three-phase generalized self-consistent method to calculate the elastic moduli of composites comprised of an isotropic, linear-elastic compliant solid hosting a spatially random monodisperse distribution of spherical liquid droplets. As opposed to similar approaches, we explicitly capture the liquid-solid interfacial stress when it is treated as an isotropic, strain-independent surface tension. Within this framework, the composite stiffness depends solely on the ratio of the elastocapillarity length L to the inclusion radius R. Independent of inclusion volume fraction, we find that the composite is stiffened by the inclusions whenever R < 3L/2. Over the same range of parameters, we compare our results with alternative approaches (dilute and Mori-Tanaka theories that include surface tension). Our framework can be easily extended to calculate the composite properties of more general soft materials where surface tension plays a role.

Surface tension and the mechanics of liquid inclusions in compliant solids.

Eshelby's theory of inclusions has wide-reaching implications across the mechanics of materials and structures including the theories of composites, fracture, and plasticity. However, it does not include the effects of surface stress, which has recently been shown to control many processes in soft materials such as gels, elastomers and biological tissue. To extend Eshelby's theory of inclusions to soft materials, we consider liquid inclusions within an isotropic, compressible, linear-elastic solid. We solve for the displacement and stress fields around individual stretched inclusions, accounting for the bulk elasticity of the solid and the surface tension (i.e. isotropic strain-independent surface stress) of the solid-liquid interface. Surface tension significantly alters the inclusion's shape and stiffness as well as its near- and far-field stress fields. These phenomena depend strongly on the ratio of the inclusion radius, R, to an elastocapillary length, L. Surface tension is significant whenever inclusions are smaller than 100L. While Eshelby theory predicts that liquid inclusions generically reduce the stiffness of an elastic solid, our results show that liquid inclusions can actually stiffen a solid when R<3L/2. Intriguingly, surface tension cloaks the far-field signature of liquid inclusions when R=3L/2. These results are have far-reaching applications from measuring local stresses in biological tissue, to determining the failure strength of soft composites.

Hertz beyond belief.

We examine the validity of Hertz's linear elastic theory for central collisions of spheres using a viscoelastic model. This model explains why Hertz's theory is accurate in predicting the collision time and maximum contact area even when 40% of the kinetic energy is lost due to viscous dissipation. The main reason is that both the collision time and maximum contact area have a very weak dependence on the impact velocity. Moreover, we show that colliding objects exhibit an apparent size-dependent yield strength, which results from larger objects dissipating less energy at a given impact velocity.

Surface tension and contact with soft elastic solids.

The Johnson-Kendall-Roberts theory is the basis of modern contact mechanics. It describes how two deformable objects adhere together, driven by adhesion energy and opposed by elasticity. Here we characterize the indentation of glass particles into soft, silicone substrates using confocal microscopy. We show that, whereas the Johnson-Kendall-Roberts theory holds for particles larger than a critical, elastocapillary lengthscale, it fails for smaller particles. Instead, adhesion of small particles mimics the adsorption of particles at a fluid interface, with a size-independent contact angle between the undeformed surface and the particle given by a generalized version of the Young's law. A simple theory quantitatively captures this behaviour and explains how solid surface tension dominates elasticity for small-scale indentation of soft materials.

Patterning droplets with durotaxis.

Numerous cell types have shown a remarkable ability to detect and move along gradients in stiffness of an underlying substrate--a process known as durotaxis. The mechanisms underlying durotaxis are still unresolved, but generally believed to involve active sensing and locomotion. Here, we show that simple liquid droplets also undergo durotaxis. By modulating substrate stiffness, we obtain fine control of droplet position on soft, flat substrates. Unlike other control mechanisms, droplet durotaxis works without imposing chemical, thermal, electrical, or topographical gradients. We show that droplet durotaxis can be used to create large-scale droplet patterns and is potentially useful for many applications, such as microfluidics, thermal control, and microfabrication.