Cluster analysis and classification of heart sounds

Acoustic heart signals, generated by the mechanical processes of the cardiac cycle, carry significant information about the underlying functioning of the cardiovascular system. We describe a computational analysis framework for identifying distinct morphologies of heart sounds and classifying them into physiological states. The analysis framework is based on hierarchical clustering, compact data representation in the feature space of cluster distances and a classification algorithm. We applied the proposed framework on two heart sound datasets, acquired during controlled alternations of the physiological conditions, and analyzed the morphological changes induced to the first heart sound (S1), and the ability to predict physiological variables from the morphology of S1. On the first dataset of 12 subjects, acquired while modulating the respiratory pressure, the algorithm achieved an average accuracy of 82 ± 7% in classifying the level of breathing resistance, and was able to estimate the instantaneous breathing pressure with an average error of 19 ± 6%. A strong correlation of 0.92 was obtained between the estimated and the actual breathing efforts. On the second dataset of 11 subjects, acquired during pharmacological stress tests, the average accuracy in classifying the stress stage was 86 ± 7%. The effects of the chosen raw signal representation, distance metrics and classification algorithm on the performance were studied on both real and simulated data. The results suggest that quantitative heart sound analysis may provide a new non-invasive technique for continuous cardiac monitoring and improved detection of mechanical dysfunctions caused by cardiovascular and cardiopulmonary diseases.     

Alternative methods of normalising EMG during cycling

We evaluated possible methods of normalisation for EMG measured during cycling. The MVC method, Sprint method and 70% Peak Power Output Method were investigated and their repeatability, reliability and sensitivity to change in workload were compared.Thirteen cyclists performed the same experimental protocol on three separate occasions. Each day, subjects firstly performed MVCs, followed by a 10 s maximal sprint on a cycle ergometer. Subjects then performed a Peak Power Output (PPO) test until exhaustion. After which they cycled at 70% of PPO for 5 min at 90 rpm. Results indicated that normalising EMG data to 70% PPO is more repeatable, the intra-class correlation (ICC) of 70% PPO (0.87) was significantly higher than for MVC (0.66) (p = 0.03) and 10 s sprint (0.65) (p = 0.04). The 70% PPO method also demonstrated the least intra-subject variability for five out of the six muscles. The Sprint and 70% PPO method highlighted greater sensitivity to changes in muscle activity than the MVC method. The MVC method showed the highest intra-subject variability for most muscles except VM.The data suggests that normalising EMG to dynamic methods is the most appropriate for examining muscle activity during cycling over different days and for once-off measurements.     

Myoglobin functions in the heart

The physiological role of myoglobin (Mb) within the heart depends on its oxygenation state. The myocardium exhibits a broad oxygen partial pressure (pO₂) spectrum with a transmural gradient from the epicardial to the subendocardial layer, ranging from arterial values to an average of 19.3 mm Hg down to 0 mm Hg. The function of Mb as an O₂ storage depot is well appreciated, especially during systolic compression. In addition, Mb controls myocardial nitric oxide (NO) homeostasis and thus modulates mitochondrial respiration under physiological and pathological conditions. We recently discovered the role of Mb as a myocardial O₂ sensor; in its oxygenated state Mb scavenges NO, protecting the heart from the deleterious effects of excessive NO. Under hypoxia, however, deoxygenated Mb changes its role from an NO scavenger to an NO producer. The NO produced protects the cell from short phases of hypoxia and from myocardial ischemia/reperfusion injury. In this review we summarize the traditional and novel aspects of Mb and its (patho)physiological role in the heart.     

Evaluation and comparison of four methods of ranking horses based on reactivity

