Most macroscopic organisms detect odors in intermittent bursts, that may be separated by extended regions with no odor. Organisms leverage this complex dynamics efficiently for diverse tasks, including locating and identifying an odor source Murlis et al., 1992; Mafra-Neto and Cardé, 1994; Vickers, 2000; Riffell et al., 2014; Ache et al., 2016; Ackels et al., 2021. However, what are the most informative features of intermittent odor cues remains largely unclear. There are two broad classes of measures that quantify the dynamics of olfactory cues: those that depend on odor intensity including e.g. odor gradients in space or time, and those that do not depend on odor intensity but only on its timing, i.e. on whether the odor is on or off regardless of its concentration. To compute quantities that depend on odor intensity, an accurate representation of the odor is needed. In contrast, measuring the timing of odor detection simply requires to mark at all times whether the odor is on or off, thus a binary switch is sufficient.

Behavioral evidence suggests that animals use both intensity and timing of odor encounters for olfactory navigation Baker et al., 2018. At close range, mammals appear to compare odor intensity either across nostrils or across sniffs Catania, 2013; Gire et al., 2016; Findley et al., 2021. On the other hand, mounting evidence suggests timing of odor detection also plays a key role for olfactory navigation Ache et al., 2016: moths respond to odor pulsed at specific frequencies Riffell et al., 2014; Vickers et al., 2001; fruit flies respond to timing since last odor detection van Breugel and Dickinson, 2014; Demir et al., 2020; lobsters and sharks compare odor arrival time across their paired olfactory organs and orient toward the side that detected the odor first Basil and Atema, 1994; Gardiner and Atema, 2010; many organisms will move upwind upon detection of an odor Kennedy and Marsh, 1974; Murlis et al., 1992; Steck et al., 2012.

Neural recordings upon stimulation with intermittent odor cues confirm that the brain of many animals is able to record information both about intensity (and its derivatives) as well as timing of odor encounters (most information comes from work on arthropods Nagel and Wilson, 2011; Vickers et al., 2001; Brown et al., 2005; Gorur-Shandilya et al., 2017; Jacob et al., 2017; Riffell et al., 2014, but see also Parabucki et al., 2019; Lewis et al., 2021). For example, when insects are presented with intermittent odor cues, information about intensity and timing is recorded in their antennal lobe (see e.g. Vickers et al., 2001; Brown et al., 2005). Odors that mimic natural intermittency elicit a response that preserves an accurate measure of timing in fruit flies and moths Gorur-Shandilya et al., 2017; Jacob et al., 2017. In lobsters, bursting olfactory neurons encode specifically for the time between successive odor encounters, see Park et al., 2014; Park et al., 2016 and references therein. Interestingly, the neural activity varies considerably with the dynamics of the odor cues Nagel and Wilson, 2011; Vickers et al., 2001; Lewis et al., 2021, but how intermittency of an odor affects its neural representation is not well understood.

This evidence suggests animals are able to identify when they detect an odor as well as how intense it is; but whether they record and rely on both kinds of information is not understood. From a physical perspective, these two measures clearly provide information about source location. Indeed, we know from theoretical Shraiman and Siggia, 2000; Falkovich et al., 2001; Celani et al., 2014 and experimental Justus et al., 2002; Moore and Crimaldi, 2004 work that turbulence causes the odor to be distributed in highly intermittent patches separated by blanks with no odor. Both intensity and timing of these intermittent bursts vary depending on the location of the source Celani et al., 2014, as early recognized by Atema, 1996, thus can be used to infer source location or navigate to it Vergassola et al., 2007; Schmuker et al., 2016; Boie et al., 2018; Leathers et al., 2020; Michaelis et al., 2020.

Here, we ask what salient features of turbulent odor signals best predict the location of the odor source and specifically compare quantities related to intensity vs timing of odor encounters. We first compose a dataset of realistic odor fields at scales of several meters using accurate state-of-the-art fluid dynamics simulations. We then develop machine learning algorithms that predict source location based on these synthetic odor fields.

We find that measures of odor temporal dynamics based on a short memory span (down to about 1 s) hold information about source location. Close to the source or close to the substrate, measures of intensity predict distance better than measures of timing; but this ranking is reversed at further distance from the source or from the substrate. Pairing the two kinds of measure improves dramatically the quality of the prediction robustly across all datasets, whereas pairing two measures of intensity or two measures of timing is either useless or detrimental.

Our results demonstrate that timing and intensity are complementary attributes of odor dynamics and are most effective in more dilute and concentrated conditions respectively. These different conditions exist in different portions of space because odor gets transported, mixed and diluted by the fluid. As a result, the spatial range of operation of a living organism constrains the solutions it may evolve to make predictions with turbulent odors.

Odor cues at several meters from the source are often turbulent. Figure 1a–c and show snapshots of the velocity field and odor cues in space, resulting from direct numerical simulations of the turbulent flow in a channel of length L, width W and height H (also see Figure 1—video 1). Air flows from left to right at a mean speed $U_b$ and hits a cylindrical obstacle that generates turbulence. The height of the obstacle is $H/4$ and tunes the intensity of turbulent fluctuations relative to the mean velocity. To characterize the flow we show in Figure 1d that the mean velocity profile, for $z^{+}\gtrsim$ 30, follows the law of the wall $\frac{⟨U⟩}{u_{\tau }}=\frac{1}{\kappa }\ln z^{+}+B$, where $\kappa$ and $B$ are constants and $z^{+}=\frac{z}{{\delta }_{\nu }}$, recovering classical statistics for channel turbulence Pope, 1984. The odor field is emitted from a concentrated source downstream from the obstacle; it develops as a meandering filament that fluctuates as it travels downstream and soon breaks into discrete pockets of odor (whiffs) separated by odor-less stretches (blanks) (Figure 1b–c and f). The Kolomogorov scaling for the spectra of odor fluctuations $k^{-5/3}$ holds for $k\eta \lesssim 0.1$ (see Figure 1e), consistent with previous experimental results in channel flow (see e.g. Saddoughi and Veeravalli, 1994). Typical time courses of the odor are shown in Figure 1f. Note that depending on the sampling location, odor may be more or less sparse (compare for example Figure 1f left and right). All parameters and methods are summarized respectively in Table 1 and in Materials and Methods.

