abstract: The pervasive use of new mobile devices has allowed a better characterization in space and time of human concentrations and mobility in general. Besides its theoretical interest, describing mobility is of great importance for a number of practical applications ranging from the forecast of disease spreading to the design of new spaces in urban environments. While classical data sources, such as surveys or census, have a limited level of geographical resolution (e.g., districts, municipalities, counties are typically used) or are restricted to generic workdays or weekends, the data coming from mobile devices can be precisely located both in time and space. Most previous works have used a single data source to study human mobility patterns. Here we perform instead a cross-check analysis by comparing results obtained with data collected from three different sources: Twitter, census and cell phones. The analysis is focused on the urban areas of Barcelona and Madrid, for which data of the three types is available. We assess the correlation between the datasets on different aspects: the spatial distribution of people concentration, the temporal evolution of people density and the mobility patterns of individuals. Our results show that the three data sources are providing comparable information. Even though the representativeness of Twitter geolocated data is lower than that of mobile phone and census data, the correlations between the population density profiles and mobility patterns detected by the three datasets are close to one in a grid with cells of $2\times 2$ and $1\times 1$ square kilometers. This level of correlation supports the feasibility of interchanging the three data sources at the spatio-temporal scales considered.
author: Maxime Lenormand; Miguel Picornell; Oliva G. Cantú-Ros; Antònia Tugores; Thomas Louail; Ricardo Herranz; Marc Barthelemy; Enrique Frías-Martinez; José J. Ramasco
title: Cross-checking different sources of mobility information

# INTRODUCTION

The strong penetration of ICT tools in the society's daily life is opening new opportunities for the research in socio-technical systems . Users' interactions with or through mobile devices get registered allowing a detailed description of social interactions and mobility patterns. The sheer size of these datasets opens the door to a systematic statistical treatment while searching for new information. Some examples include the analysis of the structure of (online) social networks , human cognitive limitations , information diffusion and social contagion , the role played by social groups , language coexistence or even how political movements raise and develop .

The analysis of human mobility is another aspect to which the wealth of new data has notably contributed . Statistical characteristics of mobility patterns have been studied, for instance, in Refs. , finding a heavy-tail decay in the distribution of displacement lengths across users. Most of the trips are short in everyday mobility, but some are extraordinarily long. Besides, the travels are not directed symmetrically in space but show a particular radius of gyration . The duration of stay in each location also shows a skewed distribution with a few preferred places clearly ranking on the top of the list, typically corresponding to home and work . All the insights gained in mobility, together with realistic data, have been used as proxies for modeling the way in which viruses spread among people or among electronic devices . Recently, geolocated data has been also used to analyze the structure of urban areas , the relation between different cities or even between countries .

Most mobility and urban studies have been performed using data coming essentially from a single data source such as: cell phone data , geolocated tweets , census-like surveys or commercial information . There is only a few recent exceptions, for instance, epidemic spreading studies . When the data has not been "generated" or gathered ad hoc to address a specific question, one fair doubt is how much the results are biased by the data source used. In this work, we compare spatial and temporal population density distributions and mobility patterns in the form of Origin-Destination (OD) matrices obtained from three different data sources for the metropolitan areas of Barcelona and Madrid. This comparison will allow to discern whether or not the results are source dependent. In the first part of the paper the datasets and the methods used to extract the OD tables are described. In the second part of the paper, we present the results. First, a comparison of the spatial distribution of users according to the hour of the day and the day of the week showing that both Twitter and cell phone data are highly correlated on this aspect. Then, we compare the temporal distribution of users by identifying where people are located according to the hour of the day, we show that the temporal distribution patterns obtained with the Twitter and the cell phone datasets are very similar. Finally, we compare the mobility networks (OD matrices) obtained from cell phone data, Twitter and census. We show that it is possible to extract similar patterns from all datasets, keeping always in mind the different resolution limits that each information source may inherently have.

# MATERIALS AND METHODS

This work is focused on two cities: the metropolitan areas of Barcelona and Madrid both in Spain and for which data from the three considered sources is available. The metropolitan area of Barcelona contains a population of $3,218,071$ (2009) within an area of $636$ $km^2$. The population of the metropolitan area of Madrid is larger, with $5,512,495$ inhabitants (2009) within an area of $1,935$ $km^2$ . In order to compare activity and intra mobility in each city, the metropolitan areas are divided into a regular grid of square cells of lateral size $l$ (Figure b). Two different sizes of grid cells ($l=1$ $km$ and $l=2$ $km$) are considered in order to evaluate the robustness of the results. Since mobility habits and population concentration may change along the week, we have divided the data into four groups: one, from Monday to Thursday representing a normal working day and three more for Friday, Saturday and Sunday.

The concentration of phone or Twitter users is quantified by defining two three dimensional matrices $T=(T_{g,w,h})$ and $P=(P_{g,w,h})$, accounting, respectively, for the number of Twitter users and the number of mobile phone users in the grid cell $g$ at the hour of the day $h$ and for the group of days $w$. The index for cells $g$ runs in the range $[1,n]$. In the following, details for the three datasets are more thoroughly described.

## Mobile phone data

The cell phone data that we are analyzing come from anonymized users' call records collected during 55 days (noted as $D$ hereafter) between September and November 2009. The call records are registered by communication towers (Base Transceiver Station or BTS), identified each by its location coordinates. The area covered by each tower can be approximated by a Voronoi tessellation of the urban areas, as shown in Figure a for Barcelona. Each call originated or received by a user and served by a BTS is thus assigned to the corresponding BTS Voronoi area. In order to estimate the number of people in different areas per period of time, we use the following criteria: each person counts only once per hour. If a user is detected in $k$ different positions within a certain 1-hour time period, each registered position will count as ($1/k$) "units of activity". From such aggregated data, activity per zone and per hour is calculated. Consider a generic grid cell $g$ for a day $d$ and hour between $h$ and $h+1$, the $m$ Voronoi areas intersecting $g$ are found and the number of mobile phone users $P_{g,d,h}$ is calculated as follows: $$P_{g,d,h}=\sum_{v=1}^m N_{v,d,h} \, \frac{\displaystyle A_{v\cap g}}{\displaystyle A_{v}},
\label{Pgdh}$$ where $N_{v,d,h}$ is the number of users in a Voronoi cell $v$ on day $d$ at time $h$, $A_{v\cap g}$ is the area of the intersection between $v$ and $g$, and $A_v$ the area of $v$. The $D$ days available in the database are then divided in four groups according to the classification explained above and the average number of mobile phone users for each day group $w$ is computed as $$P_{g,w,h}=\frac{\sum_{d\in D_{w}}P_{g,d,h}}{|D_{w}|} .
\label{Pgwh}$$

The number of mobile phone users per day for the two the metropolitan areas as a function of the time of day, and according to the day group, are displayed in Figure . The curves in Figure a show two peaks, one between noon and $3$pm and another one between $6$pm and $9$pm. They also show that the number of mobile phone users is higher during weekdays than during the weekends. The same curve is obtained for Madrid with about twice the number of users with respect to Barcelona. Further details about the data pre-processing are given in the Appendix (Section *Mobile phone data pre-processing*, Figure S1 and Figure S2).

In order to extract OD matrices from the cell phone calls a subset of users, with a mobility reliably recoverable, was selected. For this analysis we only consider commuting patterns in workdays. The users' home and work are identified as the Voronoi cell most frequently visited on weekdays by each user between $8$ pm and $7$ am (home) and between $9$ am and $5$ pm (work). We assume that there must be a daily travel between home and work location of each individual. Users with calls in more than $40\%$ of the days under study at home or work are considered valid. Aggregating the complete flow over users, an OD commuting matrix is obtained containing in each element the flow of people traveling between a Voronoi cell of residence and another of work. Since the Voronoi areas do not exactly match the grid cells, a transition matrix to change the scale is employed (see Appendix for details).

## Twitter data

The dataset comprehends geolocated tweets of $27,707$ users in Barcelona and $50,272$ in Madrid in the time period going from September 2012 to December 2013. These users were selected because it was detected from the general data streaming with the Twitter API that they have emitted at least a geolocated tweet from one of the two cities. Later, as a way to increase the quality of our database, a specific search over their most recent tweets was carried out . As for the cell phone data, the number of Twitter users $T_{g,w,h}$ in each grid cell $g$ per hour $h$ were computed for each day group $w$. The number of Twitter users per day for the metropolitan area of Barcelona according to the hour of the day and the day group is plotted on Figure b. Analogous to the mobile phone data, this figure shows two peaks, one between noon and 3pm and another one between 6pm and 9pm. It is worth noting that the mobile phone users represents on average $2\%$ of the total population against $0.1\%$ for the Twitter data. Furthermore, in contrast with the phone users profile curve, the Twitter users' profile curve shows that the number of users does not vary much from weekdays to weekend days. Moreover, we can observe that the number of Twitter users is higher during the second peak than during the first one.

The identification of the OD commuting matrices using Twitter is similar to the one explained for the mobile phones except for two aspects. Since the number of geolocated tweets is much lower than the equivalent in calls per user, the threshold for considering a user valid is set at 100 tweets on weekdays in all the dataset. The other difference is that since the tweets are geolocated with latitude and longitude coordinates, the assignment to the grid cells is done directly without the need of intermediate steps through the Voronoi cells. As for the phone, we keep only users working and living within the metropolitan areas.

## Census data

The Spanish census survey of 2011 included a question referring to the municipality of work of each interviewed individual. This survey has been conducted among one fifth of the population. This information, along with the municipality of the household where the interview was carried out, allows for the definition of OD flow matrices at the municipal level . For privacy reasons, flows with a number of commuters lower than 10 have been removed. The metropolitan area of Barcelona is composed of $36$ municipalities, while the one of Madrid contains $27$ municipalities. In addition to the flows, we have obtained the GIS files with the border of each municipality from the census office. This information is used to map the OD matrices from Twitter or the cell phone data to this more coarse-grained spatial scale to compare mobility patterns across datasets.