Four methods of ranking horses on reactivity were evaluated and compared: isolation from conspecifics, presentation of a static novel stimulus, traversing a novel stimulus in a runway (isolation, novel stimulus and runway tests, respectively) and assigning subjective emotionality scores. In all tests, horses’ heart rates were recorded and behaviour was videotaped. To be considered a valid test of reactivity, at least one heart rate and one behavioural measurement in the test had to change significantly between treatments (tranquilizer administation versus sham tranquilizer administration), and behavioural measures had to be displayed in at least 75% of the trials. Forty horses performed each of the three tests daily on three different days in a switchback design. Horses were assigned randomly to a daily test sequence, which was maintained throughout the study. In the runway test, no significant difference in heart rate values in tranquilized and non-tranquilized horses was found, and no behavioural attribute was displayed in more than 52% of the trials; therefore it was rejected as a valid test of reactivity. Both isolation and novel stimulus tests produced valid measurements. Mean heart rate was the most precise physiological measure for these tests, and walking and defecation frequency were the most precise behavioural measures for novel stimulus and isolation tests, respectively. Mean heart rates on the novel stimulus and isolation tests were correlated (r_s = 0.79, P < 0.01) indicating that these tests produced similar rankings based on physiological responses. However, behavioural measures ranked horses differently (r_s = 0.27, P < 0.10) on the tests. Rank correlations between mean heart rates and behavioural measures were higher in the novel stimulus (r_s = 0.66, P < 0.01) than the isolation test (r_s = 0.55, P < 0.01), indicating that the novel stimulus test ranked horses based on either physiological or behavioural responses more similarly than did the isolation test. Therefore, the novel stimulus test was considered the more accurate evaluation of reactivity. Subjective emotionality scores were correlated moderately with mean heart rates (r_s > 0.33, P < 0.01) from the novel stimulus and isolation tests and with walking scores (r_s = 0.47, P < 0.01) from the novel stimulus test. Assignment of subjective emotionality scores was not as accurate as the novel stimulus or isolation tests in ranking horses for reactivity. Using physiological data alone, combining physiological and behavioural measurements or using more than one behavioural measurement in reactivity tests may reflect the reactivity of the horse better than a single behavioural measurement.     

Heart rate and arterial pressure variability and baroreflex sensitivity in ovariectomized spontaneously hypertensive rats

AimsThe present study evaluated the effects of ovariectomy on heart rate and arterial pressure variability and cardiac baroreflex sensitivity (BRS) in female spontaneously hypertensive (SHR) and Wistar–Kyoto rats (WKY).Main methodsSham-surgery animals were used as control. Sixteen weeks after ovariectomy or sham-surgery, animals were recorded. Time series of pulse interval (PI) and systolic AP (SAP) were analyzed by means of autoregressive spectral analysis, which quantifies the power of very low (VLF = 0.01–0.25 Hz), low (LF = 0.25–0.75 Hz) and high frequency (HF = 0.75–2.5 Hz) bands. BRS was assessed by means of linear regression between changes of PI and SAP induced by vasoactive drugs or calculation of α-index, a spontaneous BRS index.Key findingsThere was no difference in baseline PI or SAP between ovariectomized and sham SHR. Spectral analysis of heart rate variability suggested a shift of sympatho-vagal balance toward sympathetic predominance in ovariectomized SHR (LF/HF = 1.8 ± 0.2 versus 0.7 ± 0.2 in sham SHR, p < 0.05). Ovariectomy increased total variance and VLF power of SAP in SHR (29.1 ± 9.6 mmHg² and 18.6 ± 6.3 mmHg² versus 9.1 ± 2.1 mmHg² and 4.3 ± 1.4 mmHg², respectively, in sham SHR, p < 0.05). In addition, ovariectomy reduced reflex bradycardia in SHR (0.18 ± 0.03 ms/mmHg versus 0.34 ± 0.06 ms/mmHg in sham SHR, p < 0.05). Ovariectomy did not affect heart rate and SAP variability or BRS in WKY.SignificanceThese data showed that ovarian hormones deprivation induced marked changes on cardiovascular control, increasing SAP variability and cardiac sympatho-vagal balance and blunting BRS in female hypertensive animals, which reinforce the possible protective role of ovarian hormones on the cardiovascular system.     

Changes in plasma orexin A and orexin B concentrations during the estrous cycle of the pig