Figure 1 with 1 supplement see all

Download asset Open asset

Do odor cues bear information about source location meters away from the source? To answer this question, we develop supervised machine learning algorithms that learn the relationship between the input (odor) and the distance from the source (output) from a large dataset of examples. In order to dissect what are the best predictors of source location and how ranking depends on the statistics of the odor, we need to detail more specifically the input and output of the algorithm.

To design the input we start with the odor concentration field ${\displaystyle c(\mathbf{z},t)}$ which varies stochastically in space and time as a result of turbulent transport. Here, $\mathbf{z}=(z_1,z_2,z_3)$ is a location in the three dimensional space and $t$ is time. We focus on a plane at a fixed height, and consider the conical region where odor can be detected, the ‘cone of detection’ (Figure 2a). We first compose time series of the odor field; each time series is indicated with ${\mathbf{c}}_i$ and consists of the odor sampled at $M$ equally spaced times with frequency $\omega$ at a discrete location ${\mathbf{z}}_i$ within the cone of detection. Thus, each time series is a vector ${\mathbf{c}}_i=(c({\mathbf{z}}_i,t_i),\mathrm{\dots },c({\mathbf{z}}_i,t_{i+M}))$, where $t_{i+M}-t_i=M/\omega$ is the temporal span of the time series, or memory. From each time series, ${\mathbf{c}}_i$ we calculate five features $x_i^1,\mathrm{\dots },x_i^5$, where $x_i^1$ is the temporal average of the concentration during whiffs in the time series ${\displaystyle {\mathbf{c}}_i}$; $x_i^2$ is its average slope (time derivative of odor upon detection, averaged across whiffs within ${\mathbf{c}}_i$); $x_i^3$ is the average duration of blanks (stretches of time when odor is below detection within ${\mathbf{c}}_i$); $x_i^4$ is the average duration of whiffs (stretches of time when odor is above threshold within ${\mathbf{c}}_i$); and $x_i^5$ is the intermittency factor (the fraction of time the time series ${\mathbf{c}}_i$ is above threshold). The detection threshold is defined adaptively as discussed in Materials and methods and Figure 2—figure supplement 1. Features $x^1$ and $x^2$ depend explicitly on odor concentration, whereas features $x^3$, $x^4$ and $x^5$ only depend on when the odor is on or off, but not on its intensity. To remark this difference, we refer to $x^1$ and $x^2$ as intensity features, and $x^3$, $x^4$ and $x^5$ as timing features. Our input ${\mathbf{x}}_i=(x_i^1,\mathrm{\dots },x_i^d)$ is composed of d-dimensional vectors of features and we will focus on $d=1,2,5$. We seek to infer distance from the source, thus our output $y$ is the coordinate of the sampling point $\mathbf{z}$ in the downwind direction, i.e. $y=z_1$, with the source placed at the origin (see sketch in Figure 2a). We refer to the figure supplements for results in the crosswind direction, $y=z_2$. We train the algorithm by providing $N$ examples of input-output pairs $({\mathbf{x}}_i,y_i)$ selected randomly from the full simulation, and obtain the function that connects input and output: $y\approx f(\mathbf{x})$.

Figure 2 with 6 supplements see all

Download asset Open asset

We propose a machine learning approach where the different odor features are ranked based on their predictive power, rather than their fitting properties. Different data-sets of odor/distance pairs are defined. The data-sets differ in the way odor measurements are represented in terms of feature vectors. For each data-set we learn a function to predict the distance to target given the corresponding odor features. The predictive power of each function, and corresponding set of features, is then assessed. More precisely, each data-set is split in a training and a test set, as custom in machine learning. Training sets are used to learn functions connecting odor to target location, whereas test sets are used to assess their prediction properties. The training/test split is crucial since the goal is to make good predictions on new, unseen points, that are not within the training sets. From a modeling perspective, a flexible nonlinear/nonparametric approach based on kernel methods is contrasted and shown to be superior to a simpler linear model (Figure 2—figure supplement 2). A careful protocol based on hold-out cross-validation is used to select the hyper-parameters of the considered learning models (we refer to Materials and methods and Figure 2—figure supplement 3 for more details).

To illustrate the results we pick the two-dimensional plane at height $H/4$ that contains the source. The first result is that individual features ($d=1$) bear useful information for two-dimensional source localization even at several meters from the source. Performance is quantified by the normalized squared error averaged over the $N_t$ points in the test set $\chi ={\sum }_{i=1}^{N_t}{[y_i-f({\mathbf{x}}_i)]}^2/{\sum }_{i=1}^{N_t}{[y_i-\overline{y}]}^2$. For this dataset, intensity features rank higher than timing features (Figure 2b–c), consistent with previous work Atema, 1996 and predictions are more accurate in the crosswind than in the downwind direction (compare with Figure 2—figure supplement 4). For reference, a random guess with flat probability within the correct lower and upper bounds yields ${\chi }_{\text{random}}=2$, whereas a target function $f_{\text{trivial}}(x)={⟨y⟩}_{\text{test}}$ that learns the average of the output over the test set yields ${\chi }_{\text{trivial}}=1$.

