<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML">
<link rel="stylesheet" type="text/css" href="../../../CSS_FULL/edps_full.css">
<body><div id="contenu_olm">

<!-- DOI: 10.1051/0004-6361/200811351 -->

<h2 class="sec">Online Material</h2>

<p>

<p>

<h2 class="sec"><a name="SECTION000100000000000000000"></a>&#160; <A NAME="sect:techcal"></A>
Appendix A: The principles of surface brightness observations with ISOPHOT: data
reduction  and calibration
</h2>

<p>
The most detailed description of the ISOPHOT instrument, its observing modes (so called Astronomical Observation Templates, AOTs) and the corresponding data analysis and calibration steps is given in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>). In the following we describe recent calibration techniques which are beyond the scope of the Handbook and which are essential
for the determination of the EBL&nbsp;surface brightness.

<p>

<h3 class="sec2"><a name="SECTION000101000000000000000"></a>
<A NAME="sect:abssurfacebrightnesscal"></A>
A.1 Absolute surface brightness calibration of ISOPHOT observations
</h3>

<p>
ISOPHOT was absolutely calibrated against a flux grid of celestial point source standards consisting of stars, asteroids and planets, thus covering a fair fraction of the entire dynamic flux range from <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img32.png"
 ALT="$\approx$">100&nbsp;mJy up to about 1000&nbsp;Jy. Each detector aperture/pixel was individually calibrated 
against these standards. Therefore, the basic ISOPHOT calibration is in&nbsp;Jy&nbsp;pixel<sup>-1</sup>.

<p>
In order to derive proper surface brightness values in MJy&nbsp;sr<sup>-1</sup>, the solid angles of each detector aperture/pixel must be accurately known

<br><p></p>
<DIV ALIGN="CENTER">

<!-- MATH: \begin{equation}
%
B_{\rm\lambda} = \frac{f_{\rm psf}^{\rm aper}(0,0)}{ \Omega^{\rm aper}_{\rm eff}}
\cdot F_{\rm\lambda},

\end{equation} -->

<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><IMG
 WIDTH="131" HEIGHT="78"
 SRC="img33.png"
 ALT="\begin{displaymath}%
B_{\rm\lambda} = \frac{f_{\rm psf}^{\rm aper}(0,0)}{ \Omega^{\rm aper}_{\rm eff}}
\cdot F_{\rm\lambda},
\end{displaymath}"></td>
<TD WIDTH=10 ALIGN="RIGHT">
(A.1)</td></tr>
</TABLE>
</DIV><BR CLEAR="ALL"><p></p>
with 
<!-- MATH: $B_{\rm\lambda}$ -->
<IMG
 WIDTH="17" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img34.png"
 ALT="$B_{\rm\lambda}$">
being the surface brightness, 
<!-- MATH: $\Omega^{\rm aper}_{\rm eff}$ -->
<IMG
 WIDTH="32" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img35.png"
 ALT="$\Omega^{\rm aper}_{\rm eff}$">
the effective solid angle of the pixel/aperture, 
<!-- MATH: $F_{\rm\lambda}$ -->
<IMG
 WIDTH="18" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img36.png"
 ALT="$F_{\rm\lambda}$">&nbsp;the total flux of a celestial standard and 
<!-- MATH: $f_{\rm psf}^{\rm aper}(0,0)$ -->
<IMG
 WIDTH="58" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img37.png"
 ALT="$f_{\rm psf}^{\rm aper}(0,0)$">
the fraction of the Point Spread Function contained in the pixel/aperture (i.e.&nbsp;the convolution of the PSF with the aperture response) when being centred at
position&nbsp;(0,&nbsp;0). Hence, 
<!-- MATH: $F_{\rm\lambda} \cdot f_{\rm psf}^{\rm aper}(0,0)$ -->
<IMG
 WIDTH="82" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img38.png"
 ALT="$F_{\rm\lambda} \cdot f_{\rm psf}^{\rm aper}(0,0)$">
is the 
flux per pixel.

<p>
ISOPHOT's effective solid angles have been determined by 2D-scanning of a point source over the pixel/aperture in fine steps&nbsp;d<I>x</I>' and&nbsp;d<I>y</I>' and measuring the resulting intensity at each measurement point&nbsp;
<!-- MATH: $(x'_{\rm i},y'_{\rm j})$ -->
<IMG
 WIDTH="41" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img39.png"
 ALT="$(x'_{\rm i},y'_{\rm j})$">,
the footprint, taking into account a non-flat aperture/pixel response:
<br><p></p>
<DIV ALIGN="CENTER">

<!-- MATH: \begin{equation}
%
\Omega^{\rm aper}_{\rm eff} =
\sum_{i}^{A}\sum_{j}~f_{\rm psf}^{\rm aper}(x'_{\rm i},y'_{\rm j})~{\rm d}x'~{\rm d}y'.

\end{equation} -->

<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><IMG
 WIDTH="219" HEIGHT="82"
 SRC="img40.png"
 ALT="\begin{displaymath}%
\Omega^{\rm aper}_{\rm eff} =
\sum_{i}^{A}\sum_{j}~f_{\rm psf}^{\rm aper}(x'_{\rm i},y'_{\rm j})~{\rm d}x'~{\rm d}y'.
\end{displaymath}"></td>
<TD WIDTH=10 ALIGN="RIGHT">
(A.2)</td></tr>
</TABLE>
</DIV><BR CLEAR="ALL"><p></p>
If the peak of the point source was located outside the aperture by&nbsp;<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">1/2 of the aperture size, the S/N of the resulting intensity dropped so much that at this border the summation was complemented by a model of the broad band telescope PSF and adding up the corresponding PSF&nbsp;fractions 
<!-- MATH: $f_{\rm psf}^{\rm aper}(x'_{\rm i},y'_{\rm j})$ -->
<IMG
 WIDTH="66" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img41.png"
 ALT="$f_{\rm psf}^{\rm aper}(x'_{\rm i},y'_{\rm j})$">
out to <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img10.png"
 ALT="$\pm$">10&nbsp;arcmin assuming a flat response, but taking into account a cut by ISO's pyramidal central mirror feeding the 4&nbsp;instrument beams. An example of such a synthetic footprint is shown in&nbsp;Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100_60umsyntheticfootprint">A.1</a>.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:c100_60umsyntheticfootprint">&#160;</A><!-- end Label--><A NAME="959"></A><A NAME="figure565"
 HREF="img42.png"><IMG
 WIDTH="95" HEIGHT="91" SRC="Timg42.png"
 ALT="\begin{figure}
\par\includegraphics[angle=90,width=8.5cm,clip]{1351fg04.eps}
\end{figure}"></A><!-- HTML Figure number: 2 --></td>
<td class="img-txt"><span class="bold">Figure A.1:</span><p>
Synthetic (outer part, i.e. green and blue coloured areas, modelled) footprints (convolution of the ISO telescope PSF with the pixel aperture response) of the 3&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
3&nbsp;pixels of ISOPHOT's C100 array for the 60&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m broad band filter. The solid angles of each pixel are obtained by integration over the footprint area.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=2&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>
The values of the solid angles used in PIA V11.3 are listed in Tables&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c100omega">A.1</a> and&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c200omega">A.2</a>.

<p>
<A NAME="tab:c100omega"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableA.1.html"><span class="bold">Table A.1:</span></a>&#160;&#160;
Effective solid angles for the 3&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
3&nbsp;pixels of ISOPHOT's C100 array for the 6&nbsp;filters with central wavelengths&nbsp;
<!-- MATH: $\lambda_{\rm c}$ -->
<IMG
 WIDTH="15" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img4.png"
 ALT="$\lambda _{\rm c}$">.</p>
