Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 43862.
Published online 2017 March 8. doi:  10.1038/srep43862
PMCID: PMC5341060

Emergence of linguistic laws in human voice


Linguistic laws constitute one of the quantitative cornerstones of modern cognitive sciences and have been routinely investigated in written corpora, or in the equivalent transcription of oral corpora. This means that inferences of statistical patterns of language in acoustics are biased by the arbitrary, language-dependent segmentation of the signal, and virtually precludes the possibility of making comparative studies between human voice and other animal communication systems. Here we bridge this gap by proposing a method that allows to measure such patterns in acoustic signals of arbitrary origin, without needs to have access to the language corpus underneath. The method has been applied to sixteen different human languages, recovering successfully some well-known laws of human communication at timescales even below the phoneme and finding yet another link between complexity and criticality in a biological system. These methods further pave the way for new comparative studies in animal communication or the analysis of signals of unknown code.

The main objective of quantitative linguistics is to explore the emergence of statistical patterns (often called linguistic laws) in language and general communication systems (see refs 1,2 for a review). The most celebrated of such regularities is Zipf’s law describing the uneven abundance of word frequencies3,4. This law presents many variations in human language5,6,7,8 but also shows ubiquity9 in many linguistic scales10, has been claimed to be universal11,12 and has consequences for syntax and symbolic reference13. On the other hand Heaps’ law, also called Herdan’s law14,15 states that the vocabulary of a text grows allometrically with the text length16,17, and is mathematically connected with Zipf’s law18,19,20, the scaling exponent being dependent on both Zipf law and the vocabulary size. Finally the Zipf’s Law of abbreviation (or brevity law for short) is the statistical tendency of more frequent elements in communication systems to be shorter or smaller4,21 and has been recently claimed as an universal trend derived from fundamental principles of information processing22. As such, this statistical regularity holds also phonetically23,24 and implies that the higher the frequency of a word, the shorter its length or duration, probably caused by a principle of compression25,26, although this is a general pattern that can change depending on other acoustical factors like noise27, pressure to communicate at long distances calls28 or communicative efficiency4 and energetic constraints29.

Linguistic laws extend beyond written language and have been shown to hold for different biological data7,30. According to some authors31, the presence of scaling laws in communication is indicative of the existence of processes taking place across different cognitive scales. Interpreting linguistic laws as scaling laws which emerged in communication systems32,33 actually opens the door for speculating on the existence of underlying scale-invariant physical laws operating underneath34. Of course, in order to explore the presence or absence of such patterns one needs to directly study the acoustic corpus, i.e. human voice, as every linguistic phenomenon or candidate for language law can be camouflaged or diluted after the transcription process. Notwithstanding the deep relevance of linguistic laws reported in written texts, we still wonder which of these laws found in written corpora are related or derive from more fundamental properties of the acoustics of language -and are thus candidates for full-fledged linguistic laws-, or emerge as an artifact of scripture codification.

Linguistic Laws in Acoustics?

Acoustic communication is fully determined by three physical magnitudes extracted from the signals: frequency, energy and time35,36. The configuration space induced by these magnitudes results from an intrinsic evolutionary relationship between the production and perception of sound systems37 that further shapes the range of hearing and producing sounds for a variety of life forms35,38. Since animals use their acoustic abilities both to monitor their environment and to communicate we should expect that natural selection has in some sense optimized these sensorial capabilities, but, despite the great differences that evolution has involved for different animals that communicate acoustically, there are many similarities between their mechanisms of sound production and hearing39. Focusing on primates, traditionally language has been distinguished from vocalizations of nonhuman primates only at a qualitative, semantic level40. Interestingly, it is well known that children use statistical cues to segment the input41,42,43 and probably share with non-human primates some of these mechanisms, albeit with some differences44. The discovery of these statistical learning abilities has boosted a new approach to the study of language45,46 and suggests that statistical learning of language could be based on patterns or more generally linguistic laws22: research on language acquisition shows that higher frequency facilitates learning, and Zipf’s law tell us that vocabulary learning is easier than expected a priori given the skewness of the distribution, for instance. We advance that, as will be shown below, these patterns are already present in the physical sound waves produced by human voice even at levels below those generally considered linguistically significant, i.e., below the phoneme timescale.

Empirical evidence of robust linguistic laws holding in written texts across different human languages has been reported many times (see refs 2, 7 and 19 and references therein), and it has been shown that these laws are not fully observed in random texts48. Studies with oral corpus are however much less abundant, and they systematically imply a transcription of the acoustical waves into words in the case of human speech or some ill-defined analog of words in the case of animal communication, as the main segments to analyze49,50. A few current efforts take a different road and consider other possible written units such as lemmas51 or compare written and oral production for some linguistic patterns52,53, in general showing that frequencies of elements in written corpora can be taken as a rough estimate for their frequency in spoken language54.

All in all, the exploration of linguistic laws in oral corpora is scarce. In fact, all linguistic studies in oral and written corpora are influenced by our segmentation decisions and our definition of ‘word’, intimately biased by an inherently anthropocentric perspective and, of course, by our linguistic tradition. The idea that speech is produced like writing, as a linear sequence of vowels and consonants may indeed be a relic of our scripture technology. As a matter of fact, it is well-known in linguistics that both vowels and consonants are produced linearly but also depend on their surrounding elements: this is the traditional and well-studied concept of coarticulation55. The boundaries between acoustical elements are therefore difficult to identify if we are not native speakers of a language, and that’s yet a crucial problem of phonotactics and speech segmentation and recognition56 with differences across languages57.

