The intensity of exercise above which oxygen uptake (Vo2) does not account for all of the required energy to perform work has been associated with lactate accumulation in the blood (lactate threshold, LT) and elevated carbon dioxide output (gas exchange threshold). An increase in hydrogen ion concentration ([H+]) is approximately concurrent with elevation of blood lactate and CO2 output during exercise. Near-infrared spectra (NIRS) and invasive interstitial fluid pH (pHm) were measured in the flexor digitorum profundus during handgrip exercise to produce a mathematical model relating the two measures with an estimated error of 0.035 pH units. This NIRS pHm model was subsequently applied to spectra collected from the vastus lateralis of 10 subjects performing an incremental-intensity cycle protocol. Muscle oxygen saturation (SmO2) was also calculated from spectra. We hypothesized that a H+ threshold could be identified for these subjects and that it would be different from but correlated with the LT. Lactate, gas exchange, SmO2, and H+ thresholds were determined as a function of Vo2 using bilinear regression. LT was significantly different from both the gas exchange threshold (Delta = 0.27 +/- 0.29 l/min) and H+ threshold (Delta = 0.29 +/- 0.23 l/min), but the gas exchange threshold was not significantly different from the H+ threshold (Delta = 0.00 +/- 0.38 l/min). The H+ threshold was strongly correlated with LT (R2 = 0.95) and the gas exchange threshold (R2 = 0.85). This initial study demonstrates the feasibility of noninvasive pHm estimations, the determination of H+ threshold, and the relationship between H+ and classical metabolic thresholds during incremental exercise.

A method to non-invasively and quantitatively measure muscle oxygen saturation (SmO(2)) using broadband continuous-wave diffuse reflectance near infrared (NIR) spectroscopy is presented. The method obtained SmO(2) by first correcting NIR spectra for absorption and scattering of skin pigment and fat, then fitting to a Taylor expansion attenuation model. A non-linear least squares optimization algorithm with set boundary constraints on the fitting parameters was used to fit the model to the acquired spectra. A data preprocessing/optimization scheme for accurately determining the initial values needed for the optimization was also employed. The method was evaluated on simulated muscle spectra with 4 different scattering properties, as well as on in vivo forearm spectra from 5 healthy volunteer subjects during arterial occlusion. Measurement repeatability was assessed on 24 healthy volunteers with 5 repeated measurements, each separated by at least 48 hours.