Effects of cricothyroid and thyroarytenoid interaction on voice control: Muscle activity, vocal fold biomechanics, flow, and acoustics

Sensitivity of Phonation Onset Pressure to Vocal Fold Stiffness Distribution

Phonation onset is characterized by the unstable growth of vocal fold (VF) vibrations that ultimately results in self-sustained oscillation and the production of modal voice. Motivated by histological studies, much research has focused on the role of the layered structure of the vocal folds in influencing phonation onset, wherein the outer "cover" layer is relatively soft and the inner "body" layer is relatively stiff. Recent research, however, suggests that the body-cover (BC) structure over-simplifies actual stiffness distributions by neglecting important spatial variations, such as inferior-superior (IS) and anterior-posterior gradients and smooth transitions in stiffness from one histological layer to another. Herein, we explore sensitivity of phonation onset to stiffness gradients and smoothness. By assuming no a priori stiffness distribution and considering a second-order Taylor series sensitivity analysis of phonation onset pressure with respect to stiffness, we find two general smooth stiffness distributions most strongly influence onset pressure: a smooth stiffness containing aspects of BC differences and IS gradients in the cover, which plays a role in minimizing onset pressure, and uniform increases in stiffness, which raise onset pressure and frequency. While the smooth stiffness change contains aspects qualitatively similar to layered BC distributions used in computational studies, smooth transitions in stiffness result in higher sensitivity of onset pressure than discrete layering. These two general stiffness distributions also provide a simple, low-dimensional, interpretation of how complex variations in VF stiffness affect onset pressure, enabling refined exploration of the effects of stiffness distributions on phonation onset.

The effect of swelling on vocal fold kinematics and dynamics

Swelling in the vocal folds is caused by the local accumulation of fluid, and has been implicated as a phase in the development of phonotraumatic vocal hyperfunction and related structural pathologies, such as vocal fold nodules. It has been posited that small degrees of swelling may be protective, but large amounts may lead to a vicious cycle wherein the engorged folds lead to conditions that promote further swelling, leading to pathologies. As a first effort to explore the mechanics of vocal fold swelling and its potential role in the etiology of voice disorders, this study employs a finite-element model with swelling confined to the superficial lamina propria, which changes the volume, mass, and stiffness of the cover layer. The impacts of swelling on a number of vocal fold kinematic and damage measures, including von Mises stress, internal viscous dissipation, and collision pressure, are presented. Swelling has small but consistent effects on voice outputs, including a reduction in fundamental frequency with increasing swelling (10 Hz at 30 % swelling). Average von Mises stress decreases slightly for small degrees of swelling but increases at large magnitudes, consistent with expectations for a vicious cycle. Both viscous dissipation and collision pressure consistently increase with the magnitude of swelling. This first effort at modeling the impact of swelling on vocal fold kinematics, kinetics, and damage measures highlights the complexity with which phonotrauma can influence performance metrics. Further identification and exploration of salient candidate measures of damage and refined studies coupling swelling with local phonotrauma are expected to shed further light on the etiological pathways of phonotraumatic vocal hyperfunction.

Modeling the influence of the extrinsic musculature on phonation

Neck muscles play important roles in various physiological tasks, including swallowing, head stabilization, and phonation. The mechanisms by which neck muscles influence phonation are not well understood, with conflicting reports on the change in fundamental frequency for ostensibly the same neck muscle activation scenarios. In this work, we introduce a reduced-order muscle-controlled vocal fold model, comprising both intrinsic muscle control and extrinsic muscle effects. The model predicts that when the neck muscles pull the thyroid cartilage in the superior-anterior direction (with a sufficiently large anterior component), inferior direction, or inferior-anterior direction, tension in the vocal folds increases, leading to fundamental frequency rise during sustained phonation. On the other hand, pulling in the superior direction, superior-posterior direction, or inferior-posterior direction (with a sufficiently large posterior component) tends to decrease vocal fold tension and phonation fundamental frequency. Varying the pulling force location alters the posture and phonation biomechanics, depending on the force direction. These findings suggest potential roles of particular neck muscles in modulating phonation fundamental frequency, with implications for vocal hyperfunction.

Predicting 3D soft tissue dynamics from 2D imaging using physics informed neural networks