Precisely because of this, classical signal segmentation based on the concept of “word” inherited from writing has lead some researchers to search how to transform artificially written corpora into phoneme or syllabic chains with different objectives58, at the time it involves two major problems in communication studies, namely (i) the impossibility of performing fully objective comparative studies between human and non-human signals59, where signals could be physical events, behaviours or structures to which receivers respond60. This problem leads researchers sometimes to manually segment acoustic signals guided only by their expertise, and prevents to explore signals of unknown origin including for instance the search for possible extraterrestrial intelligence61. And (ii) a rather arbitrary definition of the units of study guided by ortographic conventions already produces non-negligible epistemological problems at the core of Linguistics62,63.

In this work we explore the acoustic analog of classical linguistic laws in the statistics of speech waveform from extensive real databases64 that for the first time extend all the way into the intraphoneme range (t < 10−2 s)65. In order to do so in a systematic way, in what follows we will present a methodology that enables the direct linguistic analysis of any acoustical waveform without needs to invoke to ad hoc codes or assume any particular formal communication system underneath. This method only makes use of physical magnitudes of the signal such as energy or time duration which therefore allow for a non-ambiguous definition of our units. Speech and voice are indeed physical phenomena and as such in this work we interpret them as a concatenation of events of energy release, and propose a mathematically well-defined way for segmenting in time this orchestrated suite of energy release avalanches47. We will show clear evidence of robust Zipf, Heaps and brevity laws emerging in this context and speculate that this might be due to the fact that human voice seems to be operating close to a critical state, hence finding an example of a biological system that, driven by evolution, has linked complexity and criticality. We expect that this methodology can open a fresh line of research in communication systems where a direct exploration of underlying statistical patterns in acoustic signals is possible without needs to predetermine any of the aforementioned non-physical concepts, and hope that this will allow researchers to develop comparative studies between human language and other acoustical communication systems or even to unravel whether if a generic signal shares these patterns2.


For this work we have used two databases: the main one is a TV broadcast speech database named KALAKA-264, and additionally a complementary database named LRE has been used to confirm results on a larger set of languages. Originally designed for language recognition evaluation purposes, the former consists of wide-band TV broadcast speech recordings (4 hours per language sampled using 2 bytes at a rate of 16 kHz) ranging six different languages: Basque, Catalan, Galician, Spanish, Portuguese and English and encompassing both planned and spontaneous speech throughout diverse environmental conditions, such as studio or outside journalist reports but excluding telephonic channels. The second one complements KALAKA with an additional set of 12 languages, taken from the NIST Language Recognition Evaluation (LRE′96) corpus and including Japanese, Vietnamese, Mandarin, Korean, Arabic, Hindi or Tamil to cite some (see SI for details).

The Method

The objects under study, speech waveforms or otherwise any generic acoustic signal, are fully described by an amplitude time series A(t) (see Fig. 1 for an illustration of the method). In order to unambiguously extract a sequence of elements -the equivalent to words and phrases- from such signal without the need to perform ad hoc segmentation, we start by considering the semi-definite positive magnitude ε(t) = |A(t)|2 which, dropping irrelevant constants, has physical units of energy per time (SI for additional details on speech waveform statistics). By defining an energy threshold Θ in this signal47 we will unambiguously separate voice events (above threshold) from silence events (below threshold). More concretely, Θ is defined as a relative percentage of the signal and its actual value in energy units depends on the signal variability range: for example Θ = 80% means that 20% of the data falls under this energy level. It has been shown that Θ decimates the signal similarly to a real-space Renormalization Group (RG) transformation47,66, in such a way that increasing Θ induces a flow in RG space. Systems operating close to a critical state lie in an unstable fixed point of this RG flow and its associated signal statistics are therefore shown to be threshold-invariant. Now, Θ not only works as an energy threshold that filters out background or environmental noise (noise filtering being a key aspect that species have learned to perform27,67) but, as previously stated, enable us to unambiguously define what we call a token or voice event, that is, a sequence of consecutive measurements ε(t) > Θ, from a silence event of duration τ. Each token is in turn characterized by a duple (Ev, Tv) where Tv is the duration of the event and Ev corresponds to the total energy released during that event An external file that holds a picture, illustration, etc.
Object name is srep43862-m1.jpg obtained summing up the instantaneous energy over the duration of the event. Accordingly, the signal is effectively transformed to an ordered sequence of tokens {(Ev(i), Tv(i))}, each of these being separated by silence events of highly heterogeneous durations τ which, incidentally, are known to be power law distributed47. Finally, by linearly binning the scale of integrated energies we can assign an energy label (the bin) to each token, hence mapping the initial acoustic signal into a symbolic sequence of fundamental units which we call types. Note that two tokens whose integrated energy fall within the range of the same energy bin are mapped to the same type even if their duration can be different, so in principle several tokens could map into the same type (see SI for a table of Type/Token ratios). The linear binning of the energy scale has a bin size b = 0.01, in such a way that the first bin agglutinates tokens with energies between E0 and E0 + b, the second bin between E0 + b and E0 + 2b, and so on (note that results reported later will be robust against changes in b). The set of bins can be understood as an abstraction of a universal language vocabulary, and accordingly some bins might be empty and in general each bin will occur with uneven frequencies. As such, types can be understood as acoustically-based universal abstractions of a fundamental unit, an abstract version of words or phonemes that appear intertwined in a signal with characteristic patterns.