In order to prove that the algorithm captures the dynamics of the scalar and not of the underlying velocity field, we realize three computational fluid dynamics simulations with the odor source at different locations downstream of the obstacle. Each source is placed at a different distance from the obstacle and thus feels different velocity fields, because the flow is inhomogeneous in the downstream direction. To perform a fair comparison we train the algorithm on points sampled from a conical domain based on the most downstream source; for the other two sources, we use an identical domain shifted upstream so that the vertex of the cone is located at the respective source location (see Figure 2—figure supplement 5 left). Performance at source height varies little over the three locations, demonstrating that the algorithm is learning the dynamics of the scalar and not of the carrying flow (see Figure 2—figure supplement 5 center). Additionally we demonstrate that, irrespective of source position, timing features acquire predictive power with height, whereas the opposite is observed for intensity features (see Figure 2—figure supplement 5 right). We discuss this property in further detail in the following. . Finally, we train the algorithm with the odor fields emanating from one source location, and test its performance over odor fields generated by the two other sources. Figure 2—figure supplement 6 shows that performance of individual features varies little when test and training are performed over different dataset. Learning from pairs of features appears somewhat more sensitive to the details of the dataset. In the aggregate, the analysis corroborates that the algorithm captures odor dynamics and is rather insensitive to details of the underlying flow. In the rest of the manuscript we focus on the leftmost source, so as to exploit the full spatial range available from the simulations.

Next we analyze whether and how the sampling strategy affects performance and ranking of the features. Most results are shown for a memory of $100{\tau }_{\eta }\approx 15s$. Performance improves with longer memory (Figure 3a), because this allows to better average out noise and obtain more stable estimates of the features. But improvement follows a slow power law so that waiting for example 20 times longer yields predictions only about twice as precise. On the other hand, waiting as little as $10{\tau }_{\eta }\approx 1.5$ seconds still allows to make predictions, albeit less precise. We then verify whether performance may improve with a larger training set. Because we infer distance from an individual (scalar) feature, the problem is one dimensional and we find that a small number of training points, which we indicate with $N$, is sufficient to reach a plateau in prediction performance (Figure 3b). We choose $N=5000$ training points, which is also robust to the case with more than one feature (Figure 3—figure supplement 1). Finally, sampling more frequently than once per Kolmogorov time does not essentially affect the results nor ranking (Figure 3c). Similar results hold for the crosswind direction (Figure 3—figure supplement 2).

Figure 3 with 2 supplements see all

Download asset Open asset

Pairing two observables improves performance in some cases, but not always. In fact, pairing two features of the same category results in little to no improvement (Figure 4 and similarly for the crosswind direction, Figure 4—figure supplement 1). In contrast, combining one intensity and one timing feature improves performance considerably, up to 65%. This result can be understood by mapping the error done by individual features in space (Figure 5—figure supplement 1), showing that intensity and timing features are complementary, that is intensity features perform well in locations where timing features perform poorly.

Figure 4 with 1 supplement see all

Download asset Open asset

We next seek to clarify whether the results depend on space. To this end we compose five different dataset, a to e, obtained by extracting odor snapshots from horizontal planes at source height (b), above the source (c to e), and below the source (a) (Figure 5a). From a to e, sparsity increases and intensity decreases (Figure 5b) simply because closer to the boundary, the air slows down and the odor accumulates. By analyzing performance across these dataset, we find that ranking of individual features shifts considerably. The two intensity features outperform all timing features when the dataset is not very sparse (dataset a-b, Figure 5c and d left). In contrast, two timing features (intermittency factor and blank duration) outperform all others for the more sparse and less intense dataset d-e (Figure 5c and d right). Whiff duration performs poorly in d-e because intermittency is too severe and whiffs are short in duration thus bear little information (the average whiff duration is 1–7 time steps in over 90% of the time series). Although the ranking of individual features shifts with height, pairing one intensity and one timing feature remains the most successful strategy across all heights (Figure 5c and d). In contrast, combining all five features contributes little improvement (Figure 5—figure supplement 2).

Figure 5 with 3 supplements see all

Download asset Open asset

Let us now focus on the plane at source height and separate locations based on their distance from the source. We assemble a distal dataset and a proximal dataset, composed of points that are further and closer than 2330η from the source respectively (Figure 6a). The odor is more intense and more sparse closer to the source and it becomes more dilute and less sparse with distance from the source (Figure 6b). Performance of individual features degrades with distance (Figure 6d). Intensity features clearly outperform timing features at close range, as seen both from various percentiles of the test error (Figure 6d, left) as well as the full distribution (Figure 6c, left). The disparity between timing and intensity features disappears in the distal problem: the error distribution for all individual features is essentially superimposed except for the tails (Figure 6c, right and inset), which cause small differences in the median and other percentiles of the error (Figure 6d, right). Remarkably, mixed pairs outperform all individual features in both the distal and proximal problems (Figure 6c–d). In the aggregate, results demonstrate that, even within a single turbulent flow, ranking shifts considerably. Namely, measuring timing of odor encounters is most useful in regions where the odor is dilute, that is far from the source and from the substrate, whereas measuring intensity is most useful in concentrated conditions, that is close to the source or the substrate.