Tissue dynamics play critical roles in many physiological functions and provide important metrics for clinical diagnosis. Capturing real-time high-resolution 3D images of tissue dynamics, however, remains a challenge. This study presents a hybrid physics-informed neural network algorithm that infers 3D flow-induced tissue dynamics and other physical quantities from sparse 2D images. The algorithm combines a recurrent neural network model of soft tissue with a differentiable fluid solver, leveraging prior knowledge in solid mechanics to project the governing equation on a discrete eigen space. The algorithm uses a Long-short-term memory-based recurrent encoder-decoder connected with a fully connected neural network to capture the temporal dependence of flow-structure-interaction. The effectiveness and merit of the proposed algorithm is demonstrated on synthetic data from a canine vocal fold model and experimental data from excised pigeon syringes. The results showed that the algorithm accurately reconstructs 3D vocal dynamics, aerodynamics, and acoustics from sparse 2D vibration profiles.

Examining the influence of epithelium layer modeling approaches on vocal fold kinematics and kinetics

Grouping the thin epithelium and thicker superficial lamina propria layers into a single cover layer has been widely adopted in finite element vocal fold models. Recent silicone vocal fold studies have suggested, however, that inclusion of a distinct epithelial layer leads to more physiologically representative motion. This study systematically explores the ramifications of incorporating an epithelial layer into a cover grouping for finite element vocal fold modeling. A membrane model for the epithelium is introduced to facilitate parametric investigation by reducing the mesh density requirement of the epithelium into a single infinitesimally thin layer. Excluding the epithelium entirely leads to increased energy in higher order modes and larger inferior-superior excursion of the folds. Integrating the epithelium into a cover layer with volume-weighted average stiffness results in similar kinematics to that of a model treating the epithelium as a distinct layer. However, the internal stress/strain and contact pressure during collision are higher, and viscous dissipation is lower, when the epithelium is integrated into the cover. Thus, careful treatment of the epithelium is recommended for finite element studies, particularly when employing the models for estimating measures dependent upon internal stress/strain and/or collision pressure, such as vocal dose.

Exploring the mechanics of fundamental frequency variation during phonation onset

Fundamental frequency patterns during phonation onset have received renewed interest due to their promising application in objective classification of normal and pathological voices. However, the associated underlying mechanisms producing the wide array of patterns observed in different phonetic contexts are not yet fully understood. Herein, we employ theoretical and numerical analyses in an effort to elucidate the potential mechanisms driving opposing frequency patterns for initial/isolated vowels versus vowels preceded by voiceless consonants. Utilizing deterministic lumped-mass oscillator models of the vocal folds, we systematically explore the roles of collision and muscle activation in the dynamics of phonation onset. We find that an increasing trend in fundamental frequency, as observed for initial/isolated vowels, arises naturally through a progressive increase in system stiffness as collision intensifies as onset progresses, without the need for time-varying vocal fold tension or changes in aerodynamic loading. In contrast, reduction in cricothyroid muscle activation during onset is required to generate the decrease in fundamental frequency observed for vowels preceded by voiceless consonants. For such phonetic contexts, our analysis shows that the magnitude of reduction in the cricothyroid muscle activation and the activation level of the thyroarytenoid muscle are potential factors underlying observed differences in (relative) fundamental frequency between speakers with healthy and hyperfunctional voices. This work highlights the roles of sometimes competing laryngeal factors in producing the complex array of observed fundamental frequency patterns during phonation onset.

Control of Pre-phonatory Glottal Shape by Intrinsic Laryngeal Muscles

OBJECTIVES

Surgical manipulations to treat glottic insufficiency aim to restore the physiologic pre-phonatory glottal shape. However, the physiologic pre-phonatory glottal shape as a function of interactions between all intrinsic laryngeal muscles (ILMs) has not been described. Vocal fold posture and medial surface shape were investigated across concurrent activation and interactions of thyroarytenoid (TA), cricothyroid (CT), and lateral cricoarytenoid/interarytenoid (LCA/IA) muscles.

STUDY DESIGN

In vivo canine hemilarynx model.

METHODS

The ILMs were stimulated across combinations of four graded levels each from low-to-high activation. A total of 64 distinct medial surface postures (4 TA × 4 CT × 4 LCA/IA levels) were captured using high-speed video. Using a custom 3D interpolation algorithm, the medial surface shape was reconstructed.