Figure 1
This figure illustrates the methodology to extract a sequence of types from an acoustic signal.

To synthesize, with this methodology we are able to map an arbitrary acoustic signal into a sequence of types separated by silence events (Fig. 1). Standard linguistic laws can then be directly explored in acoustic signals without needs to have an a priori knowledge neither of the signal code nor of the adequate segmentation process or the particular syntax of the language underlying the signal. This protocol is thus independent of the communication system and can be used to make unbiased comparisons across different systems and signals. Needless to say, results could in principle depend on the particular value of Θ, as this scans the signals at different energy thresholds. However human voice has been recently shown47 to be invariant under changes in Θ -an evidence of self-organized criticality (SOC)68 in this system- and, accordingly, parameter-free laws can be extracted using a proper collapse theory as it will be shown in the results section. Also, in order to guarantee that the emergence of linguistic laws is only due to the structure and correlations of the signal and not due to the process of symbolization we will compare the results obtained from speech signals to properly defined null models which randomize the signal ε(t). These null models thus maintain the marginal instantaneous energy distribution and remove any other correlation structure, yielding non-Gaussian white noise with a fat-tailed marginal distribution.

It should be highlighted that, according to the abovementioned protocol, the types extracted do not share in principle a direct relation with real words or with any other linguistic units, but only a formal one as these types constitute a base of fundamental acoustic units of any oral communication system. As a matter of fact, we consider this possibility to be unlikely, as our basic units of study, while ranging several scales, are in most of the cases smaller than the typical timescales of words. Further work should explore to what extent there is a connection between tokens, or sequences of them, and phonemes, words or other linguistic units. Notwithstanding, this a priori lack of matching is indeed required and desired if one aims at developing a method which is universally applicable in communication signals of arbitrary origin, not necessarily oral communication based on words. Such is the case, for instance, in the realm of comparative analysis across different species. In any case, note that our method relies on retrieving voice events via energy thresholding. As many of the detected events typically belong to the intraphoneme range, one could argue that the method focus on the rich composition of voice and silence events within a given linguistic unit (phonemes or words).

Despite not being the purpose of this work, this method may be seen as an unsupervised Voice Activity Detection (VAD) algorithm. It is well known that VAD algorithms based on temporal domain analysis or on energy features are highly affected by noise but note that VAD systems are typically optimized through reference transcriptions, manually annotated at the word level69,70. The purpose of our approach is however not to perform word or phoneme level segmentation, but rather to explore organisational patterns of language-independent acoustic units taking place at different scales in oral communication, and quantify them via classical linguistic analysis. That being said, assessing how informative are those cues in (unsupervised) word segmentation or speech recognition tasks71,72,73 is an interesting question that will deserve further investigation.


Gutenberg-Richter law

The energy E released during voice events is a direct measure of the vocal fold response function under air pressure perturbations, and its distribution PΘ(E) could in principle depend both on the threshold Θ and on the language under study. In the inset of Fig. 2 we observe that PΘ(E) is power law distributed over about six decases, saturated by an exponential cut-off. This distribution has been interpreted before as the analogue of a Gutenberg-Richter law in voice, as the precise shape of energy release fluctuations during voice production parallel those occurring in earthquakes47. As increasing Θ induces a flow in RG space, systems which lie close to a critical point (unstable fixed point in RG space) show scale invariance under Θ and hence the distributions can be collapsed into a Θ-independent shape, thereby eliminating the trivial dependence on Θ. This has been shown to be the case for human voice and accordingly (technical details can be found in the SI) we can express the collapsed energy distribution as P(E) = EϕAn external file that holds a picture, illustration, etc.
Object name is srep43862-m2.jpg(E/Eξ) for An external file that holds a picture, illustration, etc.
Object name is srep43862-m3.jpg, where El is the lower limit beyond which this law is fulfilled, An external file that holds a picture, illustration, etc.
Object name is srep43862-m4.jpg is a scaling function and the relevant variable is ϕ, the scaling exponent. In order to collapse every curve PΘ(E) the theory predicts to rescale E  Eleft angle bracketEright angle bracket/left angle bracketE2right angle bracket, PΘ(E)  PΘleft angle bracketE2right angle bracket2/left angle bracketEright angle bracket3. In the outset of Fig. 2 we show the result of this analysis for the case of Spanish language, where we find that ϕ  1.15, this exponent being approximately language-independent (see Table 1 for other languages and SI for additional details). Interestingly, these exponents are compatible with those found in rainfall, another natural system that has been shown to be compatible with SOC dynamics74, and cannot be explained by simple null models47. While the precise values of the exponents are not extremely relevant, one should take these with caution as thresholding is known to introduce in some cases finite size biases in the form of a second scaling region for short sizes which is only removed when the series is long enough75. In what follows we explore the emergence of classical linguistic laws in these acoustic signals.

Figure 2
Log-log plot of the collapsed shape (see the text) and threshold-independent energy release distribution P(E) in the case of Spanish language (KALAKA database) for several thresholds, after logarithmic binning.
Table 1
Summary of scaling exponents associated to the energy release distribution (ϕ), Zipf’s law (ζ), Heaps’ law (α) and Brevity law (β) for six different languages (KALAKA database).

Zipf’s Law