# RESULTS

## Spatial distribution

A first question to address is how much the human activity level is similar or not when estimated from Twitter, $T$, or from cell phone data $P$ across the urban space in grid cells of $2$ by $2$ $km$. To quantify similarity, we start by depicting in Figure a scatter plot composed by each pair $(T_{g,w,h},P_{g,w,h})$ for every grid cell of the metropolitan area of Barcelona taking $w$ as the weekdays (aggregation from Monday to Thursday). The hour $h$ is set from midday to 1pm. A first visual inspection tells us that the agreement between the activity inferred from each dataset is quite good. In fact, the Pearson correlation coefficient between the two estimators of activity is of $\rho = 0.96$. Furthermore, the portion of activity can be depicted on two maps as in Figure b and c. The similarity of the areas of concentration of the activity is patent.

More systematically, we plot in Figure a, the box-plots of the Pearson correlation coefficients for each day group and both case studies as observed for different hours. We obtain in average a correlation of $0.93$ for Barcelona and $0.89$ for Madrid. Globally, the correlation coefficients have higher value for Barcelona than for Madrid probably because the metropolitan area of Madrid is about four times larger than the one of Barcelona. It is interesting to note that the average correlation remains high even if we increase the resolution by using a value of $l$ equal to $1$ $km$. Indeed, we obtain in average a correlation of $0.85$ for Barcelona and $0.83$ for Madrid at that new scale (Figure b).

## Temporal distribution

After the spatial distribution of activity, we investigate the correlation between the temporal activity patterns as observed from each grid cell. We start by normalizing $T$ and $P$ such that the total number of users at a given time on a given day is equal to $1$ $$\begin{aligned}
\hat{T}_{g_0,w,h} & =  \frac{\displaystyle T_{g_0,w,h}}{\sum_{g=1}^n \displaystyle T_{g,w,h}}, \\
\hat{P}_{g_0,w,h} & =  \frac{\displaystyle P_{g_0,w,h}}{\sum_{g=1}^n \displaystyle P_{g,w,h}} .
\label{chap}
\end{aligned}$$ This normalization allows for a direct comparison between sources with different absolute user's activity. For a given grid cell $g=g_0$, we defined the temporal distribution of users $\hat{P}_{g_0}$ as the concatenation of the temporal distribution of users associated with each day group. For each grid cell we obtained a temporal distribution of users represented by a vector of length $96$ corresponding to the $4 \times 24$ hours.

After removing cells with zero temporal distribution, cells of common temporal profies were found using the ascending hierarchical clustering (AHC) method. The average linkage clustering and the Pearson correlation coefficient were taken as agglomeration method and similarity metric, respectively . We have also implemented the k-means algorithm for extracting clusters but better silhouette index values were obtained with the AHC algorithm (see details in Figure S3 in Appendix). To choose the number of clusters, we used the average silhouette index $\bar{S}$ . For each cell $g$, we can compute $a(g)$ the average dissimilarity of $g$ (based on the Pearson correlation coefficient in our case) with all the other cells in the cluster to which $g$ belongs. In the same way, we can compute the average dissimilarities of $g$ to the other clusters and define $b(g)$ as the lowest average dissimilarity among them. Using these two quantities, we compute the silhouette index $s(g)$ defined as $$s(g)=\frac{b(g)-a(g)}{max\{a(g),b(g)\}} ,
\label{sg}$$ which measures how well clustered $g$ is. This measure is comprised between $-1$ for a very poor clustering quality and $1$ for an appropriately clustered $g$. We choose the number of clusters that maximize the average silhouette index over all the grid cells $\bar{S}=\sum_{g=1}^n s(g)/n$.

For the mobile phone data, three clusters were found with an average silhouette index equal to $0.38$ for Barcelona and to $0.43$ for Madrid. The three temporal distribution patterns of mobile phone users are shown in Figure for Barcelona. These three clusters can be associated with the following land uses:

-   **Business:** this cluster is characterized by a higher activity during the weekdays than the weekend days. In Figure a, we observe that the activity takes place between $6$ am and $3$ pm with a higher activity during the morning.

-   **Residential:** this cluster is characterized by a higher activity during the weekend days than during the weekdays. Figure c shows that the activity is almost constant from $9$ am during the weekend days. During the weekdays we observe two peaks, the first one between $7$ am and $8$ am and the second one during the evening.

-   **Nightlife:** this cluster is characterized by a high activity during the night especially the weekend (Figure e).

It is remarkable to note that we obtain the same three patterns for Madrid and that these patterns are robust for different values of the scale parameter $l$ (see details in Figure S4, S5 and S6 in Appendix).

For Twitter data, considering a number of clusters smaller than 10, silhouette index values lower than 0.1 are obtained for both case studies. These low values mean that no clusters have been detected in the data probably because the Twitter data are too noisy. A way to bypass this limitation is to check if, for both data sources, the same patterns are obtained considering the different clusters obtained with the mobile phone data. To do so the temporal distribution patterns of Twitter users associated with the three clusters obtained with the mobile phone data are computed. We note in Figure that for Barcelona the temporal distribution patterns obtained with the Twitter data are very similar to those obtained with the mobile phone data. We obtain the same correlation for Madrid and for different values of the scale $l$ (see details in Figure S4, S5 and S6 in Appendix).

## Users' mobility

In this section, we study the similarity between the OD matrices extracted from Twitter and cell phone data. As it involves a change of spatial resolution needing extra attention, the comparison with the census is relegated to a coming section. We are able to infer for the metropolitan areas of Barcelona and Madrid the number of individuals living in the cell $i$ and working in the cell $j$. Figure shows a scattered plot with the comparison between the flows obtained in the OD matrices for links present in both networks. In order to compare the two networks, the values have been normalized by the total number of commuters.

The overall agreement is good, the Pearson correlation coefficient is around $\rho \approx 0.9$. This coefficient measures the strength of the linear relationship between the normalized flows extracted from both networks, including the zero flows (i.e. flows with zero commuters). However, a high correlation value is not sufficient to assess the goodness of fit. Since we are estimating the fraction of commuters on each link, the values obtained from Twitter and the cell phone data should be ideally not only linearly related but the same. That is, if $y$ if the estimated fraction of mobile phone users on a connection and $x$ the estimated Twitter users on the same link, there should be not only a linear relation, which involves a high Pearson correlation, but also $y = x$. It is, therefore, important to verify that the slope of the relationship is equal to one. To do so, the coefficients of determination $R^2$ are computed to measure how well the scatterplot if fitted by the curve $y = x$. Since there is no particular preference for any set of data as $x$ or $y$, two coefficients $R^2$ can be measured, one using Twitter data as the independent variable $x$ and another using cell phone data. Note that if the slope of the relationship is strictly equal to one the two $R^2$ must be equal to the square of the correlation coefficient, we obtain a value around $R^2 = 0.85$ for Barcelona and around $0.81$ for Madrid. The slope of the best fit is in both cases very close to one.

The dispersion in the points is higher in low flow links. This can be explained by the stronger role played by the statistical fluctuations in low traffic numbers. Moreover, if we increase the resolution by using a value of $l$ equal to $1$ $km$, the Pearson correlation coefficient remains high with a value around $0.8$ (see details in Figure S7 in Appendix). The extreme situation of these fluctuations occurs when a link is present in one network and it has zero flow in the other (missing links). On average $90 \%$ of these links have a number of commuters equal to one in the network in which they are present. This shows that the two networks are not only inferring the same mobility patterns, but that the information left outside in the cross-check corresponds to the weakest links in the system. In order to assess the relevance of the missing links, the weight distributions of these links is displayed in Figure for all the networks and case studies. As a comparison line, the weight distribution of all the links are also shown in the different panels. In all cases, the missing links have flows at least one order of magnitude, sometimes two orders, lower than the strongest links in the corresponding networks. To be more precise, the strongest flow of the missing links is, depending on the case, between $25$ and $464$ times lower than the highest weight of all the links. Furthermore, the average weight of the missing links is between $4$ and $9$ times lower than that obtained over all the links. Most of the missing links are therefore negligible in the general network picture.

With the aim of going a little further, we analyze and compare next the distance distribution for the trips obtained from both datasets. The geographical distance along each link in the OD matrices is calculated and the number of people traveling in the links is taken into account to evaluate the travel-length distribution. Figure shows these distributions for each network. Strong similarity between the two distributions can be observed in the two cities considered.

## Census, Twitter and cell phone

As a final cross-validation, we compare the OD matrices estimated in workdays from Twitter and cell phone data to those extracted from the $2011$ census in Barcelona and Madrid. The census data is at the municipal level, which implies that to be able to perform the comparative analysis the geographical scale of both Twitter and phone data must be modified. To this end, the GIS files with the border of each municipality were used, instead of the grid, to compute the OD matrices from Twitter and cell phone data. Figure shows a scattered plot with the comparison between the flows obtained with the three networks. A good agreement between the three datasets is obtained with a Pearson correlation coefficient around $\rho \approx 0.99$. As mentioned previously, the correlation coefficient is not sufficient to assess the goodness of fit between the two networks. Thus, we have also computed two coefficients of determination $R^2$ for each one of the three relationships to measure how well the line $x=y$ approximates the scatter plots. For the two first relationships, the comparison between the Twitter and the mobile phone and the comparison between the mobile phone and the census OD tables, we obtain $R^2$ values higher than $0.95$. For the last relationship (Twitter vs census), two different $R^2$ values are obtained because the best fit slope of the scatter plot is not strictly equal to one (0.85). The first $R^2$ value, which measure how well the normalized flows obtained in the Twitter's OD matrix approximate the normalized flows obtained in the census's OD matrix, is equal to $0.8$ and the second value, which assess the quality of the opposite relationship, is equal to $0.9$. A better result is instead obtained for Madrid with a Pearson correlation coefficient around $0.99$ and coefficients of determination higher than $0.97$ (see details in Figure S8 in Appendix).