<p>
<A NAME="tab:c200omega"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableA.2.html"><span class="bold">Table A.2:</span></a>&#160;&#160;
Effective solid angles for the 2&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
2&nbsp;pixels of ISOPHOT's C200&nbsp;array for the 5&nbsp;filters with central wavelengths&nbsp;
<!-- MATH: $\lambda_{\rm c}$ -->
<IMG
 WIDTH="15" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img4.png"
 ALT="$\lambda _{\rm c}$">.</p>
<p>
It should be noted that an absolute surface brightness calibration is more accurate than an absolute calibration of a compact source of similar brightness, since no background subtraction has to be performed, which introduces an additional uncertainty. The accuracies quoted in the ISOPHOT
Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>), Table&nbsp;9.1 for extended emission take COBE/DIRBE photometry as the reference. By not referring to COBE/DIRBE photometry, the absolute surface brightness calibration for ISOPHOT's C100 and C200 array is as good as that for bright compact sources, i.e. better than&nbsp;15%.

<p>

<h3 class="sec2"><a name="SECTION000102000000000000000"></a>
A.2 New calibration products and strategies for PIA V11.3
</h3>

<p>
For the very sensitive analysis needed for the EBL&nbsp;determination and, in particular, an absolute surface brightness calibration that is as accurate as possible, a number of calibration upgrades and new calibration features have been developed and implemented in PIA&nbsp;V11.3. For the ones which are not described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>), we provide a description and examples for the C100&nbsp;detector in the following. An overview of the individual calibration steps associated with different instrument components is shown in Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:phtcalibrationschema">A.2</a>. By
application of all these steps, instrumental artefacts are minimized, the resulting detector signals are homogenized and a high calibration reproducibility and accuracy is achieved.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:phtcalibrationschema">&#160;</A><!-- end Label--><A NAME="966"></A><A NAME="figure607"
 HREF="img46.png"><IMG
 WIDTH="95" HEIGHT="59" SRC="Timg46.png"
 ALT="\begin{figure}
\par\includegraphics[width=8.5cm,clip]{1351fg05.eps}
\end{figure}"></A><!-- HTML Figure number: 3 --></td>
<td class="img-txt"><span class="bold">Figure A.2:</span><p>
Scheme of the ISOPHOT calibration steps associated with the different instrument components. The meaning of the abbreviations is the following: BSL&nbsp;= Bypassing Sky Light correction, DS&nbsp;= detector Dark Signal, RL&nbsp;= Ramp Linearisation, TC&nbsp;= signal Transient Correction, and RIC&nbsp;= Reset Interval Correction.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=3&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h4 class="sec3"><a name="SECTION000102100000000000000"></a>
A.2.1 Detector responsivity calibration
</h4>

<p>
The absolute photometric calibration of an individual measurement is performed via a transfer calibration using the internal calibration sources. This measures the actual responsivity of the detector and provides the absolute signal-to-flux conversion. It is a separate measurement of each 
observation mode by deflecting the chopper mirror to the field of view of the internal calibrator (Fine Calibration Source, FCS). The illumination level of the internal calibrator was not fixed but adjusted as much as possible to the expected brightness level of the sky. This was achieved by selecting an appropriate heating power for the internal source. There exists a calibration 
relation between this heating power and the optical power received by each detector pixel which is established from measurements on celestial standards.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:fcscalschema">&#160;</A><!-- end Label--><A NAME="969"></A><A NAME="figure613"
 HREF="img47.png"><IMG
 WIDTH="192" HEIGHT="113" SRC="Timg47.png"
 ALT="\begin{figure}
\par\mbox{\hspace*{2cm}\includegraphics[clip,height=5cm,angle=0]{...
...}\includegraphics[clip=true,width=5cm,angle=0]{1351fg10.eps} }
\par
\end{figure}"></A><!-- HTML Figure number: 4 --></td>
<td class="img-txt"><span class="bold">Figure A.3:</span><p>
Steps in the generation of a homogeneous and most complete calibration of ISOPHOT's long wavelength internal calibration sources (FCS). This is illustrated for the central pixel&nbsp;(#5) of the C100 array camera. <I> Upper left</I>: measured relation between optical power received on the detector and the heating power applied to the internal source. Dots indicate the discrete measurements, the solid line is a fit. <I> Upper right</I>: display of the input curves for all C100&nbsp;filters within the reliable heating power range. <I> Middle left</I>: for a selected heating power (here: 1.0&nbsp;mW) monochromatic and colour corrected fluxes of all filters are fitted by a modified BB&nbsp;curve. <I> Middle right</I>: by repeating the fits with the same modified BB&nbsp;type for the whole heating power range covered the relation between heating power and temperature of the internal source is established. <I> Lower centre</I>: by applying the FCS&nbsp;model the relation between optical power and heating power is homogenized and extended to the maximum heating power range covered by at least one measurement in any of the C100 or C200&nbsp;filters.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=4&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>
Therefore, for reliable and accurate transfer calibrations, the following requirements are put on the&nbsp;FCS:
<DL COMPACT>
<DT>1)
<DD>High reproducibility. This was better than 1%, since the monitoring
of the flux of faint standards was reproducible within a few percent,
          and this uncertainty was dominated by the signal noise (Klaas et&nbsp;al. <A NAME="aaref23"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Klaas01">2001</a>).

<p>
<DT>2)
<DD>A very detailed characterization of the illuminated power depending
          on the heating power applied to the source. This is illustrated in
          Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a>. It involves the following steps:
          <DL COMPACT>
<DT>a)
<DD>For each C100 and C200 array filter all measurements of
celestial standards done in raster map mode were evaluated
                    such that for each pixel the background signal was 
                    properly subtracted and the resulting source signal was 
                    associated with the celestial standard flux. The ratio of 
                    the source signal and the simultaneously obtained FCS 
                    signal gave the illumination power by the FCS for the 
                    selected heating power. The discrete results were fitted 
                    and the reliable lower and upper heating power limits 
                    covered by measurements were identified 
                    (Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> upper left). The heating power
                    ranges were not identical or equally large for each filter
                    (Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> upper right). In general they
                    were shifted to smaller heating power values for longer
                    wavelengths.

<p>
<DT>b)
<DD>For fine discrete steps in heating power the inband powers
                    were read from the relations and were converted to 
                    monochromatic surface brightnesses by applying the 
                    bandpass conversions derived from the relative system 
                    response profiles, see ISOPHOT Handbook 
                    (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>), Sect.&nbsp;A.2,
                    and the solid angles of Tables&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c100omega">A.1</a> 
                    and&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c200omega">A.2</a>. These fluxes were fitted with
                    a modified BB&nbsp;curve after appropriate colour correction
                    (Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> middle left). If for a 
                    certain filter the selected heating power was outside
                    the reliable limits, the value of this filter was excluded
                    from the fit. The fit gave the temperature of the FCS for
                    the selected heating power. An additional constraint was
                    that the temperature had to be the same for the fit curves
                    of all pixels. C100 and C200&nbsp;filter values had to be 
                    fitted independently because of the different detector 
                    areas and hence illumination factors, however, the fits 
                    were checked for consistent temperatures, because the 
                    illuminating FCS was the same for both detectors.

<p>
<DT>c)
<DD>This was achieved for the heating power range from 0.07
                    up to 6.5&nbsp;mW adopting an emissivity of the source 
                    