Orexin A (OXA) and orexin B (OXB) are neuropeptides synthesized mainly in the lateral hypothalamus, which are involved in the control of various physiological functions such as energy homeostasis, sleep, wakefulness and feeding behavior. The present study analyzes orexins A and B levels in the porcine plasma during the estrous cycle. The highest plasma concentrations of orexin A were observed on days 2–3 of the estrous cycle (p < 0.05 relative to days 10–12 and 14–16) and the lowest (p < 0.05) on days 14–16. The highest orexin B levels in the blood plasma were noted on days 17–19 (p < 0.05 vs. days 14–16). We demonstrated the presence of OXA and OXB in porcine blood plasma and the impact of the phase of the estrous cycle on the observed changes in plasma orexin levels.     

Genetic studies in the recently divergent Eligmodontia puerulus and E. moreni (Rodentia,Cricetidae, Sigmodontinae) from Puna and Monte deserts of South America

Eligmodontia is a genus of phyllotine rodents adapted to arid environments with seven recognized species. The sister species E. puerulus and E. moreni are distributed in the adjacent highland Puna and lowland Monte deserts respectively, and show remarkable morphological and chromosomal differences. However, analyses of the cytochrome b gene showed important variability, without reciprocal monophyly between them. In order to study the evolutionary processes involved in the diversification of both taxa, we analyzed 1161 bp of the mitochondrial control region and flanking sequences (N = 60), as well as 759 bp of the first exon of the nuclear gene IRBP (N = 14). Individuals of both species from Jujuy, Catamarca and Mendoza Provinces of Argentina were previously karyotyped. Results showed that the mitochondrial sequences present high haplotype and nucleotide diversity within all population, and no haplotype was shared between both species. F_ST indicated that populations of both species were moderately structured. The network was constituted by two major haplogroups, one composed by E. puerulus samples from Jujuy, and the other composed of sequences of all studied populations. The Bayesian analysis showed three clusters, matching the network. Phylogenetic analysis recovered two clades with high support, in coincidence with the network groups. There was only one close join between sequences of both species, corresponding to samples from Catamarca. Thus, mitochondrial data suggested hybridization between both species in Catamarca, with asymmetric introgression. The IRBP showed low variability and, in the phylogenetic analysis, the sequences of E. puerulus form a monophyletic group with intermediate support, whereas those of E. moreni collapse into a basal polytomy. Our data indicated a recent divergence and absence of introgression in the nuclear genomes. The results at the population level with mitochondrial sequences, together with integrative taxonomy at the species level in a biogeographic context, suggest that climatic and geologic changes could have had an important role in the determination of genetic variability patterns observed in these rodents.     

Heart position variability during voluntary moderate deep inspiration breath-hold radiotherapy for breast cancer determined by repeat CBCT scans

Voluntary moderate deep inspiration breath hold (vmDIBH) in left-sided breast cancer radiotherapy reduces cardiac dose. The aim of this study was to investigate heart position variability in vmDIBH using CBCT and to compare this variability with differences in heart position between vmDIBH and free breathing (FB).For 50 patients initial heart position with respect to the field edge (HP-FE) was measured on a vmDIBH planning CT scan. Breath-hold was monitored using an in-house developed vertical plastic stick. On pre-treatment CBCT scans, heart position variability with respect to the field edge (Δ_HP-FE) was measured, reflecting heart position variability when using an offline correction protocol. After registering the CBCT scan to the planning CT, heart position variability with respect to the chest wall (Δ_HP-CW) was measured, reflecting heart position variability when using an online correction protocol. As a control group, vmDIBH and FB computed tomography (CT) scans were acquired for 30 patients and registering both scans on the chest wall.For 34 out of 50 patients, the average HP-FE and HP-CW increased over the treatment course in comparison to the planning CT. Averaged over all patients and all treatment fractions, the Δ_HP-FE and the Δ_HP-CW was 0.8 ± 4.2 mm (range −9.4–+10.6 mm) and 1.0 ± 4.4 mm (range −8.3–+10.4 mm) respectively. The average gain in heart to chest wall distance was 11.8 ± 4.6 mm when using vmDIBH instead of FB. In conclusion, substantial variability in heart position using vmDIBH was observed during the treatment course.     