# DISCUSSION

In summary, we have analyzed mobility in urban areas extracted from different sources: cell phones, Twitter and census. The nature of the three data sources is very different, as also is the resolution scales in which the mobility information is recovered. For this reason, the aim of this work has been to run a thorough comparison between the information collected at different spatial and temporal scales. The first aspect considered refers to the population concentration in different parts of the cities. This point is of great importance in the analysis and planning of urban environments, including the design of new services or of contingency plans in case of disasters. Our results show that both Twitter and cell phone data produce similar density patterns both in space and time, with a Pearson correlation close to $0.9$ in the two cities analyzed. The second aspect considered has been the temporal distribution of individuals which allow us to determine the type of activity that are most common in specific urban areas. We show that similar temporal distribution patterns can be extracted from both Twitter and cell phone datasets. The last question studied has been the extraction of mobility networks in the shape of Origin-Destination commuting matrices. We observe that at high spatial resolution, in grid cells with sides of $1$ or $2$ $km$, the networks obtained with both cell phones and Twitter are comparable. Of course, the integration time needed for Twitter is higher in order to obtain similar results. Twitter data can run in serious problems too if instead of recurrent mobility the focus is on shorter term mobility, but this point falls beyond the scope of this work. Finally, the comparison with census data is also acceptable: both Twitter and cell phone data reproduce the commuting networks at the municipal scale from an overall perspective. Still and although good on average, the agreement between the three different datasets is broken in some particular connections that deviate from the diagonal in our scatterplots. This can be explained by the fact that the datasets come from different sources, were collected in different years and may have different biases and level of representativeness. For example, Twitter is supposed to be used more by younger people. The explanation of these deviations and whether they are just stochastic fluctuations or follow some rationale could be an interesting avenue for further research.

These results set a basis for the reliability of previous works basing their analysis on single datasets. Similarly, the door to extract conclusions from data coming from a single data source (due to convenience of facility of access) is open as long as the spatio-temporal scales tested here are respected.

# ACKNOWLEDGEMENTS

Partial financial support has been received from the Spanish Ministry of Economy (MINECO) and FEDER (EU) under projects MODASS (FIS2011-24785) and INTENSECOSYP (FIS2012-30634), and from the EU Commission through projects EUNOIA, LASAGNE and INSIGHT. ML acknowledges funding from the Conselleria d'Educació, Cultura i Universitats of the Government of the Balearic Islands and JJR from the Ramón y Cajal program of MINECO.

# APPENDIX

## Mobile phone data pre-processing

### Outliers detection

For both datasets we need to identify the outlier days to remove them from the data base. There are two types of outlier days, the special days (for example the National day) and the day for which we do not have the data for few hours. For example, for the metropolitan area of Barcelona, we can observe in Figure Sa eight days (from Monday to Monday) without outliers and in Figure Sb eight days with two outliers, Sunday, October 11$^{\mbox{th}}$ 2009 for which we do not have the data from 5PM to 11PM and Monday, October 12$^{\mbox{th}}$ 2009 the Spain's National Day.

### Voronoi cells

We remove the BTSs with zero mobile phone users and we compute the Voronoi cells associated with each BTSs of the metropolitan area (hereafter called MA). We remark in Figure Sa that there are four types of Voronoi cells:

1.  The Voronoi cells contained in the MA.

2.  The Voronoi cells between the MA and the territory outside the metropolitan area.

3.  The Voronoi cells between the MA and the sea (noted S).

4.  The Voronoi cells between the MA, the territory outside the metropolitan area and the sea.

To compute the number of users associated with the intersections between the Voronoi cells and the MA we have to take into account these different types of Voronoi cells. Let $m$ be the number of Voronoi cells, $N_{v}$ the number of mobile phone users in the Voronoi cell $v$ and $A_{v}$ the area of the Voronoi cell $v$, $v \in |[1,m]|$. The number of users $N_{v\cap MA}$ in the intersection between $v$ and MA is given by the following equation:

$$N_{v\cap MA}=N_v \left(\frac{\displaystyle A_{v\cap MA}}{A_v - A_{v\cap S}}\right)
    \label{vMA}$$

We note in Equation that we remove the intersection of the Voronoi area with the sea, indeed, we assume that the number of users calling from the sea are negligeable. Now we consider the number of mobile phone users $N_v$ and the associated area $A_v$ of the Voronoi cells intersecting the MA (Figure Sb).

## Origin-Destination matrices

As mentioned in the section *Extraction of commuting matrices* unlike the Twitter data we cannot directly extract an OD matrix between the grid cells with the mobile phone data because each users' home and work locations are identified by the Voronoi cells. Thus, we need a transition matrix $P$ to transform the BTS OD matrix $B$ into a grid OD matrix $G$.

Let $m$ be the number of Voronoi cells and $n$ be the number of grid cells. Let $B$ be the OD matrix between BTSs where $B_{ij}$ is the number of commuters between the BTS $i$ and the BTS $j$. To transform the matrix $B$ into an OD matrix between grid cells $G$ we define the transition matrix $P$ where $P_{ij}$ is the area of the intersection between the grid cell $i$ and the BTS $j$. Then we normalize $P$ by column in order to consider a proportion of the BTSs areas instead of an absolut value, thus we obtain a new matrix $\hat{P}$ (Equation S).

$$\hat{P}_{ij}=\frac{\displaystyle P_{ij}}{\sum_{k=1}^m \displaystyle P_{kj}}
    \label{pchap}$$

The OD matrix between the grid cells $G$ is given by a matrices multiplication given in the following equation:

$$G=P B P^t
    \label{OD3}$$


abstract: Here we describe the design and performance of the Spider instrument. Spider is a balloon-borne cosmic microwave background polarization imager that will map part of the sky at 90, 145, and 280 GHz with sub-degree resolution and high sensitivity. This paper discusses the general design principles of the instrument inserts, mechanical structures, optics, focal plane architecture, thermal architecture, and magnetic shielding of the TES sensors and SQUID multiplexer. We also describe the optical, noise, and magnetic shielding performance of the 145 GHz prototype instrument insert.
author: M. C. Runyan, P.A.R. Ade, M. Amiri, S. Benton, R. Bihary, J.J. Bock, J.R. Bond, J.A. Bonetti, S.A. Bryan, H.C. Chiang, C.R. Contaldi, B.P. Crill, O. Dore, D. O'Dea, M. Farhang, J.P. Filippini, L. Fissel, N. Gandilo, S.R. Golwala, J.E. Gudmundsson, M. Hasselfield, M. Halpern, G. Hilton, W. Holmes, V.V. Hristov, K.D. Irwin, W.C. Jones, C.L. Kuo, C.J. MacTavish, P.V. Mason, T.A. Morford, T.E. Montroy, C.B. Netterfield, A.S. Rahlin, C.D. Reintsema, J.E. Ruhl, M.C. Runyan, M.A. Schenker, J. Shariff, J.D. Soler, A. Trangsrud, R.S. Tucker, C. Tucker, and A. Turner Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, CA, USA;
School of Physics and Astronomy, Cardiff University, Cardiff, UK;
Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada;
Department of Physics, University of Toronto, Toronto, ON, Canada;
Department of Physics, Case Western Reserve University, Cleveland, OH, USA;
Jet Propulsion Laboratory, Pasadena, CA, USA;
Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, ON, Canada;
Department of Physics, Princeton University, Princeton, NJ, USA;
Department of Physics, Imperial College, University of London, London, UK;
National Institute of Standards and Technology, Boulder, CO, USA;
Department of Physics, Stanford University, Stanford, CA, USA;
Kavli Institute for Cosmology, University of Cambridge, Cambridge, UK
bibliography: mcrspie.bib
title: Design and performance of the Spider instrument

# INTRODUCTION

## The Spider Project

The Spider instrument in a balloon-borne millimeter-wave polarimeter designed to make very high sensitivity measurements of the polarization in the cosmic microwave background on mid- to large-angular scales with the goal of detecting primordial gravity waves from the inflationary epoch of the early universe. Spider will observe at 90, 145, and 280 GHz with beam sizes of $51'$, $31'$, and $17'$, respectively. An Antarctic flight of 24 days will allow Spider to map 8% of the sky to a depth of 0.25, 0.21, and 0.74 $\mu K_{CMB}$ per square degree at 90, 145, and 280 GHz, respectively. A more detailed description of the Spider project can be found in Filippini *et al.*, in these proceedings, as well as Crill *et al.*.

## The Spider Instrument Payload and Cryostat

At the heart of the Spider balloon payload is a large liquid helium cryostat. The cryostat houses up to six instrument inserts (described below) bolted to a $\sim1000\ell$ helium tank with all of the inserts pointed in the same direction. The cryostat cryogenics are briefly described in section below and a more thorough discussion can be found in Gudmundsson *et al.* in this volume. The cryostat measures 2.2 m in length and 2.0 m in diameter. The cryostat is mounted to a carbon fiber gondola via two pillow blocks that allow the cryostat to tilt in elevation. The gondola will be suspended from a 8 million cubic foot helium stratospheric balloon and is steerable in both elevation and azimuth. Power to the instrument is provided by solar arrays mounted to the back side of the payload and the instrument always points no closer than $45\deg$ to the sun.

# Instrument Inserts

## Introduction and Design Considerations

