stars form in dense and cool molecular clouds .
when the local density is high enough , the matter can gravitationally collapse and form a young stellar object ( yso ) . in the early phases , the thick envelope dominates the emission from the yso and hides what is going on within ( class  i ) .
eventually , the envelope flattens out to a circumstellar accretion disk .
this disk still causes an infrared ( ir ) excess above the level of a stellar photosphere ( class  ii or classical t tauri star - ctts ) , which can be used to distinguish those objects from main - sequence stars , for example using the _ spitzer _ space telescope @xcite infrared array camera ( irac , * ? ? ?
* ) . when the disk is cleared , the ir colors of the yso match those of main - sequence stars ( class  iii or weak - lined t tauri star - wtts ) .
in addition to the circumstellar absorption , many ysos are embedded in the molecular cloud , so that even class  iii objects can appear reddened .
the accretion disk does not reach down to the central star . instead
, the inner edge of the gas disk is truncated by the stellar magnetic field .
the inner radius of the optically thick dust in the disk is larger than the inner radius of the gas disk and mostly given by the dust - sublimation temperature .
some of the mass in the circumstellar disk condenses into planets , some is blown out by accretion - driven disk and stellar winds , and is accreted onto the central star .
this accretion can happen via magnetically confined accretion funnels ( e.g. , * ? ? ? * ) or via some magneto - hydrodynamical instability ( e.g. , * ? ? ?
t  tauri stars ( tts ) were originally identified by their variability @xcite  long before anybody realized that tts are indeed pre - main sequence stars .
the dominant timescale in the optical is the stellar rotation period , typically a few days to a week or more .
ysos can have cool spots caused by magnetic activity similar to our sun and also hot spots which mark the impact points of the accretion funnels onto the stellar surface ( see , e.g. , review by * ? ? ?
this impact happens at free - fall velocities up to 500  km  s@xmath1 ; thus , the accretion shock heats the accreted mass to x - ray emitting temperatures . in the optical
, the accretion region appears as emission that often is approximated as a blackbody with temperature @xmath2  k ( @xcite , but see also @xcite who argue that line emission contributes to the veiling in addition to a continuum ) .
variability in the mass accretion rate can lead to changes in the hot spot signatures .
the dynamical timescale that controls the accretion is the keplerian period of the inner disk where the accretion funnels start .
the inner disk radius is found close to the co - rotation radius leading to a typical timescale of a few hours for typical masses and rotation periods of ysos . indeed optical variability with amplitudes around 0.1  mag
is often observed in ctts on this timescale @xcite .
another source of variability related to the accretion could be oscillations of the accretion shock on timescales of seconds .
this has been predicted theoretically ( e.g. , * ? ? ?
* ) , but is not observed so far , possibly because the accretion spot separates into many small funnels that oscillate independently at different phases and frequencies .
however , @xcite find indications that strong accretion in v1647  ori could excite radial pulsations of the star itself .
one of the largest classes of short timescale ( @xmath3 days ) optical and ir variability in ysos is that due to variable extinction events @xcite .
these come in three categories - aa tau - type variables ( stars with broad , periodic flux dips , whose amplitudes can be up to a magnitude or more in the optical ) , presumably due to our line of sight passing through a warp in the inner circumstellar disk ; stars with similar or narrower flux dips that have no obvious periodicity - presumably due to stars where our line of sight passes close to the disk and where disk instabilities can levitate dust high enough above the plane to intersect our line of sight briefly ; and stars with narrow , periodic flux dips - perhaps where our line of sight is being intersected by dust entrained in material accreting onto the star in a funnel flow . about 20% of the ysos in ngc  2264 fall into one of these categories in the sample of @xcite .
ysos can also vary on much longer timescales .
variability on the timescales of years could be caused by changing circumstellar extinction for a keplarian disk around a solar - mass yso this timescale translates to a radius of a few au ) or by massive accretion events when a significant fraction of the disk mass drains onto the yso @xcite . in this case
the accretion luminosity can outshine the yso by orders of magnitudes and it takes months to years ( in the case of exor outbursts , * ? ? ? * ) or even centuries until the accretion decays back down to the original level .
in any lightcurve , several of the processes dicussed above can contibute to the observed variability at the same time and it depends on the properties of each object which one dominates and if secondary effects can be detected in the lightcurve .
for example , cool spots , hot spots , absorption , and massive accretion events can all influence the same optical light curve .
another case are x - rays , where the flux and the spectrum can change due to periodic absorption , variability in the accretion rate ( tw  hya , * ? ? ?
* ) or coronal activity similar to what is seen on the sun . in most ysos ,
the last point is dominant and x - ray lightcurves often show the fast rise in flux and temperature and a slower decay characteristic of coronal activity ( see , e.g. , the _ chandra _ monitoring of the orion nebular cluster , * ? ? ?
* ; * ? ? ?
* ) .    the spectral energy distribution ( sed ) of ysos in the optical is dominated by the stellar photosphere and the accretion spot .
thus , optical monitoring is very effective for understanding the stellar rotation and the accreting spot .
however , the disk radiates mostly at longer wavelengths , which are probed in the ir observations presented in this article .
depending on the mass of the disk and the size of the inner hole , the disk will start to dominate the sed at the @xmath4 band or in the irac bands at @xmath5 m and @xmath6 m .
simple disk models still treat the disk as a static and axisymmetric structure , but observationally it now seems that the disk is in fact `` a bubbling , boiling , wrinkled , dented , warped mass of gas and dust '' @xcite , see also @xcite .
this paper is part of the ysovar ( young stellar object variability ) project , which has monitored the ( onc ) and eleven smaller star forming regions with irac in @xmath5 m and @xmath6 m to understand the mid - ir variability of ysos .
first results on the onc are published in @xcite .
more details of the observing strategy and an overview of the data can be found in @xcite ( from now on `` paper  i '' ) . comparing data from all clusters , paper
i defines certain cut - off values for the data reduction , e.g.  how much variability in a lightcurve is required to reliably identify an object as variable . in the analysis ,
paper  i concentrates on variability in the ir on timescales of years . in this article
, we present a _ spitzer_/irac monitoring campaign of the star forming region l1688 in the mid - ir to characterize the variability timescales and amplitudes as well as the color changes in the mid - ir in much more detail for the objects in l1688 than paper
i on timescales up to two years .
the structure of this paper is as follows : first , we introduce l1688 , the star forming region targeted by these observations ( section  [ sect : thestarformingregionl1688 ] ) . in section  [
sect : observationsdatareductionandauxiliarydata ] we introduce the data reduction and discuss source lists and stellar properties obtained from the literature .
section  [ sect : midirvariability ] classifies all sources according to their variability .
section  [ sect : resultsanddiscussion ] presents our results and discusses physical models to explain the observed features in the lightcurves .
we end with a summary and some conclusions in section  [ sect : summary ] .
lynds  1688 ( l1688 ) is a sub - cloud of the @xmath7  ophiuchus star forming region , one of the best - studied young clusters in the sky ( see e.g.  review by * ? ? ?
the central region of l1688 is very dense and deeply embedded ( @xmath8  mag , see figure  [ fig : map ] ) .
thus , all surveys of the regions necessarily miss some cluster members .
an extinction limited , spectroscopic survey @xcite finds an average age of 3.1  myr for a 6.8  pc@xmath9 region centered on l1688 and no significant deviation from the initial mass function .
earlier studies concentrated on the deeper embedded core and found a much younger age of 0.3  myr .
some , but not all , of this difference is due to the specific reddening laws or pre - main sequence evolutionary tracks used in these studies ( see discussion in * ? ? ?
@xcite present photometric and spectroscopic monitoring for five ysos in l1688 nearly simultaneous with the _ spitzer _ observations discussed here .
they do not see any correlation between the hydrogen emission lines that are usually considered accretion indicators and the features in the ir lightcurves of their targets , indicating that the relatively modest variability they observed is not caused by changes in the accretion rate .
additional notable objects with well - sampled nir lightcurves are @xcite and , which show eclipses with periods of 131 and 93  days , respectively . these sources can be interpreted as multiple systems , where one or more components are eclipsed by a warped circumstellar disk .
thus far , the most comprehensive study of near - ir variability of ysos in l1688 is @xcite ( see references therein for other ir variability studies ) , who make use of a two micron all sky survey ( 2mass * ? ? ?
* ) calibration field that overlaps l1688 so each source has up to 1584 datapoints in @xmath10 , @xmath11 , and @xmath4 spanning 2.5  years .
they find 79% of the known ysos to be variable . in total ,
32 sources are periodic ( including cool starspots , hot accretion spots and 6 systems with eclipses ) , 31 sources show a long - term trend and 40 sources vary aperiodically on shorter timescales .
the new data presented in this article complements the @xcite study with observations at longer wavelengths .
in this section we briefly describe the data reduction for the _ spitzer _ and _ chandra _ observations .
a detailed account of the observations , data processing and the source extraction is given in paper  i. we also give an overview of auxiliary data on the stellar properties retrieved from the literature , which we need to test if the variability characteristics depend on the central star .
we then assess the cluster membership and sed class for every source with a usable lightcurve .
ccc 29267200 & 2010 - 04 - 12 11:00:41 & 2010 - 04 - 12 11:18:44 + 29266688 & 2010 - 04 - 12 17:02:37 & 2010 - 04 - 12 17:28:01 + 29266176 & 2010 - 04 - 12 22:48:49 & 2010 - 04 - 12 23:09:28 + 29265664 & 2010 - 04 - 13 08:44:31 & 2010 - 04 - 13 09:05:35 + 29265408 & 2010 - 04 - 13 20:09:17 & 2010 - 04 - 13 20:26:33 +    three fields in l1688 were observed with _
spitzer _ in four observing windows from 2010 - 04 - 12 to 2010 - 05 - 16 ( visibility window 1 ) , 2010 - 09 - 22 to 2010 - 10 - 27 ( visibility window 2 ) , 2011 - 04 - 20 to 2011 - 05 - 23 ( visibility window 3 ) , and 2011 - 10 - 01 to 2011 - 11 - 06 ( visibility window 4 ) .
these windows are consecutive visibility periods dictated by the _ spitzer _
orbit @xcite .
they are shown visually in figure  [ fig : cadence ] . in the first visibility window ,
the sampling is much denser in time than in the later visibility windows . for most sources ,
about 70 observations with irregular time intervals to aid period detection were obtained in visibility window 1 . in the first visibility window
, we use a repeating pattern of 8 observations every 3.5 days . within the 3.5 day period the time - step increases from 2 to 16 hours . in the last three visibility windows , the time steps increase linearly with step lengths of roughly 1 , 2 , 3 , ... days to again sample multiple variability frequencies equally .
this irregular sampling minimizes aliasing with the observation frequency .
less than ten datapoints per source are taken in each of the later three visibility windows . in total , there are 108 observations with a total mapping time of 30.7 hours
. table  [ tab : aor ] lists the time of each observation .
they can be found under program identification number ( pid ) 61024 in the spitzer heritage archive .
each observation consists of six dithers in irac mapping mode using the high - dynamic - range ( hdr ) data acquisition mode which obtains a 0.4  s and a 10.4  s exposure for each pointing .
the three fields chosen were observed with the irac  1 and irac  2 channels ( effective wavelengths 3.6@xmath12 m and 4.5@xmath12 m ) .
both channels operate simultaneously , but their fields - of - view are non - overlapping . thus ,
each target field is observed in two consecutive pointings , one for each channel .
a secondary field is observed in the secondary channel while the primary channel is observing the target field . in visibility windows 1 and 3 ,
sources south of the main fields have only irac  1 data , while those to the north only have irac  2 data . in visibility windows 2 and 4 ,
the situation is reversed .
not all sources in the central fields have usable data in both bands , because they might be too bright or too faint in one channel , or fall on the edge of the map .
additional sources with two band coverage are found where the northern side field of one target field overlaps with the southern side field of another target field .
here we summarize the main data reduction and processing steps described in detail in paper  i. basic calibrated data ( bcd ) are obtained from the _ spitzer _ archive .
further data reduction is performed with the idl package cluster grinder @xcite , that treats each bcd image for bright source artifacts .
aperture photometry is performed on individual bcds with an aperture radius of @xmath13 . to increase the signal - to - noise ratio and to reject cosmic rays , the photometry from all bcds in each observation
is combined .
the reported value is the average brightness of all bcds within that observation that contain the source in question , after rejecting outliers .
the photometric uncertainties obtained from the aperture photometry are , particularly for faint sources , only lower limits to the total uncertainty , since distributed nebulosity often found in star forming regions can contribute to the noise .
to improve these estimates , paper  i introduces an error floor value that is added in quadrature to the uncertainties of individual photometric points .
the value of the error floor is 0.01  mag for irac  1 and 0.007  mag for irac  2 .
we cross - match sources from individual observations with a matching radius of @xmath14 with each other and with the 2mass catalog , which is used as a coordinate reference .
all photometric measurements performed in the context of the ysovar project are collected in a central database , which we intend to deliver to the infrared science archive ( irsa ) for general distribution .
data for this article were retrieved from the ysovar database on 2013 - 10 - 31 and further processed using custom routines in python available at https://github.com/ysovar .
we visually checked all frames for lightcurves that are classified as variable in section  [ sect : midirvariability ] and removed datapoints visibly affected by instrumental artifacts ( cosmic rays , read - out streaks for bright neighbors ) .
figure  [ fig : map ] shows the distribution of the sources with lightcurves in our input catalog overlayed on a larger irac  1 map observed during the cryogenic mission .    in this article
, we consider only objects that have at least five datapoints in our irac  1 or irac  2 lightcurves .
a stricter definition is employed in paper  i , where only sources with more than five datapoints _ in the fast - cadence data ( first visibility window ) _ are used . in l1688 , 822 of the total list of 882 sources fullfill this stricter condition . in table
[ tab : tab2 ] they are marked in the column ` standardset ` .      despite the careful data reduction described above
, some residual artifacts remain in the _ spitzer _ lightcurves . in this section
, we search for artifacts that are related to the position on the detector . compared with other clusters in the ysovar project
, l1688 is particularly suited to discover these kinds of effects because the fields observed in l1688 have almost no rotation within one visibility window .
the spacecraft orientation flips between visibility windows , so that most instrumental artifacts produce lightcurves that have one level in visibility window 1 and 3 and another level in visibility window 2 and 4 .
figure  [ fig : instrumentallcs ] shows sources with a magnitude between 8 and 15 in @xmath6 m , that might fall in this category .
the relative difference in magnitude is smaller for brighter sources , but since the photometric errors are also smaller these instrumental effects can still be significant .
we visually inspected every frame for a sample of sources with lightcurves similar to those in figure  [ fig : instrumentallcs ] to identify the cause for the artifacts .
we found that several different effects can cause these steps in the lightcurve .
some sources are close to a detector edge , such that the background is not well - determined , some are in the wings of the point - spread function ( psf ) of a bright neighbor , and some show residuals from hot pixels .
proximity to a detector edge or hot pixel affects only one of the spacecraft orientations ; even the psf wings change with the orientation since the psf is not circular .
sufficiently strong intrinsic variability can mask this offset between visibility windows and for any individual source , this instrumental effect can not be distinguished from intrinsic variability with a period of one year .
thus , we use a statistical approach to quantify the number of lightcurves that suffer from this problem . for each channel
, we have about 200 lightcurves with datapoints in all four observing visibility windows ( the total number of lightcurves is larger , but only for sources in the three target fields do we have data for all four visibility windows ) .
if the mean magnitude in a visibility window depends on the detector position of the source , this source will on average be brighter in visibility window 1 and 3 and fainter in visibility window 2 and 4 ( or vice - versa ) .
we calculate the mean magnitude in each visibility window and test whether the two brightest mean magnitudes belong to visibility window [ 1,3 ] or [ 2,4 ] or any other combination ( [ 1,2 ] , [ 1,4 ] , [ 2,3 ] , [ 3,4 ] ) .
we use the following abbreviated notation : when we calculate the mean magnitudes in each visibility window for a sample of @xmath15 sources , then @xmath16:@xmath17:@xmath18 means that @xmath16 sources have the brightest mean magnitudes in visibility windows [ 1,3 ] , @xmath17 sources in visibility windows [ 2,4 ] and @xmath18 sources have their two brightest mean magnitudes in any other combination of visibility windows .
if the difference in the mean magnitude between visibility windows is unrelated to the position on the detector , we expect the ratio @xmath1917%:17%:67% .
in contrast , the observed lightcurves have the ratio 47:66:75=25%:35%:40% for irac  1 and 52:59:95=25%:29%:46% for irac  2 .
this is incompatible with the expected multinomial distribution ( the probability to observe a distribution at least as far from the expected @xmath20 by chance is @xmath21 ) . in each case , about 50 lightcurves , a quarter of the sample , need to be shifted from the first two bins to the last bin to make the observed distribution compatible with the expected distribution .
this implies that about a quarter of all lightcurves suffer from the artifacts described above .
however , both irac channels use independent detectors and thus the chance that both channels are affected for the same source and that the effect goes the same way ( bright - faint , vs. faint - bright ) is low .    below in section
[ sect : midirvariability ] we identify variable sources using a stetson and a @xmath22 test .
the limits in those tests are designed to be conservative and indeed we find a distribution of 5:7:24=14%:19%:67% for irac  1 and 5:5:27=14%:14%:73% for irac  2 , when we restrict the sample to those lightcurves that will be classified as variable below .
both ratios are fully compatible with the expected multinomial distribution @xmath20 .
this shows that the limits we apply are conservative enough that the sample of stars we identify as variable has no or only few sources where the variability is not due to intrinsic source variability .
therefore it is not necessary to remove any source based on magnitude jumps between visibility windows .
only one of the four lightcurves shown in figure  [ fig : instrumentallcs ] ( sstysv j162727.53 - 242611.2 ) will be identified as variable below .
we visually inspected all lightcurves that are marked as variable and , apart form the example shown above , we did not see lightcurves where the variability seems to be due to the pattern discussed in this section .    in summary ,
about a quarter of all lightcurves are affected by detector position dependent artifacts . in extreme cases ,
the associated jumps reach 0.5  mag for a 14  mag source , and up to 0.05  mag in bright sources around 8  mag .
however , we show statistically that the definition of variability we use is so conservative that this effect does not contribute a significant number of objects to our sample of variable sources .
disk - bearing ysos can be identified from _ spitzer _
data alone , but information in other spectral bands is required to find the other cluster members .
x - ray observations are one way to identify diskless ( class  iii ) ysos . in the ir
the seds of those sources are indistinguishable from a main - sequence field star , but due to their rapid rotation , ysos are much brighter in x - rays than field stars .
l1688 was observed by _
chandra _ on 2000 - 04 - 13 for 100  ks exposure time in the ` faint ` mode with the acis instrument ( ) .
we reprocessed this exposure with the anchors pipeline @xcite using a recent calibration ; see discussion in paper  i. these data has been analyzed in detail to study the distribution of x - ray properties in ctts and to identify brown dwarfs in this star forming region @xcite .    during the observation ,
five acis chips were operational , four from the central acis - i imaging array , as well as one acis - s chip .
the point - spread - function ( psf ) degrades significantly for sources located off - axis , and thus the coordinates of the outer sources are less reliable .
to cross - match x - ray sources with _
sources we used a matching radius of 1 for x - ray sources within 3 of the optical axis of _ chandra _ , 1.5 for sources between 3 and 6 away from the optical axis and 2 for all sources located further than 6 from the optical axis .
the observed _ chandra _ field overlaps about two thirds of the area covered in the _ spitzer _ monitoring .
however , the variable psf leads to a sensitivity that varies over the observed field , thus the absence of an x - ray detection for _ spitzer _ sources does not necessarily imply the absence of x - ray emission .
the acis detector has an intrinsic energy resolution , and we use the net flux to characterize the x - ray properties . for sources with more than 20 counts , we also fit an absorbed single - temperature apec model @xcite with abundances fixed at 0.3 times the solar value from @xcite .
sources are extracted down to a very low significance . in total , there are 315 detected x - ray sources , but only 31 of them match an object with a _ spitzer _
we disregard all unmatched sources ; for sources with a lightcurve and an x - ray counterpart , the x - ray properties are given in table  [ tab : tab2 ] . to estimate the number of spurious matches , we multiply the fraction of the total survey area that is included in the positional error circles of the x - ray sources with the number of _ spitzer _ sources with lightcurves .
the result is the average number of spurious matches .
we expect at most 2 - 4 _ spitzer _ sources to be matched to a spurious x - ray source .
x - ray sources that are cross - matched successfully are marked in figure  [ fig : map ] .
the star forming region l1688 has been the target of intense study over the past decades and a wealth of additional information exists in the literature . in particular
, we refer the reader to two reviews @xcite .
the latter review compiles a list of objects with a high probability of membership from a variety of published sources .
the membership criteria employed are ( i )  x - ray emission , which at the distance of l1688 is detectable only from young , and thus active stars ; ( ii )  optical spectroscopy , with h@xmath23 in emission or li in absorption ; ( iii )  a location above the main - sequence in the hr diagram ; or ( iv )  ir emission that is indicative of a circumstellar disk .
l1688 was also observed with _
spitzer _ in the cryogenic mission phase with all four irac channels and the 24@xmath24 m channel of the multiband imaging photometer for spitzer ( mips , * ? ? ?
objects classified as ysos from these data @xcite are already contained in the membership list of @xcite .
we augment our own _ spitzer _
data reduction with values from the catalog published by the c2d project ( c2d =  from cores to disks  ; * ? ? ?
if we did not obtain a photometric value for an irac band or the @xmath25 m mips , but a value with the quality specifier a , b or c is present in the c2d catalog then we use that value . the data are given in table  [ tab : tab2 ] , which specifies if a datapoint is taken from our own data reduction (  g09  : using the pipeline from * ? ? ? * ) or the c2d database .    near - ir data
is taken from 2mass @xcite and cross - matched by the cluster grinder pipeline @xcite .
additionally , we take detections from the ukirt infrared deep sky survey ( ukidss ) galactic cluster survey , data release 9 .
the ukidss project is defined in @xcite .
ukidss uses the united kingdom infrared telescope ( ukirt ) wide field camera and a photometric system described in @xcite .
the pipeline processing and science archive are described in @xcite .
we only retain detections with a ` mergedclass ` flag between -3 and 0 .
sources brighter than @xmath26  mag can be saturated and we discard the ukidss data for any source that lies above this threshold either in ukidss or 2mass .
this leaves us with 2mass data in the range @xmath27  mag and ukidss @xmath28  mag .
the luminosity function for both dataset is almost identical from @xmath26 to @xmath29 . for fainter sources ,
2mass is incomplete ; ukidss is incomplete for @xmath30  mag .
the ysovar data is also cross - matched with data from the simbad service to provide an identification with known objects from the literature . in all cases the matching radius is set to @xmath14 . if a catalog contains multiple entries within @xmath14 of a ysovar source , we match it to the closest catalog entry . in some cases
the best cross - match is not obvious .
those sources are discussed in appendix  [ sect : remarksaboutcrossmatchingindividualsources ] .
ukidss has a better spatial resolution than our irac data .
there are nine sources where more than one ukidss source is found within the size of the aperture we use for irac photometry . in six cases ( wsb 52 , iso - oph 152 , iso - oph 131 , sstysv j162728.13 - 243719.6 , sstysv j162718.11 - 244814.1 , sstysv j162718.25 - 244955.8 )
the second source is visible in the @xmath31 band , the ukidss band that is closest in wavelength to the irac data , so it is likely that both sources contribute to the observed irac emission . in the remaining three cases ( iso - oph 28 , iso - oph 57 , sstysv j162741.14 - 242038.3 ) the second source is not visible at @xmath31 band .
ccccc 1 & ra & deg &  & j2000.0 right ascension + 2 & dec & deg &  & j2000.0 declination + 3 & name &  &  & identifier for object + 4 & iau_name &  &  & j2000.0 iau designation within the ysovar program + 5 & other_names &  &  & alternative identifiers for object + 6 & c2d_id & &  &  + 7 & wil08_id &  &  & wilking et al .
( 2008 ) + 8 & adoc08_aoc &  &  & j2000.0 iau designation ( jhhmmss.ss+ddmmss.s ) + 9 & ukidss_sourceid &  &  &  + 10 & jmag & mag & @xmath10 &  + 11 & e_jmag & mag & @xmath10 & observational uncertainty + 12 & r_jmag &  & @xmath10 & data source + 13 & hmag & mag & @xmath11 &  + 14 & e_hmag & mag & @xmath11 & observational uncertainty + 15 & r_hmag &  & @xmath11 & data source + 16 & kmag & mag & @xmath31 &  + 17 & e_kmag & mag & @xmath31 & observational uncertainty + 18 & r_kmag &  & @xmath31 & data source + 19 & 3.6mag & mag & @xmath32 m &  + 20 & e_3.6mag & mag & @xmath32 m & observational uncertainty + 21 & r_3.6mag &  & @xmath32 m & data source + 22 & 4.5mag & mag & @xmath33 m &  + 23 & e_4.5mag & mag & @xmath33 m & observational uncertainty + 24 & r_4.5mag &  & @xmath33 m & data source + 25 & 5.8mag & mag & @xmath34 m &  + 26 & e_5.8mag & mag & @xmath34 m & observational uncertainty + 27 & r_5.8mag &  & @xmath34 m & data source + 28 & 8.0mag & mag & @xmath35 m &  + 29 & e_8.0mag & mag & @xmath35 m & observational uncertainty + 30 & r_8.0mag &  & @xmath35 m & data source + 31 & 24mag & mag & @xmath36 m &  + 32 & e_24mag & mag & @xmath36 m & observational uncertainty + 33 & r_24mag &  & @xmath36 m & data source + 34 & sedclass &  &  & ir class according to sed slope + 35 & s1_sedclass &  &  & ir class according to sed slope ( visibility window 1 ) + 36 & member(ysovar ) &  &  & cluster membership according to ysovar standard + 37 & standardset &  &  & source in ysovar standard set ? + 38 & ns1_36 & ct & @xmath32 m & number of datapoints ( visibility window 1 ) + 39 & ns1_45 & ct & @xmath33 m & number of datapoints ( visibility window 1 ) + 40 & maxs1_36 & mag & @xmath32 m & maximum magnitude in lightcurve ( visibility window 1 ) + 41 & mins1_36 & mag & @xmath32 m & minimum magnitude in lightcurve ( visibility window 1 ) + 42 & maxs1_45 & mag & @xmath33 m & maximum magnitude in lightcurve ( visibility window 1 ) + 43 & mins1_45 & mag & @xmath33 m & minimum magnitude in lightcurve ( visibility window 1 ) + 44 & means1_36 & mag & @xmath32 m & mean magnitude ( visibility window 1 ) + 45 & stddevs1_36 & mag & @xmath32 m & standard deviation calculated from non - biased variance ( visibility window 1 ) + 46 & deltas1_36 & mag & @xmath32 m & width of distribution from 10% to 90% ( visibility window 1 ) + 47 & means1_45 & mag & @xmath33 m & mean magnitude ( visibility window 1 ) + 48 & stddevs1_45 & mag & @xmath33 m & standard deviation calculated from non - biased variance ( visibility window 1 ) + 49 & deltas1_45 & mag & @xmath33 m & width of distribution from 10% to 90% ( visibility window 1 ) + 50 & redchi2tomeans1_36 &  & @xmath32 m & reduced @xmath22 to mean ( visibility window 1 ) + 51 & redchi2tomeans1_45 &  & @xmath33 m & reduced @xmath22 to mean ( visibility window 1 ) + 52 & coherence_time_36 & d & @xmath32 m & decay time of acf ( visibility window 1 ) + 53 & coherence_time_45 & d & @xmath33 m & decay time of acf ( visibility window 1 ) + 54 & s1_stetson_36_45 &  & @xmath32 m , @xmath33 m & stetson index for a two - band lightcurve .
( visibility window 1 ) + 55 & s1_cmd_alpha_36_45 & rad & @xmath32 m , @xmath33 m & angle of best - fit line in cmd ( visibility window 1 ) + 56 & s1_cmd_alpha_error_36_45 & rad & @xmath32 m , @xmath33 m & uncertainty on angle ( visibility window 1 ) + 57 & n_36 & ct & @xmath32 m & number of datapoints + 58 & n_45 & ct & @xmath33 m & number of datapoints + 59 & max_36 & mag & @xmath32 m & maximum magnitude in lightcurve + 60 & min_36 & mag & @xmath32 m & minimum magnitude in lightcurve + 61 & max_45 & mag & @xmath33 m & maximum magnitude in lightcurve + 62 & min_45 & mag & @xmath33 m & minimum magnitude in lightcurve + 63 & mean_36 & mag & @xmath32 m & mean magnitude + 64 & stddev_36 & mag & @xmath32 m & standard deviation calculated from non - biased variance + 65 & delta_36 & mag & @xmath32 m & width of distribution from 10% to 90% + 66 & mean_45 & mag & @xmath33 m & mean magnitude + 67 & stddev_45 & mag & @xmath33 m & standard deviation calculated from non - biased variance + 68 & delta_45 & mag & @xmath33 m & width of distribution from 10% to 90% + 69 & redchi2tomean_36 &  & @xmath32 m & reduced @xmath22 to mean + 70 & redchi2tomean_45 &  & @xmath33 m & reduced @xmath22 to mean + 71 & stetson_36_45 &  & @xmath32 m , @xmath33 m & stetson index for a two - band lightcurve .
+ 72 & cmd_alpha_36_45 & rad & @xmath32 m , @xmath33 m & angle of best - fit line in cmd + 73 & cmd_alpha_error_36_45 & rad & @xmath32 m , @xmath33 m & uncertainty on angle + 74 & teff & k &  & effective temperature from literature + 75 & r_teff &  &  & reference for teff + 76 & x - ray &  &  & chandra counterpart ?
+    table  [ tab : tab2 ] contains the position , the designation , the flux densities of each source and properties of their lightcurves .
the properties of the lightcurve will be discussed in detail in the remainder of this article .
most properties of the lightcurve , e.g.  mean , minimum and maximum , appear twice .
they are calculated once over the entire available lightcurve and once for the first visibility window only ; the fast - cadence sampling is available uniformly for all clusters in the ysovar project and thus values calculated over the fast - cadence only can be compared between clusters ( see , e.g.  paper  i ) .
a subset of the properties of lightcurves with a mean magnitude @xmath37 in irac  1 or irac  2 is shown in table  [ tab : tab2small ] .
cccccccc cfhtwir - oph 29 & f & yes & 0.57 & 0.55 & 368.58 & 18.61 & 3.20 + [ edj2009 ] 809 & ii & yes & 0.17 & 0.18 & 41.13 &  &  + wl 6 & i & yes & 0.58 &  & 440.74 &  & 4.70 + cfhtwir - oph 16 & ii & yes & 0.05 & 0.05 & 2.42 &  &  + iso - oph 138 & ii & yes & 0.11 & 0.12 & 15.90 &  & 1.00 + iso - oph 53 & ii & yes & 0.06 & 0.08 & 5.56 & 1.72 & 1.10 + wsb 52 & ii & yes & 0.20 & 0.24 & 69.45 & 7.81 & 7.80 + wl 4 & ii & yes & 0.46 & 0.38 & 366.61 & 17.97 & 5.80 + iso - oph 137 & i & yes & 0.20 & 0.17 & 65.12 & 5.60 & 1.40 + wl 3 & i & yes & 0.23 & 0.25 & 61.16 & 8.67 & 8.50 + iso - oph 139 & f & yes & 0.12 & 0.14 & 17.48 & 4.10 & 1.60 + iso - oph 51 & f & no & 0.40 & 0.30 & 239.84 &  & 4.20 + iso - oph 122 & f & yes & 0.19 &  & 45.78 &  & 1.10 + wsb 49 & ii & no & 0.23 & 0.18 & 84.90 &  &  + iso - oph 161 & i & yes & 0.34 & 0.34 & 241.40 & 13.16 & 5.40 + rox 25 & ii & yes &  & 0.20 &  &  &  + iso - oph 140 & ii & yes & 0.17 & 0.19 & 42.95 & 5.99 & 4.70 + iso - oph 120 & f & yes & 0.18 & 0.24 & 69.78 & 6.44 & 2.20 + sstysv j162636.08 - 242404.2 & i & no & 0.09 & 0.07 & 6.14 & 1.28 & 0.90 + iso - oph 152 & ii & yes & 0.08 & 0.07 & 5.45 & 0.73 & 1.00 + iso - oph 21 & i & yes & 0.52 & 0.35 & 329.65 & 13.48 & 7.70 + roxn 44 & ii & no & 0.06 & 0.05 & 5.13 & 1.86 & 3.20 + ylw 15 & i & yes & 0.13 &  & 34.38 &  & 7.20 + [ gy92 ] 30 & i & yes & 0.13 & 0.10 & 12.90 & 2.03 & 3.60 + sstysv j162721.82 - 241842.4 & ii & no & 0.03 & 0.06 & 0.76 &  & 1.30 + wl 11 & ii & yes & 0.21 & 0.25 & 28.38 &  &  + ylw 47 & ii & yes & 0.16 & 0.13 & 49.72 &  & 3.10 + iso - oph 35 & ii & yes & 0.06 & 0.06 & 3.68 &  &  + [ gy92 ] 264 & ii & yes & 0.20 & 0.26 & 62.16 & 7.17 & 2.10 + 2mass j16271881 - 2448523 & iii & no & 0.10 & 0.09 & 4.79 &  &  + iso - oph 153 & ii & no & 0.29 & 0.23 & 118.17 & 8.42 & 0.70 + sstysv j162622.19 - 242352.2 & iii & no & 0.02 & 0.02 & 0.68 & 0.32 & 0.70 + cfhtwir - oph 74 & ii & no & 0.21 & 0.38 & 3.56 & 1.08 & 5.00 + iso - oph 34 & f & yes & 0.05 & 0.05 & 2.26 & 1.03 & 0.50 + crbr 2322.3 - 1143 & ii & yes & 0.07 & 0.12 & 4.53 &  & 2.00 + iso - oph 33 & f & yes & 0.10 & 0.12 & 6.25 & 1.90 & 4.70 + iso - oph 145 & f & yes & 0.36 & 0.41 & 216.30 & 13.03 & 1.80 + roxn 41 & ii & no & 0.05 & 0.03 & 3.02 & 1.04 & 1.30 + iso - oph 144 & f & yes & 0.13 & 0.08 & 21.70 & 3.32 & 1.50 + sstysv j162617.46 - 242314.3 & ii & no & 0.35 & 0.11 & 21.64 & 6.37 &  + iso - oph 50 & i & yes & 0.99 & 1.90 & 821.89 &  & 5.60 + [ gmm2009 ] oph l1688 30 & i & yes & 0.24 & 0.20 & 42.64 & 5.27 & 3.70 + iso - oph 165 & i & yes & 0.24 & 0.26 & 73.32 & 6.23 & 6.90 + [ edj2009 ] 892 & f & yes & 0.95 & 1.13 & 1084.07 & 33.56 & 7.60 + iso - oph 26 & f & yes & 0.11 & 0.12 & 16.35 & 2.65 & 3.00 + iso - oph 154 & ii & yes & 0.15 & 0.24 & 22.39 & 4.74 & 1.30 + cfhtwir - oph 21 & f & no & 0.06 & 0.05 & 5.10 &  & 3.00 + iso - oph 124 & i & yes & 0.17 & 0.19 & 33.31 & 4.76 & 1.20 + wl 13 & ii & yes & 0.11 &  & 14.34 &  & 1.70 + rox 26 & i & yes & 0.29 &  & 120.65 &  & 8.40 + iso - oph 37 & i & yes & 0.46 & 0.36 & 228.62 & 11.54 & 6.40 + [ edj2009 ] 824 & f & yes & 0.08 & 0.07 & 5.94 & 2.19 & 3.70 + iso - oph 118 & f & yes & 0.27 & 0.29 & 52.56 & 7.06 &  + sstysv j162727.53 - 242611.2 & i & no & 0.31 & 0.27 & 5.03 & 1.46 & 1.10 + iso - oph 52 & f & yes & 0.30 & 0.27 & 110.30 & 9.03 & 4.90 + sstysv j162621.66 - 241820.1 & f & no & 0.23 & 0.28 & 7.15 &  & 0.70 + iso - oph 19 & ii & yes & 0.09 & 0.12 & 13.24 & 3.00 &  + 2mass j16263046 - 2422571 & f & yes & 0.40 & 0.44 & 185.94 & 12.69 & 4.20 + sstysv j162728.30 - 244029.5 & i & no & 0.52 & 0.54 & 4.72 & 1.97 & 6.20 +      we build two l1688 membership lists based on different criteria . the first is defined in paper  i and is applied uniformly for all clusters in the ysovar project .
sources are treated as cluster members if they fulfill at least one of the following criteria : ( i ) they are classified as ysos by @xcite based on their ir excess in cryogenic mission _
spitzer _ data or ( ii ) they are detected as x - ray sources in _ chandra _ imaging and have a spectral slope compatible with a stellar photosphere ( sed class  iii , see section  [ sect : spectralslope ] ) . at the distance of l1688
( we use 120  pc from , but see also discussion in paper  i ) , cluster members that are young and thus still magnetically active stars can be detected in x - rays .
a total of 57 sources fulfill one or both conditions . in table
[ tab : tab2 ] , these objects are marked as `` member ( ysovar ) '' .
the main biases in this sample are that the ir criterion selects only those members with disks and not class  iii sources , while the _ chandra _ criterion suffers from incomplete spatial coverage and it may include late - type foreground stars .
the different biases are a common problem in multiwavelength studies of star forming regions ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the second membership list is taken from @xcite . on this list , 74 of our 884 sources with _
lightcurves are cluster members ; 51 of those 74 are also included in the ysovar membership list due to their x - ray or ir emission .
when we compare properties of members and non - members below without referring to a specific set , then we mean membership as defined by the standard ysovar criteria .
six sources in the ysovar standard member set are not part of the @xcite list .
more details on those source are given in appendix  [ sect : g09notwil ] .
on the other hand , all of the sources in @xcite that are not part of the ysovar standard set were selected based on an x - ray detection from observations other than _ chandra _ ( for the specific references see * ? ? ?
figure  [ fig : hist36 ] shows a histogram of the observed mean magnitudes in @xmath32 m .
almost all bright sources are cluster members .
the fraction of members drops below 11th  mag and no cluster member is found below 14th  mag .
we expect background sources to be fainter than cluster members because of their larger distance and because they are seen through the cloud , but also fainter sources can not be reliably classified .
the thickness of the cloud is inhomogenous .
the three primary fields cover roughly @xmath38  mag in the extinction maps from the c2d project @xcite .
this correspondes to @xmath39}=1.3 - 2.2 $ ]  mag @xcite . at the distance of l1688 , a k5 star with little or no extinction
has @xmath40 } = 12.2 $ ]  mag according to the evolutionary tracks of and the yso colors from @xcite .
thus , even late - k or early - m star members in the cloud are bright enough to be detected .
just over 10% of the 884 sources with lightcurves are classified by @xcite at all , because this classification scheme is conservative , representing a set of sources that can be dereddened and classified to a high degree of confidence .
however , for many sources , we lack the required spectral coverage . of those sources that can be classified in this way , 16
are class  i sources , 36 are class  ii sources and 59 are class  iii candidates , i.e. , sources with a weak or absent ir excess in their sed ( in this case the ir sed does not provide sufficient information to decide if a star is a yso or a field star ) .
ten of those 59 have an x - ray counterpart in our _ chandra _ data .    in order to classify all sources , including those not classified in @xcite , paper  i defines a simpler approach , which uses the observed colors only . in this scheme
, we fit the spectral slope @xmath41 , where @xmath42 is the wavelength and @xmath43 the flux density per unit wavelength interval at that wavelength .
we make use of all measured flux densities ( no upper limits ) in the 2 - 24@xmath12 m range , which corresponds to the range from the @xmath4 or @xmath31 filter , which we take from 2mass and/or ukidss , respectively , to the mips 24@xmath12 m channel .
we use the flux densities for all _ spitzer _ bands from the observations presented in @xcite ( and we include all sources extracted using this pipeline including those that are not published in @xcite because they are not classified as yso in that work ) or c2d and summarize the new lightcurves by calculating the mean and standard deviation for the new irac1 and irac2 data .
thus , sources detected both in the cryogenic mission and in the new dataset will have two datapoints for irac1 and irac2 .
if a source is detected in both 2mass and ukidss we use both values for the fit , since they are independent measurements .
we perform a least - squares fit and call sources with @xmath44 class  i , @xmath45 flat - spectrum ( f ) , @xmath46 class  ii sources and @xmath47 class  iii candidates .
the classifications are given in table  [ tab : tab2 ] in two columns . `
s1_sedclass ` presents the derived sed class using all available literature data as described above and the mean of the lightcurves from visibility window  1 for comparison with paper  i ; ` sedclass ` uses the same literature data but the mean flux density for the lightcurves calculated for all visibility windows . in this scheme
a significantly reddened main - sequence star , which has an intrinsic slope @xmath48 , may appear as class  ii object , so this observational classification can not be translated directly into the evolutionary stage of an object . of the remaining sources we find 110 class  i , 78 flat - spectrum and 455 class  ii sources and 221 class  iii candidates .
we can not classify 20 sources , because they are seen in one band only .    not only the evolutionary status but also other properties of the individual source such as inclination or stellar mass influence the value of @xmath23 .
massive stars emit more energetic radiation and can thus change the structure of the accretion disk .
however , l1688 is a region without massive ysos @xcite and we define @xmath23 using wavelengths longward of the @xmath31 band far from the peak of the stellar sed , so the shape of the photospheric sed has only negligible influence on the total sed .
@xcite simulated different ysos of low mass .
their results show that seds depend mostly on the evolutionary stage except for very extreme cases such as stars with an edge - on disk .
thus , @xmath23 provides a good proxy for the evolutionary state of the ysos in l1688 .
comparing with @xcite , the resulting classification is very similar , particularly for the class  i and flat - spectrum sources .
most sources with a spectral slope @xmath49 according to our slope - fitting are class  i sources in both classification schemes . with two exceptions
our flat - spectrum sources are either class  i or ii in @xcite . of the 52 sources with a spectral slope between -1.6 and -0.3 , 21
are also called class  ii by @xcite , but 30 are class  iii in that paper , indicating that a significant fraction of what we call class  ii might indeed be reddened background stars . the general agreement between the more complex classification scheme and the observed spectral slope is also found for other star forming regions ( paper  i ) .
we use three different methods to detect variability in all lightcurves , independent of their membership status or sed slope .
sources are considered variable , if they fullfill at least one of the following conditions : ( i ) a two - band lightcurve exisits and their stetson index is larger than 0.9 ( sect .  [ sect : stetsonindex ] ) ; ( ii ) only a one - band lightcurve exisits and @xmath50 ( sect .
[ sect : chi2 ] ) ; ( iii ) the lightcurve is periodic ( sect .
[ sect : periodicity ] ) .
sources in the target fields are observed in irac  1 and irac  2 almost simultaneously ( within a few minutes ) . for those sources ,
we calculate the stetson index @xmath51 with points weighted evenly @xcite : @xmath52 where the sum is taken over all @xmath53 pairs of observations in irac1 and irac2 with observed magnitudes @xmath54 and @xmath55 and uncertainties @xmath56 and @xmath57 ; @xmath58 denotes the sign of @xmath59 .
@xmath59 is the product of the normalized residuals in both bands : @xmath60 here @xmath61 denotes the error - weighted mean of all irac  1 magnitudes and @xmath62 the error - weighted mean of all irac  2 magnitudes .
the stetson index is very robust to observational errors since those are unlikely to affect both bands in the same way .
following paper  i , we define a source as variable if @xmath63 . in paper
i , this is calculated over the first visibility window only ; in this article we use the lightcurve from all four visibility windows . the stetson index calculated for the first visibility window only can be found in table  [ tab : tab2 ] in column ` s1_stetson_36_45 ` and for the entire lightcurve in column ` stetson_36_45 ` .
figure [ fig : stetsonlcs ] shows examples of lightcurves that are classified as variable according to the stetson index . in the first visibility window , we find 34 sources to be variable , that is 18% of all lightcurves where we can calculate the stetson index .
the number is similar ( 38 sources , 18% ) if we consider all observing visibility windows .
this is not surprising , given that the stetson index is designed to be a `` robust '' statistic that is little influenced by a small number of datapoints , and compared to the first visibility window , the other visibility windows only contribute a small number of datapoints .
we caution that equation  [ eqn : stetson ] will introduce a systematic bias in such a way that a source that is variable in only one band will not be recognized as variable .
[ sect : chi2 ]    for sources outside the primary target fields or for sources which are too bright or too faint in one irac channel , the stetson index can not be calculated . for those sources we rely on the reduced chi - squared value to detect variability : @xmath64 instrumental uncertainties lead to a non - gaussian error distribution ,
so we use a conservative cut - off and mark sources as variable only if @xmath65 ( paper  i ) . using this metric
, we find 22 sources that exhibit variability in irac1 and 18 that exhibit variability in irac2 .
for comparison , we note that of the 38 sources classified as variable according to their stetson index in the last section , 34 are also variable according to the @xmath22 test for their irac1 lightcurve and 37 according to their irac2 lightcurve .
the @xmath66 calculated for the first visibility window only can be found in table  [ tab : tab2 ] in column ` redchi2tomeans1_36 ` and ` redchi2tomeans1_45 ` for irac  i and irac  2 , respectively , and for the entire lightcurve in columns ` redchi2tomean_36 ` and ` redchi2tomean_45 ` , again for irac  1 and irac  2 , respectively .
examples are shown in figure  [ fig : chi2lc ] . if we fit a linear slope instead of comparing to the mean ( i.e.  a constant ) , the results are the same .
many sources have significantly larger values of @xmath66 for the whole observation series than within one visibility window , indicating a constant luminosity over 40  days , but a change over timescales between six months and two years .
ccccc cfhtwir - oph 29 & [ 3.6]-[4.5 ] & 2.6 & 0.03 & f + iso - oph 138 & [ 3.6 ] & 3.3 & 0.00 & ii + sstysv j162636.08 - 242404.2 & [ 3.6 ] & 4.1 & 0.00 & i + iso - oph 152 & [ 3.6 ] & 4.7 & 0.03 & ii + [ gy92 ] 30 & [ 3.6 ] & 14.5 & 0.00 & i + wl 11 & [ 4.5 ] & 3.0 & 0.01 & ii + 2mass j16271881 - 2448523 & [ 4.5 ] & 6.4 & 0.00 & iii + sstysv j162622.19 - 242352.2 & [ 4.5 ] & 6.1 & 0.00 & iii + iso - oph 34 & [ 3.6 ] & 2.2 & 0.00 & f + iso - oph 33 & [ 4.5 ] & 2.4 & 0.02 & f + roxn 41 & [ 3.6 ] & 6.5 & 0.00 & ii + iso - oph 154 & [ 3.6 ] & 5.6 & 0.00 & ii + iso - oph 124 & [ 3.6 ] & 3.4 & 0.00 & i + wl 13 & [ 3.6 ] & 10.7 & 0.00 & ii +    lastly , we search for periodicity in the fast - cadence ( visibility window 1 ) lightcurves using a lomb - scargle periodogram ; again , see paper  i for details . we require a false alarm probability @xmath67 and to further reduce the number of false positives , we additionally run a period detection of all lightcurves on the nasa exoplanet archive periodogram service .
this service employs several algorithms because each algorithm is particularly suited to a different signal shape ( see paper  i or the website for details on the other algorithms ) .
we also calculate the autocorrelation function of each signal .
we find that the lomb - scargle periodogram is advantageous for lightcurves with @xmath68 points as is the case for most of our data .
we require that at least one of three supplementary algorithms that we run retrieves a periodicity with a similar timescale ( see paper  i for details ) .
we only search for periods between 0.1 and 14.5 days , so that at least three periods fit in each visibility window .
allowing longer periods leads to the detection of many long - term trends , where the data do not show that these trends are actually periodic .
finally , all algorithmically detected periodicity is vetted by eye .
if periods are found in multiple bands , we report the period in irac  1 , which is generally the most reliable due to the lower measurement uncertainties .
if irac  1 does not reveal a periodicity , we report the value for the irac  2 lightcurve , and , as a last resort the period found in the [ 3.6]-[4.5 ] color to include objects where a periodic signal is overlayed by a long - term trend so that it is undetectable in each individual channel , but might be visible in the color .
the final list of adopted periods is given in table  [ tab : period ] .
it contains sources from all sed classes .
we find 14 sources that show a periodic behavior in the first visibility window , but this periodicity is not stable over more than one visibility window .
the datapoints taken in the first visibility window follow the folded period much better than the data from the later visibility windows , indicating that the period is not stable for longer than a few months . in roxn  41 ,
the datapoints of the later visibility windows , folded with the same period , seem to follow a different , yet clearly defined lightcurve with a similar period .
ten periodic sources have information in both irac bands .
eight of those are already classified as variable by the stetson index .
three out of four sources with data in one band only are variable according to the @xmath66 test .
the remaining two sources with information in both irac bands show periodicty , but fail the stetson index test because a larger variability amplitude is required for a significant detection in stetson index , which does not make any assumption about the form of the variability compared with the lomb - scargle periodigram , that only detects periodic signals .
equally , the remaining source with data in only one band fails the @xmath66 test , because the @xmath66 test also requires a larger variability amplitude for a significant detection than the lomb - scarge periodogram .
the number of datapoints in visibility windows 2 , 3 , and 4 is too low to search for periodicity in those visibility windows alone .
figure  [ fig : periodlc ] shows the phase - folded lightcurve in one band for six periodic sources where the periodicity is significant in a single band .
figure  [ fig : periodcolorlc ] shows the phased lightcurve for , where the period is seen only in the color term [ 3.6]-[4.5 ] .
most of the periods found in table  [ tab : period ] are in the range 3 - 7 days .
only two sources have longer periods of 11 and 14 days .
the largest amplitude is around 0.3  mag , and the smallest around 0.05  mag .
all periodic sources except roxn  41 , sstysv j162636.08 - 242404.2 , 2mass j16271881 - 2448523 , and sstysv j162622.19 - 242352.2 are cluster members according to our membership criteria .
roxn  41 and sstysv j162636.08 - 242404.2 have class i and ii seds , respectively .
all but cfhtwir - oph  29 ( see appendix  [ sect : g09notwil ] ) , sstysv j162636.08 - 242404.2 ( see appendix  [ sect : g09notwil ] ) , and 2mass j16271881 - 2448523 ( no information beyond the @xmath69 magnitudes is available in the literature ) are cluster members according to @xcite .
ten of the 56 cluster members are periodic ; in total there are 14 periodic sources out of 60 sources that are variable in the first visibility window . in comparison
, @xcite find 18% of all variable stars to be periodic in @xmath70 monitoring of the orion  a molecular cloud with a similar time coverage as we have for l1688 .
@xcite find a third of all variable stars in l1688 to be periodic ; due to their longer time baseline , @xcite are sensitive to different periods , but only 4 out of their 32 periods have values that are outside of the range we could detect
see their table .
we are sensitive to periods of up to 14.5  days here .
the largest group of lightcurves in the sample of l1688 does not have any significant variability .
the lower limit where variability is detected depends on the brightness of the object - fainter sources need a stronger relative variability due to the larger measurement uncertainties .
paper  i presents monte - carlo simulations to show that the stetson test finds variability when the amplitude is a few times larger than the average uncertainty ; the exact number depends on the signal shape .
for example , the variability in a source that switches between two discrete levels will be detected when the amplitude is at least twice the uncertainty ( @xmath71  mag for a star with magnitude 14 ) .
if data from only one band are available , the step size must be more than four times the uncertainty to be found by the @xmath22-test ( @xmath72  mag for a star with magnitude 14 . ) .
given the observing cadence , the monte - carlo simulations show that we are sensitive to periods between 0.1 and 14.5 days ( paper  i ) .
in summary , we call a source variable if it is either periodic with a low false alarm probability ( 14 sources ) or fulfills one of the following conditions : sources with simultaneous data in two bands need to have a stetson index @xmath73 ( 34 sources in the ysovar standard set fast cadence data and 38 sources in total ) and sources without simultaneous data need to have @xmath65 ( 22 sources in the ysovar standard set fast cadence data and 29 sources in total ) .
we consider the lightcurves for all 882 distinct sources with at least five datapoints in at least one irac band . of those lightcurves ,
70 are classified as variable ; 56 sources are cluster members according to the ysovar criteria ( 44 of them are variable ) and 73 sources are cluster members according to @xcite ( 47 of them are variable ) .
both membership samples have considerable overlap .    in the following subsections ,
we compare properties of the lightcurves between different sed classes . in most cases , the lightcurves of irac1 ( 3.6@xmath12 m ) and irac2 ( @xmath6 m ) have very similar properties .
figure  [ fig : probvarsed ] contains the variability fraction sorted by sed class ( see section  [ sect : spectralslope ] ) . within each class ,
the population is sub - divided into previously identified members , bright ( @xmath74  mag in @xmath32 m or @xmath33 m ) stars not previously identified as members , and faint ( @xmath75  mag ) stars not previously identified as members .
sources are considered variable if they have a low false alarm probability for periodicity or a stetson index @xmath73 ( if simultaneous data in irac  1 and irac  2 exist ) or @xmath76 ( if no simultaneous data exist ) .
note that the classes shown here are based on the observed sed and thus background non - member stars seen through the cloud might appear red ( like a class ii source ) even if their intrinsic sed slope is compatible with sed class  iii .
no star classified as a member is fainter than 14  mag
.    the observed fraction of variable stars as a function of sed class is the best guess for the probability to find variability for a star of a given sed class .
we calculate the uncertainty for the probability based on the observed number of variable stars @xmath77 in a bin with @xmath53 total sources as follows : if the true probability for a source in the bin to be variable is @xmath78 , then the probability to observe @xmath77 variables in @xmath53 stars is given by @xmath79 , which follows a binomial distribution : @xmath80 the peak of this distribution is @xmath81 .
we then calculate the boundaries of the confidence interval @xmath82 for @xmath78 such that @xmath83 and the area under the curve includes 68% of the probability .
those uncertainty ranges are shown in figure  [ fig : probvarsed ] .
almost all known cluster members are variable ; the fraction of variables is decreasing slightly from class  i to ii .
there are five class  iii members and none of them shows variability .
this is compatible with the variability fraction observed for class  iii non - members . by definition
, sources with a class  iii sed do not have a large ir excess over a stellar photosphere , thus they are not expected to harbor a substantial disk or show disk - related variability .
however , even field stars can show significant variability , at least in the near - ir .
@xcite find 1.6% of all field stars to be variable in @xmath69 .
this includes eclipsing systems , stars with unusually large photospheric spots and other unidentified variables .
all effects that influence the @xmath69 lightcurve are likely also visible in the mid - ir .
our fraction of variable objects with a class  iii sed is similar ( 1.8% ) to what wolk et al .  found . in the optical and near - ir ,
the fraction of variable field stars is lower than that value .    _
spitzer _ observations of other young clusters also find that younger stars are more variable , but with a lower variability fraction .
@xcite find only 70% of all stars with disks in their onc sample to have detectable variability .
@xcite find 60% in ic  348 and @xcite find @xmath84% of all members in ngc  2264 to be variable in the irac bands .
all three clusters are located at larger distances and thus observations are less sensitive to small variations , the observations of @xcite are more densely sampled and thus provide a better signal . in l1688 , a k5 star with little extinction will have a magnitude of 12.2 in @xmath5 m using the evolutionary tracks of and the yso colors from @xcite . in this magnitude range ,
variability down to 0.05  mag can be detected ( paper  i ) .
in contrast , for the same star in ic  348 , the variability has to be about twice as large to be detected and three times as large to the detected in the onc .
of course , even our observations miss the faint end of the distribution .
thus , it seems likely that essentially all ysos show substantial mid - ir variability on timescales of days to weeks .
sources with a class  iii sed are not contaminated by disk - bearing stars , since reddening by the cloud can only increase the spectral slope @xmath23 but never hide an existing ir excess
. the larger fraction of non - member variable objects in the other classes in figure  [ fig : probvarsed ] thus shows that the membership lists are incomplete .
one caveat here is that the class  iii sources in the sample of members are selected in a different way than the other classes . if true class  iii cluster members with x - ray emission ( that are included in our member sample )
are systematically less variable than class  iii members where we do not detect x - ray emission ( that are therefore not included in the member sample ) , that would also lead to a lower observed variability fraction in class  iii member sources .
however , given the size of the observed effect in variability such a bias seems to be unlikely to be the sole reason for the observed distribution .    the observed probability that a source is variable for all classes of bright non - members ( blue bars ) is consistent with around 5%
the fact that this barely changes with the sed class , quite unlike the distribution for the cluster members , indicates that unidentified cluster members can not make up a large fraction of the non - member sample since they would bias the observed variability fraction to higher values for the earlier classes .
the vastly different number of sources in the different bins makes it highly unlikely that each bin is contaminated by the same fraction of class  i - ii cluster members .
the sample of faint stars has almost the same fraction of variable sources in every sed class .
the low fraction can be explained by two effects .
first , variability of faint sources can not be reliably detected unless the amplitude is exceptionally large .
second , as most cluster members are brighter than 14  mag , the faint sample contains fewer unrecognized cluster members than the sample of brighter sources .
figure  [ fig : meanampsed ] shows that there is a wide distribution of amplitudes in ysos of class i , f and ii .
the error bars in this figure represent the standard deviation of the mean amplitudes within one class .
amplitude and standard deviation in the figure are measured in magnitudes ; a median amplitude of 0.1  mag means that the flux in @xmath33 m varies by 10% for a typical source . for members ,
the variability amplitude is larger in class  i sources than in flat - spectrum and class  ii sources , where the amplitude is still larger than in class  iii sources .
also , the spread of the amplitudes in class i and f is much larger than for class ii and iii . for non - members ,
the observed amplitudes seem to follow a similar , but less pronounced trend .
this suggests that the non - member category still includes not only reddened background objects with an intrinsic class  iii sed , but also some unidentified cluster members .
indeed , we propose that all variable class  i , f and ii sources are members and table  [ tab : newmembers ] lists those new members that are neither contained in our ysovar standard membership set ( see table  [ tab : tab2 ] ) nor in @xcite .
however , for consistency , we do not modify our sample of members at this point and continue to treat the sources in the table as non - members for the remainder of the analysis .
two of the objects in the table , and were suggested as a @xmath7  oph substellar candidate members by , the remaining objects have not been studied before .    in summary ,
the probability for any given source to be observed as a mid - ir variable decreases little between class  i and ii sources , but the mean amplitude as well as the differences within a class are much larger for sources which have more circumstellar material .
cccc 246.88244 & -24.79978 & sstysv j162731.78 - 244759.2 & ii + 246.64060 & -24.31667 & sstysv j162633.74 - 241900.0 & ii + 246.84236 & -24.78363 & sstysv j162722.16 - 244701.0 & ii + 246.59001 & -24.49297 & sstysv j162621.60 - 242934.7 & ii + 246.65036 & -24.40117 & sstysv j162636.08 - 242404.2 & i + 246.84092 & -24.31180 & sstysv j162721.82 - 241842.4 & ii + 246.84135 & -24.78713 & sstysv j162721.92 - 244713.6 & ii + 246.84336 & -24.64383 & cfhtwir - oph 74 & ii + 246.90225 & -24.64854 & sstysv j162736.53 - 243854.7 & ii + 246.57278 & -24.38732 & sstysv j162617.46 - 242314.3 & ii + 246.67446 & -24.25616 & sstysv j162641.86 - 241522.1 & ii + 246.91134 & -24.55516 & sstysv j162738.72 - 243318.5 & ii + 246.60121 & -24.26383 & cfhtwir - oph 21 & f + 246.88281 & -24.34855 & sstysv j162731.87 - 242054.7 & f + 246.86474 & -24.43647 & sstysv j162727.53 - 242611.2 & i + 246.59027 & -24.30560 & sstysv j162621.66 - 241820.1 & f + 246.86794 & -24.67487 & sstysv j162728.30 - 244029.5 & i +      next , we will look at the timescales of the lightcurves .
the list of periodic sources contains objects of all evolutionary stages , but is too small to recognize any trends in the period with evolutionary stage .
we calculate the discrete auto - correlation function @xmath85 for each fast - cadence ( visibility window 1 ) lightcurve to characterise the time - scale in all sources , not just the periodic ones .
one complication here is that our data are unevenly sampled . in order to calculate the @xmath86 we linearly interpolate the lightcurve on a grid with time steps of 0.1  days .
this process can change the properties of the lightcurves on short timescales .
however , for most lightcurves the relevant timescales are longer than the distance between two observations ( 0.1 to 0.8 days ) .
we use the following definition of the @xmath86 @xmath87 where @xmath53 is the total number of points in the lightcurve @xmath88 with a mean of @xmath61 and a standard deviation of @xmath89 . in the discrete @xmath86 @xmath90
is the number of timesteps .    by definition @xmath91 and
the @xmath86 decays for longer time - lags @xmath90 .
we take the first value of @xmath90 with @xmath92 as the characteristic time - scale for a lightcurve .
a more common definition would be to use the position of the first local maximum in the @xmath86 , but due to the low number of datapoints in our lightcurves , the noise in the @xmath86 is large and this value is often not well defined . while our definition might not give us the timescale of the physical processes at work , the @xmath86 timescale still provides relative comparisons .
we find that the average @xmath86 timescale decreases from 4  days for class i sources to one day for class  iii sources ( figure  [ fig : acf ] ) .
a more detailed analysis of stochastic and quasi - periodic properties of the lightcurves is beyond the scope of the current paper ( for a discussion of these properties in @xmath70 lightcurves , see * ? ? ?
the majority of the stars included in our data show no significant variability , and we exclude them from further discussion .
the remaining stars are definitely variable , with light curve shapes that display a variety of morphologies .
visual inspection , and some quantitative analysis , allows us to group these stars into sets with similar light curve morphologies - which can be the first step in attaching physical mechanisms as the cause of the observed variability .
based on a near - ir ( @xmath69 ) monitoring campaign of similar cadence and duration to ours , @xcite identified four light curve morphological classes : ( a ) periodic ; ( b ) quasi - periodic - which they defined as stars with cyclic brightening and fading , but where the frequency and amplitude of the variations varied from cycle to cycle ; ( c ) long - duration - stars with relatively long term monotonic changes in brightness over weeks or months , with eventual changes in sign of the variability ; ( d ) stochastic - which indicated all other variability types , where no obvious pattern to the variability was present .    in another recent paper ,
@xcite analysed a @xmath93 day monitoring campaign for the star - forming region ngc 2264 , using optical data from corot and ir data from irac , to assign light curve morphology classes to their young star set .
their proposed light curve morphologies were : ( i ) periodic - stars whose light curves show periodic waveforms whose amplitude and shape are unchanging or only change in minor ways over the 30 day observing window .
these light curves were ascribed generally to cold spots on the stellar photosphere ; ( ii ) dippers - stars showing a well - defined maximum brightness , upon which are superposed flux dips of variable shape and amplitude .
these are sub - divided into stars where the flux dips occur at an approximately constant period - designated as quasi - periodic systems - and stars with no obvious periodicity to the dips - designated aperiodic . in previous literature , these two sub - classes , or portions of them , have been referred to as aa tau systems and ux ori systems , respectively ; ( iii ) short duration bursters - stars with relatively well - defined minimum light curve levels , superposed on which are brief ( hours to day ) flux increases , attributed to accretion bursts .
see @xcite for further discussion of this set ; ( iv ) quasi - periodic variables - stars lacking a well - defined maximum brightness but showing periodic variability
whose waveform changes shape and amplitude from cycle to cycle ; ( v ) stochastic - stars with prominent luminosity changes on a variety of timescales , with no preference for `` up '' or `` down '' .
( vi ) long timescale variability - stars with slow ( weeks to months ) changes in brightness .    while there are some clear similarities in the two schemes , their usage of the terms quasi - periodic and stochastic are not the same , and light curves described as belonging to those classes in one scheme would not necessarily be so classified in the other scheme
. it will be important to resolve these nomenclature issues in the future , perhaps in the same manner as was done for sorting out similar issues on classifying pre - main sequence disk sed morphologies ( evans et al . 2009 ) . in the meantime
, we must choose which scheme to adopt and be clear that is what we have done .
for this paper , we adopt the @xcite scheme .
periodic lightcurves are described in section  [ sect : midirvariability ] and examples can be seen in figure  [ fig : periodlc ] .
examples of quasi - periodic lightcurves are shown in figure  [ fig : morelcs ] .
shows variability with a timescale @xmath94  days in the first half of the first visibility window , but around mjd 55316 , the flux increases significantly and this rise masks out any periodicity . the dip around mjd  55325
again has a similar duration as those observed in the beginning of the visibility window .
the other two sources shown might be similar , but the signal - to - noise is not as good .
other sources are variable , but no preferred timescale for the variability is discernible . and 51 ( figure  [ fig : morelcs2 ] ) are examples of this aperiodic behavior . in yet other cases , the timescale of the variability is so long that we can not decide if a feature is a singular event or part of a recurring pattern .
good examples of this are wl  4 and , which are shown in figure  [ fig : cmd ] .    in some cases ,
we see short duration features in the lightcurve in addition to a longer trend .
this can either be a brightening ( e.g.  mjd 55303 and 55311 in , figure  [ fig : morelcs2 ] ) or short dips in the lightcurve that last between one and ten days ( in the same figure ) .
like @xcite , we find that even strongly periodic sources are not perfect clocks in the mid - ir in that the amplitude can vary from cycle to cycle .
there might also be phase shifts from visibility window to visibility window . in a very similar study to the one we present here , @xcite observed ic 348 over a 40 day window with irac1 and 2 and find 25% of the variable stars are likely periodic .
here we find a similar division ; 10 of 52 members are strongly periodic .
we present the first mid - ir monitoring of l1688 , but the field has been monitored before in the near - ir .
figure  [ fig : amplitudekvs36 ] compares the @xmath4 band amplitude found in 2mass data taken between 1997 and 2001 @xcite and the @xmath31 band amplitude observed at ukirt between 2005 and 2006 with the amplitude of our @xmath32 m lightcurves .
given the different bands , it is not surprising that the absolute value of the amplitude differs , but there is a good correlation such that the sources with the largest @xmath31 or @xmath4 band amplitudes in earlier observations also have the largest amplitudes in the irac bands five to ten years later . in most objects ,
color changes are small or happen in parallel with luminosity changes ( section  [ sect : colorchangesandreddening ] ) , so the @xmath31 or @xmath4 band and the [ 3.6 ] variability should be strongly correlated .
figure  [ fig : amplitudekvs36 ] shows that the amplitude of the variability is relatively stable over at least one decade .
@xcite used a larger number of observations than and observe a larger spread between a typical yso s brightest and faintest @xmath31-band magnitude , indicating that the longest timescale of variability is longer than the time span of the observations .
@xcite noted that roughly 6% of the stars in ic  348 had data which indicated significant changes in the source s mid - ir flux over the three year interval since the last irac observation .
@xcite also displayed numerous examples of stars which showed continuous change over the course of at least one observing visibility window .
all this shows that the amplitude of variability of a given source may evolve over very long timescales , but is consistent over at least one decade .
ysos apparently do not switch between highly variable and much less variable states .       of time - variable extinction .
the largest @xmath95 observed is 1 , but for clearity only the region to @xmath96  mag is shown .
_ bottom _ : in this panel the length of each line gives the sed class of that object .
only cluster members have slopes where the source bluens as it fades .
[ fig : cmdslope],title="fig:",scaledwidth=40.0% ]   of time - variable extinction .
the largest @xmath95 observed is 1 , but for clearity only the region to @xmath96  mag is shown .
_ bottom _ : in this panel the length of each line gives the sed class of that object .
only cluster members have slopes where the source bluens as it fades .
[ fig : cmdslope],title="fig:",scaledwidth=40.0% ]    for sources in the primary target fields , the observations in @xmath5 m and @xmath6 m are separated by only a few minutes . for those sources ,
we compare the color and the magnitude in the @xmath5 m band in a color - magnitude diagram ( cmd ) ( figures  [ fig : morelcs ] and [ fig : cmd ] ) . in some cases
, the color is fairly stable within one visibility window , but changes between visibility windows ; in others the timescale of variability is much shorter and color changes are seen within one visibility window as well . for non - variable sources , the cmd forms a point cloud with the size set by the photometric uncertainties , but for variable sources , the shape of the cmd can reveal the physical cause of the variability  for example , if a disk warp or accretion funnel passes in front of a yso , we expect its color to become redder as it becomes fainter . if the absorber has the same gas and dust properties as the inter - stellar medium ( ism ) and absorbs star and inner disk at the same time , then the datapoints in the cmd follow a line with the slope of the interstellar reddening law ( e.g. * ? ? ?
* ) .    in each cmd ,
the observational data ( i1 , i1-i2 ) are fit to a line segment using an orthogonal distance regression method that takes the errors in both the x and y directions into account .
we calculate the length in magnitudes of the line segment excluding the 10% of the data that are outliers on either side .
the cmd slopes of the 37 sources with a well - defined slope are shown in figure  [ fig : cmdslope ] .
the x and the y axis in the cmd are correlated , because they both depend on the @xmath5 m magnitude so sources which are noise - dominated would show a slope of @xmath97 .
no source is seen in this region .    from a total of 42 variable sources with cmds with at least ten datapoints ,
37 can be fitted with a formal uncertainty on the best fit slope below @xmath98 . in four of the five sources with a statistical error on the slope in the cmd @xmath99 ,
the color variability is dominated by observational uncertainties ; the remaining source iso - oph  140 is discussed below .    as can be seen in figure  [ fig : cmd ] , even in those sources with a well defined slope in the cmd ,
the scatter around the best linear fit is larger than the measurement uncertainties ; individual visibility windows are systematically above or below the fitted line .
the slopes shown in figure  [ fig : cmdslope ] can be separated ( somewhat arbitrarily ) in two groups .
one group becomes redder when the sources are fainter with slopes comparable to the ism reddening .
this group contains almost all variables that are not classified as cluster members in section  [ sect : l1688membership ] .
most sources in this group have class  i or flat - spectrum seds ( bottom panel ) .
sources that bluen ( we use the term `` bluen '' as a verb to mean that a source comes bluer similar to the common expressions `` redden '' or `` reddening '' ) when they are fainter are mostly class  ii and flat spectrum sources . figure  [ fig : morelcs ] shows the lightcurves and cmds for two of the bluening sources , and .
the top panel in figure  [ fig : cmdslope ] quantifies the magnitude of the reddening or bluening observed in each source . for a source
where the slope in the cmd is compatible with ism reddening , this can be interpreted as the @xmath95 of the intervening material .
a range of values is observed , but in most sources it is @xmath49  mag .
no significant difference in vector length is seen between sources that bluen or redden when they dim .
we discuss the physical mechanisms that could cause these slopes in section  [ sect : scenariosforcolorchanges ] .    comparing the amplitudes of all sources in figure  [ fig : cmdslope ]
that either redden or bluen over the full time span of the observations , we do not find significant differences in value or distribution of the variability amplitudes .
similar to our results , find reddening and bluening slopes in the cmd in their @xmath100 monitoring , but compared to figure  [ fig : cmdslope ] they see a higher fraction of stars that bluen as they dim ( 33/49 ) .
only six sources have well defined slopes in both their study and ours . in the data of ,
all of them belong to the group that bluen with lower fluxes . in five cases ,
our data agree .
in contrast , roxn  44 reddens in our observations , so its behavior changed over the time interval of five years , or it is behaving differently in the near and mid - ir .      is the only source with a statistical error on the slope in the cmd @xmath99 where the color variability is not dominated by observational uncertainties , but where the slope varies substantially over time .
its lightcurve and cmd for the first visibility window are shown in figure  [ fig : cmdsquare ] .
it has a class  ii sed slope and is a low mass yso ( spectral type m1 ; * ? ? ?
there is no single slope in the cmd ; instead the source seems to switch between different modes .
initially , the luminosity drops sharply , while the color becomes slightly bluer . the behavior changes around mjd 55305 . for the next 15 days , the source becomes redder and dimmer , but the slope is less steep than the ism reddening law .
the evolution is not monotonic , but a brightening ( e.g.  around mjd 55313 ) corresponds to a bluer color and the source reddens as it dims . at the end of this period ,
the brightness increases again sharply and the slope in the cmd is comparable to the period mjd 55295 to 55305 : the source becomes noticeably redder with higher luminosity . in the last 10 days of the monitoring , the source again has a slope close to the ism reddening law .
its luminosity continues to increase and the color bluens until it has a similar color and luminosity as it did at the beginning of the monitoring .
this closes the circle in the cmd .
this is the only source in our sample where we observe multiple changes of the reddening slope over time . in this case ,
variability similar to ism reddening corresponds to slow changes in luminosity , while the other direction of the slope in the cmd is associated with faster luminosity changes .
now we discuss the lightcurves with the largest changes in magnitude in our sample .
wl  4 , shown in figure  [ fig : cmd ] , has an outburst event consisting of a brightness increase by about 0.6  mag and bluening at the same time . the rise and
fall take only a few days and the entire outburst lasts about a month .
[ edj2009 ]  892 ( shown in the same figure ) presents a slow rise by about 0.2  mag in @xmath33 m and smooth decay over two years where the source gradually bluens ( although more activity between visibility windows can not be excluded ) .
there is an indication of a slow and smooth increase in brightness again in the last observing visibility window , pointing to a recurring phenomenon .
these timescales are much longer than the dynamical time in the inner disk , so they are presumably driven by disk phenomena that originate at larger radii .
@xcite identified a few such bursting or fading events in their analysis of ysovar data for the orion nebula cluster ( onc ) and they have also been observed in the optical @xcite . the physical cause is not known yet .
@xcite find a few lightcurves that match theoretical predictions for short mass accretion events , but without simultaneous spectroscopy this is hard to prove .
the best example presented here , wl  4 , becomes bluer during the burst , but the opposite happens for wl  3 , which otherwise has a similar lightcurve .
thus , it is likely that different mechanisms cause these burst events .
the reddening sources in figure  [ fig : cmdslope ] have slopes that are roughly compatible with the slope of an interstellar reddening law .
small deviations can be explained by modifications of the dust - to - gas ratio or the grain size and composition .
this is commonly observed in individual young stars and star forming regions as a whole @xcite .
the most prominent example of a class  ii source where the extinction changes on timescales of days is . in this case ,
the absorption dip is periodic and caused by an inner disk warp partially occulting the star .
many examples of comparable lightcurves have been found in the onc by @xcite and in ngc  2264 by @xcite .
these authors call objects with absorption events in the lightcuve `` dippers '' .
we observe reddening on timescales of months to years in about half of all sources with a well - defined slope in the cmd .
an example of that is shown in the top right panel of figure  [ fig : cmd ] .
we do not know the inclination of the objects in our sample , but if it is close to edge - on , this might be caused by an asymmetric disk similar to aa  tau .
the required warp or local change in scale height that lets some gas and dust protrude above the average disk height , could be caused by a low - mass star , brown dwarf or planetary mass object orbiting in the disk ( e.g. * ? ? ?
* ; * ? ? ?
alternatively , vortices in the disk or dust traps can also change the local scale height .
changes in the optical brightness and color of a yso are often attributed to spots , either hot accretion spots or cold magnetic spots , rotating in and out of view or to time - variable extinction ( e.g. * ? ? ?
* ; * ? ? ?
if the reddening increases and , at the same time , the visible fraction of the hot spot decreases , then the resulting slope in the cmd would be intermediate between ism reddening and colorless .
since spots on the star can not explain sources that bluen as they dim , this observed phenomenon must be related to the disk .
several parameters of the inner disk could be time variable .
the most obvious one is the accretion rate .
however , the direct effect of increased accretion is a larger accretion spot , which would make the star bluer and brighter , not fainter .
also , @xcite find no relation between mid - ir lightcurves and spectroscopic accretion tracers in their limited sample .
this indicates that the relevant parameter can change without affecting the accretion rate . in the mid - ir
, we see the optically thick dust at the inner disk edge .
because the disk is heated over its entire vertical height at the inner edge while only the disk surface absorbs radiation at larger radii , it can form a puffed - up inner rim that reaches above the usual disk scale height and casts a shadow on the remaining disk , reducing the luminosity at longer wavelengths ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
if this rim grows , it intercepts more stellar light .
thus , the emission from the inner wall increases and the emission from the outer disk decreases , which results in bluer colors .
since the total emitting area also decreases , the source can become bluer and fainter .
@xcite also present radiative transfer models with hot spots and warped disks that reproduce features of the observed lightcurves .
alternatively , @xcite suggested a non - axisymmetric model . if the emission from the stellar photosphere is stronger at some longitudes due to a strong accretion spot or a spot from magnetic activity on the star , it irradiates the inner disk like a searchlight beam and causes the dust at this longitude to retreat .
the result is a disk where the inner rim is broken at one ( or several ) positions .
this would reduce the emission in @xmath5 m , but cause more energy to be emitted at a larger radius , where the dust is cooler on average and thus emits at longer wavelengths such as @xmath6 m . in this way
, the source can become redder as it brightens .
there is one group of ysos , the so - called ux ori ( uxor ) variables where several members have been observed to bluen when they fade in optical observations .
this can be explained by a larger fraction of scattered light that contributes to the observed sed in their faint state .
a similar scenario for the ysos discussed here requires , first , scattering by relatively large grains ( @xmath101  mm ) , because only for large grains radiation at @xmath32 m will be scattered more than radiation at @xmath33 m .
grains this large will be present in the disk , but not at all radii and all disk heights .
second , for this mechanism to work , a large fraction of the ir radiation must be scattered light and not intrinsic emission . in contrast
, many bluening sources are class  i to flat - spectrum sources , where the ir luminosity is comparable to the total stellar irradiation .
none of the simple models presented so far offers a convincing explanation of the observed lightcurves and cmds .
however , detailed simulations that take into account turbulent transport processes in the disk and different dust species with a complete chemical network have not yet been performed .
the real situation is likely to be much more complex than sketched above . since different dust species form and sublimate with different speeds and at different temperatures ,
time variable irradiation can actually cause a very complex mixture of species with different opacities , which will not react linearly to increased irradiation . such a complex network might cause a delayed and non - linear response of the disk , which masks the underlying relation between accretion and disk emission .
the characteristics of the time variability in the lightcurves could be related to the spectral type of the central star as hotter stars irradiate the disk with harder spectra .
we compared the 90% quantiles of the observed @xmath6 m magnitudes for those sources in our sample that have spectroscopically determined spectral types from @xcite .
we use the relation between spectral type and @xmath102 from @xcite .
this gives values for ten objects in our sample and we do not find a correlation in this set .
similarly , @xcite do not identify any trends in infarared variability versus effective temperature among few myr old ngc 2264 stars .
in apparent contrast , @xcite find that the variability in [ 3.6 ] and [ 4.5 ] increases with increasing @xmath102 , but they show that this trend probably does not reflect a change in disk properties , but is instead due to the lower relative contribution of the photosphere compared to the disk for hotter stars .      in this section ,
we analyze the subsample of sources with x - ray counterparts .
since the subsample is not selected for its ir properties , it provides a clean sample to calculate the variability fraction .
there are 31 sources with detected x - ray emission , of which 20 are variable in the ir .
again , we find high variability fractions in class  i ( 6/7 ) , f ( 3/4 ) and ii ( 11/15 ) , while none of the four class  iii sources is variable .
the sample size is small , but this is consistent with our analysis of the ir selected sample in section  [ sect : evolutionarytrendsofthevariability ] .
we searched for correlations between the parameters of the x - ray emission ( median energy , fitted temperature and absorbing column density ) and the amplitude of the mid - ir variability .
no such correlation is apparent within each sed class . for a given absorbing column density ,
soft x - ray emission is more strongly absorbed than hard x - ray emission .
thus , the observed x - ray spectrum from more embedded sources always appears harder .
since class  ii sources have more circumstellar matter than class  iii sources , it is not surprising that they appear harder . @xcite already discuss this for the l1688 x - ray data ; figure  [ fig : probvarsed ] and [ fig : meanampsed ] show that class  ii sources are more variable than class  iii sources .
we see no additional effect of the x - ray emission on the variability characteristics in the mid - ir .
we present _ spitzer _ observations of ysos in the star forming region l1688 .
observations were taken in four visibility windows in spring and fall of 2010 and 2011 , with about 70 observations in spring 2010 and about ten observations in the remaining visibility windows .
the cadence of the observations is non - uniform to avoid bias in the period detection .
our sample consists of 882 sources with lightcurves in irac1 and irac2 with at least 5 datapoints . of those 882 sources , we classify 70 sources as variable using the stetson test , the @xmath22 test and the lomb - scagle periodogram .
the faintest sources in the sample have @xmath103  mag , but naturally the measurement uncertainties are larger for fainter sources .
the algorithms detect variability if the amplitude is larger than @xmath104  mag for a source of @xmath105  mag .
we define a sample of cluster members , including sources with an ir excess due to a circumstellar disk or with x - ray emission . for cluster members
, there is a clear correlation between evolutionary status and ir variability .
more embedded sources are more often detected to be variable , and they have on average larger variability amplitudes .
overall , the data are consistent with the idea that all ysos are variable at @xmath5 m and @xmath6 m , and we thus propose that all variable sources in our sample are members of the l1688 star forming region .    qualitatively different morphological types of lightcurves can be distinguished : 14 lightcurves are detected to be periodic ; beyond that , we find quasi - periodic lightcurves , where the variability has an apparent timescale , but is not regular enough to be detected as periodic ; aperiodic lightcurves without a preferred timescale ; and long - term variable lightcurves where variability is apparent , but no periodicity or time - scale can be determined within the observational window .
in addition , there are lightcurves with short , non - repeating bursts or dips .
roughly half of all sources become redder when they are fainter ; the other half becomes bluer .
the reddening values of the first group are compatible with ism reddenig .
the color changes in the second group require variability in the inner disk structure as proposed by @xcite or @xcite .
this work is based on observations made with the spitzer space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa .
support for this work was provided by nasa through an award issued by jpl / caltech .
this research made use of astropy , a community - developed core python package for astronomy @xcite .
this research has made use of the simbad database and the vizier catalogue access tool , both operated at cds , strasbourg , france and of data products from the two micron all sky survey , which is a joint project of the university of massachusetts and the infrared processing and analysis center / california institute of technology , funded by the national aeronautics and space administration and the national science foundation .
acknowledges spitzer grant 1490851 .
h.y.a.m .  and p.p .
acknowledge support by the ipac visiting graduate fellowship program at caltech / ipac .
p.p .  also acknowledges the jpl research and technology development and exoplanet exploration programs .
rag gratefully acknowledges funding support from nasa adap grants nnx11ad14 g and nnx13af08 g and caltech / jpl awards 1373081 , 1424329 , and 1440160 in support of spitzer space telescope observing programs .
s.j.w . was supported by nasa contract nas8 - 03060 .
_ facilities : _
: :    we identify 2mass j16272802 - 2439335 with because they match within the    positional accuracy and the observed flux densities fit nicely    together in a sed . : :    after a visual inspection of the 2mass images , which form the basis of    @xcite , we also identify 2mass j16273288 - 2428116 with . : :    based on the position and the expected shape of the sed for a yso , we    also identify [ gpj2008 ] source 3 with .
: :    we propose that [ amd2002 ] j162724 - 242850 and are the same source ,    because their positions match and the fluxes form an sed that is fully    consistent with a yso in l1688 .   and
: :    [ sect : misidentified ] in comparison with @xcite , we find that is    erroneously identified with source oph  l1688  3 in simbad .
the distance    between the positions on the sky is 13 , much larger than the    positional uncertainty .
therefore , the name is not used in our source    table .
ylw 16 : :    in one case , the simbad database does contain a reference to a    multiple source ( ) , that is resolved in the ir .
we associate our    lightcurves with , which is the dominant component .
: :    the ukidss @xmath31-band image shows at least two sources    within 1 of iso - oph  152 with partially overlapping point - spread    functions ( sourceids 442426789990 and 442426789990 ) .
the position of    the second source is uncertain and fitted differently in the other    bands ( where it is called source i d 442426789991 ) .
all of this is not    resolved in the irac data ; thus , in each band , we assign the brightest    magnitude of either of those sources to iso - oph  152 .
[ sect : g09notwil ]    : :    this source detected in several near - ir observations and is known to    be variable in that wavelength range .
this source is part of the triple system that is well resolved in    sub - mm observations . based on its sed stretching out to 6  cm ,
this is    a class  i source , which matches the classification we derive from the    ir data .
most likely , this souce contributes to the seizable outflows    observed from the triple system , which would further confirm its youth    and thus its status as a cluster member @xcite .
: :    this source is detected at @xmath106 m , but not at    @xmath107 m by @xcite but no further information is    available in the literature .
: :    in @xcite , this objects is listed as a yso and potential member of    l1688 based on the g09 sed classification and a match to a scuba    source at @xmath108 m . however , the next peak of the scuba    flux is 29 from the position of [ gmm2009 ] oph l1688 30 , just below the    maximum distance for a match that is accepted in that work .
: :    this source is also listed as a yso with a scuba flux peak 20 away    @xcite .
: :    we find that this source is mismatched in simbad ( see    section  [ sect : misidentified ] ) .
there are several radio observations in    molecular lines and in the continuum which indicate substantial    circumstellar material and thus confirm the status as a yso , but if
sstc2d j162621.7 - 242250 and [ gmm2009 ] oph l1688 3 are different    sources either of them could be the source of the radio signatures .                                                                                ,
a. , whitney , b. , wood , k. , plavchan , p. , terebey , s. , stauffer , j.  r. , morales - calderon , m. , & ysovar .
2013 , in american astronomical society meeting abstracts , vol .
221 , american astronomical society meeting abstracts , # 256.10