Respiratory cycle related EEG changes: Modified respiratory cycle segmentation

Respiratory cycle related EEG change (RCREC) is characterized by significant relative EEG power changes within different stages of respiration during sleep. RCREC has been demonstrated to predict sleepiness in patients with obstructive sleep apnoea and is hypothesized to represent microarousals. As such RCREC may provide a sensitive marker of respiratory arousals. A key step in quantification of RCREC is respiratory signal segmentation which is conventionally based on local maxima and minima of the nasal flow signal. We have investigated an alternative respiratory cycle segmentation method based on inspiratory/expiratory transitions. Sixty two healthy paediatric participants were recruited through staff of local universities in Bolivia. Subjects underwent attended polysomnography on a single night (Compumedics PS2 system). Studies were sleep staged according to standard criteria. C3/A2 EEG channel and time-locked nasal flow (thermistor) were used in RCREC quantification. Forty Seven subjects aged 7–17 (11.4 ± 3) years (24M:23F) were found to have usable polysomnographs for the purpose of RCREC calculation. Respiratory cycles were segmented using both the conventional and novel (transition) methods and differences in RCREC derived from the two methods were compared in each frequency band. Significance of transition RCREC as measured by Fisher's F value through analysis of variance (ANOVA) was found to be significantly higher than the conventional RCREC in all frequency bands (P < 0.05) but beta. This increase in statistical significance of RCREC as demonstrated with the novel transition segmentation approach suggests better alignment of the respiratory cycle segments with the underlying physiology driving RCREC.     

The influence of hydration status on heart rate variability after exercise heat stress

While exercise heat stress and hydration status are known to independently influence heart rate variability (HRV), the combined effect of these physiological stressors is unknown. Thus, heat-acclimated subjects (n=5) performed exercise heat trials (40 °C, 20% relative humidity) in the euhydrated and hypohydrated state (3.9±0.7% body weight loss). During each trial, cardiac cycle R–R interval data were collected for 45 min at rest (pre-) and after (post-) completing 90 min of cycle ergometer exercise. Pre- and post-exercise RRI data were analyzed by Fast Fourier Power Spectral analysis to determine the high-frequency (HF), low-frequency (LF), very low-frequency (VLF), and total power (TP) components of HRV. Overall HRV was decreased by both hypohydration and exercise heat stress. Hypohydration reduced TP, LF, VLF, and LF:HF ratio $(P<0.05)$ while HF was significantly higher. The change in both LF and HF power (pre- vs. post-exercise) were blunted during hypohydration compared to euhydration. These data suggest that dehydration alone positively influences the parasympathetic (HF) control of HRV, but the reduction in overall HRV and the blunted oscillations in LF and HF power following exercise heat stress support an overall deleterious effect of dehydration on autonomic cardiac stability.     

An investigation of the simultaneously recorded occlusal contact and surface electromyographic activity of jaw-closing muscles for patients with temporomandibular disorders and a scissors-bite relationship

Surface electromyographic (SEMG) activity of the masseter and anterior temporalis (TA) muscles has been reported to be associated with occlusion and orofacial pain. However, our recent report did not reveal an association between the side of orofacial pain and the side showing higher or lower level of SEMG activity of masseter or TA. The present purpose was to re-test this association in patients who had unilateral scissors-bite relationship. Thirty-two unilateral scissors-bite femalepatients complaining of unilateral orofacial pain (n = 15) or TMJ sounds (n = 17) were enrolled to simultaneously record contacts, force distribution of occlusion, and SEMG activity of masseter and TA during centric maximal voluntary clenching (MVC). The results indicated that neither orofacial pain nor the TMJ sounds had an association with the masseter’s SEMG values, while scissors-bite had (P < 0.05). A lower SEMG value for masseter was found on the scissors-bite side where there was a smaller number of contacts and a lower biting force distribution (P < 0.05). No such association was revealed in TA. In conclusion, in patients with unilateral TMD symptom(s) and scissors-bite, the jawclosing muscles’ SEMG activity during centric MVC was associated with the scissors-bite rather than the symptoms of orofacial pain or TMJ sounds.     