<!-- MATH: ${\propto}\lambda^{1.25}$ -->
<IMG
 WIDTH="37" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img48.png"
 ALT="${\propto}\lambda^{1.25}$">
yielding the temperature&nbsp;vs. 
                    heating power relation as shown in Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> middle right.

<p>
<DT>d)
<DD>By applying this FCS temperature model and the established
                    illumination factors for each pixel it was possible to
                    establish homogeneous calibration curves of the internal
                    reference source, thus polishing out measurement outliers
                    affecting the initial empirical curves. The multi-filter
                    approach connecting all curves and not treating them
                    individually enabled a large extension and a common range
                    for all filters: compare Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> lower
                    centre with Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:fcscalschema">A.3</a> upper right. 
          </DL></DL> 

<p>

<h4 class="sec3"><a name="SECTION000102200000000000000"></a>
A.2.2 Bypassing sky light correction of FCS signal
</h4>

<p>
As a safety design against single point failures, ISOPHOT was not equipped with any cold shutter to suppress straylight when performing internal calibration measurements. Therefore, when deflecting the chopper onto the illuminated internal calibration (FCS) sources, some fraction of the power 
received on the detector did not come from the FCS but from sky light bypassing along non nominal light paths. Since this depends on the sky brightness it is subtracted in the transfer calibration measurements on celestial standards and hence has to be subtracted for any FCS&nbsp;measurement
in order to get a reproducible zero point. This was achieved by performing a number of measurements on the switched-off, i.e. cold FCS, so that only the bypassing sky light contribution was measured.

<p>
The result for one C100 array pixel is shown in Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100pix4bypassskylightcorr">A.4</a> which demonstrates a linear dependence of the bypassing sky light contribution to the FCS&nbsp;signal on the 
sky background. This correction was established for all C100 and C200&nbsp;array pixels. The bypassing sky light contribution contains the detector dark signal contribution, cf.&nbsp;Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:darksignal">A.2.5</a>. 

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:c100pix4bypassskylightcorr">&#160;</A><!-- end Label--><A NAME="970"></A><A NAME="figure649"
 HREF="img50.png"><IMG
 WIDTH="93" HEIGHT="64" SRC="Timg50.png"
 ALT="\begin{figure}
\par\includegraphics[clip=true,width=8.5cm,angle=0]{1351fg11.eps}
\end{figure}"></A><!-- HTML Figure number: 5 --></td>
<td class="img-txt"><span class="bold">Figure A.4:</span><p>
Bypassing sky light contribution to the FCS&nbsp;signal depending on the sky background.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=5&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h4 class="sec3"><a name="SECTION000102300000000000000"></a>
A.2.3 Effective pixel/aperture solid angles
</h4>

<p>
These are described in the previous Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:abssurfacebrightnesscal">A.1</a> and their values are compiled in Tables&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c100omega">A.1</a> and&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#tab:c200omega">A.2</a>.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:c100pix5darksignal">&#160;</A><!-- end Label--><A NAME="971"></A><A NAME="figure658"
 HREF="img51.png"><IMG
 WIDTH="94" HEIGHT="69" SRC="Timg51.png"
 ALT="\begin{figure}
\par\includegraphics[clip=true,width=8.5cm,angle=0]{1351fg12.eps}
\end{figure}"></A><!-- HTML Figure number: 6 --></td>
<td class="img-txt"><span class="bold">Figure A.5:</span><p>
Orbit dependent dark signal determination for the central pixel&nbsp;5 of ISOPHOT's C100 array. Dots represent individual measurements obtained during the entire ISO mission, filled and open signals identify a different reset interval in the integration of the dark signal. The solid line is the fit to the measurements providing the so-called default dark level. The dotted line is the default dark level of an older calibration version used before&nbsp;2001.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=6&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h4 class="sec3"><a name="SECTION000102400000000000000"></a>
A.2.4 Filter profiles
</h4>

<p>
The bandpass system responses and the conversion factors from inband power to a monochromatic flux, as well as colour correction factors are described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>).

<p>

<h4 class="sec3"><a name="SECTION000102500000000000000"></a>
<A NAME="sect:darksignal"></A>
A.2.5 Detector dark signal
</h4>

<p>
The detector dark signals were re-analyzed as described in del Burgo (<A NAME="aaref3"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#delBurgo02">2002</a>). In this analysis special care was given to exclude those dark measurements suffering from memory effects by preceding bright illuminations, thus not representing the true dark level. An example of the new results is shown in Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100pix5darksignal">A.5</a> for the central pixel&nbsp;5 of the C100&nbsp;array. A slight orbital dependence is visible with an increase of the dark signal towards the beginning and the end of the observational window. It can also be noticed that there is a scatter of the dark signals at the same orbit position and there are occasional large outliers. These are not due to signal determination uncertainties, but are real variations due to space weather effects on different revolutions over the ISO&nbsp;mission.

<p>

<h4 class="sec3"><a name="SECTION000102600000000000000"></a>
A.2.6 Ramp linearisation
</h4>

<p>
This was performed as described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>). For ISOPHOT's far-infrared detectors two types of effects cause non-linearities of the integration ramps:
<DL COMPACT>
<DT>1)
<DD>De-biasing effects of the photoconductors operated with low bias 
caused by feed-back from the integration capacitor.

<p>
<DT>2)
<DD>Non-linearities generated in the cold read-out electronics.
</DL>
<p>

<h4 class="sec3"><a name="SECTION000102700000000000000"></a>
A.2.7 Signal dependence on reset interval correction
</h4>

<p>
Despite the ramp linearisation step, signals obtained under constant illumination, but with different reset intervals show a systematic difference, see Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100resetintervalcorr">A.6</a> upper panel. In order to have a consistent signal handling of measurements with different reset interval settings applied&nbsp;- to optimize the dynamic range of the signal&nbsp;- all signals were converted as if they were taken with a&nbsp;1/4&nbsp;s reset interval. The correction relations were established from special calibration measurements applying the full suite of reset intervals under constant illumination and this for different illumination levels. In this way signal corrections were established for all reset intervals in the range 1/32&nbsp;s to 8&nbsp;s (Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100resetintervalcorr">A.6</a> middle and lower panel). While previously, as still described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>) a linear correlation with offset was used, a re-analysis (del Burgo et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#delBurgo02">2002</a>) yielded non-linear relations as shown in
Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100resetintervalcorr">A.6</a>.  This latter analysis also found a bi-modal behaviour for C100
array pixels, such that the pixels on the main diagonal, #1, 5 and&nbsp;9, behaved differently from the rest of the pixels. For the C200&nbsp;array all pixels behaved in the same way.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:c100resetintervalcorr">&#160;</A><!-- end Label--><A NAME="972"></A><A NAME="figure681"
 HREF="img52.png"><IMG
 WIDTH="97" HEIGHT="227" SRC="Timg52.png"
 ALT="\begin{figure}
\par\includegraphics[clip=true,width=8.7cm,angle=0]{1351fg13.eps}...
...s}\par\includegraphics[clip=true,width=8.7cm,angle=0]{1351fg15.eps}
\end{figure}"></A><!-- HTML Figure number: 7 --></td>
<td class="img-txt"><span class="bold">Figure A.6:</span><p>
Correction of the signal dependence on the selected reset interval. <I> Upper panel</I>: demonstration of the effect, showing the resulting signal versus
the selected reset interval over the range from 1/32&nbsp;s up to 8&nbsp;s (reset intervals were commanded in powers of&nbsp;2) under constant illumination. <I> Middle panel</I>: solid line: correction relation 
for a reset interval of 8&nbsp;s w.r.t. the reference reset interval of 1/4&nbsp;s for all C100 array pixels, except the ones on the main diagonal. Dotted line: old linear correlation used before the 
re-analysis. <I> Lower panel</I>: solid line: Correction relation for a reset interval of 8&nbsp;s w.r.t. the reference reset interval of 1/4&nbsp;s for all C100&nbsp;array pixels on the main diagonal (pixels&nbsp;#1, 
5, and&nbsp;9). Dotted line: old linear correlation used before the re-analysis (same as for <I> middle panel</I>).</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=7&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h4 class="sec3"><a name="SECTION000102800000000000000"></a>
A.2.8 Signal transient correction
</h4>

