Optical and near-infrared variability is a well-known property of young stellar objects. However, a growing number of recent studies claim that a considerable fraction of them also exhibit mid-infrared flux changes. With the aim of studying -2and interpreting variability on a decadal timescale, here we present a mid-infrared spectral atlas containing observations of 68 low-and intermediate mass young stellar objects. The atlas consists of 2.5-11.6 µm low-resolution spectra obtained with the ISOPHOT-S instrument on-board the Infrared Space Observatory (ISO) between 1996 and 1998, as well as 5.2-14.5 µm low-resolution spectra obtained with the IRS instrument on-board the Spitzer Space Telescope between 2004 and 2007. The observations were retrieved from the ISO and Spitzer archives and were post-processed interactively by our own routines. For those 47 objects where multi-epoch spectra were available, we analyze mid-infrared spectral variability on annual and/or decadal timescales. We identify 37 variable candidate sources. Many stars show wavelength-independent flux changes, possibly due to variable accretion rate. In several systems, all exhibiting 10 µm silicate emission, the variability of the 6−8 µm continuum and the silicate feature exhibit different amplitudes. A possible explanation is variable shadowing of the silicate emitting region by an inner disk structure of changing height or extra silicate emission from dust clouds in the disk atmosphere. Our results suggest that mid-infrared variability, in particular the wavelength-dependent changes, are more ubiquitous than was known before. Interpreting this variability is a new possibility to explore the structure of the disk and its dynamical processes.The Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope had two low resolution and two high resolution modules, covering the 5.2-38 µm wavelength range. Since our aim is to look for mid-IR variability, here we utilize the Short-Low channel that covers the 5.2-14.5 µm wavelength range, thus has the largest wavelength overlap with ISOPHOT-S. The Short-Low channel has a spectral resolution of R ≈ 60-127, comparable to that of ISOPHOT-S. The width and length of the IRS Short-Low slit is 3. ′′ 6 and 57 ′′ , respectively. According to the IRS instrument handbook (Version 4.0), the uncertainty of the absolute flux calibration is 5-10%, while the repeatability is 2-5%.We downloaded basic calibrated data (BCD) files from the Spitzer archive, processed with the pipeline version 18.18.0. Measurements were taken either in staring mode (two nod positions), or in spectral mapping mode (small 2×3 maps). In the latter case, we used the central 2 positions, and considered them as normal staring measurements. We first subtracted the two nod positions from each other. Then we used the Spitzer IRS Custom Extraction software (SPICE) to extract the positive signal from the two-dimensional dispersed images. We extracted the spectra from a wavelength-dependent, tapered aperture around the target, and we averaged the two s...