The illustrious George Kingsley Zipf formulated a statistical observation which is popularly known as Zipf’s law76. In its original formulation3, it establishes that in a sizable sample of language the number of different words (vocabulary) N(n) which occur exactly n times decays as N(n)~nζ, where the exponent ζ varies from text to text77 but is usually close to 2. An alternative and perhaps more common formulation of this law4 is defined in terms of the rank, such that if words are ranked in decreasing order by their frequency of appearance then the number of occurrences of words with a given rank r goes like n(r)~rz, where it is easy to see that both exponents are related via An external file that holds a picture, illustration, etc.
Object name is srep43862-m5.jpg and thus z approximates to 116,78. Here for convenience we make use of the former and explore NΘ(n) applied to the statistics of types. Interestingly, while Zipf’s (and Heaps’s) laws were originally proposed to quantify order in the word arrangement on written texts, these laws indeed describe the arrangement of informational units that go beyond words in texts (see for instance7,11) and there exist several theoretical justifications for the presence of these laws which do not require resorting to linguistics but to more fundamental concepts such as principles of compression in information theory12,26 which don’t necessarily require the units to adhere to any particular type of communication.

As a matter of fact, a mathematical relation between the energy release statistics and Zipf’s law can be found. First, note that the frequency of type r is proportional to

An external file that holds a picture, illustration, etc.
Object name is srep43862-m6.jpg

where the asymptotic approximation is easily found expanding An external file that holds a picture, illustration, etc.
Object name is srep43862-m7.jpg up to second order. Now, since P(E) is monotonically decreasing, then there is a formal equivalence between type r and rank r: the most frequent type corresponds to type ‘1’ and thus has rank 1, the second most frequent type is ‘2’, and so on: type with label r will have rank r, as frequencies of types are naturally ordered in a monotonically decreasing way by construction. Therefore, the frequency of rank r is proportional to eq. 1. Accordingly An external file that holds a picture, illustration, etc.
Object name is srep43862-m8.jpg so we predict z  ϕ and thus An external file that holds a picture, illustration, etc.
Object name is srep43862-m9.jpg.We now consider experimental results.

Again NΘ(N) could in principle depend on the threshold but assuming that the signal complies with the scale-invariance mentioned above, one can collapse all threshold-dependent curves into a universal shape and thus remove any dependence on this parameter by rescaling n  n/LV, NΘ(n)  NΘ(n)VL where V is the total number of different types present in the signal and L as the total number of tokens (see SI for technical details). Results are shown in the case of Basque language in Fig. 3, where a clear threshold-independent decaying power law emerges with a scaling exponent ζ  1.77. Analogous results with compatible exponents for other languages can be found in SI and Table 1. The relation between exponents ζ and ϕ is on good agreement with our theoretical prediction. Null models systematically deviate from these results, and neither display the characteristic power law decay nor any invariance under variation of the energy threshold (SI).

Figure 3
Log-log plot of Zipf’s law for the case of Basque language (KALAKA database, see SI for equivalent results in other languages), for different thresholds Θ.

Heaps’ law

Together with Zipf’s law and connected mathematically (see ref. 2 and references therein), the second classical linguistic law is Heaps’ law, the sublinear growth of the number of different words V in a text with text size L (measured in total number of words): V ~ Lα, α < 114,79 (a constant rate for appearance of new words leads to α = 1). Here the vocabulary V is defined as the total number of different types that appear in the signal, whereas L is defined as the total number of tokens found for a given threshold. Results are shown for a specific language in Fig. 4 (see SI for the rest). In the outset panel we present the collapsed, threshold-independent curves, where again we find a scaling law with an effective exponent α′ related to the original exponent α = α/(1 + α). In this case equivalent computation on the null model yield a Heaps law with the trivial exponent α   1 (SI).

Figure 4
Log-log plot of the Heaps’ law for the Portuguese sample (KALAKA database) and several thresholds.

These results are quantitatively consistent with previous results on written texts16. In particular, several authors18 point out that, at least asymptotically, the relation ζ = 1 + α holds with good approximation, and this is on reasonably good agreement with our findings in human voice as well. Interestingly, a recent work80 has found that, as opposed to Indoeuropean (alphabetically based) languages, Zipf’s law breaks down and Heaps’ law reduces to the trivial case for written texts in Chinese, Japanese, Korean and other logosyllabic languages. Applying our methodology in the database of such logosyllabic languages (Japanese, Mandarin, Korean) our results are at odds with80 as we still find evidence of Zipf and Heaps laws holding for these languages in terms of the energetic voice fluctuations. We conclude that those differences found when written texts are analyzed do not arise when we study directly human voice.

Brevity Law

The tendency of more frequent words to be shorter3,4,21 can be generalized as the tendency of more frequent elements to be shorter or smaller, and its origin has been suggested to be related to optimization and information compression arguments in connection with other linguistic laws81. In acoustics, spontaneous speech indeed tends to obey this law after text segmentation82, and has been found also in other non-human primates28,81. Here we can test brevity law in essentially two different ways. First, note that voice events (tokens) map into types according to the logarithmic binning of their associated energy, hence voice events with different duration might yield the same type as previously noted. Thus for each type we can compute its mean duration t averaging over all voice events that fall within that type, and then plot the histogram MΘ(t) that describes frequency of each type versus its mean duration. Brevity law would require MΘ(t) to be a monotonically decreasing function. These results are shown for a particular language in log-log scales in Fig. 5, finding initially a power law decaying relation which is indicative of a brevity law (results are again found to be language independent, see SI for additional results). The inset provide the threshold-dependent distributions and the outset panel provides the collapsed, threshold-independent shape M(t) (see Table 1 for scaling exponents). Again in this case results in null models deviate from such behavior (SI) and are clearly different from the random typing12. Alternatively, one can also directly observe the duration frequency at the level of voice events, finding similar results (see SI).