One of the biggest concerns for Spider is instrument payload mass. The science payload lifting capacity of the balloon is limited to $\sim5000$ pounds and the duration of our flight will depend, in part, upon payload mass. To that end we have tried to balance the constraint on mass with the desire to build a stiff structure that will not deflect under the periodic acceleration of the gondola scan pattern, nor break during launch (and hopefully landing). Materials with high thermal conductivity (such as copper) and good magnetic shielding properties (niobium, lead, and high-permeability materials) tend to be dense. So we have tried to minimize their use where ever possible. We had the benefit of starting the Spider insert design after much of the design of our sister experiment BICEP2 was completed. BICEP2 was designed for operation from the ground where weight is not as serious of a concern. We identified many areas where we could reduce the weight of the insert design with either a change in materials, strategic light weighting, or completely modifying a component (such as changing large metal tubes for carbon fiber trusses sheathed in light-weight copper-clad fiberglass).

## Optics Truss Structure and Cold Plate

A gold-plated $1/2''$ thick $\varnothing17.35''$ 1100-H14 aluminum cold plate forms the base of each Spider insert and is the interface to which the inserts mount to the cryostat helium tank. This aluminum plate has been light-weighted significantly. On top of this cold plate is mounted a ${}^3$He sorption fridge and the truss structure that supports the optics and focal plane. A gold-plated C10100 copper bus bar is mounted to the bottom side of the cold plate and the sorption fridge is bolted to this. The other side of this thermal bus is connected to copper tabs embedded in the liquid helium tank of the cryostat. A gold-plated copper bar penetrates the cold plate; the bottom side is connected to the pumped helium auxiliary tank of the cryostat and the other provides a 1.6 K point for the condensation point of the sorption fridge, a thermal intercept ring in the focal plane unit (FPU), a high-$\mu$ magnetic shield, and a cooled blackened light baffle, all of which are described below.

The focal plane unit is supported above the cold plate by a hexapod carbon fiber truss with $\varnothing1/2''$ rods called the camera truss. The two lenses are then supported from the top of the camera truss by two octopod carbon fiber trusses with $\varnothing3/8''$ rods. The individual carbon fiber legs have aluminum end caps epoxied to each end with Stycast 2850 and we have used fishing line to center the carbon fiber rods in the end caps to avoid possible galvanic corrosion from electrical contact. The mounting holes in the end caps have all been reamed to have a slip fit tolerance over $\varnothing1/4''$ shoulder bolts. The mounting pads on the feet that mount the rod assemblies to the aluminum rings are all coplanar for each leg. This allows the length of any truss to be changed simply by swapping out legs of different length. The insert volume in the single-insert test cryostat used for development and characterization is longer than that in the flight cryostat and we will compensate by shrinking the length of each camera truss by $4''$ before installing the inserts into the flight cryostat.

All of the carbon fiber trusses are wrapped in thin copper-clad G-10 sheets to form a light and RF tight sleeve around the inserts. We have the option of coating the inside of these wraps with a flexible blackening material made from silicone RTV, carbon lampblack, and 316 stainless powder. We measured the mass loss during cure and heating of many commercially available silicones and found Dow Corning 748 to have very low mass loss. We have run our cryostat with $\sim1.5$ m${}^2$ of this material and had no issues with out gassing. Outside of the wraps run four C10100 copper bars to cool the optics trusses. Each insert is capped with a 1100 aluminum snout that holds a set of filters in front of the objective lens. These filters are described in section below.

## Optics

### Design Description

Each Spider insert is a cooled two-lens refracting telescope. The lens design is identical to that of the BICEP2 experiment and is described in detail in Aikin *et al.* (in these proceedings). Both sides of the two lenses are simple conics and are separated by 550 mm with an effective focal length of 583.5 mm. This yields a plate scale of 0.98 deg/cm on the focal plane. The optical system is telecentric to accommodate the flat focal plane geometry. The field of view of each insert is $20\deg$ across the diagonal. The entire telescope optics are cooled to 4.2 K to reduce the in-band loading from loss in the thick dielectric lenses.

### Lenses, AR Coats, and Flexures

The lenses used in Spider are machined from $2''$ thick cast HDPE slabs. The lens fabrication process involves a number of steps designed to reduce the stress in the lens and achieve the desired lens figure. Before machining, the HDPE slabs are annealed in a temperature-controlled oven between metal plates and then rough cut to shape. The rough cut lenses are then clamped in an annealing jig and re-annealed before finish machining on a CNC mill. The finished lenses are then measured on a CMM to assure the lens figure and surface finish have been archived. The maximum allowable deviation in the lens figure is $0.005''$ and the lens surface finish is $<63\mu ''$. The final lenses are $12.4''$ in diameter with an optically-active diameter of $11.4''$.

Because the index of refraction of HDPE in the millimeter is $\sim1.52$ , the instrument would suffer reflections of order $[(n_1-n_2)/(n_1+n_2)]^2 \sim 4\%$ at each lens surface. To reduce these reflections we bond a tuned layer of porous PTFE sheet to each lens surface manufactured by Porex. We have chosen a material with a density that yields an index of refraction close to the ideal of $n_{AR} = \sqrt{n_{HDPE}}$ and we order the sheets in thicknesses of $\lambda/4n_{AR}$. The toleracing of the AR material from the manufacturer is only good to a few mils. But modeling of the transimission of HDPE (and nylon) slabs with AR coats of various thickness and indicies of refraction shows that $>99\%$ transmission is achieved with thickness and density variations within manufacturing tolerances.

The plastic lenses contract during cooling substantially more than the aluminum mounts from which they are supported ($\sim 2\%$ for HDPE and only $\sim0.4\%$ for aluminum). The radial differential contraction between the lens and support ring is $\sim0.1''$ at the edge. To allow for this contraction while keeping the lens well centered and rigid we use eight $1/32''$ thick copper flexures spaced equally around the lenes. The copper flexures also provide the cooling path for the lenses. However, the cooldown time of the lenses is limited by the internal thermal time constant of the HDPE. The length of the flexures are tuned so that the spacing of the lenses is correct after thermal contraction as the insert cools. Most of the support structure is carbon fiber, so the contraction in the lenses themselves dominates this correction.

### Optical Filtering

Spider employs multiple filter types to reduce IR loading on the various thermal stages. We use a stack of four metal-patterned mylar "IR shaders" on both of the vapor cooled shield (VCS) stages to reduce the IR loading on the helium stage . We also use hot-pressed resonant metal-mesh filters on the cold side of the 20 K VCS stage and the entrance to the insert, as well as just above the focal plane at 1.6 K. The entrance to the insert has an AR-coated sheet of $3/32''$ nylon on the cold side to help absorb additional IR.

### Baffling and Optical Stop

The space between the two lenses is filled with a blackened sleeve that is capped with an annulus that defines the optical stop of the system $9.5''$ in diameter just behind the objective lens. The stop limits the spillover of the beam onto warmer stages upstream along the optical chain, including the waveplate, filters, and window. The square phased-array antennas produce a 2D sinc beam pattern and $\sim 25\%$ of the beam power falls outside of the stop. The blackened sleeve absorbs this sidelobe power and prevents it from escaping out the front of the cryostat. Because a significant fraction of power is incident upon this blackened sleeve, we susped is from the front of the optics tube using a carbon fiber truss and cool the entire tube to 1.6 K using the pumped helium auxilliary tank (described below). The sleeve is formed by soldering a tube of copper-clad G-10 and then lining the sleeve with a mixture of Stycast 2850, carbon lampblack, and 316 stainless powder.

### Waveplate

To increase polarized angle coverage and mitigate the effect of beam asymmetries, polarization modulation in Spider is achieved via an AR-coated sapphire half-wave plate (HWP) mounted at 4 K skyward of the primary optic of each insert . The HWP will be periodically rotated to different angles throughout the flight with a cryogenic stepper motor. Each HWP will be optimized for the single frequency band of the 90, 145, or 280 GHz insert in which it is mounted. A fused-quartz AR coat will be applied to each side. We have measured the millimeter-wave transmission spectra of birefringent sapphire at room and liquid helium temperatures in the lab, which are consistent with our physical optics model. We have also taken preliminary spectra of an AR-coated 145 GHz HWP integrated with the Spider optics and detector system in the prototype Spider receiver. Preliminary results show good performance.

## Focal Plane

The development of the Spider focal plane was a combined effort with the BICEP2 project. Orlando *et al.* (in these proceedings) describes that effort to develop these platter-style focal planes in detail (these cover revisions A through D). During the course of testing we found that we could not adequately shield the system from magnetic fields at a level where we would not have to worry about signals induced by spinning the cryostat in earth's field. This was a fundamental limitation of the focal plane architecture that places the detector tiles and SQUID multiplexer chips on the same plane without the ability to surround the SQUIDs with magnetic shielding. We have opted to redesign the focal plane and move all of the SQUIDs inside of a closed superconducting niobium box with secondary high-permeability and superconducting shields within. This focal plane architecture is described in detail below and is referred to as RevX.

### Architecture

The Spider RevX focal plane architecture is shown in figure . The detector tiles are mounted onto a square gold-plated C10100 copper detector plate $8.55''$ across. The focal plane structure is cooled to 300 mK with the ${}^3$He sorption fridge. Pins in the detector plate and slots in the silicon detector tiles register the detector tiles and allow for differential thermal contraction between the silicon and copper. The detector tiles are held in place with beryllium copper clips. The detectors are patterned onto the back side of the silicon tiles with respect to the incident light and NSG-N quartz anti-reflection tiles are placed on the front side of the detector tiles to minimize surface reflection from the silicon. The detectors are spaced $\lambda/4$ from a niobium backshort plate. The spacing between the detectors and backshort is set using lapped Macor ceramic spacers around the perimeter of the focal plane with a single spacer in the center. Custom made niobium fasteners attach the copper detector plate to the niobium backshort and the holes in the copper plate are over sized to allow for differential thermal contraction between the copper and niobium.

