Dependence of phonation threshold pressure on vocal tract acoustics and vocal fold tissue mechanics

Analytical and computer simulation studies have shown that the acoustic impedance of the vocal tract as well as the viscoelastic properties of vocal fold tissues are critical for determining the dynamics and the energy transfer mechanism of vocal fold oscillation. In the present study, a linear, small-amplitude oscillation theory was revised by taking into account the propagation of a mucosal wave and the inertive reactance (inertance) of the supraglottal vocal tract as the major energy transfer mechanisms for flow-induced self-oscillation of the vocal fold. Specifically, analytical results predicted that phonation threshold pressure (Pth) increases with the viscous shear properties of the vocal fold, but decreases with vocal tract inertance. This theory was empirically tested using a physical model of the larynx, where biological materials (fat, hyaluronic acid, and fibronectin) were implanted into the vocal fold cover to investigate the effect of vocal fold tissue viscoelasticity on Pth. A uniform-tube supraglottal vocal tract was also introduced to examine the effect of vocal tract inertance on Pth. Results showed that Pth decreased with the inertive impedance of the vocal tract and increased with the viscous shear modulus (G") or dynamic viscosity (eta') of the vocal fold cover, consistent with theoretical predictions. These findings supported the potential biomechanical benefits of hyaluronic acid as a surgical bioimplant for repairing voice disorders involving the superficial layer of the lamina propria, such as scarring, sulcus vocalis, atrophy, and Reinke's edema.     

Investigating acoustic correlates of human vocal fold vibratory phase asymmetry through modeling and laryngeal high-speed videoendoscopy

Vocal fold vibratory asymmetry is often associated with inefficient sound production through its impact on source spectral tilt. This association is investigated in both a computational voice production model and a group of 47 human subjects. The model provides indirect control over the degree of left-right phase asymmetry within a nonlinear source-filter framework, and high-speed videoendoscopy provides in vivo measures of vocal fold vibratory asymmetry. Source spectral tilt measures are estimated from the inverse-filtered spectrum of the simulated and recorded radiated acoustic pressure. As expected, model simulations indicate that increasing left-right phase asymmetry induces steeper spectral tilt. Subject data, however, reveal that none of the vibratory asymmetry measures correlates with spectral tilt measures. Probing further into physiological correlates of spectral tilt that might be affected by asymmetry, the glottal area waveform is parameterized to obtain measures of the open phase (open/plateau quotient) and closing phase (speed/closing quotient). Subjects' left-right phase asymmetry exhibits low, but statistically significant, correlations with speed quotient (r=0.45) and closing quotient (r=-0.39). Results call for future studies into the effect of asymmetric vocal fold vibration on glottal airflow and the associated impact on voice source spectral properties and vocal efficiency.     

Medial surface dynamics of an in vivo canine vocal fold during phonation

Quantitative measurement of the medial surface dynamics of the vocal folds is important for understanding how sound is generated within the larynx. Building upon previous excised hemilarynx studies, the present study extended the hemilarynx methodology to the in vivo canine larynx. Through use of an in vivo model, the medial surface dynamics of the vocal fold were examined as a function of active thyroarytenoid muscle contraction. Data were collected using high-speed digital imaging at a sampling frequency of 2000 Hz, and a spatial resolution of 1024 x 1024 pixels. Chest-like and fry-like vibrations were observed, but could not be distinguished based on the input stimulation current to the recurrent laryngeal nerve. The subglottal pressure did distinguish the registers, as did an estimate of the thyroarytenoid muscle activity. Upon quantification of the three-dimensional motion, the method of Empirical Eigenfunctions was used to extract the underlying modes of vibration, and to investigate mechanisms of sustained oscillation. Results were compared with previous findings from excised larynx experiments and theoretical models.     

Vocal fold and ventricular fold vibration in period-doubling phonation: physiological description and aerodynamic modeling

Occurrences of period-doubling are found in human phonation, in particular for pathological and some singing phonations such as Sardinian A Tenore Bassu vocal performance. The combined vibration of the vocal folds and the ventricular folds has been observed during the production of such low pitch bass-type sound. The present study aims to characterize the physiological correlates of this acoustical production and to provide a better understanding of the physical interaction between ventricular fold vibration and vocal fold self-sustained oscillation. The vibratory properties of the vocal folds and the ventricular folds during phonation produced by a professional singer are analyzed by means of acoustical and electroglottographic signals and by synchronized glottal images obtained by high-speed cinematography. The periodic variation in glottal cycle duration and the effect of ventricular fold closing on glottal closing time are demonstrated. Using the detected glottal and ventricular areas, the aerodynamic behavior of the laryngeal system is simulated using a simplified physical modeling previously validated in vitro using a larynx replica. An estimate of the ventricular aperture extracted from the in vivo data allows a theoretical prediction of the glottal aperture. The in vivo measurements of the glottal aperture are then compared to the simulated estimations.     

A model for vocal fold vibratory motion, contact area, and the electroglottogram

The electroglottogram (EGG) has been conjectured to be related to the area of contact between the vocal folds. This hypothesis has been substantiated only partially via direct and indirect observations. In this paper, a simple model of vocal fold vibratory motion is used to estimate the vocal fold contact area as a function of time. This model employs a limited number of vocal fold vibratory features extracted from ultra high-speed laryngeal films. These characteristics include the opening and closing vocal fold angles and the lag (phase difference) between the upper and lower vocal fold margins. The electroglottogram is simulated using the contact area, and the EGG waveforms are compared to measured EGGs for normal male voices producing both modal and pulse register tones. The model also predicts EGG waveforms for vocal fold vibration associated with a nodule or polyp.     

Spatio-temporal analysis of irregular vocal fold oscillations: biphonation due to desynchronization of spatial modes.

This report is on direct observation and modal analysis of irregular spatio-temporal vibration patterns of vocal fold pathologies in vivo. The observed oscillation patterns are described quantitatively with multiline kymograms, spectral analysis, and spatio-temporal plots. The complex spatio-temporal vibration patterns are decomposed by empirical orthogonal functions into independent vibratory modes. It is shown quantitatively that biphonation can be induced either by left-right asymmetry or by desynchronized anterior-posterior vibratory modes, and the term "AP (anterior-posterior) biphonation" is introduced. The presented phonation examples show that for normal phonation the first two modes sufficiently explain the glottal dynamics. The spatio-temporal oscillation pattern associated with biphonation due to left-right asymmetry can be explained by the first three modes. Higher-order modes are required to describe the pattern for biphonation induced by anterior-posterior vibrations. Spatial irregularity is quantified by an entropy measure, which is significantly higher for irregular phonation than for normal phonation. Two asymmetry measures are introduced: the left-right asymmetry and the anterior-posterior asymmetry, as the ratios of the fundamental frequencies of left and right vocal fold and of anterior-posterior modes, respectively. These quantities clearly differentiate between left-right biphonation and anterior-posterior biphonation. This paper proposes methods to analyze quantitatively irregular vocal fold contour patterns in vivo and complements previous findings of desynchronization of vibration modes in computer modes and in in vitro experiments.     

Sensitivity of a continuum vocal fold model to geometric parameters, constraints, and boundary conditions

The influence of key dimensional parameters, motion constraints, and boundary conditions on the modal properties of an idealized, continuum model of the vocal folds was investigated. The Ritz method and the finite element method were used for the analysis. The model's vibratory modes were determined to be most sensitive to changes in the anterior-posterior length of the vocal fold model, due to the influence of three-dimensional stress components acting in the transverse plane. Anterior/ posterior boundary conditions were found to have a significant influence on the vibratory response. Overestimation of modal frequencies resulted when vibration of the structure was restricted to the transverse plane. The overestimation of each modal frequency was proportional to the ratio of longitudinal to transverse Young's modulus, and was significant for ratio values less than 20.     

Effect of inferior surface angle on the self-oscillation of a computational vocal fold model

Geometry of the human vocal folds strongly influences their oscillatory motion. While the effect of intraglottal geometry on phonation has been widely investigated, the study of the geometry of the inferior surface of the vocal folds has been limited. In this study the way in which the inferior vocal fold surface angle affects vocal fold vibration was explored using a two-dimensional, self-oscillating finite element vocal fold model. The geometry was parameterized to create models with five different inferior surface angles. Four of the five models exhibited self-sustained oscillations. Comparisons of model motion showed increased vertical displacement and decreased glottal width amplitude with decreasing inferior surface angle. In addition, glottal width and air flow rate waveforms changed as the inferior surface angle was varied. Structural, rather than aerodynamic, effects are shown to be the cause of the changes in model response as the inferior surface angle was varied. Supporting data including glottal pressure distribution, average intraglottal pressure, energy transfer, and flow separation point locations are discussed, and suggestions for future research are given.     

Mammalian laryngseal air sacs add variability to the vocal tract impedance: physical and computational modeling

Cavities branching off the main vocal tract are ubiquitous in nonhumans. Mammalian air sacs exist in human relatives, including all four great apes, but only a substantially reduced version exists in humans. The present paper focuses on acoustical functions of the air sacs. The hypotheses are investigated on whether the air sacs affect amplitude of utterances and/or position of formants. A multilayer synthetic model of the vocal folds coupled with a vocal tract model was utilized. As an air sac model, four configurations were considered: open and closed uniform tube-like side branches, a rigid cavity, and an inflatable cavity. Results suggest that some air sac configurations can enhance the sound level. Furthermore, an air sac model introduces one or more additional resonance frequencies, shifting formants of the main vocal tract to some extent but not as strongly as previously suggested. In addition, dynamic range of vocalization can be extended by the air sacs. A new finding is also an increased variability of the vocal tract impedance, leading to strong nonlinear source-filter interaction effects. The experiments demonstrated that air-sac-like structures can destabilize the sound source. The results were validated by a transmission line computational model.     

Identification of geometric parameters influencing the flow-induced vibration of a two-layer self-oscillating computational vocal fold model

Simplified models have been used to simulate and study the flow-induced vibrations of the human vocal folds. While it is clear that the models' responses are sensitive to geometry, it is not clear how and to what extent specific geometric features influence model motion. In this study geometric features that played significant roles in governing the motion of a two-layer (body-cover), two-dimensional, finite element vocal fold model were identified. The model was defined using a flow solver based on the viscous, unsteady, Navier-Stokes equations and a solid solver that allowed for large strain and deformation. A screening-type design-of-experiments approach was used to identify the relative importance of 13 geometric parameters. Five output measures were analyzed to assess the magnitude of each geometric parameter's effect on the model's motion. The measures related to frequency, glottal width, flow rate, intraglottal angle, and intraglottal phase delay. The most significant geometric parameters were those associated with the cover--primarily the pre-phonatory intraglottal angle--as well as the body inferior angle. Some models exhibited evidence of improved model motion, including mucosal wave-like motion and alternating convergent-divergent glottal profiles, although further improvements are still needed to more closely mimic human vocal fold motion.     

Characteristics of vocal fold adduction related to voice onset

This study examined whether vocal fold kinematics prior to phonation differed between hard (glottal), normal, or breathy onsets in men and women. Glottal landmarks were identified and digitized from videotape recorded with a rigid laryngoscope during different voice onset types. Significant linear relationships (p 0.0055) were found among onset types on measures of (a) gesture duration when moving from 80% to 20% of maximum distance during adduction, (b) maximum velocity, (c) duration between the completion of adduction and phonation onset, and (d) ratios of maximum velocity to maximum distance between the vocal processes, an estimate of stiffness. The gesture duration was greatest for breathy onsets and least for hard onsets, while the maximum velocity, latency between adduction and phonation onset, and estimated stiffness were greatest for hard onsets and least for breathy onsets. The results suggest that one trajectory seems to be used with increases in gesture duration being accompanied by decreases in articulator stiffness when moving from hard to normal to breathy voice onset types.     

Infraglottic aspect of canine vocal fold vibration: Effect of increase of mean airflow rate and lengthening of vocal fold

The mucosal upheaval (MU), where the mucosal wave starts and propagates upward, appears only when the vocal fold vibrates. The location of the MU histologically and the effect of changes in mean air flow rate (MFR) and vocal fold length on occurrence of the MU were studied in twelve excised canine larynges. The lower surface of the vocal fold was marked to serve as a landmark for subsequent study. Cricothyroid approximation was performed to lengthen the vocal fold. After taking high-speed pictures or recording stroboscopic images from the tracheal side, a small cut wound was made at the mark. This wound served to compare the position of the MU with the histologically identified location of the mark. The larynx was then sectioned in the frontal plane. Before lengthening the vocal fold, the MU occurred on the area where the lamina propria became thinner and where the muscular layer neared the epithelial layer. After lengthening the vocal fold, the MU actually shifted medially compared with its original position. The subglottic area surrounded by the bilateral MUs became longer and thinner. Whether or not complete glottal closure during a vibratory cycle was achieved did not alter these findings. In contrast, with a fixed vocal fold length the MU appeared more laterally as MFR increased, but, based on the relation with the mark, its location on the vocal fold did not change from its original position before increase of MFR.     

Videokymography: High-speed line scanning of vocal fold vibration

A digital technique for high-speed visualization of vibration, called videokymography, was developed and applied to the vocal folds. The system uses a modified video camera able to work in two modes: high-speed (nearly 8,000 images/s) and standard (50 images/s in CCIR norm). In the high-speed mode, the camera selects one active horizontal line (transverse to the glottis) from the whole laryngeal image. The successive line images are presented in real time on a commercial TV monitor, filling each video frame from top to bottom. The system makes it possible to observe left-right asymmetries, open quotient, propagation of mucosal waves, movement of the upper and, in the closing phase, the lower margins of the vocal folds, etc. The technique is suitable for further processing and quantification of recorded vibration.     

Autologous fat implantation for vocal fold scar: A preliminary report

New insights into the anatomy and physiology of phonation, along with technological advances in voice assessment and quantification, have led to dramatic improvements in medical voice care. Techniques to prevent vocal fold scar have been among the most important, especially scarring and hoarseness associated with voice surgery. Nevertheless, dysphonia due to vocal fold scar is still encountered all too frequently. Although it is not generally possible to restore such injured voices to normal, patients with scar-induced dysphonia can usually be helped. Voice improvement is optimized through a team approach. Treatment may include sophisticated voice therapy and vocal fold surgery. Although experience with collagen injection has been encouraging in selected cases (particularly in those involving limited areas of vocal fold scar), there is no consistently successful surgical technique. Attempts to treat massive vocal fold scar, such as may be seen following vocal fold stripping, have been particularly unsuccessful. This paper reports preliminary experience with the implantation of autologous fat into the vibratory margin of the vocal fold of patients with severe, extensive scarring. Using this technique, it appears possible to recreate a mucosal wave and improve voice quality. Additional research is needed.     

Electroglottographic wavegrams: a technique for visualizing vocal fold dynamics noninvasively

A method for analyzing and displaying electroglottographic (EGG) signals (and their first derivative, DEGG) is introduced: the electroglottographic wavegram ("wavegram" hereafter). To construct a wavegram, the time-varying fundamental frequency is measured and consecutive individual glottal cycles are identified. Each cycle is locally normalized in duration and amplitude, the signal values are encoded by color intensity and the cycles are concatenated to display the entire voice sample in a single image, similar as in sound spectrography. The wavegram provides an intuitive means for quickly assessing vocal fold contact phenomena and their variation over time. Variations in vocal fold contact appear here as a sequence of events rather than single phenomena, taking place over a certain period of time, and changing with pitch, loudness and register. Multiple DEGG peaks are revealed in wavegrams to behave systematically, indicating subtle changes of vocal fold oscillatory regime. As such, EGG wavegrams promise to reveal more information on vocal fold contacting and de-contacting events than previous methods.     

Cyclicity of laryngeal cavity resonance due to vocal fold vibration

Acoustic effects of the time-varying glottal area due to vocal fold vibration on the laryngeal cavity resonance were investigated based on vocal tract area functions and acoustic analysis. The laryngeal cavity consists of the vestibular and ventricular parts of the larynx, and gives rise to a regional acoustic resonance within the vocal tract, with this resonance imparting an extra formant to the vocal tract resonance pattern. Vocal tract transfer functions of the five Japanese vowels uttered by three male subjects were calculated under open- and closed-glottis conditions. The results revealed that the resonance appears at the frequency region from 3.0 to 3.7 kHz when the glottis is closed and disappears when it is open. Real spectra estimated from open- and closed-glottis periods of vowel sounds also showed the on-off pattern of the resonance within a pitch period. Furthermore, a time-domain acoustic analysis of vowels indicated that the resonance component could be observed as a pitch-synchronized rise-and-fall pattern of the bandpass amplitude. The cyclic nature of the resonance can be explained as the laryngeal cavity acting as a closed tube that generates the resonance during a closed-glottis period, but damps the resonance off during an open-glottis period.     

Nonlinear behavior of vocal fold vibration: The role of coupling between the vocal folds

Coupling between the vocal folds is one of the nonlinear mechanisms allowing regulation and synchronization of mucosal vibration. The purpose of this study was to establish that modulations such as diplophonia and abnormalities observed in vocal signals that may be observed in some cases of laryngeal pathology can be considered as nonlinear behavior due to the persistence of some physical interaction (coupling). An experimental model using excised porcine larynx was designed to create tension asymmetry between the vocal folds and to obtain vocal signals with modulations. Signals were analyzed by spectral analysis and the phase portrait method. Results were compared with computer-generated synthetic signals corresponding to nonlinear combinations of sinusoid signals. Under these conditions, evidence of nonlinear behavior was detected in 85% of experimental signals. These findings were interpreted as a demonstration of vocal fold interaction. Based on these findings, the authors conclude that (1) coupling must be taken into account in physical models of laryngeal physiology, and that (2) methods of nonlinear dynamics may be used for objective voice analysis.