Tomoe Note 021: Ensemble measurements in 20191120: statistical properties and moving sources

Michael Richmond
Jun 12, 2020

Executive Summary

I've examined the ensemble photometry outputs of quad 1 of 20191120. This note shows some statistical properties of the objects contained in the ensemble output. It also discusses briefly how one can find moving objects by analyzing the ensemble-level output of the pipeline.


The dataset for 20191120

You can find a description of the pipeline processing of this night in

The weather was fair that night, and the observations lasted about 41 minutes. All 84 chips produced good results.

As the FITS headers of one of these images states, this data was taken in 2Hz mode:

OBJECT  = 'J0337+1648_dith1'   / object name                                    
EXPTIME =            59.988480 / [s] total exposure time                        
TELAPSE =            60.500020 / [s] elapsed time                               
EXPTIME1=             0.499904 / [s] exposure time per frame                    
TFRAME  =             0.500000 / [s] frame interval in seconds                  
DATA-FPS=             2.000000 / [Hz] frame per second                          
DATA-TYP= 'OBJECT'             / data type (OBJECT,FLAT,DARK)                   
OBS-MOD = 'Imaging'            / observation mode                               
FILTER  = 'BLANK'              / filter name                                    
PROJECT = 'Earth Shadow 2Hz Survey' / project name                              
OBSERVER= 'Noriaki Arima'      / observer name                                  
PIPELINE= 'wcs,stack,raw'      / reduction pipeline templete                    

Now, the pipeline produces output of two kinds:

This note examines some properties of objects in the ensemble output data.


Statistical properties of ensemble objects

During this night, the camera pointed at a few locations and acquired multiple images in each chip at each location. Each set of 120 images were packaged together into a "chunk", and the objects in that chunk were collected into an ensemble. We expect bright objects should appear in all chunks at one particular location -- since they should appear in all images -- while very faint objects might not appear in some chunks.

The graph below shows the median number of chunks in which objects appeared, as a function of their magnitudes. And, indeed, we see the expected pattern: objects brighter than mag 17 are detected in many chunks, while very faint objects are detected in only a few.

Next, let's examine the measurements of brightness of the objects. Each chunk of 120 images yields a single, average magnitude for an object, reported in the ensemble output. If a non-variable object is detected in multiple chunks, we expect that its magnitude ought to be the same in each one. For objects detected in multiple chunks, I computed the mean and standard deviation of all the reported magnitudes. The graph below shows the results.

At the bright end, magnitude measurements suffer from saturation at V < 10; and at the faint end, the uncertainty grows quickly fainter than V = 16.5 or so.

You can see the value of the uncertainty for stars between these limits more clearly in the graph below, which uses a logarithmic scale on the y-axis.

Most objects do not move appreciably during the course of a single night. Therefore, our repeated measurements of their positions should yield the same value. The scatter around the mean position is an indication of the precision of our measurements of position. The graph below shows that scatter around the mean position of objects measured in multiple chunks.

Just as the photometric values are bad for stars brighter than V = 10, the astrometric positions are also poor for these saturated stars. The same is true for the very faint stars: the uncertainty in both the magnitude and position increase after V = 16.5 or so.

The same information is shown in the graph below, but I have zoomed in on the y-axis in order to show the precision for stars of intermediate brightness. For these "good" stars, the positions repeat to within about 0.05 arcseconds.


Moving objects -- based on ensemble measurements

One can search for moving objects in a large dataset in many ways. I chose a method which is probably not a great one, but which was easy to do with the ensemble output. Remember, the ensemble method produces just one position for each chunk of 120 images.

The method was as follows:

  1. start with "ensemble objects", which have one position and magnitude per chunk
  2. collect all "ensemble objects" with similar positions (within a radius of 5 arcseconds) to create a "source." A "source" represents a real, physical object in space. Each "source" has a number of "ensemble objects", which are the detections of it within each chunk.
  3. Make a linear fit to both RA and Dec as functions of time:
    
               RA  =   a  +  b*JD
    
              Dec  =   c  +  d*JD
       
  4. Select as "significant" any motion in which the measured value of coefficient "b" or "d" is larger than 3 times its uncertainty.

Remember, the timespan covered by these measurements is at most 41 minutes. Therefore, this technique can only detect objects with considerable speeds, of at least several arcseconds per hour.

In the quadrant 1 dataset, there were a total of 59 sources identified as having "significant" motion. The graph below shows the motions of each one.

There are two main groups: one centered around an RA motion of -35 arcseconds per hour, and a second group centered around RA motion of 0 arcseconds per hour. In the following text file, the two groups are separated by a set of blank lines:

The first group has motions consistent with those of main-belt asteroids:

In order to identify objects, one can use the

For example, entering into this webpage the time of the observations (Date = 2019 11 20.61 UT), the observatory code for Kiso (381), and the position for the first item in the list (source ID 8238 RA = 03 29 59.0 Dec = 19 26 59.0), then clicking on "Produce list", yields the identification


 The following objects, brighter than V = 19.0, were found in the 5.0-arcminute region around R.A. = 03 29 59.0, Decl. = +19 26 59.0 (J2000.0) on 2019 11 20.61 UT:

 Object designation         R.A.      Decl.     V       Offsets     Motion/hr   Orbit  Further observations?
                           h  m  s     °  '  "        R.A.   Decl.  R.A.  Decl.        Comment (Elong/Decl/V at date 1)

  (4024) Ronan           03 29 59.3 +19 26 59  15.3   0.1E   0.0N    44-     1+   28o  None needed at this time.
 

Note that the motion computed from the Tomo-e measurements is a good match to the prediction in this case:



              measured                        prediction
            (arcsec/hour)                   (arcsec/hour)
------------------------------------------------------------------

RA           -44.70  +/-  1.20                  -44

Dec            1.55  +/-  2.28                    1

------------------------------------------------------------------

The second group of moving sources is centered near the origin:

In most cases, I believe that these motion are NOT real, but simply chance combinations of a small number of positions which are really random noise around a fixed position, but happen to produce an apparently linear trend. Note how the uncertainties in most of these objects are very large.

I'm not sure if the graph below is useful, but in trying to figure out which motions were real, and which spurious, I tried plotting the significance (measured motion / uncertainty in motion) in Dec vs. the significance in RA. Can one use this diagram to distinguish real from spurious motions?