Figure 5
Log-log plot of the Brevity law in the case of English (KALAKA database), for several thresholds. In the inner panel we plot, for different thresholds, the histogram MΘ(t) that describes the relative frequency of a type of mean duration t.


In this work we have explored the equivalent of linguistic laws directly in acoustic signals. We have found that human voice manifests the analog of classical linguistic laws found in written texts (Zipf’s law, Heaps’ law and the brevity law or law of abbreviation). These laws are found to be invariant under variation of the energy threshold Θ, and can be collapsed under universal functions accordingly. As Θ is the only free parameter of the method, this invariance determines that the results are not afflicted by ambiguities associated to arbitrarily defining unit boundaries. Results are robust across a list of 16 different languages (indoeuropean and non-indoeuropean, including some logosyllabic), across timescales (extending all the way into the intraphoneme range, where no cognitive effects operate, and invariant under energy threshold variation) and across conversational modes (one or more speakers, both planned and spontaneous speech throughout diverse environment conditions, such as studio, outside journalist reports, and telephonic channel). Interestingly, an equivalent analysis performed on null models defined by randomizing the signal ε(t) (yielding white noise with the same instantaneous energy distribution of the original signal) fail to reproduce this phenomenology (SI). The concrete range of exponents found for both Zipf and Heaps laws are compatible between each other and somewhat similar -but not identical- to the typical ones observed in the literature for written texts2,7,76,80, whereas to the best of our knowledge this is the first observation of scaling behavior with a clear exponent in the case of brevity law in speech. Actually, our finding of a power law in brevity law differs from the case of random typing where a power law doesn’t conform12.

The specific and complex alternation of air stops (silences) intertwined with voice production are at the core of the microscopic voice fluctuations. During voice production, acoustic communication is governed by the so-called biphasic cycle (breath and glottal cycle, see ref. 83 for a review) that together with some other acoustic considerations (pitch period, voice onset time, the relation between duration, stress and syllabic structure82) determines the microscopic structure of human voice, including silence stops. However, these timescales are in general very large: as previously stated, this current study focus and scans voice properties even at intraphonemic timescales, where the statistical laws of language emerge directly from the physical magnitudes that govern acoustic communication. Our results therefore open the possibility of speculating whether the fact that these laws have been found in upper levels of human communication might be a result of a scaling process and a byproduct of the physics rather than derived from the choice of the typical units of study on the analysis of written corpus (phonemes, syllabus, words, …), like differences between analysis of Indoeuropean and logosyllabic languages demonstrates80. As a matter of fact, in a previous work human voice has been framed within self-organized criticality (SOC), speculating that the fractal structure of lungs drives human voice close to a critical state47, this mechanism being ultimately responsible for the microscopic self-similar fluctuations of the signal. This constitutes a new example of the emergence of SOC in a physiological system, different in principle from the classical one found in neuronal activity84. In the same line, note that the set of critical exponents characterizing the SOC nature of human voice have been shown to be very similar to the ones found in another SOC system: rainfall74. The celebrated theory of critical phenomena thoroughly explains that many physical processes, when poised close to critical point, essentially reduce to the same phenomenon despite of having different physical origins, and can be described by a set of critical exponents that classifies the core phenomenon in universality classes. As rainfall and voice production share the same critical exponent for energy release, this suggests that in essence a similar mechanism could be triggering both types of threshold-like dynamics.

In any case, one could thus speculate that the emergence of linguistic scaling laws is just a consequence of the presence of SOC, what in turn would support the physical origin of linguistic laws. From an evolutionary viewpoint, under this latter perspective human voice, understood as a communication system which has been optimized under evolutionary pressures, would constitute an example where complexity (described in terms of robust linguistic laws) emerges when a system is driven close to criticality, something reminiscent of the celebrated edge of chaos hypothesis85. Furthermore, under this interpretation the onset of complexity in language (as reported by linguistic laws in written texts) could be a direct byproduct of similar patterns already emanating at purely physiological grounds (the physics of energy fluctuation in voice production): in this respect the role played by cognitive processes would be unclear, as the timescales involved in the physiological mechanisms of voice production (intraphonemic) are typically shorter than those associated to cognitive ones. More succinctly, the classical explanation for the origin of scaling laws in language in terms of scaling phenomena in cognitive processes could be challenged by this new evidence. Certainly, we are not claiming that statistical laws of language are only a consequence of human physiology, among other reasons because the clearest evidence leads us to consider the presence of these laws in texts2,7, but the plausible causal link between scaling laws in acoustics and the analogous ones at the language level is suggestive and should be studied in depth.

On more practical grounds, the method used and proposed here also addresses the longstanding problem of signal segmentation. It has been acknowledged that there is no such thing as a ‘correct’ segmentation strategy86. In written corpus white space is usually taken as an easy marker of the separation between words, however this is far from evident in continuous speech where separation between words or concepts is technologically harder to detect, conceptually vague and probably ill-defined. Few exceptions that used oral corpus for animal communication still require to define ad hoc segmentation algorithms28, or manual segmentation strategies which usually give an arbitrary or overestimated segmenting times or windows81, what might even raise epistemological questions. As such, this segmentation problem unfortunately has prevented wider, comparative studies in areas such as animal communication or the search for signs of possible extraterrestrial intelligence from radio signals (in this line only few proposals have been made61). By varying the energy threshold the method presented here automatically partitions and symbolizes the signal at various energy scales, providing a recipe to establish an automatic, general and systematic way of segmenting and thus enabling comparison of across acoustic signals of arbitrary origin for which we may lack the syntax, code or exact size of its constituents.