Improving and evaluating the soil cover indicator for agricultural land in Canada

A soil cover days (SCD) model has been developed by Agriculture and Agri-Food Canada for use as an agri-environmental indicator to monitor the relationship between agricultural production activities and agri-environmental quality. The SCD indicator integrates information on crops, soils, climate, and field activities to estimate the total equivalent number of days that agricultural soils are covered by crop canopy, crop residue and snow in a given year. Daily cover fractions of plant and residue for a given crop in an ecoregion are simulated using typical crop calendar and field management practices, and the equivalent number of days that soil is covered by snow in winter is derived from long term climate normals. The equivalent SCD for a spatial unit is then derived as the area-weighted sum of different crops and different management practices within the unit. This paper presents the SCD framework, details an assessment of the accuracy of the model and outlines future improvements. Annual snow days derived from 30-year climate normals as used in the model was strongly correlated (excluding mountain areas) with that derived from satellite data (R² = 0.45, n = 48), even though the remote sensing product showed significant temporal and spatial variability. Crop residue fraction estimated by the model was strongly correlated with field data collected over major crop areas and crop types (R² = 0.74, n = 55), and modelled plant cover fraction was well correlated with that derived from remote sensing data (R² = 0.57, n = 57). Large discrepancies were observed for some samples due to deviation of the actual crop calendar from that estimated using climate normals. National map showing the change in the indicator from 1981 to 2011 reveals changes in crop and residue management practices.     

Quantified proarrhythmic potential of selected human embryonic stem cell-derived cardiomyocytes

To improve proarrhythmic predictability of preclinical models, we assessed whether human ventricular-like embryonic stem cell-derived cardiomyocytes (hESC-CMs) can be selected following a standardized protocol. Also, we quantified their arrhythmogenic response and compared this to a contemporary used rabbit Purkinje fiber (PF) model. Multiple transmembrane action potentials (AP) were recorded from 164 hESC-CM clusters (9 different batches), and 12 isolated PFs from New Zealand White rabbits. AP duration (APD), early afterdepolarizations (EADs), triangulation (T), and short-term variability of repolarization (STV) were determined on application of the I_Kr blocker E-4031 (0.03/0.1/0.3/1 μM). Isoproterenol (0.1 μM) was used to assess adrenergic response. To validate the phenotype, RNA isolated from atrial- and ventricular-like clusters (n = 8) was analyzed using low-density Taqman arrays. Based on initial experiments, slow beating rate (< 50 bpm) and long APD (> 200 ms) were used to select 31 ventricular-like clusters. E-4031 (1 μM) prolonged APD (31/31) and induced EADs only in clusters with APD90 > 300 ms (11/16). EADs were associated with increased T (1.6 ± 0.2 vs 2.0 ± 0.3?) and STV (2.7 ± 1.5 vs 6.9 ± 1.9?). Rabbit PF reacted in a similar way with regards to EADs (5/12), increased T (1.3 ± 0.1 vs 1.9 ± 0.4?), and STV (1.2 ± 0.9 vs 7.1 ± 5.6?). According to ROC values, hESC-CMs (STV 0.91) could predict EADs at least equivalent to PF (STV 0.69). Isoproterenol shortened APD and completely suppressed EADs. Gene expression analysis revealed that HCN1/2, KCNA5, and GJA5 were higher? in atrial/nodal-like cells, whereas KCNJ2 and SCN1B were higher? in ventricular-like cells (?P < 0.05). Selection of hESC-CM clusters with a ventricular-like phenotype can be standardized. The proarrhythmic results are qualitatively and quantitatively comparable between hESC-CMs and rabbit PF. Our results indicate that additional validation of this new safety pharmacology model is warranted.     

Effect of seat positions on discomfort,muscle activation,pressure distribution and pedal force during cycling