<p>
The ISOPHOT detectors were photoconductors operated under low background conditions provided by a cryogenically cooled spacecraft. Under these conditions they showed the behaviour that the output signal was not instantaneously adjusted to a flux change but rather, following an initial
jump by a certain fraction of the flux step, the signal adjusted with some time constant to the final level, see e.g. Acosta et&nbsp;al. (<A NAME="aaref1"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Acosta-Pulido00">2000</a>). In particular the ISOPHOT C100&nbsp;detector showed significant transient behaviour. This time constant depended on the detector
material (doping of the semi-conductor and its contacts), the flux step, the direction of the flux step (dark to bright versus bright to dark) and the illumination history. Attempts had been made to model this behaviour (Acosta et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Acosta-Pulido00">2000</a>), but no unique description could be found for the FIR&nbsp;detectors. To overcome this effect at least partly the method of transient recognition was implemented in the ISOPHOT analysis software as described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>) using the most stable part of the measurement for signal determination. Finally, another approach was to use a data base of long measurements with signals stabilising and to determine the deviation from the end level for shorter intermediate times (del Burgo et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#delBurgo02">2002</a>), see Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100transientcorr">A.7</a> for an illustration. A measurement time of 128&nbsp;s was used as reference, because
<DL COMPACT>
<DT>1)
<DD>Most calibration measurements in staring mode were performed with this basic measurement time.
<p>
<DT>2)
<DD>In most cases the signals stabilised within this time.
</DL>For the C200 array pixels the signal transient effect is considerably smaller and faster and therefore it is sufficient to apply the transient recognition as described in the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>).

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fi:c100transientcorr">&#160;</A><!-- end Label--><A NAME="973"></A><A NAME="figure702"
 HREF="img53.png"><IMG
 WIDTH="95" HEIGHT="119" SRC="Timg53.png"
 ALT="\begin{figure}
\par\includegraphics[clip=true,width=8.5cm,angle=0]{1351fg16.eps}
\end{figure}"></A><!-- HTML Figure number: 8 --></td>
<td class="img-txt"><span class="bold">Figure A.7:</span><p>
Empirical signal transient correction for ISOPHOT's C100 array. The left column shows the signal loss for integration times of&nbsp;4, 8, 16, 32, and 64&nbsp;s (commendable integration times of ISOPHOT&nbsp;detectors) with regard to the reference time of 128&nbsp;s. The red line is a fit through the measured points over the covered signal range and is used as the correction relation. The right column shows the residuals after applying this correction.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=8&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h2 class="sec"><a name="SECTION000110000000000000000"></a>&#160; <A NAME="app:obs"></A>
Appendix B: Observations and data reduction  for the EBL fields
</h2>

<p>

<h3 class="sec2"><a name="SECTION000111000000000000000"></a>
B.1 ISOPHOT observations
</h3> 

<p>
The following tables give details of the ISOPHOT observations used in the paper. Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:tdt">B.1</a> lists the raster maps and absolute photometry measurements that were made at&nbsp;90, 150, and 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m. Correspondingly, Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_observations">B.2</a> lists observations used for the determination of the zodiacal light levels. These include both mid-infrared measurements carried out with the ISOPHOT-P&nbsp;detector and longer wavelength absolute photometry measurements carried out with the C100 and C200&nbsp;cameras.

<p>
<A NAME="table:tdt"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableB.1.html"><span class="bold">Table B.1:</span></a>&#160;&#160;
List of ISOPHOT observations of EBL fields carried out in the
PHT-22 and PHT-25 observation modes.</p>
<p>
<A NAME="table:ZL_observations"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableB.2.html"><span class="bold">Table B.2:</span></a>&#160;&#160;
Observations used for the determination of the zodiacal light emission. The columns are: (1)&nbsp;name of the field; (2),&nbsp;(3)&nbsp;position; (4)&nbsp;wavelength; (5)&nbsp;the ISO identifier number (TDT) of the observation; and (6)&nbsp;time difference between the listed observation and the observation of the EBL raster maps of Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:tdt">B.1</a>. A time difference is quoted only when the observations were not performed within the same day. Observations at wavelengths below 60&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m are made with the ISOPHOT-P&nbsp;detector.</p>
<p>
Each field was mapped in the PHT22 staring raster map mode (ISOPHOT Handbook, Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>) using filters&nbsp;C_90, C_135, and C_180. The corresponding reference wavelengths of the filters are 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, 150&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, and 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m. The 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m observations were made with the C100&nbsp;detector consisting of 3&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
3&nbsp;pixels, with 43.5
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
43.5
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">
each. The longer wavelength observations were made with the C200&nbsp;detector which has 2&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
2&nbsp;detector
pixels, with 89.4
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
89.4
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">
each. The same raster maps were used in Juvela et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Juvela00">2000</a>). Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:fields">B.3</a> lists the coordinates and the sizes of the maps.
Additionally, we make use of PHT25 absolute photometry measurements (see ISOPHOT Handbook, Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>) made at the same three wavelengths. Two positions in NGP, two positions in EBL26, and one position in EBL22 were observed in this mode.

<p>
<A NAME="table:fields"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableB.3.html"><span class="bold">Table B.3:</span></a>&#160;&#160;
The positions and sizes of the observed fields. Columns are: (1)&nbsp;name
of the field; (2),&nbsp;(3)&nbsp;equatorial coordinates of the centre of the field; (4),&nbsp;(5)&nbsp;galactic coordinates; (6),&nbsp;(7)&nbsp;ecliptic coordinates; (8)&nbsp;number of raster points; (9)&nbsp;area in square degrees; and (10)&nbsp;additional remarks. All areas were observed at 90, 150, and 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m. In NGP(N) an additional square map was observed at 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m only. Details of the individual measurements are listed in Appendix, in Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:tdt">B.1</a>.</p>
<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fig:allsky">&#160;</A><!-- end Label--><A NAME="979"></A><A NAME="figure737"
 HREF="img57.png"><IMG
 WIDTH="203" HEIGHT="132" SRC="Timg57.png"
 ALT="\begin{figure}
\par\resizebox{18cm}{!}{\includegraphics{1351fg17.eps}}
\end{figure}"></A><!-- HTML Figure number: 9 --></td>
<td class="img-txt"><span class="bold">Figure B.1:</span><p>
The ISOPHOT EBL fields. The three frames show the 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m ISOPHOT maps. The coordinates correspond to the 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps. The  90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps cover the same area but, in the figure, the 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps have been plotted south of the 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps. The small yellow circles indicate the positions observed with the ISOPHOT P-detector for the determination of the zodiacal light levels. To indicate the locations of the fields with respect to the galactic and ecliptica planes, the positions are shown on an all-sky map that is combined from DIRBE observations between the wavelengths of 12&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and&nbsp;240&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=9&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h3 class="sec2"><a name="SECTION000112000000000000000"></a>&#160; <A NAME="sect:isophot_reduction"></A>
B.2 Reduction of EBL field observations
</h3>