Behind the niobium backshort is a $1/32''$ thick welded niobium box that encloses all of the SQUIDs. The signals from the detector tiles are brought around from the detector plate to the inside of this box using flexible aluminum superconducting circuits (flexi circuit). The aluminum flexi circuit has a measured superconducting transition temperature of 775 mK. The other side of the flexi circuits are connected to the input stages of the SQUID multiplexing system, which is described below. Each detector tile has an individual flexi circuit that can carry 128 TES signals. The flexi circuits are made in two parts which overlap at the center to minimize the size of the slot in the niobium box required for their passage; the effectiveness of the niobium box is limited by the width of these slots. The niobium box is fastened to the niobium backshort with $40\times$ 7075 aluminum fasteners with a measured superconducting transition temperature of 900 mK. The total mass of the 300 mK focal plane structure is 8 kg.

### Detectors and band definition

The antennas and detectors used in Spider are described in Kuo *et al.* as well as Orlando *et al.* (in these proceedings). A focal plane consists of four detector tiles, each of which has an 8$\times$`<!-- -->`{=html}8 array of dual-polarization pixels (128 detectors) for 145 and 280 GHz, and a 6$\times$`<!-- -->`{=html}6 array (72 detectors) at 90 GHz. Each pixel consists of two arrays of polarized slot antennas (called A and B and are rotated 90 degrees w.r.t. each other), each of which are summed in phase on superconducting microstrip. This microstrip passes through a resonant filter that defines the millimeter-wave frequency band and then the signal is fed into a resistive gold meander located on a suspended bolometer island (see figure ). Also on each island are both a titanium and aluminum transition edge sensor (TES) in series. The millimeter power dissipated in the gold meander heats up the bolometer island and the voltage-biased TES responds to that change in temperature with a change in resistance, thus producing a measurable change in current.

The thermal conductivity ($G$) of the Spider bolometers is low (20 pW/K at $T_C$) to take advantage of the decreased optical loading from a stratospheric balloon platform. The $G$ of the bolometers are tuned to put the titanium TES ($T_C\sim500$ mK) on transition at float with a margin of safety of $\sim3$. These low-$G$ devices are saturated on the titanium transition when doing testing in the laboratory, so the aluminum TES ($T_C\sim1.3$ K) in series is used. This allows characterization of beams and frequency bands in the lab before flight. The TES bolometer islands have some response to high-frequency "blue leaks" since the microstrip running to the island acts like an antenna within the cutout of the niobium ground plane. In an effort to push that frequency response above our filter cutoff we have shrunk the window cutout of the islands. The lower $G$s of the Spider devices required that the island legs that determine the device $G$ be meandered to obtain enough length for the required thermal conductivity.

### SQUID multiplexer and cabling

At the heart of Spider's signal readout chain is a three-stage SQUID time-domain multiplexer made by NIST . SQUIDs are remarkably sensitive to changes in magnetic flux, so coupling an input coil into a SQUID loop forms the basis of a very sensitive current readout. Each bolometer signal is fed through a low-pass $L/R$ filter that reduces aliased noise and then into the input coil of a clover leaf SQUID (SQ1). A single mux chip handles the signals from 32 TESs and a "dark" SQUID (used to monitor noise and magnetic field pickup) and combines these signals into a summing coil that is fed into a second stage SQUID (SQ2). The individual SQ1 are biased in sequence, so the summing coil (and hence, SQ2) only has the signal from one of the TESs at any time. The SQ1 bias lines are referred to as "row selects" because they select which SQ1 is active. The SQ2 signal is then fed into a SQUID series array amplifier (SSA) that has 100 SQUIDs in series per channel. The timing, biasing, and readout of the SQUIDs is controlled by the Multi-Channel Electronics (MCE) crate made by UBC . Each MCE crate handles the 512 detector signals from an entire insert with only a power cable and fiber optic pair leaving the crate.

As mentioned above, the signals from the detectors are routed to the SQUID chips on superconducting flexible circuits. The multiplexer chips are mounted on circuit boards with copper traces onto which superconducting PbSn solder has been electroplated between the mux chips and the bond edge. Each of these "MUX boards" has four sets of SQUID chips which are enough to handle the signals from one detector tile (up to 128 devices). The typical impedance of the TESs on the titanium transition is $\sim30-40$ m$\Omega$ and the TES bias shunt resistors are only 3 m$\Omega$. The superconducting aluminum and PbSn circuit between the TESs and mux chips reduce the parasitic impedance along this path to a level well below the shunt resistance.

The signals from the SQ2s, SQ1 row select lines, detector biases, and SQUID feedback lines are feed to Samtec headers at the edge of the MUX boards and the four MUX boards plug into one breakout board which routes the signals to four 100-way micro-D connectors. All of the outputs of the SQ2s are run through one of these 100-way connectors up to a separate circuit board that houses the two SSA modules. The other three 100-way connectors have all of the signals, biases, and feedback lines and run to the cold plate on NbTi ribbon cables. There are 300 wires leaving the focal plane and the detector bias and SQ2 bias lines can carry currents up to a couple mA each. The joule power dissipation by manganin cables would have been too much for our ${}^3$He fridge to handle. So we have chosen to make all of the wires running to the focal plane superconducting (with the exception of some thermometry). The cables running from the 4 K cold plate to the 300 K vacuum feedthrough to the MCE crate are all teflon-jacketed shielded twisted-pair manganin.

### Magnetic shielding

A serious concern to the Spider project is magnetic field pickup in the system that cannot be separated from astrophysical signals. One of the strengths of the Spider instrument is the ability to make long scans across the sky and measure signals on large angular scales. The earth's magnetic field forms a dipole that we will scan through and magnetic pickup in the various SQUID stages and the TESs themselves will be difficult to discriminate from the desired signal. Changing magnetic flux through the SQUIDs produces a shift in the $V-\phi$ curve. To reduce this effect, NIST has fabricated the SQUIDs in a counter-wound cloverleaf pattern so that they are nominally only sensitive to 2nd order gradients in magnetic field. Changing magnetic flux will also induce a current in the superconducting summing loop that combines the signals from all of the SQ1s and feeds it into SQ2. This loop area has been minimized but is still finite. The exact value of $T_C$ of the TESs depends upon the magnetic field and so a change in magnetic environment will shift $R(T)$ and produce a change in current.

To that end we have tried to shield the focal plane from magnetic fields so that pickup from the earth's field will be below the expect signal from the CMB. Our shielding scheme consists of multiple layers of both high-permeability and superconducting materials. The TESs are spaced $\lambda/4$ away from the large superconducting niobium backshort and the entire focal plane is surrounded by a large high-permeability "spittoon" made of Amuneal A4K material to reduce changes in $T_C$ of the TESs. As mentioned above, the TES signals are routed inside of a superconducting niobium box that is closed with superconducting aluminum fasteners. We have tried to minimize the size of all slots penetrating the niobium box by overlapping the flexi circuits to halve their width as well as truncate the two thermal straps where they enter the box. Immediately inside the niobium box at the level of the slots is another high-permeability A4K plate that wicks entering magnetic field away before it can get to the SQUIDs. There is a reentrant superconducting 1100 aluminum open-top box that the flexi circuits run up and over where the SQUID mux chips are located. All of the SQ1 and SQ2 chips (which include the superconducting summing coils) are further shielded with squat, open-ended sleeves that surround the mux boards (each has four sets of mux chips). The SSA boxes are suspended in the aluminum box and are housed in compact open-ended niobium boxes with Cryoperm sheets inside. These niobium boxes are further wrapped in ten layers of 0.6 mil Metglas 2714A. Outside of the insert we have wrapped a $20''$ long, $0.006''$ thick superconducting lead sleeve that is centered on the focal plane structure. Lastly, each of the insert tubes is lined with a two-layer $0.040''$ thick A4K shield that runs the length of the helium tank.

We have modeled the magnetic shielding configuration using COMSOL Multiphysics. This model does not have the outer-most Cryoperm-10 shield to keep the mesh size manageable. The results of the model for a 1 T incident magnetic field are shown in figure where we display the $log_{10}$ of the resulting field amplitude. COMSOL does not handle multiply-connected superconducting regions properly. So we model the effects of the narrow $0.10''\times1.35''$ slots in the niobium box through which the flexi circuits pass one at a time. The direction of the magnetic field is into the page which is along the long axis of the slots (the direction of maximum magnetic field penetration). This model is idealized but suggests that the shielding configuration of the RevX focal plane is capable of yielding shielding factors in excess of $10^8$ at the SQUIDs, not including the additional A4K tubes that run the length of the helium tank and are expected to provide an additional shielding factor of $\sim50$. We have estimated the shielding factor to earth's magnetic field necessary for the flight to be $\sim10^7$.

## Cryogenics

### Overview

The Spider Long Duration Balloon (LDB) cryostat and cryogenic systems are described in detail in Gudmundsson *et al.* (in these proceedings). The cryostat holds $\sim1000\ell$ of liquid helium in a tank shaped like a 6-shot revolver cylinder. There are seven holes in the tank; the six around the outside house the six instrument inserts and the central hole is used to pass wiring up to the waveplates and thermometry in front of the inserts. The six instrument cold plates (described above) bolt to the liquid helium tank via gold-plated aluminum interface plates. In addition to the main helium tank, a smaller $20\ell$ auxiliary tank (Aux tank) of helium is fed via a capillary tube from the main tank and is pumped to 1.6 K. On the ground, this pumping is achieved with a vacuum pump but in flight the reduced atmospheric pressure at float provides the pumping. As described below, the Aux tank provides an intermediate temperature stage between our sub-K fridge and the main liquid helium bath.

### Truss structure and thermal management

The focal plane unit (FPU) is supported on eight 316 stainless "heat capacity blocks". The cooling path from the FPU to the ${}^3$He fridge still runs through these blocks and we have measured them to have a combined thermal impedance of 2.3 mK/$\mu$W, in excellent agreement with COMSOL simulations. They are gold plated on the ends to minimize thermal boundary impedance. Stainless steel has a relatively high heat capacity and these blocks act as a passive thermal filter to temperature fluctuations generated below in the truss structure, thermal straps, or the fridge. We have measured the thermal transfer function of the passive thermal filter to have a 3dB point of 2 mHz.