The aim of this study was to measure and analyse discomfort and biomechanics of cycling, i.e., muscle activation, centre of pressure of seat pressure profiles and pedal forces as a function of seat position. Twenty-one recreationally active individuals cycled for 10 min at 100 W on an ergometer cycle using five different seat positions. The neutral position was considered as basic seat position and was compared with upward, downward, forward and backward seat positions. The initial bout was repeated at the end of the recording session. Discomfort increased for upward and backward condition compared with neutral (P < 0.05). Normalized surface electromyography from gastrocnemius decreased in the downward and forward position but increased in the upward and backward position. The minimum force became less negative for forward position compared with neutral seat position (P < 0.05). The degree of variability of centre of pressure increased in the upward and backward position and the entropy of the centre of pressure of sitting posture for backward position decreased compared with neutral seat position (P < 0.05). The present study revealed that consecutive changes of seat position over time lead to increase in discomfort as well as alterations of the biomechanics of cycling.     

EMG normalization to study muscle activation in cycling

The value of electromyography (EMG) is sensitive to many physiological and non-physiological factors. The purpose of the present study was to determine if the torque–velocity test (T–V) can be used to normalize EMG signals into a framework of biological significance. Peak EMG amplitude of gluteus maximus (GMAX), vastus lateralis (VL), rectus femoris (RF), biceps femoris long head (BF), gastrocnemius medialis (GAS) and soleus (SOL) was calculated for nine subjects during isometric maximal voluntary contractions (IMVC) and torque–velocity bicycling tests (T–V). Then, the reference EMG signals obtained from IMVC and T–V bicycling tests were used to normalize the amplitude of the EMG signals collected for 15 different submaximal pedaling conditions. The results of this study showed that the repeatability of the measurements between IMVC (from 10% to 23%) and T–V (from 8% to 20%) was comparable. The amplitude of the peak EMG of VL was 99 ± 43% higher (p < 0.001) when measured during T–V. Moreover, the inter-individual variability of the EMG patterns calculated for submaximal cycling exercises differed significantly when using T–V bicycling normalization method (GMAX: 0.33 ± 0.16 vs. 1.09 ± 0.04, VL: 0.07 ± 0.02 vs. 0.64 ± 0.14, SOL: 0.07 ± 0.03 vs. 1.00 ± 0.07, RF: 1.21 ± 0.20 vs. 0.92 ± 0.13, BF: 1.47 ± 0.47 vs. 0.84 ± 0.11). It was concluded that T–V bicycling test offers the advantage to be less time and energy-consuming and to be as repeatable as IMVC tests to measure peak EMG amplitude. Furthermore, this normalization method avoids the impact of non-physiological factors on the amplitude of the EMG signals so that it allows quantifying better the activation level of lower limb muscles and the variability of the EMG patterns during submaximal bicycling exercises.     

Enfermedad cardiovascular tras infección por SARS-CoV-2 en pacientes ancianos. Resultados del seguimiento anual de una cohorte de supervivientes

IntroductionAlthough the effects of SARS-CoV-2 infection on the cardiovascular system is well known in the acute phase, the cardiovascular impact of the elderly population surviving COVID-19 respiratory infection after 1 year of follow-up has not been sufficiently studied.MethodsObservational registry of 240 elderly patients (75 years or older), consecutively admitted for COVID-19 respiratory infection and survivors of the same, between March 1 and April 30, 2020, at the Hospital General Universitario de Ciudad Real. The incidence of major cardiovascular events [MACE] (cardiovascular death [CD], acute coronary syndrome [ACS], cerebrovascular disease [CVD], venous thromboembolic disease [VTE] and heart failure [HF]) was prospectively analysed.ResultsThe mean age was 83.75 ± 5.75 years. After a mean follow-up of 352.2 ± 70.4 days, 13.8% of patients died and 9.6% had MACE, the most frequent being heart failure, with no differences in severity or overall course of acute illness. In the multivariate Cox regression model, the risk of developing MACE was higher in patients with chronic obstructive pulmonary disease and (HR 4.29; 95%CI 1.62-11.39; P = .003) and loop diuretic (HR 2.99; 95%CI 1.27-7.07; P = .01).ConclusionsIn elderly COVID-19 survivors, the incidence of MACE after one year of follow-up is high, the main manifestation being heart failure.     