To round off, we hope that this work paves the way for new research avenues in comparative studies. Open questions that deserve further work abound; just to name a few: in the light of this new method, what can we say about the acoustic structure in other animal communication systems? Can we find evidence of universal traits in communication that do not depend on a particular species but are only physically and physiologically constrained, or on the other hand are linguistic universals a myth87? How these laws evolve with aging6,88? Are they affected by cognitive or pulmonary diseases? What is the precise relation between SOC and linguistic laws in this context? And in particular, can we find mathematical evidence of a minimal, analytically tractable SOC model that produce these patterns? These and other questions are interesting avenues for future work.

Additional Information

How to cite this article: Torre, I. G. et al. Emergence of linguistic laws in human voice. Sci. Rep. 7, 43862; doi: 10.1038/srep43862 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:


We acknowledge financial support from ONCE foundation (IGT), grant 2014SGR890 (MACDA) from AGAUR, Generalitat de Catalunya and the APCOM Project TIN2014-57226-P (AHF), grant FIS2013-41057P from Ministerio de Economia y Competitividad (BL and IGT) and an EPSRC fellowship EP/P01660X/1 (LL). BL acknowledges the hospitality and support of Queen Mary University of London, where part of this research was developed, and a Salvador de Madariaga fellowship. We also thank Ramon Ferrer-i-Cancho and Alvaro Corral for helpful comments, as well as the editor and anonymous referees.


The authors declare no competing financial interests.