The heat capacity blocks connect the four sides of the FPU to a copper ring via copper supports. This whole structure is cooled to 300 mK and is referred to as the sub-Kelvin stage (sub-K stage). The sub-K stage is supported off of another copper ring cooled to 1.6 K with the pumped helium auxiliary tank with carbon fiber rods. This "aux ring" is itself supported from a 4.2 K aluminum ring with carbon fiber rods. We measured the thermal conductivity of many polymeric and composite materials during the course of the development of the Spider truss structure and found carbon fiber rods to have the highest ratio of elastic modulus to thermal conductivity of non-brittle materials over the range of temperatures of interest . The truss structure appears gossamer but is very stiff. We have modeled the deflection with COMSOL and expect the 8 kg focal plane to deflect a few mils over a $90\deg$ tilt, corresponding to much less than an arcminute of beam deflection on the sky.

There is a concern that the vibration of rigid metal thermal straps attached to the cold stages might vibrate at frequencies within the science band and couple power into the detectors. The stainless heat capacity blocks described above are part of the strategy to mitigate that effect. We have opted instead to use flexible thermal straps formed with many layers of $0.001''$ C11000 copper foil that have been electron beam welded into copper mounting tabs. Despite the foil alloy being C11000, we have measured its RRR to be $\sim200$ and confirmed that they have correspondingly high thermal conductivity.

We use a combination of DT-series silicon diodes and Cernox resistance thermometers (both made by Lakeshore Cryotronics) to monitor the temperatures within the cryostat and inserts. We exclusively use resistance thermometers at all sub-K stages and diodes for all stages with temperatures greater than 1 K.

# Spider System Performance

## Optical Performance

Spider will observe in three frequency bands at 90, 145, and 280 GHz. Most of the insert development effort has focused on the 145 GHz system, and consequently, it is the most well characterized of our frequencies. Here we present some of the optical performance results from the 145 GHz system.

### Bands

As mentioned above, Spider's frequency bands are defined with resonant on-tile filters between the detector antenna and the resistive meander on the TES island . We have measured the spectal response of one of the inserts configured for 145 GHz operation. The average bandwidth for the 145 GHz channels is 34 GHz (using the convention of Runyan *et al.*), or 25%. Even at float we are concerned about residual atmospheric emission and have worked to fit the frequency bands within the available atmospheric windows, paying particular attention to avoiding water lines.

Spider's observing strategy will involve differencing the two polarizations in a given pixel and a mismatch in spectra could limit the effectiveness of this strategy. Figure shows examples of measured spectral bands at 145 GHz. In the figure, the spectral bands have been normalized to have the same integrated value over the band. The oscillatory nature of the difference band suggests that common mode emission sources with smoothly-varying spectra will be well-differenced. The Spider filter strategy also effectively removes high-frequency spectral leaks that could couple high-frequency power onto the detectors. We have measured the high-frequency response of our 145 GHz band using a chopped thermal source and a thick-grill filter with a cut off of 185 GHz and found the response to be less than a few tenths of a percent of the in-band signal (consistent with the noise floor of the measurement).

### Optical efficiency

We needed a way to characterize the Spider instrument inserts under flight-like loading conditions. To do this we developed a helium cold load cryostat that bolts to the front of either the insert test cryostat or the flight cryostat and presents a cold, millimeter-black, beam-filling surface to the instrument. This cold load cryostat is a $23\ell$ helium cryostat with a $15.5''$ cold plate. Using a carbon fiber truss, we stand off a copper plate covered in millimeter-black pyramidal tiles manufactured by Thomas Keating, Ltd. Embedded within these tiles are multiple Lakeshore SD diodes to monitor the load temperature. We have mounted a resistive heater on the back side of the black load and can elevate it to any temperature above 5 K. The cold load cryostat has radiation shields attached to the helium bath and its vapor cooled shield and we have mounted IR-blocking filters to each of these stages to reduce the optical loading on the black absorber.

A convenient feature of being able to elevate the cold load temperature is that we can measure the end-to-end optical efficiency of our system by measuring the optical loading on the detectors, $P_{opt}$, as a function of cold load temperature, $T_{CL}$. The only optical element missing in this configuration is the cryostat vacuum window (the cold load and main cryostats share a common vacuum) and we have added a handful of IR filters on the cold load side which will not be there in flight. We ramp the cold load temperature from 5 K to $\sim20$ K and take load curves to measure the loading at each temperature. We convert the actual cold load temperature into an equivalent Rayleigh-Jeans temperature, $T_{RJ}$, and then take the slope $dP_{opt}/dT_{RJ} = k\int\eta_\nu~d\nu$, where $k$ is Boltzmann's constant, $\eta_\nu$ is the measured spectral response normalized by the optical efficiency, and we have assumed single-moded performance ($A\Omega=\lambda^2$). We measure $dP_{opt}/dT_{RJ}$ to be $\sim0.17$ pW/K$_{RJ}$ at 145 GHz. We can then calculate a band-average optical efficiency $\overline{\eta} = (\int{\eta_{\nu}~d\nu})/\Delta\nu$. The 145 GHz Spider inserts have a typical end-to-end optical efficiency of 36%.

### Internal loading

To fully take advantage of the stratospheric balloon platform we need to reduce the internal loading to a level where the additional photon noise contribution does not significantly increase the overall system noise. We have estimated the loading from the CMB and atmosphere at float to be $\sim0.28$ pW per polarization at 145 GHz. This corresponds to a RJ temperature of $\sim1.7$ K for a system with 36% optical efficiency and a 34 GHz bandwidth. The loading from a 300 K vacuum window with 0.5% loss is $\sim0.2$ pW. So we would like the loading from within the cryostat to be a small fraction of a pW, or only a few K${}_{RJ}$.

The cold load provides a convenient way to measure the internal loading from within the cryostat. We extrapolate the $P_{opt}$ versus $T_{RJ}$ curve described in the previous section back to zero and read off the internal loading as the y-intercept. Note that this technique does not include the loading from the warm vacuum window, which is estimated to be $\sim0.25$ pW. Before installing the cooled black optics sleeve between the lenses we measured the internal loading from within the cryostat to be $\sim0.6$ pW, or $\sim3.5$ K${}_{RJ}$ at 145 GHz. We have since installed the cooled optics sleeve and by measuring the coupling of the detectors to the sleeve we expect the internal loading to decrease by 0.3 to 0.4 pW.

### Beams

We have also done some preliminary beam characterization at 145 GHz by looking into the laboratory through a $6''$ thick Zotefoam PPA-30 window (this will not be the flight window configuration). We have made near-field beam maps a few inches away from the window aperture using a chopped thermal source on an X-Y translation stage as well as far-field slices using a large chopped thermal source on a linear stage that can be rotated about the bore sight. The near field beam maps are in qualitative agreement with physical optics modeling using Zemax and the far-field beam slices for the few pixels measured show well-matched gaussian beam pairs for the two polarization antennae within a pixel with a FWHM of $31'$ at 145 GHz. The difference in the A and B polarizations within an individual pixel are measured to be $<2\%$ of the average beam. See Aikin *et al.* in these proceedings for more details on beams from the BICEP2 optical system, which is very similar to Spider's.

## Detector performance

We have measured detector parameters for a handful of prototype detector tiles. The Spider-specific engineering tiles have two different values of $G$ so that we can determine the optimum value for the flight.

### Noise

We have measured the instrument noise at 145 GHz versus $R/R_n$ on the TES transition while staring into the cold load cryostat at 5.5 K to simulate flight loading. We calibrate the noise using bolometer load curves and the measured optical efficiency of the system, $dP/dT_{RJ}$, and then convert the noise PSDs into $NET_{CMB}$. We also confirm this calibration by varying the temperature of the blackened cold load plate and measure the response on the detectors. Some example noise spectra from an A/B pixel pair (as well as the difference PSD) is shown in figure .

The angular scales on which Spider will be able to measure the CMB power spectrum will be affected by the combination of scan speed and $1/f$-noise. We have measured the $1/f$-knee frequency of undifferenced bolometers to be typically 0.1 to 0.4 Hz. Taking the difference of the two detectors in an A/B pixel pair removes the common optical and electrical components and lowers the $1/f$ frequency to $\sim60$ mHz. This corresponds to $\ell$ of a few for the likely scan speeds of Spider. As seen in the figure, pixel pair differencing produces three decades of clean signal bandwidth.

## Magnetic shielding performance

As mentioned above, the response of the system to magnetic fields is a significant concern to Spider. To measure the magnetic field response of the system we use two Helmholtz coils 79 cm in diameter. These coils can be spaced in true Helmholtz configuration for the direction along the optic axis of the test cryostat, but must be spaced further than $D/2$ for the two orthogonal axes due to the diameter of the cryostat. In all configurations we center the coils about the focal plane and calculate the nominal field strength at the center in the absence of magnetic materials. The coils are driven with a AE Techron current amplifier with a sine wave frequency generator and can produce a field strength of many 10s times earth's field.

The principle motivation for converting from the flat focal plane architecture of RevA-C to the shielded box architecture of RevX was to reduce the magnetic field pickup in the system. We have reduced the magnetic field response by at least an order of magnitude for most channels by converting to the RevX design. In fact, for the majority of channels we can only place upper limits on the magnetic response.