<p>
The ISOPHOT data were processed with PIA (PHT Interactive Analysis) program version&nbsp;11.3. For details of the analysis steps, see the ISOPHOT Handbook (Laureijs et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Handbook">2003</a>) and Appendix, Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:techcal">A</a>. For C100 a method of signal transient correction was introduced in PIA&nbsp;11.3. This procedure was used for all C100&nbsp;measurements. Nevertheless, some of the internal calibrator (FCS) measurements show residual drifts. In those cases we applied transient recognition which removes the initial, unstabilised part of the measurements. The flux density calibration
was made using the internal calibrator measurements (FCS1) performed immediately before and after each map for actual detector response assessment. The calibration was applied using the average response of the two FCS&nbsp;measurements. 

<p>
The reduced data contained a few artifacts. These include short time scale detector drifts at the beginning of some C100&nbsp;observations, temporary signal variations caused by cosmic ray glitches, and occasional drifting of some detector pixels that may also be connected with cosmic ray hits. The time ordered data were examined by eye. For rasters and detector pixels affected by clear anomalies (glitches or drifting) the corresponding PIA&nbsp;error estimates were scaled upwards, typically by a factor of a&nbsp;few. For each detector pixel the signal values were scaled so that their average
value over a map became equal to the overall average over all detector pixels. The scaling takes into account the already manually adjusted error estimates. The flat fielding would actually not be necessary, because FIR&nbsp;fluxes are compared only with observed HI&nbsp;21&nbsp;cm lines and, therefore, averaged over areas that are large compared with the size of the ISOPHOT&nbsp;rasters.

<p>
Long term detector response drifts are not taken out by a simple averaging of the FCS&nbsp;measurements, nor is an initial non-linear drift corrected for by linear interpolation between the two FCS&nbsp;measurements. Both could introduce an artificial gradient in the time ordered data and, because of the systematic scan pattern, also in the maps. The maps were compared with IRAS data in order
to see if there were any gradients uncorrelated with the IRAS 100&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m signal. The only significant difference was found in the C200&nbsp;observations of the southern NGP&nbsp;field. The gradient was removed while keeping the average surface brightness unchanged. The correction has little effect on the subsequent analysis. Apart from the EBL22&nbsp;field, all maps contain four detector scans that run alternatively in opposite directions along the longer map dimension. When data are correlated with the lower resolution HI&nbsp;observations, the subsequent scan legs tend to cancel out any long term
drifts.

<p>
The raster map observations themselves do not contain any direct measurement of the dark current. In such cases one usually relies on the orbit dependent ``default'' dark current estimates included in the PIA.  However, absolute photometry PHT-25&nbsp;measurements were carried out within a couple of hours
before or after each raster map. The data reduction was carried out also using the dark current and cold FCS&nbsp;values obtained from those measurements. In the subsequent analysis, we use maps that are averages of those obtained using default dark current values and those obtained using PHT-25&nbsp;dark current measurements.

<p>
When absolute photometry points were inside the mapped area they were compared with the surface brightness of the raster maps. The maps were re-scaled so that the final surface brightness corresponds to the average of the original FCS calibrated maps and the values given by the absolute photometry measurements. This causes systematic lowering of the surface brightness values
of the original maps. For EBL26, NGP(N), and NGP(S) the change is typically <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">4%, for both C100 and C200&nbsp;observations. In the case of EBL22 the correction is larger, some 20%, for the C200&nbsp;detector.

<p>
In the region NGP there are separate northern and southern fields that overlap by a few arcminutes. The maps, each containing 32&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
4&nbsp;raster points, were fitted together using the overlapping area, where the final map is at a level equal to the average of the northern and the southern maps.  The resulting change in the surface brightness levels of individual maps was <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">5% or less. In the north there is yet another 15&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
15&nbsp;raster map that was observed only at 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m. Because that measurement includes only very short FCS&nbsp;measurements, it was scaled to fit the already combined long 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m map. This required scaling of the surface brightness values by a factor of&nbsp;1.05.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fig:HI_positions">&#160;</A><!-- end Label--><A NAME="983"></A><A NAME="figure748"
 HREF="img58.png"><IMG
 WIDTH="89" HEIGHT="197" SRC="Timg58.png"
 ALT="\begin{figure}
\par\resizebox{8cm}{!}{\includegraphics{1351fg18.eps}}\par\resize...
...51fg19.eps}}\par\resizebox{8cm}{!}{\includegraphics{1351fg20.eps}}
\end{figure}"></A><!-- HTML Figure number: 10 --></td>
<td class="img-txt"><span class="bold">Figure B.2:</span><p>
The figures shows as black rectangles the areas mapped with ISOPHOT&nbsp;(90, 150, and 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m) and as circles the pointings used in the Effelsberg HI&nbsp;observations. The diameter of the circles, 9
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">,
is equal to the <I> FWHM</I> of the Effelsberg beam.  The frames&nbsp;<B> a)</B>-<B> c)</B> correspond to regions&nbsp;EBL22, EBL26, and NGP. In the case of NGP, the dashed red line indicates the area that was mapped at 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m only.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=10&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>
The main maps of the field EBL22 cover an area of low cirrus emission. There are additional one-directional scans that extend to a region of higher surface brightness in the west. In the absence of scans in the opposite direction, it is not possible to directly determine the presence of detector response drifts. However, these observations were reduced using the average of the responsivities given by the two FCS&nbsp;measurements and the error bars reflect also the difference in the responsivity before and after the measurement. Using the overlapping area, the 32&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
1&nbsp;raster strips were scaled to the same level with the 32&nbsp;<IMG
 WIDTH="12" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img2.png"
 ALT="$\times $">
3&nbsp;raster maps. The scalings applied were&nbsp;0.97, 1.02, and 0.84 at 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, 150&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, and 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, respectively.

<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fig:HI_comparison">&#160;</A><!-- end Label--><A NAME="985"></A><A NAME="figure764"
 HREF="img59.png"><IMG
 WIDTH="100" HEIGHT="78" SRC="Timg59.png"
 ALT="\begin{figure}
