On the relation between subglottal pressure and fundamental frequency in phonation

The change in fundamental frequency with subglottal pressure in phonation is quantified on the basis of the ratio between vibrational amplitude and vocal fold length. This ratio is typically very small in stringed instruments, but becomes quite appreciable in vocal fold vibration. Tension in vocal fold tissues is, therefore, not constant over the vibratory cycle, and a dynamic tension gives rise to amplitude-frequency dependence. It is shown that the typical 2-6 Hz/cm H2O rise in fundamental frequency with subglottal pressure observed in human and canine larynges is a direct and predictable consequence of this amplitude-frequency dependence. Results are presently limited to phonation in the chest register.     

Mechanisms of irregular vibration in a physical model of the vocal folds

Previous investigations have shown that one mechanism of irregular vocal fold vibration may be a desynchronization of two or more vibratory modes of the vocal fold tissues. In the current investigation, mechanisms of irregular vibration were further examined using a self-oscillating, physical model of vocal fold vibration, a hemi-model methodology, and high-speed, stereoscopic, digital imaging. Using the method of empirical eigen-functions, a spatiotemporal analysis revealed mechanisms of irregular vibration in subharmonic phonation and biphonation, which were not disclosed in a standard acoustic spectrum.     

Deviant Vocal Fold Vibration as Observed During Videokymography: The Effect on Voice Quality

Videokymographic images of deviant or irregular vocal fold vibration, including diplophonia, the transition from falsetto to modal voice, irregular vibration onset and offset, and phonation following partial laryngectomy were compared with the synchronously recorded acoustic speech signals. A clear relation was shown between videokymographic image sequences and acoustic speech signals, and the effect of irregular or incomplete vocal fold vibration patterns was recognized in the amount of perceived breathiness and roughness and by the harmonics-to-noise ratio in the speech signal. Mechanisms causing roughness are the presence of mucus, phase differences between the left and right vocal fold, and short-term frequency and amplitude modulation. It can be concluded that the use of simultaneously recorded videokymographic image sequences and speech signals contributes to the understanding of the effect of irregular vocal fold vibration on voice quality.     

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.     

Model-based classification of nonstationary vocal fold vibrations

Classification of vocal fold vibrations is an essential task of the objective assessment of voice disorders. For historical reasons, the conventional clinical examination of vocal fold vibrations is done during stationary, sustained phonation. However, the conclusions drawn from a stationary phonation are restricted to the observed steady-state vocal fold vibrations and cannot be generalized to voice mechanisms during running speech. This study addresses the approach of classifying real-time recordings of vocal fold oscillations during a nonstationary phonation paradigm in the form of a pitch raise. The classification is based on asymmetry measures derived from a time-dependent biomechanical two-mass model of the vocal folds which is adapted to observed vocal fold motion curves with an optimization procedure. After verification of the algorithm performance the method was applied to clinical problems. Recordings of ten subjects with normal voice and ten dysphonic subjects have been evaluated during stationary as well as nonstationary phonation. In the case of nonstationary phonation the model-based classification into "normal" and "dysphonic" succeeds in all cases, while it fails in the case of sustained phonation. The nonstationary vocal fold vibrations contain additional information about vocal fold irregularities, which are needed for an objective interpretation and classification of voice disorders.     

Simulation of vocal fold impact pressures with a self-oscillating finite-element model

Vocal fold impact pressures were studied using a self-oscillating finite-element model capable of simulating vocal fold vibration and airflow. The calculated airflow pressure is applied on the vocal fold as the driving force. The airflow region is then adjusted according to the calculated vocal fold displacement. The interaction between airflow and the vocal folds produces a self-oscillating solution. Lung pressures between 0.2 and 2.5 kPa were used to drive this self-oscillating model. The spatial distribution of the impact pressure was studied. Studies revealed that the tissue collision during phonation produces a very large impact pressure which correlates with the lung pressure and glottal width. Larger lung pressure and a narrower glottal width increase the impact pressure. The impact pressure was found to be roughly the square root of lung pressure. In the inferior-superior direction, the maximum impact pressure is related to the narrowest glottis. In the anterior-posteriorfirection, the greatest impact pressure appears at the midpoint of the vocal fold. The match between our numerical simulations and clinical observations suggests that this self-oscillating finite-element model might be valuable for predicting mechanical trauma of the vocal folds.     

Strobophotoglottographic transillumination as a method for the analysis of vocal fold vibration patterns