RESULTS

Combined activation of ILMs yielded a range of unique pre-phonatory postures. Both LCA/IA and TA activation adducted the vocal fold but with greater contribution from TA. The transition from a convergent to a rectangular glottal shape was primarily mediated by TA muscle activation but LCA/IA and TA together resulted in a smooth rectangular glottis compared to TA alone, which caused rectangular glottis with inferomedial bulging. CT activation resulted in a lengthened but slightly abducted glottis.

CONCLUSIONS

TA was primarily responsible for the rectangular shape of the adducted glottis with synergistic contribution from the LCA/IA. CT contributed minimally to vocal fold medial shape but elongated the glottis. These findings further refine laryngeal posture goals in surgical correction of glottic insufficiency.

LEVEL OF EVIDENCE

N/A, Basic science Laryngoscope, 2022.

LaDIVA: A neurocomputational model providing laryngeal motor control for speech acquisition and production

Many voice disorders are the result of intricate neural and/or biomechanical impairments that are poorly understood. The limited knowledge of their etiological and pathophysiological mechanisms hampers effective clinical management. Behavioral studies have been used concurrently with computational models to better understand typical and pathological laryngeal motor control. Thus far, however, a unified computational framework that quantitatively integrates physiologically relevant models of phonation with the neural control of speech has not been developed. Here, we introduce LaDIVA, a novel neurocomputational model with physiologically based laryngeal motor control. We combined the DIVA model (an established neural network model of speech motor control) with the extended body-cover model (a physics-based vocal fold model). The resulting integrated model, LaDIVA, was validated by comparing its model simulations with behavioral responses to perturbations of auditory vocal fundamental frequency (fo) feedback in adults with typical speech. LaDIVA demonstrated capability to simulate different modes of laryngeal motor control, ranging from short-term (i.e., reflexive) and long-term (i.e., adaptive) auditory feedback paradigms, to generating prosodic contours in speech. Simulations showed that LaDIVA’s laryngeal motor control displays properties of motor equivalence, i.e., LaDIVA could robustly generate compensatory responses to reflexive vocal fo perturbations with varying initial laryngeal muscle activation levels leading to the same output. The model can also generate prosodic contours for studying laryngeal motor control in running speech. LaDIVA can expand the understanding of the physiology of human phonation to enable, for the first time, the investigation of causal effects of neural motor control in the fine structure of the vocal signal.

Triangular body-cover model of the vocal folds with coordinated activation of the five intrinsic laryngeal muscles

Poor laryngeal muscle coordination that results in abnormal glottal posturing is believed to be a primary etiologic factor in common voice disorders such as non-phonotraumatic vocal hyperfunction. Abnormal activity of antagonistic laryngeal muscles is hypothesized to play a key role in the alteration of normal vocal fold biomechanics that results in the dysphonia associated with such disorders. Current low-order models of the vocal folds are unsatisfactory to test this hypothesis since they do not capture the co-contraction of antagonist laryngeal muscle pairs. To address this limitation, a self-sustained triangular body-cover model with full intrinsic muscle control is introduced. The proposed scheme shows good agreement with prior studies using finite element models, excised larynges, and clinical studies in sustained and time-varying vocal gestures. Simulations of vocal fold posturing obtained with distinct antagonistic muscle activation yield clear differences in kinematic, aerodynamic, and acoustic measures. The proposed tool is deemed sufficiently accurate and flexible for future comprehensive investigations of non-phonotraumatic vocal hyperfunction and other laryngeal motor control disorders.

A computational framework for patient-specific surgical planning of type 1 thyroplasty

A computational framework is proposed for virtual optimization of implant configurations of type 1 thyroplasty based on patient-specific laryngeal structures reconstructed from MRI images. Through integration of a muscle mechanics-based laryngeal posturing model, a flow-structure-acoustics interaction voice production model, a real-coded genetic algorithm, and virtual implant insertion, the framework acquires the implant configuration that achieves the optimal acoustic objectives. The framework is showcased by successfully optimizing an implant that restores acoustic features of a diseased voice resulted from unilateral vocal fold paralysis (UVFP) in producing a sustained vowel utterance. The sound intensity is improved from 62 dB (UVFP) to 81 dB (post-correction).