Figure 6

Download asset Open asset

Thus, the predictive power of different features varies greatly in space and time. Next, we show how this spatial variation is dictated by the statistical properties of the odor plume. To this end, we provide an analytical characterization of the test error of individual features, $\chi$, that connects directly to the physics of the problem. Such a characterization is consistent with the observed variation in performance. Note that for large enough samples, the test error $\chi$ approximates the following expected error

The latter is called expected error and is the ideal version of the test error. It can be interpreted as the error summed over all possible input-output pairs, weighted by their corresponding joint probability to be sampled. Here, $f^{*}(x)=⟨y|x⟩={\int }_0^{R}yp(y|x)dy$ is the so called regression function, which minimizes the expected error among all possible functions. The regression function can be approximated using Kernel ridge regression and sufficiently rich kernels. Indeed, kernel ridge regression is known to be a so called universal estimator Hastie et al., 2001; Steinwart, 2002. In the above expression, $R$ is the length of the cone, $p(y|x)$ is the conditional probability distribution of the output given the input, $p(x,y)$ is the joint probability distribution of the input output pair ${\displaystyle (x,y)}$ is the prior distribution on the output $y$. The idea is to relate $f^{*}$, hence the expected error, to the distribution $p(x|y)$ of observing feature $x$ depending on distance $y$. Indeed, the latter is dictated by the fluid dynamics of odor plumes from a concentrated source, and hence provides a more direct connection between the expected error and the physics of the problem. Since the considered features are sample averages, in the limit of large samples, their distribution is well approximated by a Gaussian, hence fully characterized by the mean and standard deviation. Then simple estimates can be computed empirically from our data. Equation (1), provided with these estimates, reproduces the predictive power of individual features showed in Figure 5c (see Figure 5—figure supplement 3). Note that whiffs deviate considerably from a normal distribution, hence the argument needs to be revised for this feature. To move beyond empirical estimates and extend the above reasoning to other flows and generic combinations of (possibly non Gaussian) features, a generalization of the asymptotic arguments proposed in Celani et al., 2014 is needed.

Our results demonstrate that within the cone of detection, the time course of an odor bears useful information for source localization even at meters from the source. We find that the concentration and the slope of a turbulent odor signal, averaged over a memory lag, are particularly useful to predict source location at close range or near the boundary. These features quantify the intensity of the odor and its variation. The primacy of the intensity features wanes in more challenging conditions, for example moving away from the source or away from the boundary. In these portions of space, where the odor is scarcer, features that quantify timing of odor detection become as effective as intensity features, or more effective. One of the best studied example of olfactory search in dilute conditions is arguably the case of insects.

Interestingly, olfactory receptor neurons in insects appear to encode efficiently information about timing across a wide range of intensities Gorur-Shandilya et al., 2017; Martelli et al., 2013.

Note that while the statistics of an odor plume clearly depends on all details of the flow and the source, see e.g. Justus et al., 2002; Celani et al., 2014; Fackrell and Robins, 1982, here we keep all of these parameters constant and demonstrate that even within a single flow, odor dynamics and the best predictors vary considerably in space. This begs the next question: do organisms switch between different modalities depending on attributes of odor dynamics, which will vary in space? This could be the case for mice, where the neural activity in the first relay of olfactory processing does in fact depend on how sparse is the odor Lewis et al., 2021. Specifically, sparse odor cues elicit individual responses that follow closely the ups and downs of the odor in time. In contrast, continuous signals elicit intense responses which are however uncorrelated to the temporal dynamics of the odor itself Lewis et al., 2021.

We find that features within the same class are redundant whereas features from different classes are complementary. Indeed, features of the same class have similar patterns of performance in space, but each class has a distinct pattern. As a consequence, measuring both timing and intensity is beneficial, but using more than one feature to quantify either timing or intensity provides no advantage. Combining all features does not improve over the performance of mixed pairs, consistent with redundancy within each class. Note that there is no fundamental reason to expect features from the same class to be redundant, and further work with a larger library of features is needed to prove or disprove this notion.

Importantly, mixed time/intensity pairs of features outrank individual features robustly, that is in all portions of space, regardless of distance from the source and from the ground. This is in contrast with individual features and suggests relying on simultaneous timing and intensity features is advantageous when odors are sensed at various distances from the source and from the substrate. Interestingly, the coexistence of bursting olfactory neurons and canonical olfactory neurons in lobsters suggests these animals are in fact able to measure simultaneously timing and intensity Park et al., 2014; Ache et al., 2016, which is consistent with the increased predictive power of the mixed pairs of features. Similarly, in mammals, optogenetic activation of the olfactory bulb Smear et al., 2013 demonstrates that both kinds of measures guide behavior (lick vs no lick).

In this work, we have investigated the problem of predicting the location of a target from measures of the time course of a turbulent odor. Previous work explored a related question, that is how to best represent instantaneous snapshots of the odor to encode maximum information about source location Victor et al., 2019. The two approaches are not immediately comparable: first, Victor et al., 2019 consider few snapshots of the odor, rather than measures of its time course. Second, maximizing information does not guarantee good predictions (to make predictions information needs to be extracted and processed, and importantly the focus is on new data that were not previously seen). We provide two comments that are relevant if information is the limiting factor for prediction accuracy: (i) binary representations were suboptimal in all conditions considered in Victor et al., 2019; Boie et al., 2018, that is at few tens of cm from the source. This is consistent with our results in concentrated conditions, where timing features -accessible through binary representations- are suboptimal. Our evidences suggest, however, that the result may not hold in more dilute conditions, where the gap between binary and more accurate representations should become increasingly small. (ii) Individual snapshots of odor from Victor et al., 2019; Boie et al., 2018 contained 1–2 bits of information about source location, but allocating more resources to represent how the odor varies in time was found informative Victor et al., 2019; Boie et al., 2018. Our mixed pairs of features at close range achieve precisions of 5–6%, corresponding to coding for position with words of 4–4.3 bits. Our results thus confirm that memory is indeed useful, but the gain does not increase indefinitely with further memory.