An electrophysiological study of the sequential water swallowing

BackgroundSequential water swallowing (SWS) was mostly investigated by the videofluoroscopic and endoscopic methods. However some physiological features of SWS was rarely evaluated by neurophysiological techniques. Our aim was to investigate some neural and muscular changes on sequences of SWS using electromyography (EMG) methods.MethodsFifty-eight normal adults were investigated. SWS was initiated voluntarily with 50 ml and 100 ml water volumes from a cup. Submental EMG, respiratory signals, heart rate, and sympathetic skin responses (SSR) were measured during SWS.Key resultsAll parameters were increased significantly during the 100 ml SWS. During swallowing apnea period, compensatory respiration cycles occurred in 24% and 48% of participants in the 50 ml and 100 ml SWS, respectively. Heart rate increased during swallowing apnea. SSR were evoked just before and just after the SWS in more than halves of participants. A foreburst EMG in SM muscles at the initiation of SWS was recorded in 86% of normal participants. Older age was associated with a prolonged duration of the apnea period.ConclusionsAll parameters of the SWS could be recorded numerically and objectively using electrphysiological methods. These are similar to those obtained by videofluoroscopic and similar methods. The foreburst activity of the initiation of SWS may represents preparatory activity from the activation of the fast cortical descending motor pathway. Increasing heart rate and the prolonged apnea urged that older people and patients could be carefully tested for respiratory and cardiac rhythm disorders.     

Carbon and oxygen isotope evidence for high-frequency (104–105 yr) and My-scale glacio-eustasy in Middle Pennsylvanian cyclic carbonates (Gray Mesa Formation), central New Mexico

We combine cyclo- and sequence stratigraphy along with whole rock δ¹³C and conodont apatite δ¹⁸O analysis to document high-frequency (10⁴–10⁵ yr) and My-scale sea-level changes for the Middle Pennsylvanian (Desmoinesian or Moscovian) Gray Mesa Formation of central New Mexico. Approximately 75 subtidal cycles (1–8 m) are grouped into 4 1/2 My-scale depositional sequences (40–80 m). About 50% of the cycles show evidence of prolonged subaerial exposure at cycle tops with the development of calcretes, diagenetic mottling, and regolith intraclasts. High-resolution δ¹³C analysis of whole rock limestones across nine of the cycles indicates that the cycle tops were diagenetically altered by isotopically light, meteoric fluids during sea-level fall and lowstand. These δ¹³C trends support the interpretation that high-frequency sea-level changes were responsible for cycle development.Conodont apatite δ¹⁸O values from sampled cycles indicate that the high-frequency sea-level changes were driven by glacio-eustasy combined with changes in surface seawater temperature (SST). δ¹⁸O values from conodont apatite, spanning parts of three depositional sequences indicate that My-scale glacio-eustasy and/or SST changes controlled sequence development. δ¹⁸O shifts indicate that the magnitudes of 10⁴–10⁵ yr glacio-eustasy were between ～ 55 and 170⁺ m combined with tropical SST changes of ～ 1.5°–6 °C. Calculated My-scale glacio-eustatic oscillations were between ～ 60 and 140 m with SST changes of < 3.5 °C. The most plausible driver for the My-scale paleoclimate changes is long-period obliquity (～ 1.2 My) variations. These calculated high-frequency, glacio-eustatic values are similar or greater than Pleistocene values, and lie within the range estimated for other Middle Pennsylvanian successions using a variety of independent eustatic proxies. The similarity in range of magnitudes between high-frequency and My-scale sea-level changes combined with the large differences in magnitudes between individual high-frequency sea-level oscillations helps explain the lack of systematic cycle-stacking patterns within these Pennsylvanian icehouse sequences.