We have also been unable to detect significant magnetic response in the TESs themselves due to shifting $T_c$. Although we haven't seen response of the detectors themselves due to changing magnetic fields, the response of the SQUID multiplexing system will appear as a signal in the detector time streams. We can calibrate the response of the system to magnetic fields into K${}_{CMB}$ per earth's field ($B_e$) and find that we can place a limit of the magnetic field response of $<10~\mu K_{CMB}/B_e$ in the majority of channels at 145 GHz for all three axes without any attempt to remove pickup (*e.g.* with A/B pixel differencing or deprojection of the non-antenna coupled SQUID in each MUX column). In previous focal plane shielding iterations we found that pixel pair differencing reduces the magnetic signal by 1 to 2 orders of magnitude. Even at this limit, the signal from earth's field would be much less than the CMB dipole but would still be detectable in Spider's CMB maps. We hope to do longer integrations at high magnetic fields to put even tighter limits on the magnetic response.

The Spider collaboration gratefully acknowledges the support of NASA (grant number NNX07AL64G), the Gordon and Betty Moore Foundation, and NSERC. With great sadness, the Spider collaboration acknowledges the countless contributions of Andrew E. Lange, the late PI of the Spider project. His wisdom and selfless leadership will be sorely missed. WCJ acknowledges the support of the Alfred P. Sloan Foundation. The authors gratefully acknowledge our collaboration with the BICEP2 and Keck projects. The author thanks J. Lazear for his help with the design and construction of the cold load.


abstract: *Characterizing atmospheres beyond our Solar System is now within our reach.
  Kevin Heng received his education in astrophysics (M.S. and Ph.D) at JILA (the Joint Institute for Laboratory Astrophysics) and the University of Colorado at Boulder. Subsequently, he was a postdoctoral researcher at the Institute for Advanced Study, Princeton, from 2007 to 2010 (including holding the Frank & Peggy Taplin Membership from 2009 to 2010). He is currently a Zwicky Prize Fellow at the Institute for Astronomy at ETH Zürich (the Swiss Federal Institute of Technology) in the Star and Planet Formation Group, where he is involved in the Exoplanet Characterization Observatory (EChO) mission proposed to the European Space Agency. He has worked on several topics in astrophysics, including shocks, planetesimal disks and fluid dynamics. His current, main research interest is in developing a hierarchy of theoretical tools to understand the basic physics and chemistry of exoplanetary atmospheres from the perspective of an astrophysicist. He spends a fair amount of time humbly learning the lessons gleaned from studying the Earth and Solar System planets, as related to him by atmospheric, climate and planetary scientists. He received a Sigma Xi Grant-in-Aid of Research in 2006.
  Text-only version of article, edited by Fenella Saunders. Full version is available at: `www.americanscientist.org `*
address: ETH Zürich, Institute for Astronomy
  Wolfgang-Pauli-Strasse 27, CH-8093, Zürich, Switzerland
author: Kevin Heng
title: The Study of Climate on Alien Worlds

It is a distracting, inconvenient coincidence that we are living in times of paradigm-shifting astronomical discoveries overshadowed by the deepest financial crisis since the Great Depression. Amid a battery of budget cuts, the astronomical community has discovered more planets outside of our Solar System—called extrasolar planets or simply exoplanets—in the past decade than in previous millennia. In the last couple of years alone, the Kepler Space Telescope has located more than 2,000 exoplanet candidates, including Earth-sized ones potentially capable of sustaining liquid water, demonstrating the ease at which nature seems to form them and hinting that we may be uncovering the tip of an iceberg. Discovering and characterizing distant, alien worlds is an endeavor no longer confined to the realm of science fiction.

In tandem with numerous surveys of the night sky performed from the ground, the Hubble, Kepler and Spitzer Space Telescopes observe the universe from outside of Earth's atmosphere. These devices detect an exoplanet by recording the diminution of light as the body, residing in an edge-on orbit, passes in front of its host star. In the past few years, astronomers also have achieved the remarkable feat of measuring the diminution of light as the exoplanet passes behind its star, known as the secondary eclipse. In other words, astronomical techniques have advanced to the point where we can detect a star masking the light from its exoplanet, which is a demonstrably small effect—at most a few parts in a thousand in the infrared and much smaller in the optical range of wavelengths. During a secondary eclipse, the light from an exoplanetary system originates only from the star, and these data can be used to subtract out the starlight when the exoplanet is not eclipsed. All that remains is the light of the exoplanet and its atmosphere (if it exists). Such a technique has enabled astronomers to make the first detections of the light directly emitted by an exoplanet, which typically appears at its brightest in the infrared.

Measuring transits and eclipses at several different wavelengths allows one to construct a spectrum of the exoplanetary atmosphere, of which a spectral analysis yields its composition and elemental abundances. (A spectrum describes the range of colors of the photons emanating from the exoplanet, but it generally extends beyond what our eyes can see toward both shorter and longer wavelengths.) In some cases, astronomers were able to record the ebb and rise of the brightness of the exoplanet as it orbits its parent star, otherwise known as the phase curve. An inversion technique, developed by Nick Cowan of Northwestern University and Eric Agol of the University of Washington, allows one to convert the phase curve into a "brightness map," which is the latitudinally averaged brightness of the exoplanet across longitude. Recent work by the same researchers has yielded two-dimensional information on the brightness of the exoplanet HD 189733b as a function of both latitude and longitude. In other words, we have started to do cartography on exoplanets!

**Tidal Locks**

The first studies of exoplanetary atmospheres were performed on a class of objects known as hot Jupiters. A combination of the transit technique with a measurement of the radial velocity (which is the gravitational wobble of a star as its exoplanet orbits around their common center of mass) yields the radius and mass of a hot Jupiter, respectively, and reveals that they are similar in these aspects to our own Jupiter. The startling difference is that hot Jupiters are found about a hundred times closer to their parent stars than Jupiter, which raises their surface temperatures to between 1,000 and 3,000 degrees Kelvin. With spatial separations of a hundredth to a tenth of an astronomical unit (the average distance from the Earth to the Sun) from their stars, the discovery of hot Jupiters caught the astronomical community by surprise, because their existence was neither predicted from astrophysical theory nor subsequently explained by it.

Their large sizes render hot Jupiters easier to observe and thus the most obvious laboratories for extrasolar atmospheric studies. Furthermore, the belief that their atmospheres are dominated by molecular hydrogen—which is consistent with the densities of the exoplanets, inferred from the astronomical observations to be about 1 gram per cubic centimeter—offers some hope that the atmospheres are primary, reflecting the composition of the primordial nebulae from which they formed, rather than secondary and reprocessed by geological mechanisms (such as on Earth).

Given enough time, an exoplanet's position and rotation tend to relax toward a state of minimum energy—a spin synchronized state, such that one hemisphere of the exoplanet always faces its parent star with the other hemisphere shrouded in perpetual darkness. The characteristic time scale associated with this process is typically 1,000 times less than the age of the star. (As a more familiar example, the Moon is in a spin synchronized state with respect to the Earth, notwithstanding its tiny rotational corrections called librations.) In other words, one hot Jovian day is equal to one hot Jovian year. The unfamiliar configuration of permanent day- and night-side hemispheres on hot Jupiters opens up an unexplored regime of atmospheric circulation with no precedent in the Solar System and motivates theoreticians to test their tools in unfamiliar territory.

Understanding these hot Jovian atmospheres requires clarifying the complex interplay between irradiation, atmospheric dynamics, chemistry and possibly magnetic fields. On the most irradiated hot Jupiters, the exoplanet viewed from the poles resembles a sphere painted half white and half black—the phase curve is a sinusoidal function that peaks at secondary eclipse and becomes dimmest at transit. Any shift of this peak from its reference point at secondary eclipse may be interpreted as being due to the presence of horizontal winds in the atmosphere, which act to transport heat from the day- to the night-side hemisphere. This angular shift was first measured for an exoplanet, the hot Jupiter HD 189733b, by Heather Knutson of the California Institute of Technology and her collaborators, who reported a peak shift of about 30 degrees east—in the direction of rotation. This angular shift was also measured for the hot Jupiters Ups And b (by Ian Crossfield of the University of California at Los Angeles and his collaborators) and WASP-12b (by Cowan and his collaborators).

Other astronomers continue to push the envelope. Ignas Snellen of Leiden University and his colleagues, using the ground-based European Very Large Telescope (VLT), used a technique called absorption spectroscopy to measure the speed of the horizontal winds on the hot Jupiter HD 209458b. The technique compares the relative size of the exoplanet across a range of wavelengths. At a wavelength where an atmospheric atom or molecule is the most absorbent, the exoplanet appears larger. By monitoring the shift in wavelength of an absorption line of carbon monoxide, the group determined that HD 209458b's winds clock in at about 2 kilometers per second, roughly 100 times faster than those on Earth. More attempts to measure atmospheric wind speeds are in the works, and these measurements remain at the cutting—if not the bleeding—edge of what astronomers can achieve.

**New Languages**

The importance of these discoveries to astronomy cannot be overstated—they signal the dawn of exoplanetary meteorology, or at least legitimize its study in the eyes of astronomers and astrophysicists.

Astronomers now possess a tool kit not only to measure the masses and sizes of exoplanets but also to characterize their atmospheric dynamics and chemistry. Besides galvanizing the astronomical community, this newfound field is starting to exert a profound sociological impact on related fields of study: atmospheric and climate science, geophysics and planetary science. It marks the first great confluence of these fields with astrophysics, a gathering of scientists with different scientific and modeling philosophies, which is especially evident at interdisciplinary conferences where we struggle to understand one another's jargon. Atmospheric and climate scientists, as well as geophysicists, are firmly grounded in a data-rich regime, living within the system they study. Awash in an abundance of data from the terrestrial atmosphere and the geological record, no single model is capable of accounting for all of the observed phenomena. Instead, a hierarchy of models with different degrees of sophistication is utilized, with each model isolating key pieces of physics. The strategy is to first divide and conquer, then to unify and rule.