The literature on olfactory navigation is vast. Although a complete review of available algorithms is beyond the scope of the present work, we remark that recent results investigated gradient descent algorithms using either concentration alone Gire et al., 2016, or various measures of timing and intensity Park et al., 2016; Michaelis et al., 2020; Leathers et al., 2020. Overall, both intensity and timing appear to have a potential to lead to an odor source, consistent with our results on individual features. A combination of the two kinds of features was found beneficial in Leathers et al., 2020, consistent with our results on mixed pairs. Note that predicting source location and navigating to reach it are distinct tasks. Although they are often assumed to be intimately connected, whether good predictors may be good variables for navigation in more general contexts remains to be understood.

Here, we have analyzed the features that enable the most accurate prediction of source location. We add a few observations about the significance of the results for animal behavior. First: whether animals rely on features from either class will depend on what features best support behavior. It is often implicitly assumed that features that bear reliable information on source location are also the most useful for navigation. However, this connection between prediction and navigation is far from straightforward and more work is needed to establish whether accurate predictions imply efficient navigation. Second: animals are unlikely to have prior information on the details of the odor source, for example its intensity. Timing features are more robust than intensity features with respect to the intensity of the source and may thus be favored regardless of their performance, which was argued in Schmuker et al., 2016. In our work, timing features are precisely invariant with source intensity because we define the detection threshold adaptively (see Materials and methods). More realistic conditions will need to be evaluated, where dependence on source intensity emerges as a result of non-linearities that we did not model in this work. These effects emerge for example, close to a boundary which partially absorbs the odor Gorur-Shandilya et al., 2019, or in the case of fixed thresholds, although this dependence is weak in the far field where timing features are most useful Celani et al., 2014. Third: we have focused on predicting source location from within the cone of detection, where an agent will detect the odor quite often. However, a crucial difficulty of turbulent navigation is to find the cone itself. We cannot address the problem of predicting source location from outside the cone because detections are so rare that we lack statistics. The distinction between inside and outside the cone of detection is key for navigation with sparse cues (see Reddy et al., 2021) and deserves further attention.