The purpose of this exploratory study was to determine if laryngeal transillumination in combination with stroboscopy (strobophotoglottography; SPGG) is useful for (1) the visualization of vocal fold vibration (VFV) opening patterns, (2) the localization of initial vocal fold opening in horizontal glottal thirds (anterior, midmembranous, and posterior), (3) determination of the temporal correspondence of the so-called electroglottography (EGG)-knee and initial vocal fold separation, and, finally, (4) automatized quantitative measurements of glottal area function within endoscopic images. With stroboscopic transillumination, initial inferior vocal fold separation was detectable during the "closed" phase, where the vocal folds were still closed in the upper portion and therefore initial inferior vocal fold separation could not be visualized with usual laryngoscopy techniques. In the horizontal plane within similar fundamental frequencies in modal voice registers in two male subjects, localization of initial glottal opening depended on the voice types used (soft, normal, or pressed phonation). We found zipperlike posterior-to-anterior openings, initial midmembranous openings, initial anterior openings, as well as simultaneous initial opening of all three portions in the two healthy male adults examined. This technique proved to add temporal and spatial information to vocal fold opening patterns and extends our examination techniques to the very beginning of vocal fold opening at the inferior portion. Simultaneous electroglottogram tracking and comparison with bidirectionally illuminated stroboscopic images revealed a time-locked correspondence of the EGG-knee with the aforementioned initial inferior vocal fold separation. Bidirectional illumination combined with digital color extraction techniques allowed for image separation of subglottally and supraglottally illuminated structures. This facilitated vocal fold contour detection and automatized image processing, for example, for determination of glottal area function, and is considered to be a further step to objective automatized quantitative measurements within endoscopic images.     

Three-dimensional biomechanical properties of human vocal folds: parameter optimization of a numerical model to match in vitro dynamics

The human voice signal originates from the vibrations of the two vocal folds within the larynx. The interactions of several intrinsic laryngeal muscles adduct and shape the vocal folds to facilitate vibration in response to airflow. Three-dimensional vocal fold dynamics are extracted from in vitro hemilarynx experiments and fitted by a numerical three-dimensional-multi-mass-model (3DM) using an optimization procedure. In this work, the 3DM dynamics are optimized over 24 experimental data sets to estimate biomechanical vocal fold properties during phonation. Accuracy of the optimization is verified by low normalized error (0.13 ± 0.02), high correlation (83% ± 2%), and reproducible subglottal pressure values. The optimized, 3DM parameters yielded biomechanical variations in tissue properties along the vocal fold surface, including variations in both the local mass and stiffness of vocal folds. That is, both mass and stiffness increased along the superior-to-inferior direction. These variations were statistically analyzed under different experimental conditions (e.g., an increase in tension as a function of vocal fold elongation and an increase in stiffness and a decrease in mass as a function of glottal airflow). The study showed that physiologically relevant vocal fold tissue properties, which cannot be directly measured during in vivo human phonation, can be captured using this 3D-modeling technique.     

Frequency and amplitude perturbation analysis of electroglottograph during sustained phonation