\par\resizebox{9cm}{!}{\includegraphics{1351fg21.eps}}
\end{figure}"></A><!-- HTML Figure number: 11 --></td>
<td class="img-txt"><span class="bold">Figure B.3:</span><p>
Comparison of HI spectra from the Leiden/Dwingeloo survey (Kalberla et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kalberla2005">2005</a>; dashed lines) and our Effelsberg data convolved with a beam of&nbsp;36
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">
(solid lines). The spectra correspond to positions at the southern and northern end of the NGP map (13<IMG
 WIDTH="9" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
 SRC="img5.png"
 ALT="$^{\rm h}$">42<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img6.png"
 ALT="$^{\rm m}$">0<IMG
 WIDTH="8" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img7.png"
 ALT="$^{\rm s}$">
+40
<!-- MATH: $\hbox{$^\circ$ }$ -->
<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img8.png"
 ALT="$\hbox {$^\circ $ }$">30
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">0
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">
and 13<IMG
 WIDTH="9" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
 SRC="img5.png"
 ALT="$^{\rm h}$">52<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img6.png"
 ALT="$^{\rm m}$">0<IMG
 WIDTH="8" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img7.png"
 ALT="$^{\rm s}$">
+38
<!-- MATH: $\hbox{$^\circ$ }$ -->
<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img8.png"
 ALT="$\hbox {$^\circ $ }$">40
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">0
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">)
and one position in the field&nbsp;EBL26 (1<IMG
 WIDTH="9" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
 SRC="img5.png"
 ALT="$^{\rm h}$">17<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img6.png"
 ALT="$^{\rm m}$">0<IMG
 WIDTH="8" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img7.png"
 ALT="$^{\rm s}$">
+2
<!-- MATH: $\hbox{$^\circ$ }$ -->
<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img8.png"
 ALT="$\hbox {$^\circ $ }$">20
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">0
<!-- MATH: $\hbox{$^{\prime\prime}$ }$ -->
<IMG
 WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img9.png"
 ALT="$\hbox {$^{\prime \prime }$ }$">). The EBL26&nbsp;spectra have been scaled by a&nbsp;factor&nbsp;0.5.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=11&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>
The final FIR errorbars show the uncertainty for the weighted means over the Effelsberg beam. The noise of each HI&nbsp;spectrum was estimated separately using the velocity channels outside the line. The uncertainty of the line area was calculated assuming the same, uncorrelated noise for the integrated velocity interval. This might underestimate the total uncertainty, because it ignores the uncertainties in the stray radiation subtraction that do not affect the signal in the line wings. However, for a small field the stray radiation causes a constant systematic error rather than statistical uncertainty and does not affect the weighting of the observations when the linear fit is&nbsp;made.

<p>
For selected positions there exist mid-infrared observations made with the ISOPHOT P-detectors as well as further absolute photometry measurements with the C100 and C200&nbsp;cameras (see Appendix, Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_observations">B.2</a>). These observations were performed for the purpose of estimating the zodiacal light. The data reduction of P-detector data is similar to that of the C100 and C200&nbsp;cameras, except that also signal linearisation is included.

<p>

<h3 class="sec2"><a name="SECTION000113000000000000000"></a>&#160; <A NAME="sect:HImeas"></A>
B.3 HI measurements
</h3>

<p>
The observations of the hydrogen 21&nbsp;cm line were made with the Effelsberg radio telescope in May&nbsp;2002. The observed positions, 580 in number, are indicated in Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:HI_positions">B.2</a>. The integration times were 30&nbsp;s in EBL22 and 62&nbsp;s in EBL26. In the field NGP the observations were done with 62&nbsp;s integrations except for the northern part where the integration time was 94&nbsp;s. The average noise estimated from the velocity intervals outside the HI&nbsp;line is 0.15&nbsp;K&nbsp;per channel of 1&nbsp;km&nbsp;s<sup>-1</sup>. This corresponds to a typical uncertainty of 1.7&nbsp;K&nbsp;km&nbsp;s<sup>-1</sup> in the integrated line area. 

<p>
For calibration purposes and for precise subtraction of the stray radiation, regular observations of the standard region&nbsp;S7 were made. The stray radiation subtraction is crucial because it affects the zero point of the estimated HI&nbsp;column densities. The observed fields, NGP in particular, have some of the lowest line-of-sight column densities over the whole sky. Under these conditions the stray radiation received by the telescope side lobes becomes a significant fraction of the total signal. The stray radiation was removed with a program developed by Kalberla (Kalberla et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kalberla2005">2005</a>; see Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:obs">3</a>).

<p>
In Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:HI_comparison">B.3</a> we compare our data with spectra from the Leiden/Dwingeloo survey (Hartmann &amp; Burton <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Hartmann1997">1997</a>; Kalberla et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kalberla2005">2005</a>). For this comparison, in order to match the resolution of the Leiden/Dwingeloo survey, the Effelsberg data were convolved
with a Gaussian with <I>FWHM</I> equal to&nbsp;36
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">.
The HI&nbsp;profiles agree very well. Part of the differences may be caused by the fact that our HI&nbsp;maps do not cover the whole area of the 36
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">&nbsp;beams. Nevertheless, the figure shows that the observations and the stray radiation subtraction (see Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:obs">3</a>) are consistent with the Kalberla (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kalberla2005">2005</a>)
results.

<p>

<h2 class="sec"><a name="SECTION000120000000000000000"></a>&#160; <A NAME="sect:calibration"></A>
Appendix C: Calibration accuracy
</h2>

<p>
The error estimates listed in Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:fit">1</a> are based on the statistical uncertainties in the fits between FIR and HI&nbsp;data.  The scatter of data points around the fitted lines is usually larger than their estimated uncertainty. This could be a sign of underestimated measurement uncertainties but is more likely caused by true scatter in the relation. If the formal uncertainties of the line parameters were estimated based on the error estimates of the individual points, the uncertainties could be severely underestimated. Therefore, instead of relying only on the measurement uncertainties, the uncertainty of the fit parameters was estimated separately with the bootstrap method so that they reflect the true scatter of observed points. The error estimates corresponding to a&nbsp;67% confidence interval are given in&nbsp;Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:fit">1</a>. 

<p>
These uncertainties do not include estimates for the systematic errors introduced by the independent calibration of each map or the absolute accuracy of the overall ISOPHOT calibration. There are both multiplicative and additive sources of uncertainty. The former include, for example, uncertainties in the internal calibration source (FCS) measurements (e.g., detector drifts) that alter the estimated detector response. The uncertainties that affect the zero point of the intensity scale are more critical, because the CIRB is small compared with the observed signal and can be recovered only as the residual after the subtraction of the&nbsp;ZL. 

<p>
Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:uncertainties">C.1</a> lists an assessment of uncertainty that, using data in Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:fit">1</a>, have been converted into uncertainty of the FIR&nbsp;flux at zero HI&nbsp;column density. The quoted values are half of the difference of two values obtained in two independent ways. Thereby the quoted values are also <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">1-<IMG
 WIDTH="13" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img11.png"
 ALT="$\sigma$">
estimates for the uncertainty of the average of the two values.

<p>
In Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:uncertainties">C.1</a> Col.&nbsp;4 has been obtained by comparing the fine calibration source measurements performed before and after each map. The numbers indicate the statistical uncertainty of the detector response measurements. The FCS&nbsp;measurements are generally very consistent, particularly in the case of the C200&nbsp;detector. On the other hand, the effect of the drift affecting the first FCS&nbsp;measurement of the one-dimensional strip map of EBL22 is clearly visible at&nbsp;90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.

<p>
The dark signal subtraction is the most important correction affecting the zero point of the FIR&nbsp;intensity. Close to each of the raster map observations, we have one or two absolute photometry observations which include dark signal measurements of their own. In PIA, the default dark current calibration is based on a larger set (<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">70) of dark current measurements for which the
orbit trend has been determined. Therefore, the PIA default dark current calibration is less affected by the noise of individual measurements but may not take into account short time scale variations in the detector dark current on a specific orbit. The maps were reduced using the default dark current values and the actually measured dark current values. In Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:uncertainties">C.1</a> Col.&nbsp;5 shows the associated uncertainty in the FIR&nbsp;signal at zero HI&nbsp;column density. The observed uncertainty in the dark current values is comparable with the
variation observed in the systematic analysis of a large sample of ISOPHOT observations (del Burgo et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#delBurgo02">2002</a>; see also Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fi:c100pix5darksignal">A.5</a>).

<p>
When absolute photometry measurements existed within mapped areas, those were used to re-scale the surface brightness values of the maps (see Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:isophot_reduction">B.2</a>). The difference in the absolute photometry and mapping measurements is used to derive the values in Col.&nbsp;6 of Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:uncertainties">C.1</a>. The final column reflects the difference in the surface brightness in areas where two independently calibrated maps overlap. The numbers in Cols.&nbsp;6 and&nbsp;7 include, of course, dark current and FCS&nbsp;uncertainties as one of their components. For the C100&nbsp;observations at 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m the uncertainty is close to 1&nbsp;MJy&nbsp;sr<sup>-1</sup>, i.e., comparable with the expected EBL&nbsp;signal. On the other hand, for the C200&nbsp;detector the uncertainty of an individual map is
<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">0.3&nbsp;MJy&nbsp;sr<sup>-1</sup>. Most of this is caused by the uncertainty in the dark current values. 

<p>
<A NAME="table:uncertainties"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableC.1.html"><span class="bold">Table C.1:</span></a>&#160;&#160;
Assessment of the calibration uncertainty for the ISOPHOT&nbsp;maps. Columns are (1)&nbsp;name of the field; (2)&nbsp;wavelength; (3)&nbsp;average surface brightness of the map; (4)&nbsp;difference between calibration measurements performed before and after each map; (5)&nbsp;difference between actual dark current measurements and default dark current values; (6)&nbsp;difference between the independently calibrated absolute photometry measurements and raster maps; and (7)&nbsp;difference between partially overlapping maps. These uncertainties have been converted to correspond to the uncertainties at zero hydrogen column density using the fit parameters listed in&nbsp;Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:fit">1</a>.</p>
<p>
<div class="inset-old">
<table>
<tr><td><!-- init Label --><A NAME="fig:zlfit">&#160;</A><!-- end Label--><A NAME="989"></A><A NAME="figure822"
 HREF="img62.png"><IMG
 WIDTH="89" HEIGHT="171" SRC="Timg62.png"
 ALT="\begin{figure}
\par\resizebox{8cm}{!}{\includegraphics{1351fg22.eps}}
\end{figure}"></A><!-- HTML Figure number: 12 --></td>
<td class="img-txt"><span class="bold">Figure C.1:</span><p>
Fits used to estimate the ZL levels in the three fields EBL22, EBL26, and NGP (frames&nbsp;<B> a)</B>-<B> c)</B>, respectively). The red circles are ISOPHOT observations. The lower lines are the cirrus (blue solid line) and the ZL (red solid line) templates, the uppermost solid green line is their sum.  The figures also show DIRBE&nbsp;values for the closest DIRBE pixel, read from the DIRBE weekly maps. The solid squares correspond to observations with the same solar elongation as in the case of ISOPHOT observations, the open squares to the other measurement with identical absolute value of the solar aspect angle but opposite solar elongation. For clarity, the latter have been shifted slightly in wavelength. The dashed line shows the predictions of the Kelsall et&nbsp;al. (<A NAME="aaref18"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kelsall98">1998</a>) ZL&nbsp;model.</p></td>
</tr><tr><td colspan="2"><a href="http://dexter.edpsciences.org/applet.php?pdf_id=12&DOI=10.1051/0004-6361/200811351" target="DEXTER">Open with DEXTER</a></td></tr>

</table></div>
<p>

<h3 class="sec2"><a name="SECTION000121000000000000000"></a>
C.1 Straylight radiation
</h3>

<p>
Straylight may be another instrumental artefact affecting the zero level of the FIR&nbsp;surface brightness. By design and operation ISO's viewing direction stayed by several tens of degrees away from the brightest FIR&nbsp;emitters in the sky, the Sun, the Earth and the Moon (Kessler et&nbsp;al.
<A NAME="aaref19"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kessler2003">2003</a>). A dedicated straylight program was executed verifying by deep ``differential'' integrations that the uniform straylight level due to these sources was below ISOPHOT's detection limit, even under the most unfavourable pointing conditions close to the visibility constraints (Lemke et&nbsp;al. <A NAME="aaref32"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Lemke2001">2001</a>).

<p>
Specular straylight by the second brightest class of objects, the giant planets Jupiter and Saturn, was observed when pointing to within 15
<!-- MATH: $\hbox{$^\prime$ }$ -->
<IMG
 WIDTH="11" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img1.png"
 ALT="$\hbox {$^\prime $ }$">
to&nbsp;1
<!-- MATH: $\hbox{$^\circ$ }$ -->
<IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img8.png"
 ALT="$\hbox {$^\circ $ }$">
of the planet, expressing itself as finger-like stripes or faint ghost rings (Kessler et&nbsp;al. <a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kessler2003">2003</a>; Lemke et&nbsp;al.
<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Lemke2001">2001</a>). The NGP and EBL&nbsp;22&nbsp;fields are far away from the ecliptic and can thus not suffer from this type of straylight. For EBL&nbsp;26 we checked the positions of the planets Mars, Jupiter, Saturn, Uranus and Neptune at the time of the observations, 1997-06-26 and 1997-07-11, respectively. Mars, Jupiter, Uranus, Neptune were all far off. Saturn was at a distance of 3.25&nbsp;degrees, which is still more than a factor of 3&nbsp;off of any known straylight-critical distance. 

<p>

<h2 class="sec"><a name="SECTION000130000000000000000"></a>&#160; <A NAME="sect:zlfit"></A>
Appendix D: Determination of the ZL levels
</h2>

<p>
The ZL level was estimated by fitting ZL and cirrus templates to ISOPHOT observations in the wavelength range from 7.3&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m to 200&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m (see Sect.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#sect:ZL">4.2</a>).  Figure&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:zlfit">C.1</a> shows the results of these fits. In the field&nbsp;EBL22 we had observations of one position and in the fields&nbsp;EBL26  and NGP of two positions (see Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_observations">B.2</a>). For the latter two fields, the figures show the fit to data combined from the two positions.

<p>
Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_observations">B.2</a> lists the time difference between the listed observations and the observations of the raster maps. In the case of NGP these are relative to the 150&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m observations. The 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps were observed four days before and the 180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m one day after the 150&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m maps. According to the Kelsall et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kelsall98">1998</a>) ZL&nbsp;model the four day difference causes only <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">1.5% change in the expected&nbsp;ZL. The combined NGP&nbsp;map is almost 1.5&nbsp;degrees long. In the Kelsall model the difference in the centre positions of the southern and northern parts corresponds to about&nbsp;1% difference in the ZL. Therefore, we use only one zodiacal estimate value for both NGP(N) and NGP(S) and for all observations made during the five day interval.

<p>
In the fields EBL26 and NGP, MIR observations exist for two separate positions (see Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:allsky">B.1</a>). In both fields, the measurements at these two positions are close to each other, both in time and position. Therefore, their ZL&nbsp;values should be identical and also the cirrus levels should be very similar. Comparison of the fits performed using these independent sets of
measurements gives the first indication of the statistical uncertainty of the ZL&nbsp;values. In both fields, the ZL&nbsp;values obtained for the two positions agree within&nbsp;10%. 

<p>
The observations are fitted as a sum of ZL and cirrus components. The ZL&nbsp;template is a black body curve at the temperature obtained from Leinert et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Leinert2002">2002</a>). The cirrus template is based on the model by Li &amp; Draine (<A NAME="aaref33"></A><a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Li2001">2001</a>). Using the ISOPHOT filter profiles we calculate for both radiation components, ZL and cirrus, and for each filter the in-band power values that can be directly compared with the observed values. In the fit we have only two free parameters, the intensity of the ZL&nbsp;component and the intensity of the cirrus component. The ZL estimates should be based mainly on data between 10&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and 60&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m where the ZL is clearly the dominant
component. Therefore, in the fit, the weight of the data points in this
wavelength range is increased by a factor of&nbsp;two. 

<p>
The level of the cirrus component is determined mostly by the longer wavelength data. In reality, the component corresponds to the sum of the cirrus and CIRB&nbsp;signals. As long as the component is small in the MIR, the ZL estimates are almost independent of the exact shape of this template. We
confirmed this by replacing the Li &amp; Draine (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Li2001">2001</a>) cirrus template by a pure CIRB template, using the model curve from Dole et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Dole2006">2006</a>; Fig.&nbsp;13). The resulting change in the ZL estimates was less than one per cent.

<p>
The actual statistical errors of the ZL&nbsp;values are estimated using the standard deviation of the relative errors when observations are compared with the fitted ZL&nbsp;curve. The last column of Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_estimates">2</a> lists the corresponding error of the mean, calculated using data points between 7.3&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.  In the case of fields&nbsp;NGP and EBL26, the error estimates are calculated from the fits where we have combined the data from the two measured positions within each field. In all three fields, the obtained relative uncertainties are <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">10%. In the fields&nbsp;EBL26 and NGP the uncertainties are also consistent with the difference of the ZL&nbsp;values obtained for the two individual positions. The ZL&nbsp;fits are shown in&nbsp;Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:zlfit">C.1</a>.

<p>
In this paper we have used original ISOPHOT observations without applying colour corrections. Therefore, in the fitting procedure also the ZL and cirrus templates were converted to corresponding values using the ISOPHOT filter profiles. However, for Fig.&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#fig:zlfit">C.1</a> we have performed colour
corrections. The templates are plotted by connecting the values at the nominal wavelengths by a straight line. The template spectra used in the ZL&nbsp;fitting are colour corrected using their respective spectral shapes. In the figure, the colour correction of the observed surface brightness values is done assuming the blackbody ZL&nbsp;spectrum below 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, and a modified black body
cirrus spectrum, 
<!-- MATH: $B_{\nu}(T=18~{\rm K}) \nu^2$ -->
<IMG
 WIDTH="93" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img63.png"
 ALT="$B_{\nu}(T=18~{\rm K}) \nu^2$">,
at 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and longer wavelengths.

<p>
The plots include DIRBE values from the DIRBE weekly maps. These correspond to the DIRBE pixel closest to the centre of the corresponding ISOPHOT map. Linear interpolation was performed between the weeks in order to accurately match the solar elongation of the ISOPHOT observations. In addition to the DIRBE&nbsp;value that corresponds directly to the ISOPHOT observations (solid squares) we plot the DIRBE value for the same solar aspect angle and opposite sign of the solar elongation. Assuming that the zodiacal dust cloud is symmetric along the ecliptic, the two values should be identical.  The
predictions of the ZL&nbsp;model of Kelsall et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Kelsall98">1998</a>) are also plotted. The DIRBE&nbsp;values are colour corrected. As in the case of ISOPHOT data, colour correction of the observations assumes a blackbody ZL&nbsp;spectrum at and below 60&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m, and a modified black body cirrus spectrum, 
<!-- MATH: $B_{\nu}(T=18~{\rm K}) \nu^2$ -->
<IMG
 WIDTH="93" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
 SRC="img63.png"
 ALT="$B_{\nu}(T=18~{\rm K}) \nu^2$">,
at the longer wavelengths, 100, 140, and&nbsp;240&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.

<p>
There is a clear difference in the ISOPHOT and DIRBE surface brightness scales. The DIRBE&nbsp;values are consistently lower by some 20-30%, in the MIR range. In the FIR&nbsp;bands the extended cirrus structures combined with the much larger pixel size and noise in the DIRBE&nbsp;pixels precludes direct comparison. The determination of the CIRB&nbsp;values is not directly affected by a possible calibration difference between DIRBE and ISOPHOT because, in this paper, we use exclusively ISOPHOT&nbsp;measurements. 

<p>
Systematic uncertainties affecting all ISOPHOT bands have only little impact on the derived CIRB&nbsp;values. The relative calibration accuracy between the FIR&nbsp;cameras and the ISOPHOT-P&nbsp;photometer is more important, because the zodiacal light estimates are based on the latter. When the absolute level of the zodiacal light was estimated we calculated the scatter between the SED&nbsp;model and the observations at different wavelengths (see Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:ZL_estimates">2</a>). The scatter was typically <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">10-20%. The importance of this error source depends, of course, on the absolute level of the ZL&nbsp;emission. The field&nbsp;EBL26 is located near the ecliptic plane and at 90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m the observed signal and the ZL are both of the order of&nbsp;20&nbsp;MJy&nbsp;sr<sup>-1</sup>. Therefore, a relative uncertainty of 10% would already correspond to about twice the expected level of the CIRB. For EBL22 and especially for NGP the zodiacal light level is much lower so that more meaningful limits can be derived for the CIRB also at&nbsp;90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.

<p>
The quoted ZL error estimates reflect the uncertainty in the determined ZL&nbsp;level in the mid-infrared. If there were a systematic difference in the calibration of the mid- and FIR-bands, the ZL estimates could be wrong by the corresponding amount. Generally the relative calibration accuracy is considered to be within 15%. This uncertainty would not necessarily be reflected in the quality of the ZL&nbsp;spectrum fits, because a systematic calibration error could have been partly compensated by a change in the intensity of the cirrus component.

<p>
The ZL spectrum was assumed to be a pure black body with the temperature given by Leinert et&nbsp;al. (<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#Leinert2002">2002</a>). As far as the mid-infrared points are concerned, a wrong temperature would, at some level, be reflected also in our error estimate. However, if the ZL&nbsp;spectrum deviated from the assumed shape only in the FIR this could again be masked by a change in the fitted cirrus component without a corresponding increase in the rms&nbsp;value. Therefore, we must explicitly assume that the same ZL&nbsp;temperature is applicable both at mid-infrared and far-infrared wavelengths. However, because a&nbsp;5&nbsp;K change in the ZL&nbsp;temperature corresponds to only <IMG
 WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
 SRC="img12.png"
 ALT="$\sim$">2% relative change in the ratio of 150&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and 25&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m intensities, this source of uncertainty is unimportant compared with the uncertainty in the relative calibrations of the different detectors.

<p>

<h2 class="sec"><a name="SECTION000140000000000000000"></a>&#160; <A NAME="sect:dirbe_ebl"></A>
Appendix E: Comparison with DIRBE EBL estimates
</h2>

<p>
This present study represents the first determination of the absolute level of the FIR&nbsp;EBL that is independent of measurements of the COBE DIRBE&nbsp;instrument. In Table&nbsp;<a href="/articles/aa/full_html/2009/23/aa11351-08/aa11351-08.html#table:comparison">E.1</a> we list FIR&nbsp;EBL estimates given in seven publications based on the DIRBE&nbsp;measurements. Included are also our
2-<IMG
 WIDTH="13" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img11.png"
 ALT="$\sigma$">&nbsp;upper limit at&nbsp;90&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m and the EBL estimate for the range 150-180&nbsp;<IMG
 WIDTH="11" HEIGHT="26" ALIGN="MIDDLE" BORDER="0"
 SRC="img3.png"
 ALT="$\mu $">m.

<p>
<A NAME="table:comparison"></A><p class="inset-old"><a href="/articles/aa/full_html/2009/23/aa11351-08/tableE.1.html"><span class="bold">Table E.1:</span></a>&#160;&#160;
Comparison of existing CIRB estimates in the FIR range. The
error estimates quoted by the authors are shown in parenthesis. In our case, we include only the statistical uncertainty.</p>
<p>
<br>

</div></body></html>