The knowledge gleaned from studying Earth and the Solar System planets serves as an invaluable guide, but there is a cautionary tale to be told. As a rule of thumb, there are two characteristic length scales describing an atmosphere: the Rhines length is the typical width of zonal (east-west) jets, whereas the Rossby length is the typical size of vortices or eddies. For Solar System objects, both length scales are much smaller than planetary radii. On close-in exoplanets, the Rhines and Rossby lengths are comparable to exoplanetary radii, implying that the atmospheric features are global in extent, an expectation that is borne out in three-dimensional simulations. The atmospheres of close-in exoplanets are thus in a circulation regime that is unprecedented in the Solar System. Atmospheric circulation simulations therefore have to be global instead of local, and other physical implications—such as the mixing of atmospheric constituents and its effect on the spectral appearance of the exoplanet—remain to be fully understood.

The study of exoplanets is essentially confined to the scrutiny of point sources in the night sky. Although we may obtain detailed spectral and temporal information on these point sources, the procurement of detailed spatial information remains a grand challenge for posterity. Planetary scientists benefit from the ability to obtain photographs of the Martian surface and Jovian weather patterns, a privilege unavailable to astrophysicists. It is important to recognize that astrophysicists are therefore trapped in a data-poor regime, with its myriad restrictions on how to construct models and interpret data. When faced with multiple explanations that are consistent with a given data set, astrophysicists often apply the principle of Occam's Razor: In the absence of more and better data, the simplest explanation is taken as the best one. To put it more tongue-in-cheek, one aims to be roughly accurate rather than precisely wrong. The need to recalibrate our scientific expectations and philosophies lies at the heart of this confluence of expertise.

From studying the atmospheres of Earth and the Solar System planets, researchers have realized that atmospheres are complex entities subjected to positive and negative feedback loops, exhibiting chemical, dynamical and radiative signatures over a broad range of time scales. Isaac Held of the Geophysical Fluid Dynamics Laboratory in Princeton, New Jersey, has argued that to truly understand these complex systems, one has to construct a hierarchy of theoretical models. These simulations range from one-dimensional, pen-and-paper models that isolate a key piece of physics to full-blown, three-dimensional general circulation models (GCMs)—used for climate and weather forecasting—that incorporate a complicated soup of ingredients to capture the intricate interactions between the atmosphere, land and oceans on Earth. For instance, GCMs concurrently solve a set of equations (called the Navier-Stokes equation) that treat the atmosphere as a fluid, along with a thermodynamic equation and diverse influencing factors such as orography (mountain formation) and biology. Many of these intricacies are unwarranted in theoretical investigations of exoplanetary atmospheres, and thus one of the key challenges is to realize how and when to simplify an Earth-centric model.

Whether knowingly or unknowingly, a hierarchy of one- to three-dimensional theoretical models has emerged in the astrophysical literature. Because the treatment of hot Jupiters and brown dwarfs—substellar objects not massive enough to sustain full-blown nuclear fusion at their cores—share several similarities, many of the pioneering models (by researchers such as Adam Burrows of Princeton University and Ivan Hubeny of the University of Arizona) were carried over from the latter to the former class of objects. Furthermore, the early models focused on the spectral appearance of hot Jupiters, with the most sophisticated variants borrowing from an established technique in atmospheric and climate science known as abundance and temperature retrieval. Given the spectrum of an exoplanet, this technique obtains the atmospheric chemistry and temperature-pressure profile consistent with the data. In the case of the hot Jupiter WASP-12b, Nikku Madhusudhan of Yale University and his collaborators inferred, using the retrieval technique, that the exoplanet possesses a carbon-to-oxygen ratio at least twice that of its star. If this result is confirmed—and the carbon-to- oxygen ratio is measured for other exoplanets—it offers a valuable link between the properties of an exoplanetary atmosphere and the formation history of the exoplanet.

Astrophysicists have been quick to realize that atmospheric chemistry and dynamics intertwine in a nontrivial manner to produce the observed characteristics of a hot Jupiter. Adam Showman, a planetary scientist at the University of Arizona, became one of the first researchers to harness the power of GCMs in studying hot Jovian atmospheres. Several other researchers from the astrophysical community (including myself) followed soon after. My collaborators and I generalized a benchmark test, which solves for the basic climatology of a (exo)planet using two methods of solution, to hot Jupiters. The transport of heat from the day-side to the night-side hemisphere of a spin synchronized exoplanet is—by definition—at least a two-dimensional problem. For extrasolar gas giants, the characteristic time scale on which the atmosphere reacts to such radiative disturbances spans many orders of magnitude, thus necessitating its theoretical consideration in three dimensions, an endeavor that is only tractable using GCMs. Several groups have now successfully adapted GCMs to model exoplanetary atmospheres and are obtaining consistent results. Some outstanding technical issues remain, but it is clear that three-dimensional models are necessary if one wishes to predict not just the spectral appearance of exoplanets, but simultaneously their phase curves and temporal behavior. As the astronomical state-of-the-art advances, the exoplanets being discovered will be more Earthlike, both in size and temperature, with the implication that GCMs will become even more relevant.

**Earthlike Exoplanets**

We are only starting to understand the basic properties of hot Jupiters, including why some appear more inflated than others and why some appear to redistribute heat more efficiently from their day-side to their night-side hemispheres. In the cases of HD 189733b and HD 209458b, simulated spectra and phase curves—the latter of which constrains the efficiency of heat redistribution—that are computed using GCMs are able to match their observed counterparts fairly well. Until the next generation of space telescopes becomes operational, these examples remain the cornerstones of our understanding of exoplanetary atmospheres.

HD 189733b is noninflated, meaning that its radius and mass are well matched by standard evolutionary theories of exoplanets (which predict the size of an exoplanet as it cools down from the primordial heat of formation). It also appears to be shrouded in haze of unidentified chemistry, because its spectrum in the optical—obtained via the Hubble Space Telescope by Frédéric Pont and David Sing of Exeter University, together with their collaborators—reveals a smooth, featureless slope consistent with Rayleigh scattering (the same process that causes the color of the sky as observed from Earth by preferentially affecting bluer sunlight).

By contrast, HD 209458b is free of haze, but it is markedly larger than expected from evolutionary calculations. Theoretical ideas for radius inflation include the suggestion that partially ionized, hot Jovian atmospheres behave like giant electrical circuits, which when advected past an ambient magnetic field invoke Lenz's law on a global scale: Nature abhors a change in magnetic flux. Electric currents and opposing forces are induced to counteract the horizontal winds; the consequent conversion of mechanical energy into heat, called Ohmic dissipation, is believed to be responsible for keeping some hot Jupiters inflated. However, it remains to be proven if hot Jupiters even possess magnetic fields like those of Earth and some of the Solar System planets. This field of research remains active.

The study of hot Jupiters remains relevant because we already have the data to inform our hypotheses and modeling efforts, thereby affording us the opportunity to sharpen our theoretical tools—as much of the salient physics is identical—before applying them to Neptunelike or even Earthlike exoplanets, for which the data are currently scarce or nonexistent. For many researchers, the ultimate prize is more familiar: to detect the spectrum of an Earthlike exoplanet orbiting a Sunlike star, and thereby answer age-old questions about the existence of extraterrestrial life. More succinctly, one wishes to establish if the solitary example of an Earth twin in orbit around a solar twin is the only possible cradle for life in the Universe. At the moment, such a quest remains elusive and appears out of the reach of even the next generation of space telescopes.

Instead, astronomers such as David Charbonneau of Harvard University and Jill Tarter of the SETI Institute have argued that a promising route toward detecting potentially habitable super Earths—Earthlike exoplanets with masses and radii somewhat larger than those of Earth—is to hunt for them around M stars (also known as red dwarfs). These diminutive cousins of our Sun, with only a tenth to half of its mass, comprise about three-quarters of the stellar population in our galactic neighborhood. There are several advantages to scrutinizing M stars: They are cooler in temperature than Sunlike stars, implying that their exoplanets may reside 10 to 100 times closer and yet still be able to harbor liquid water on their surfaces. Being more proximate to their M stars renders these exoplanets more amenable to detection via current, established astronomical techniques—namely, transit and radial velocity measurements. However, the price to pay is that they are expected to be spin synchronized and possess permanent day- and night-side hemispheres, much like their hot Jovian brethren. Such an expectation has led to theoretical concerns that their atmospheres may collapse due to the main constituent molecules condensing out on the frigid night sides. The astronomical approach to this conundrum is to charge forward with making new and better observations—after all, the answer is ultimately revealed by the data. For example, for the super Earth GJ 1214b, transmission spectra have already been obtained, but interpretations about its atmospheric composition remain controversial.

**Better Telescopes**

Experimentally, the next leap is to build dedicated, space-based telescopes capable of measuring high-resolution spectra of exoplanets over protracted periods of time. Astronomers around the world are mobilizing to launch such missions as the Exoplanet Characterization Observatory (EChO) and the Fast Infrared Exoplanet Spectroscopy Survey Explorer (FINESSE), as proposed to the European Space Agency (ESA) and the National Aeronautics and Space Agency (NASA), respectively. If and when these missions—or their successors—eventually fly (in the next decade or two), they will deliver a bounty of both spectral and temporal information on hundreds of exoplanets, from which we may infer their atmospheric chemistry, dynamics and climates. With a richly sampled data set of the emitted light from point-source exoplanets over time, one may construct a power spectrum that elucidates the characteristic time scales on which an exoplanet is flickering, indicating changes in temperature. Such a power spectrum of the atmosphere has been spectacularly constructed for the Earth's surface, spanning time scales of under a day (diurnal variations) to many millennia (called Milankovitch cycles, and inferred from the geological record). Certainly, space missions are saddled with demands that will not allow for the construction of power spectra on time scales longer than a few months, but it is likely that many of the characteristic peaks in the power spectra will be compressed into a shorter time span for close-in exoplanets such as hot Jupiters and super Earths.

The onus is on the theoretical community to lay down the foundation for understanding the climates of point-source exoplanets in general, thus moving us a step closer toward making more robust statements about their habitability.