Electroglottography (EGG) was used to monitor vocal fold vibration patterns in normal subjects and patients with various laryngeal disorders. In order to evaluate the regularity of vocal fold vibration, frequency and amplitude perturbation of EGG waves during sustained phonation were measured with a laboratory computer. The data were compared to the degree of hoarseness evaluated by auditory perception and by sound spectrographic analysis. Frequency and amplitude perturbation measures showed some overlap between normal and pathological groups. However, there was a close relation between perturbation analysis of EGG waves and degree of hoarseness (Spearman's rank correlation coefficient rs = 0.73, p less than 0.0005). Amplitude perturbation was found to be a more sensitive measure of the irregularity of vocal fold vibration than frequency perturbation.     

Pressure-flow relationships during phonation in the canine larynx

Measurements of air pressure and flow were made using an in vivo canine model of the larynx. Subglottic pressures at varying flow rates were taken during phonation induced by laryngeal nerve stimulation. Results showed that during constant vocal fold stiffness, subglottic pressure rose slightly with increased air flow. The larynx in the in vivo canine model exhibited a flow-dependent decrease in laryngeal airway resistance. Increasing flow rate was associated with an increase in frequency of phonation and open quotient, as measured glottographically. Results from this experiment were compared with a theoretical two-mass model of the larynx and other theoretical models of phonation. The influence of aerodynamic forces on glottal vibration is explained by increased lateral excursion of the vocal folds during the open interval and shortening of the closed interval during the glottal cycle.     

Phonation Threshold Power in Ex Vivo Laryngeal Models

This study hypothesized that phonation threshold power is measureable and sensitive to changes in the biomechanical properties of the vocal folds. Phonation threshold power was measured in three sample populations of 10 excised canine larynges treated with variable posterior glottal gap, variable bilateral vocal fold elongation, and variable vocal fold lesioning. Posterior glottal gap varied from 0 to 4 mm in 0.5 mm intervals. Bilateral vocal fold elongation varied from 0% to 20% in 5% intervals. Vocal fold lesion treatments included unilateral and bilateral vocal fold lesion groups. Each treatment was investigated independently in a sample population of 10 excised canine larynges. Linear regression analysis indicated that phonation threshold power was sensitive to posterior glottal gap (R² = 0.298, P < 0.001) and weakly to vocal fold elongation (R² = 0.052, P = 0.003). A one-way repeated measures analysis of variance indicated that phonation threshold power was sensitive to the presence of lesions (P < 0.001). Theoretical and experimental evidence presented here suggests that phonation threshold power could be used as a broad screening parameter sensitive to certain changes in the biomechanical properties of the larynx. It has not yet been measured in humans, but because it has the potential to represent the airflow-tissue energy transfer more completely than the phonation threshold pressure or flow alone, it may be a more useful parameter than these and could be used to indicate that laryngeal health is likely abnormal.     

Visualization and Quantification of the Medial Surface Dynamics of an Excised Human Vocal Fold During Phonation

SUMMARY: The purpose of this investigation was to investigate physical mechanisms of vocal fold vibration during normal phonation through quantification of the medial surface dynamics of the fold. An excised hemilarynx setup was used. The dynamics of 30 microsutures mounted on the medial surface of a human vocal fold were analyzed across 18 phonatory conditions. The vibrations were recorded with a digital high-speed camera at a frequency of 4,000 Hz. The positions of the sutures were extracted and converted to three-dimensional coordinates using a linear approximation technique. The data were reduced to principal eigenfuctions, which captured over 90% of the variance of the data, and suggested mechanisms of sustained vocal fold oscillation. The vibrations were imaged as the following phonatory conditions were manipulated: glottal airflow, an adductory force applied to the muscular process, and an elongation force applied to the thyroid cartilage. Over the range of variables studied, only the variation in glottal airflow yielded significant changes in subglottal pressure and fundamental frequency. All recordings showed high correlation for the distribution of the dynamics across the medial surface of the vocal fold. The distribution of the different displacement directions and velocities showed the highest variations around the superior region of the medial surface. Although the computed vibration patterns of the two largest empirical eigenfunctions were consistent with previous experimental observations, the relative prominence of the two eigenfunctions changed as a function of glottal airflow, impacting theories of vocal efficiency and vocal economy.     

Physical mechanisms of phonation onset: a linear stability analysis of an aeroelastic continuum model of phonation

In an investigation of phonation onset, a linear stability analysis was performed on a two-dimensional, aeroelastic, continuum model of phonation. The model consisted of a vocal fold-shaped constriction situated in a rigid pipe coupled to a potential flow which separated at the superior edge of the vocal fold. The vocal fold constriction was modeled as a plane-strain linear elastic layer. The dominant eigenvalues and eigenmodes of the fluid-structure-interaction system were investigated as a function of glottal airflow. To investigate specific aerodynamic mechanisms of phonation onset, individual components of the glottal airflow (e.g., flow-induced stiffness, inertia, and damping) were systematically added to the driving force. The investigations suggested that flow-induced stiffness was the primary mechanism of phonation onset, involving the synchronization of two structural eigenmodes. Only under conditions of negligible structural damping and a restricted set of vocal fold geometries did flow-induced damping become the primary mechanism of phonation onset. However, for moderate to high structural damping and a more generalized set of vocal fold geometries, flow-induced stiffness remained the primary mechanism of phonation onset.     

Anterior-posterior biphonation in a finite element model of vocal fold vibration

In this paper, a finite-element model is used to simulate anterior-posterior biphonation [Neubauer et al., J. Acoust. Soc. Am. 110(6), 3179-3192 (2001)]. The anterior-posterior stiffness asymmetric factor and the anterior-posterior shape asymmetric factor describe the asymmetry properties of vocal folds. Spatiotemporal plot, spectral analysis, anterior-posterior fundamental frequency ratio, cross covariation function, and correlation length quantitatively estimate the spatial asymmetry of vocal fold oscillations. Calculation results show that the anterior-posterior stiffness asymmetry decreases the spatial coherence of vocal fold vibration. When the stiffness asymmetry reaches a certain level, the drop in spatial coherence desynchronizes the vibration modes. The anterior and posterior sides of the vocal fold oscillate with two independent fundamental frequencies (f(a) and f(p)). The complex spectral characteristics of vocal fold vibration under biphonation conditions can be explained by the linear combination of f(a) and f(p). Empirical orthogonal eigenfunctions prove the existence of higher-order anterior-posterior modes when anterior-posterior biphonation occurs. Then, it is found that the anterior-posterior shape asymmetry also decreases the spatial coherence of vocal fold vibration, and shape asymmetry is a possible reason for anterior-posterior biphonation.     

On the difference between negative damping and eigenmode synchronization as two phonation onset mechanisms

Negative damping and eigenmode synchronization as two different mechanisms of phonation onset are distinguished. Although both mechanisms lead to a favorable phase relationship between the flow pressure and the vocal fold motion as required for a net energy transfer into the vocal folds, the underlying mechanisms for this favorable phase relationship are different. The negative damping mechanism relies on glottal aerodynamics or acoustics to establish before onset and maintain after onset the favorable phase relationship, and therefore has minimum requirements on vocal fold geometry and biomechanics. A single degree-of-freedom vocal fold model is all that is needed for self-oscillation in the presence of a negative damping mechanism. In contrast, the mechanism of eigenmode synchronization critically depends on the geometrical and biomechanical properties of the vocal folds (at least 2-degrees-of-freedom are required), and has little requirement on the glottal aerodynamics other than flow separation. The favorable phase relation is established once synchronization occurs, regardless of the phase relationship imposed by glottal aerodynamics before onset. Unlike that of the negative damping mechanism, initiation of eigenmode synchronization requires neither a velocity-dependent flow pressure nor an alternating convergent-divergent glottis. The clinical implications of the distinctions between these two mechanisms are discussed.     

Vocal Warm-up Increases Phonation Threshold Pressure in Soprano Singers at High Pitch

Vocal warm-up is thought to optimize singing performance. We compared effects of short-term, submaximal, vocal warm-up exercise with those of vocal rest on the soprano voice (n = 10, ages 19-21 years). Dependent variables were the minimum subglottic air pressure required for vocal fold oscillation to occur (phonation threshold pressure, Pth), and the maximum and minimum phonation fundamental frequency. Warm-up increased Pth for high pitch phonation (p = 0.033), but not for comfortable (p = 0.297) or low (p = 0.087) pitch phonation. No significant difference in the maximum phonation frequency (p = 0.193) or minimum frequency (p = 0.222) was observed. An elevated Pth at controlled high pitch, but an unchanging maximum and minimum frequency production suggests that short-term vocal exercise may increase the viscosity of the vocal fold and thus serve to stabilize the high voice.     

Phonation threshold pressure and onset frequency in a two-layer physical model of the vocal folds

The influence of vocal fold geometry and stiffness on phonation onset was experimentally investigated using a body-cover physical model of the vocal folds. Results showed that a lower phonation threshold pressure and phonation onset frequency can be achieved by reducing body-layer or cover-layer stiffness, reducing medial surface thickness, or increasing cover-layer depth. Increasing body-layer stiffness also restricted vocal fold motion to the cover layer and reduced prephonatory glottal opening. Excitation of anterior-posterior modes was also observed, particularly for large values of the body-cover stiffness ratio. The results of this study were also discussed in relation to previous theoretical and experimental studies.     

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.     

Analysis of Voice and Quantitative Measurement of Glottal Gap After Thyroplasty Type I in the Treatment of Unilateral Vocal Paralysis

Thyroplasty type I is one of several surgical treatments in which improving the voice of unilateral vocal fold paralysis is the ultimate objective. The goal of the surgery is the medialization of the paralyzed vocal fold. The purpose of this study is to evaluate the effectiveness of thyroplasty type I through acoustical analysis, aerodynamic measures, and quantitative videostroboscopic measurements. We report on 20 patients with unilateral vocal cord paralysis who underwent thyroplasty type I. We performed preoperative and postoperative video image analysis (normalized glottal gap area) and computer-assisted voice analysis (fundamental frequency, jitter, shimmer, noise-to-harmonic ratio, mean phonation time, mean flow rate, mean subglottic pressure) in all patients. The glottal gap was significantly reduced after thyroplasty type I. Postoperative voice quality was characterized by an improved pitch and amplitude pertubation (jitter and shimmer), phonation time (mean phonation time), and subglottic pressure (mean subglottic pressure). Thyroplasty type I is an effective method for regaining glottal closure and vocal function.