Author Contributions B.L. and A.H.F. designed research, I.G.T. and J.L. performed data processing, I.G.T. performed data analysis, B.L. and L.L. contributed with theoretical analysis, all authors contributed with discussions and A.H.F. and L.L. wrote the manuscript.


  • Kohler R., Altmann G. & Piotrowski R. Quantitative Linguistics (De Gruyter Mouton, Berlin, 2008).
  • Altmann E. G. & Gerlach M. Statistical Laws in Linguistics. In Creativity and Universality in Language, Lecture Notes in Morphogenesis (eds Degli Esposti M., Altmann E. & Pachet F.) 7–26 (Springer, Cham, 2016).
  • Zipf G. K. The psycho-biology of language (Houghton, Mifflin, 1935).
  • Zipf G. K. Human behavior and the principle of least effort (Addison-Wesley Press, 1949).
  • Ferrer i Cancho R. The variation of Zipf’s law in human language. Eur. Phys. J. B 44, 249–257 (2005).
  • Baixeries J., Elvevag B. & Ferrer i Cancho R. The evolution of the exponent of Zipf’s law in language ontogeny. PloS One 8, e53227 (2013). [PMC free article] [PubMed]
  • Piantadosi S. T. Zipf’s word frequency law in natural language: A critical review and future directions. Psychon. Bull. & Rev. 21, 1112–1130 (2014). [PMC free article] [PubMed]
  • van Egmond M., van Ewijk L. & Avrutin S. Zipf’s law in non-fluent aphasia. J. Quant. Linguist. 22, 233–249 (2015).
  • Li W. Zipf’s law everywhere. Glottom. 5, 14–21 (2002).
  • Ha L., Sicilia-Garcia E., Ming J. & Smith F. Extension of Zipf’s law to words and phrases in Proc. 19th Int. Conf. Comp. Ling. 1, 1–6 (Ass. Comp. Ling., 2002).
  • Corominas-Murtra B. & Solé R. Universality of Zipf’s law. Phys. Rev. E. 82, 011102 (2010). [PubMed]
  • Ferrer i Cancho R. Compression and the origins of Zipf’s law for word frequencies. Complexity 21, 409–411 (2016).
  • Ferrer i Cancho R., Riordan O. & Bollobás B. The consequences of Zipf’s law for syntax and symbolic reference. Proc. R. Soc. Lond. B: Biol. Sci. 272, 561–565 (2005). [PMC free article] [PubMed]
  • Herdan G. Quantitative linguistics (Butterworth, 1964).
  • Heaps H. Information retrieval: Computational and theoretical aspects (Acad. Press, 1978).
  • Font-Clos F., Boleda G. & Corral A. A scaling law beyond Zipf’s law and its relation to heaps’ law. New J. Phys. 15, 093033 (2013).
  • Gerlach M. & Altmann E. Scaling laws and fluctuations in the statistics of word frequencies. New J. Phys. 16, 113010 (2014).
  • Mandelbrot B. On the theory of word frequencies and on related markovian models of discourse. Struc. Lang. Mathem. Asp. 12, 190–219 (1961).
  • Baayen H. Word frequency distributions 18 (Springer Sci. & Business Media, 2001).
  • Font-Clos F. & Corral A. Log-log convexity of type-token growth in Zipf’s systems. Phys. Rev. Let. 114, 238701 (2015). [PubMed]
  • Grzybek P. Contributions to the science of text and language: word length studies and related issues 31 (Springer Sci. & Business Media, 2006).
  • Bentz C. & Ferrer i Cancho R. Zipf’s law of abbreviation as a language universal (Lorentz Center Workshop, Leiden, 2015).
  • Aylett M. & Turk A. Language redundancy predicts syllabic duration and the spectral characteristics of vocalic syllable nuclei. J. Acoust. Soc. Am. 119, 3048–3058 (2006). [PubMed]
  • Tomaschek F., Wieling M., Arnold D. & Baayen H. Word frequency, vowel length and vowel quality in speech production: an EMA study of the importance of experience. In Interspeech 1302–1306 (2013).
  • Ferrer i Cancho R., Bentz C. & Seguin C. Compression and the origins of Zipf’s law of abbreviation. Preprint at: (2015).
  • Ferrer i Cancho R. et al. . Compression as a universal principle of animal behavior. Cognitive Sci. 37, 1565–1578 (2013). [PubMed]
  • Brumm H. Animal communication and noise (Springer, 2013).
  • Ferrer i Cancho R. & Hernández-Fernández A. The failure of the law of brevity in two new world primates. Statistical caveats. Glottotheory 4, 45–55 (2013).
  • Gillooly J. & Ophir A. The energetic basis of acoustic communication. Proc. Roy. Soc. Lon. B: Biol. Sci. 277(1686), 1325–1331 (2010).
  • Schwab D. J., Nemenman I. & Mehta P. Zipf’s law and criticality in multivariate data without fine-tuning. Phys. Rev. Let. 113, 068102 (2014). [PMC free article] [PubMed]
  • Kello C. et al. . Scaling laws in cognitive sciences. Trends Cog. Sci. 14, 223–232 (2010). [PubMed]
  • Ferrer i Cancho R. & Sole R. Least effort and the origins of scaling in human language. Proc. Natl. Acad. Sci. USA 100, 788–791 (2003). [PubMed]
  • Nowak M. & Krakauer D. The evolution of language. Proc. Natl. Acad. Sci. USA 96, 8028–8033 (1999). [PubMed]
  • Chater N. & Brown G. Scale-invariance as a unifying psychological principle. Cognition 69, B17–B24 (1999). [PubMed]
  • Sueur J. & Drosopoulos S. Insect species and their songs. In Insect Sounds and Communication: Physiology, Behaviour, Ecology, and Evolution 207–217 (Taylor & Francis, Boca Raton, 2006).
  • Saposhkov M. Electroacustica (Reverte, 1983).
  • MacNeilage P. The evolution of phonology (Oxford Univ. Press, 2011).
  • Berg R. & Stork D. The Physics of Sound (Prentice Hall, 1995).
  • Fletcher N. Animal bioacoustics. In Springer Handbook of Acoustics 821–841 (Springer, 2014).
  • Fitch T. The evolution of speech: a comparative review. Trends Cog. Sci. 4, 258–267 (2000). [PubMed]
  • Saffran J., Aslin R. & Newport E. Statistical learning by 8-month-old infants. Science 274, 1926–1928 (1996). [PubMed]
  • Kuhl P. et al. . Phonetic learning as a pathway to language: new data and native language magnet theory expanded (nlm-e). Philos. Trans. Roy. Soc. Lon. B: Biol. Sci. 363, 979–1000 (2008). [PMC free article] [PubMed]
  • Romberg A. & Saffran J. Statistical learning and language acquisition. Wiley Inter. Rev.: Cog. Sci. 1, 906–914 (2010). [PMC free article] [PubMed]
  • Saffran J. et al. . Grammatical pattern learning by human infants and cotton-top tamarin monkeys. Cognition 107, 479–500 (2008). [PMC free article] [PubMed]
  • Kuhl P. A new view of language acquisition. Proc. Natl. Acad. Sci. USA 97, 11850–11857 (2000). [PubMed]
  • Emberson L. & Rubinstein D. Statistical learning is constrained to less abstract patterns in complex sensory input (but not the least). Cognition 153, 63–78 (2016). [PMC free article] [PubMed]
  • Luque J., Luque B. & Lacasa L. Scaling and universality in the human voice. J. R. Soc. Interface 12, 20141344 (2015). [PMC free article] [PubMed]
  • Ferrer-i Cancho R. & Elvevag B. Random texts do not exhibit the real Zipf’s law-like rank distribution. PLoS One 5, e9411 (2010). [PMC free article] [PubMed]
  • McCowan B., Hanser S. & Doyle L. Quantitative tools for comparing animal communication systems: information theory applied to bottlenose dolphin whistle repertoires. Anim. behav. 57, 409–419 (1999). [PubMed]
  • Ferrer i Cancho R. & McCowan B. A law of word meaning in dolphin whistle types. Entropy 11, 688–701 (2009).
  • Corral A., Boleda G. & Ferrer i Cancho R. Zipf’s law for word frequencies: Word forms versus lemmas in long texts. PloS One 10, e0129031 (2015). [PMC free article] [PubMed]
  • Nabeshima T. & Gunji Y. Zipf’s law in phonograms and Weibull distribution in ideograms: comparison of English with Japanese. Biosystems 73, 131–139 (2004). [PubMed]
  • Baroni M. Distributions in text. In Corpus linguistics: An international handbook (eds Lüdeling A. & Kytö M.) 803–821 (Mouton de Gruyter, Berlin, 2008).
  • Samlowski B., Möbius B. & Wagner P. Comparing syllable frequencies in corpora of written and spoken language. In Proc. Interspeech 637–640 (2011).
  • Farnetani E. & Recasens D. Coarticulation and Connected Speech Processes. In The Handbook of Phonetic Sciences 316–352 (Blackwell Pub., 2010).
  • Glass J. A probabilistic framework for segment-based speech recognition. Comp. Speech & Lang. 17, 137–152 (2003).
  • Tyler M. & Cutler A. Cross-language differences in cue use for speech segmentation. J. Acoust. Soc. Am. 126, 367–376 (2009). [PubMed]
  • Taylor P. Text-to-Speech Synthesis (Cambridge Univ. Press, 2009).
  • Kuhl P. Human speech and birdsong: Communication and the social brain. Proc. Natl. Acad. Sci. USA 100, 9645–9646 (2003). [PubMed]
  • Stegmann U. Animal communication theory: information and influence (Cambridge Univ. Press, 2013).
  • Doyle L., McCowan B., Johnston S. & Hanser S. Information theory, animal communication, and the search for extraterrestrial intelligence. Acta Astronaut., 68, 406–417 (2011).
  • Bunge M. What is pseudoscience? The Skeptical Inquirer 9, 36–46 (1984).
  • Kohler R. Synergetic linguistics. In Quantitative linguistics 760774 (de Gruyter, Berlin, 2005).
  • Rodriguez-Fuentes L., Peñagarikano M., Bordel G., Varona A. & Diez M. Kalaka: A TV broadcast speech database for the evaluation of language recognition systems. In Proc. 7th Int. Conf. LREC 1678–1685 (Valleta, Malta, 2010).
  • Crystal T. & House A. Segmental durations in connected-speech signals: Current results. J. Acoust. Soc. Am. 83, 1553–1573 (1988).
  • Corral A. Point-occurrence self-similarity in crackling-noise systems and in other complex systems. J. Stat. Mech.- Theory E. 2009(01), P01022 (2009).
  • Brumm H. & Slabbekoorn H. Acoustic communication in noise. Adv. St. Behav. 35, 151–209 (2005).
  • Bak P. How Nature Works (Copernicus, 1996).
  • Juang B.-H. & Rabiner L. Automatic speech recognition-a brief history of the technology development (Georgia Ins. of Tech., Atlanta Rutgers Univ. & Univ. of California, Santa Barbara, 2005).
  • Drugman T., Stylianou Y., Kida Y. & Akamine M. Voice Activity Detection: Merging Source and Filter-based Information, IEEE Signal Process. Let. 23, 252–256 (2016).
  • Lamel F., Rabiner R., Rosenberg E. & Wilpon G. An improved endpoint detector for isolated word recognition, IEEE Trans. Acoust., Speech, Signal Process., 29, 777–785 (1981).
  • Park A. S. & Glass J. R. Unsupervised pattern discovery in speech. IEEE Transactions on Audio, Speech, and Language Processing 16(1), 186–197 (2008).
  • Räsänen O., Gabriel D. & Michael C. Frank. Unsupervised word discovery from speech using automatic segmentation into syllable-like units. Interspeech Dresden, Germany, pp., 3204–3208 (2015).
  • Peters O., Deluca A., Corral A., Neelin J. & Holloway C. Universality of rain event size distributions. J. Stat. Mech.- Theory E. 2010(11), P11030 (2010).
  • Font-Clos F., Pruessner G., Deluca A. & Moloney N. R. The perils of thresholding. New Journal of Physics 17, 043066 (2015).
  • Moreno-Sanchez I., Font-Clos F. & Corral A. Large-scale analysis of zip’s law in english texts. PLoS One 11, 1–19 (2016).
  • Zanette D. Statistical patterns in written language. Preprint at: arXiv, 1412.3336 (2014).
  • Ferrer i Cancho R. & Hernández-Fernández A. Power laws and the golden number In Problems of general, germanic and slavic linguistics 518–523 (Books-XXI, 2008).
  • Egghe L. Untangling Herdan’s law and Heaps’ law: Mathematical and informetric arguments. J. Am. Soc. Inf. Sci. & Tech. 58, 702–709 (2007).
  • Lü L., Zhang Z. & Zhou T. Deviation of Zipf’s and Heaps’ laws in human languages with limited dictionary sizes. Sci. Rep. 3, 1082 (2013). [PMC free article] [PubMed]
  • Gustison M. & Semple S., Ferrer i Cancho, R. & Bergman, T. Gelada vocal sequences follow Menzerath’s linguistic law. Proc. Natl. Acad. Sci. USA 113, E2750–E2758 (2016). [PubMed]
  • Greenberg S., Carvey H., Hitchcock L. & Chang S. Temporal properties of spontaneous speech-a syllable-centric perspective. J. Phonetics 31, 465–485 (2003).
  • MacNeilage P. The Origin of Speech (Oxford Univ. Press, 2008).
  • Kello C. Critical branching neural networks. Psychol. Rev. 120, 230 (2013). [PubMed]
  • Langton C. Computation at the edge of chaos: phase transitions and emergent computation. Physica D 42, 12–37 (1990).
  • Kvale K. Segmentation and labelling of speech (Fakultet for informasjonsteknologi, matematikk og elektroteknikk, 1993).
  • Evans N. & Levinson S. The myth of language universals: Language diversity and its importance for cognitive science. Behav. Brain Sci. 32, 429–448 (2009). [PubMed]
  • Piantadosi S., Tily H. & Gibson E. Word lengths are optimized for efficient communication. Proc. Natl. Acad. Sci. USA 108, 3526–3529 (2011). [PubMed]
  • Clauset A., Shalizi C. & Newman M. Power-law distributions in empirical data. SIAM Rev. 51, 661–703 (2009).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group