Potential Drivers of Mid-Infrared Variability in Young Stars: Testing Physical Models with Multiepoch Near-Infrared Spectra of YSOs in ρ Oph

Recent studies have identified several young stellar objects (YSOs) which exhibit significant mid-infrared (mid-IR) variability. A wide range of physical mechanisms may be responsible for these variations, including changes in a YSO's accretion rate or in the extinction or emission from the inner disk. We have obtained and analyzed multiepoch near-infrared (NIR) spectra for five actively accreting YSOs in the rho Oph star-forming region along with contemporaneous mid-IR light curves obtained as part of the YSOVAR Spitzer/IRAC survey. Four of the five YSOs exhibit mid-IR light curves with modest (similar to 0.2-0.4 mag) but statistically significant variations over our 40-day observation window. Measuring the strengths of prominent photospheric absorption lines and accretion sensitive H I and He I lines in each NIR spectrum, we derive estimates of each YSO's spectral type, effective temperature (T-eff), and H-band extinction (A(H)), and analyze the time evolution of their NIR veiling (r(H) and r(K)) and mass accretion rates ((M) over dot(acc)). Defining a YSO's evolutionary stage such that heavily veiled, high accretion rate objects are less evolved than those with lower levels of veiling and ongoing accretion, we infer that GY 314 is the most evolved YSO in our sample, with GY 308 and GY 292 at progressively earlier evolutionary stages. Leveraging our multiepoch, multiwavelength dataset, we detect significant variations in mass accretion rates over timescales of days to weeks, but find that extinction levels in these YSOs remain relatively constant. We find no correlation between these YSO mid-IR light curves and time-resolved veiling or mass accretion rates, such that we are unable to link their mid-IR variability with physical processes localized near the inner edge of the circumstellar disk or within regions which are directly responsive to mass accretion. We do find, however, that redshifted He I lambda 10830 emission, where present in our spectra, shows both quantitative and qualitative temporal correlations with accretion-sensitive H I emission lines. Blueshifted He I absorption, on the other hand, does not demonstrate a similar correlation, although the time-averaged strength of this blueshifted absorption is correlated with the time-averaged accretion rate in our sample of YSOs.