Simulations of odor transport are generated through NEK5000 (43), freely available from Argonne National Laboratory (https://nek5000.mcs.anl.gov/). Outputs from DNS presented in Figure 1 are processed to extract time series and compute the five features described in the text (average concentration, slope, blank duration, whiff duration and intermittency factor). These data are available at the online repository https://osf.io/ja9xr/. The results of kernel ridge regression performed on these data are presented in Figures 2-6. We perform kernel ridge regression on these data with the freely available code FALKON (https://github.com/LCSL/FALKON_paper, copy archived at swh:1:rev:480741cf1e7da0d1d7415309cd6f254080a6ca17).

1. 1. Ache BW
   2. Hein AM
   3. Bobkov YV
   4. Principe JC

   (2016) Smelling time: a neural basis for olfactory scene analysis

   Trends in Neurosciences 39:649–655.

   https://doi.org/10.1016/j.tins.2016.08.002

   • PubMed
   • Google Scholar
2. 1. Ackels T
   2. Erskine A
   3. Dasgupta D
   4. Marin AC
   5. Warner TPA
   6. Tootoonian S
   7. Fukunaga I
   8. Harris JJ
   9. Schaefer AT

   (2021) Fast odour dynamics are encoded in the olfactory system and guide behaviour

   Nature 593:558–563.

   https://doi.org/10.1038/s41586-021-03514-2

   • PubMed
   • Google Scholar
3. 1. Atema J

   (1996) Eddy chemotaxis and odor landscapes: exploration of nature with animal sensors

   The Biological Bulletin 191:129–138.

   https://doi.org/10.2307/1543074

   • PubMed
   • Google Scholar
4. 1. Baker KL
   2. Dickinson M
   3. Findley TM
   4. Gire DH
   5. Louis M
   6. Suver MP
   7. Verhagen JV
   8. Nagel KI
   9. Smear MC

   (2018) Algorithms for olfactory search across species

   The Journal of Neuroscience 38:9383–9389.

   https://doi.org/10.1523/JNEUROSCI.1668-18.2018

   • PubMed
   • Google Scholar
5. 1. Basil J
   2. Atema J

   (1994) Lobster orientation in turbulent odor plumes: simultaneous measurement of tracking behavior and temporal odor patterns

   The Biological Bulletin 187:272–273.

   https://doi.org/10.1086/BBLv187n2p272

   • PubMed
   • Google Scholar
6. 1. Boie SD
   2. Connor EG
   3. McHugh M
   4. Nagel KI
   5. Ermentrout GB
   6. Crimaldi JP
   7. Victor JD

   (2018) Information-theoretic analysis of realistic odor plumes: What cues are useful for determining location?

   PLOS Computational Biology 14:e1006275.

   https://doi.org/10.1371/journal.pcbi.1006275

   • PubMed
   • Google Scholar
7. 1. Brown SL
   2. Joseph J
   3. Stopfer M

   (2005) Encoding a temporally structured stimulus with a temporally structured neural representation

   Nature Neuroscience 8:1568–1576.

   https://doi.org/10.1038/nn1559

   • PubMed
   • Google Scholar
8. 1. Catania KC

   (2013) Stereo and serial sniffing guide navigation to an odour source in a mammal

   Nature Communications 4:1441.

   https://doi.org/10.1038/ncomms2444

   • PubMed
   • Google Scholar
9. 1. Celani A
   2. Villermaux E
   3. Vergassola M

   (2014) Odor landscapes in turbulent environments

   Physical Review X 4:041015.

   https://doi.org/10.1103/PhysRevX.4.041015

   • Google Scholar
10. 1. Demir M
    2. Kadakia N
    3. Anderson HD
    4. Clark DA
    5. Emonet T

    (2020) Walking Drosophila navigate complex plumes using stochastic decisions biased by the timing of odor encounters

    eLife 9:e57524.

    https://doi.org/10.7554/eLife.57524

    • PubMed
    • Google Scholar
11. 1. Duplat J
    2. Jouary A
    3. Villermaux E

    (2010) Entanglement rules for random mixtures

    Physical Review Letters 105:034504.

    https://doi.org/10.1103/PhysRevLett.105.034504

    • PubMed
    • Google Scholar
12. 1. Fackrell JE
    2. Robins AG

    (1982) Concentration fluctuations and fluxes in plumes from point sources in a turbulent boundary layer

    Journal of Fluid Mechanics 117:1–26.

    https://doi.org/10.1017/S0022112082001499

    • Google Scholar
13. 1. Falkovich G
    2. Gawȩdzki K
    3. Vergassola M

    (2001) Particles and fields in fluid turbulence

    Reviews of Modern Physics 73:913–975.

    https://doi.org/10.1103/RevModPhys.73.913

    • Google Scholar
14. 1. Findley TM
    2. Wyrick DG
    3. Cramer JL
    4. Brown MA
    5. Holcomb B
    6. Attey R
    7. Yeh D
    8. Monasevitch E
    9. Nouboussi N
    10. Cullen I
    11. Songco JO
    12. King JF
    13. Ahmadian Y
    14. Smear MC

    (2021) Sniff-synchronized, gradient-guided olfactory search by freely moving mice

    eLife 10:e58523.

    https://doi.org/10.7554/eLife.58523

    • PubMed
    • Google Scholar
15. 1. Fischer PF
    2. Loth F
    3. Lee SE
    4. Lee SW
    5. Smith DS
    6. Bassiouny HS

    (2007) Simulation of high-Reynolds number vascular flows

    Computer Methods in Applied Mechanics and Engineering 196:3049–3060.

    https://doi.org/10.1016/j.cma.2006.10.015

    • Google Scholar
16. Software

    1. Fischer PF
    2. Lottes JW
    3. Kerkemeier J

    (2008) Nek5000, version v19.0

    NEK.

    http://nek5000.mcs.anl.gov
17. 1. Gardiner JM
    2. Atema J

    (2010) The function of bilateral odor arrival time differences in olfactory orientation of sharks

    Current Biology 20:1187–1191.

    https://doi.org/10.1016/j.cub.2010.04.053

    • PubMed
    • Google Scholar
18. 1. Gire DH
    2. Kapoor V
    3. Arrighi-Allisan A
    4. Seminara A
    5. Murthy VN

    (2016) Mice develop efficient strategies for foraging and navigation using complex natural stimuli

    Current Biology 26:1261–1273.

    https://doi.org/10.1016/j.cub.2016.03.040

    • PubMed
    • Google Scholar
19. 1. Gorur-Shandilya S
    2. Demir M
    3. Long J
    4. Clark DA
    5. Emonet T

    (2017) Olfactory receptor neurons use gain control and complementary kinetics to encode intermittent odorant stimuli

    eLife 6:e27670.

    https://doi.org/10.7554/eLife.27670

    • PubMed
    • Google Scholar
20. 1. Gorur-Shandilya S
    2. Martelli C
    3. Demir M
    4. Emonet T

    (2019) Controlling and measuring dynamic odorant stimuli in the laboratory

    The Journal of Experimental Biology 222:207787.

    https://doi.org/10.1242/jeb.207787

    • PubMed
    • Google Scholar
21. Book

    1. Hastie T
    2. Friedman J
    3. Tibshirani R

    (2001) The Elements of Statistical Learning

    New York, NY: Springer.

    https://doi.org/10.1007/978-0-387-21606-5

    • Google Scholar
22. Thesis

    1. Ho L

    (1989)

    A Legendre spectral element method for simulation of incompressible unsteady viscous free-surface flows

    Massachusetts Institute of Technology, Cambridge, USA.

    • Google Scholar
23. 1. Jacob V
    2. Monsempès C
    3. Rospars J-P
    4. Masson J-B
    5. Lucas P
    6. Morozov AV

    (2017) Olfactory coding in the turbulent realm

    PLOS Computational Biology 13:e1005870.

    https://doi.org/10.1371/journal.pcbi.1005870

    • Google Scholar
24. 1. Justus KA
    2. Murlis J
    3. Jones C
    4. Cardé RT

    (2002) Measurement of odor-plume structure in a wind tunnel using a photoionization detector and a tracer gas

    Environmental Fluid Mechanics 2:115–142.

    https://doi.org/10.1023/A:1016227601019

    • Google Scholar
25. 1. Kennedy JS
    2. Marsh D

    (1974) Pheromone-regulated anemotaxis in flying moths

    Science 184:999–1001.

    https://doi.org/10.1126/science.184.4140.999

    • PubMed
    • Google Scholar
26. 1. Leathers KW
    2. Michaelis BT
    3. Reidenbach MA

    (2020) Interpreting the spatial-temporal structure of turbulent chemical plumes utilized in odor tracking by lobsters

    Fluids 5:82.

    https://doi.org/10.3390/fluids5020082

    • Google Scholar
27. 1. Lewis SM
    2. Xu L
    3. Rigolli N
    4. Tariq MF
    5. Suarez LM
    6. Stern M
    7. Seminara A
    8. Gire DH

    (2021) Plume dynamics structure the spatiotemporal activity of mitral/tufted cell networks in the mouse olfactory bulb

    Frontiers in Cellular Neuroscience 15:633757.

    https://doi.org/10.3389/fncel.2021.633757

    • PubMed
    • Google Scholar
28. 1. Mafra-Neto A
    2. Cardé RT

    (1994) Fine-scale structure of pheromone plumes modulates upwind orientation of flying moths

    Nature 369:142–144.

    https://doi.org/10.1038/369142a0

    • Google Scholar
29. 1. Martelli C
    2. Carlson JR
    3. Emonet T

    (2013) Intensity invariant dynamics and odor-specific latencies in olfactory receptor neuron response

    The Journal of Neuroscience 33:6285–6297.

    https://doi.org/10.1523/JNEUROSCI.0426-12.2013

    • PubMed
    • Google Scholar
30. 1. Michaelis BT
    2. Leathers KW
    3. Bobkov YV
    4. Ache BW
    5. Principe JC
    6. Baharloo R
    7. Park IM
    8. Reidenbach MA

    (2020) Odor tracking in aquatic organisms: the importance of temporal and spatial intermittency of the turbulent plume

    Scientific Reports 10:7961.

    https://doi.org/10.1038/s41598-020-64766-y

    • PubMed
    • Google Scholar
31. 1. Moore P
    2. Crimaldi J

    (2004) Odor landscapes and animal behavior: tracking odor plumes in different physical worlds

    Journal of Marine Systems 49:55–64.

    https://doi.org/10.1016/j.jmarsys.2003.05.005

    • Google Scholar
32. 1. Murlis J
    2. Elkinton JS
    3. Cardé RT

    (1992) Odor plumes and how insects use them

    Annual Review of Entomology 37:505–532.

    https://doi.org/10.1146/annurev.en.37.010192.002445

    • Google Scholar
33. 1. Nagel KI
    2. Wilson RI

    (2011) Biophysical mechanisms underlying olfactory receptor neuron dynamics

    Nature Neuroscience 14:208–216.

    https://doi.org/10.1038/nn.2725

    • PubMed
    • Google Scholar
34. 1. Orszag SA

    (1980) Spectral methods for problems in complex geometries

    Journal of Computational Physics 37:70–92.

    https://doi.org/10.1016/0021-9991(80)90005-4

    • Google Scholar
35. 1. Parabucki A
    2. Bizer A
    3. Morris G
    4. Munoz AE
    5. Bala ADS
    6. Smear M
    7. Shusterman R

    (2019) Odor concentration change coding in the olfactory bulb

    ENeuro 6:ENEURO.0396-18.2019.

    https://doi.org/10.1523/ENEURO.0396-18.2019

    • PubMed
    • Google Scholar
36. 1. Park IM
    2. Bobkov YV
    3. Ache BW
    4. Príncipe JC

    (2014) Intermittency coding in the primary olfactory system: A neural substrate for olfactory scene analysis

    The Journal of Neuroscience 34:941–952.

    https://doi.org/10.1523/JNEUROSCI.2204-13.2014

    • PubMed
    • Google Scholar
37. 1. Park IJ
    2. Hein AM
    3. Bobkov YV
    4. Reidenbach MA
    5. Ache BW
    6. Principe JC

    (2016) Neurally encoding time for olfactory navigation

    PLOS Computational Biology 12:e1004682.

    https://doi.org/10.1371/journal.pcbi.1004682

    • PubMed
    • Google Scholar
38. 1. Patera AT

    (1984) A spectral element method for fluid dynamics: Laminar flow in A channel expansion

    Journal of Computational Physics 54:468–488.

    https://doi.org/10.1016/0021-9991(84)90128-1

    • Google Scholar
39. Book

    1. Pope SB

    (1984) Turbulent Flows

    Cambridge: Cambridge University Press.

    https://doi.org/10.1017/CBO9780511840531

    • Google Scholar
40. Preprint

    1. Reddy G
    2. Shraiman BI
    3. Vergassola M

    (2021) Sector search strategies for odor trail tracking

    bioRxiv.

    https://doi.org/10.1101/2021.03.03.433838

    • Google Scholar
41. 1. Riffell JA
    2. Shlizerman E
    3. Sanders E
    4. Abrell L
    5. Medina B
    6. Hinterwirth AJ
    7. Kutz JN

    (2014) Sensory biology flower discrimination by pollinators in a dynamic chemical environment

    Science 344:1515–1518.

    https://doi.org/10.1126/science.1251041

    • PubMed
    • Google Scholar
42. Software

    1. Rigolli N

    (2022) FALKON_paper, version swh:1:rev:480741cf1e7da0d1d7415309cd6f254080a6ca17

    Software Heritage.

    https://archive.softwareheritage.org/swh:1:dir:b152d94d711b88f99c206601afb1235de15321eb;origin=https://github.com/LCSL/FALKON_paper;visit=swh:1:snp:5032915c66d96288fedec074afe8d025600fca3b;anchor=swh:1:rev:480741cf1e7da0d1d7415309cd6f254080a6ca17
43. Conference

    1. Rudi A
    2. Carratino L
    3. Rosasco L

    (2018) Neural information processing systems

    Proceedings of the 31st International Conference on Neural Information Processing Systems.

    https://dl.acm.org/doi/10.5555/3294996.3295145

    • Google Scholar
44. 1. Saddoughi SG
    2. Veeravalli SV

    (1994) Local isotropy in turbulent boundary layers at high Reynolds number

    Journal of Fluid Mechanics 268:333–372.

    https://doi.org/10.1017/S0022112094001370

    • Google Scholar
45. 1. Schmuker M
    2. Bahr V
    3. Huerta R

    (2016) Exploiting plume structure to decode gas source distance using metal-oxide gas sensors

    Sensors and Actuators B 235:636–646.

    https://doi.org/10.1016/j.snb.2016.05.098

    • Google Scholar
46. Book

    1. Schölkopf B
    2. Smola A

    (2002)

    Learning with Kernels

    Cambridge, MA, USA: MIT Press.

    • Google Scholar
47. 1. Shraiman B
    2. Siggia E

    (2000) Scalar turbulence

    Nature 405:639–646.

    https://doi.org/10.1038/35015000

    • PubMed
    • Google Scholar
48. 1. Smear M
    2. Resulaj A
    3. Zhang J
    4. Bozza T
    5. Rinberg D

    (2013) Multiple perceptible signals from a single olfactory glomerulus

    Nature Neuroscience 16:1687–1691.

    https://doi.org/10.1038/nn.3519

    • PubMed
    • Google Scholar
49. 1. Steck K
    2. Veit D
    3. Grandy R
    4. Badia SBI
    5. Badia SBI
    6. Mathews Z
    7. Verschure P
    8. Hansson BS
    9. Knaden M

    (2012) A high-throughput behavioral paradigm for Drosophila olfaction - The Flywalk

    Scientific Reports 2:361.

    https://doi.org/10.1038/srep00361

    • PubMed
    • Google Scholar
50. 1. Steinwart I

    (2002) Support vector machines are universally consistent

    Journal of Complexity 18:768–791.

    https://doi.org/10.1006/jcom.2002.0642

    • Google Scholar
51. 1. van Breugel F
    2. Dickinson MH

    (2014) Plume-tracking behavior of flying Drosophila emerges from a set of distinct sensory-motor reflexes

    Current Biology 24:274–286.

    https://doi.org/10.1016/j.cub.2013.12.023

    • PubMed
    • Google Scholar
52. 1. Vergassola M
    2. Villermaux E
    3. Shraiman BI

    (2007) “Infotaxis” as a strategy for searching without gradients

    Nature 445:406–409.

    https://doi.org/10.1038/nature05464

    • PubMed
    • Google Scholar
53. 1. Vickers NJ

    (2000) Mechanisms of animal navigation in odor plumes

    The Biological Bulletin 198:203–212.

    https://doi.org/10.2307/1542524

    • PubMed
    • Google Scholar
54. 1. Vickers NJ
    2. Christensen TA
    3. Baker TC
    4. Hildebrand JG

    (2001) Odour-plume dynamics influence the brain’s olfactory code

    Nature 410:466–470.

    https://doi.org/10.1038/35068559

    • PubMed
    • Google Scholar
55. 1. Victor JD
    2. Boie SD
    3. Connor EG
    4. Crimaldi JP
    5. Ermentrout GB
    6. Nagel KI

    (2019) Olfactory navigation and the receptor nonlinearity

    The Journal of Neuroscience 39:3713–3727.

    https://doi.org/10.1523/JNEUROSCI.2512-18.2019

    • PubMed
    • Google Scholar

Author details

1. Nicola Rigolli

   1. Department of Physics, University of Genova, Genova, Italy
   2. Institut de Physique de Nice, Université Côte d’Azur, Centre National de la Recherche Scientifique, Nice, France
   3. National Institute of Nuclear Physics, Genova, Italy
   4. MalGa, Department of Civil, Chemical and Environmental Engineering, University of Genoa, Genoa, Italy

   Contribution

   Conceptualization, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

   For correspondence

   nicola.rigolli@edu.unige.it

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-0734-2105

2. Nicodemo Magnoli

   1. Department of Physics, University of Genova, Genova, Italy
   2. National Institute of Nuclear Physics, Genova, Italy

   Contribution

   Supervision, Investigation, Project administration, Writing – review and editing

   Competing interests

   No competing interests declared

3. Lorenzo Rosasco

   MaLGa, Department of computer science, bioengineering, robotics and systems engineering, University of Genova, Genova, Italy

   Contribution

   Software, Supervision, Funding acquisition, Methodology, Project administration, Writing – review and editing

   Competing interests

   No competing interests declared

4. Agnese Seminara

   1. Institut de Physique de Nice, Université Côte d’Azur, Centre National de la Recherche Scientifique, Nice, France
   2. MalGa, Department of Civil, Chemical and Environmental Engineering, University of Genoa, Genoa, Italy

   Contribution

   Conceptualization, Software, Formal analysis, Supervision, Funding acquisition, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   agnese.seminara@unige.it

   Competing interests

   Reviewing editor, eLife

   "This ORCID iD identifies the author of this article:" 0000-0001-5633-8180

• 848

  views

• 159

  downloads

• 17

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.