On the Magnetic Field Strength of Active Region Filaments

Abstract: We study the vector magnetic field of a filament observed over a compact Active Region Neutral Line. Spectropolarimetric data acquired with TIP-II (VTT, Tenerife, Spain) of the 10830 Åspectral region provide full Stokes vectors which were analyzed using three different methods: magnetograph analysis, Milne-Eddington inversions and PCA-based atomic polarization inversions. The inferred magnetic field strengths in the filament are of the order of 600 - 700 G by all these three methods. Longitudinal fields are found in the range of 100 - 200 G whereas the transverse components become dominant, with fields as large as 500 - 600 G. We find strong transverse fields near the Neutral Line also at photospheric levels. Our analysis indicates that strong (higher than 500 G, but below kG) transverse magnetic fields are present in Active Region filaments. This corresponds to the highest field strengths reliably measured in these structures. The profiles of the Helium 10830 Ålines observed in this Active Region filament are dominated by the Zeeman effect.

\section{Introduction}

The magnetic field strength of Active Region (AR) filaments has long
remained poorly known or understood.  The situation for quiescent
filaments is notably more satisfactory since the early
measurements back in the 70s (see, e.g., the review by
\cite{ariste&aulanier07}). Field strengths measured in quiescent
structures, mostly using the Hanle effect on the linear polarization
of He {\sc i} D$_3$ at 5876\,\AA\ (e.g., \cite{SB77}), were found
to be in the
range of 3--15\,G (\cite{leroy83}; see also the review by
\cite{anzer&heinzel07}). However, more recent measurements, which also took
into account circular polarization, consistently show a tendency
towards larger field strengths. For example, the He {\sc i} 10830
\AA~investigation of \cite{trujillo02} and the He {\sc i} D$_3$ of
\cite{casini03} found, respectively, field strengths of the order of 20--40\,G and
10--20\,G (with field strengths up to 80\,G in the latter set of observations). With
similar techniques, using carefully inverted observations of
the \mbox{He {\sc i}} 10830\,\AA\ lines, \cite{merenda06} found field strengths of
about 30\,G. Other recent works that also obtained relatively strong filament
magnetic fields as compared with older measurements are \cite{paletou01}
and \cite{wiehr&bianda03}.
Of particular interest are the high field strength areas revealed by
Zeeman-dominated Stokes $V$ profiles in conjunction with scattering
dominated $Q$ and $U$ linear polarization signals obtained by
\cite{casini03}, showing how important it is to use all four Stokes
parameters for a proper determination of the field strength in prominences.

Well established properties of AR filaments are their systematically
stronger field strengths as compared to their quiescent counterparts and that
they lie lower in the atmosphere (see, e.g., \cite{aulanier&demoulin03}).
Zeeman observations, using coronagraphs, provided field strengths for
AR prominences in the range of 50--200\,G
(\cite{tandberghanssen&mckimmalville74}, and references therein).
However, the intrinsic difficulties of observing such low-lying structures
near the limb and the use of H$\alpha$ magnetography render these results somewhat
questionable (see \cite{lopezariste05} for the problems associated with using a
Zeeman-based formulation of the measurements made in H$\alpha$).
Perhaps, the most relevant modern estimate of the field strength in
AR prominences was provided by \cite{wiehr&stellmacher91}, who
measured the longitudinal magnetic field using the Stokes $I$ and
$V$ profiles of the \mbox{Ca {\sc ii}} IR triplet, finding a value of 150\,G,
compatible with previous measurements.
It is important for the discussion presented in our paper
to note that, in that observational study, linear polarization
signals ($Q$ and $U$) were not included. The only measurements of strong magnetic fields using \mbox{He {\sc i}} 10830\,\AA\ full Stokes vector that we are aware of, belong to a multicomponent flaring active region and reached values of $\approx$ 380\,G (see \cite{sasso07}).

Because measurements are scarce, we are still far from having a
complete picture of the field strengths that permeate solar filaments.
Some progress has come from using photospheric distributions of (mostly
longitudinal) fields that are extrapolated into the corona (usually using models of
constant-$\alpha$ force-free field). The study by Aulanier \& D{\'e}moulin (2003) is of
particular relevance to the present work. According to these authors, the
field strength of filaments found from the extrapolations is about 3\,G
for the quiescent case, 40\,G for active filaments (called ``plage
case'' in their work) and 15\,G for intermediate cases.
The typical height in the atmosphere of a filament base ranges from 20 to 10\,Mm, or
even lower, as one moves from quiescent to AR structures. These authors also
find positive field gradients with height that result in values of
0.1--1\,G\,Mm$^{-1}$, i.e., stronger fields are higher up in the atmosphere.
If these gradients are used to extrapolate down into the photosphere, the
fields there would be very similar to those found some Mm up in the
corona. For example a 10\,Mm-height active filament with 100\,G in the
corona would have no less than 90\,G close to the photosphere. The
argument could also be turned around to start from the fields measured
near the photosphere, close to AR Neutral Lines (NLs), and
infer what the fields could be high in the corona. In the case of kG-strong
plage fields, with high density areal filling factors (that can
reach up to 50\%; \cite{martinezpillet97}), one could expect fields
of several hundred Gauss within the filaments. Indeed,
Aulanier \& D{\'e}moulin (1998, 2003) provide arguments supporting a
relationship between the photospheric fields at
the base of the filament and the filament fields up in the corona. They
conclude that the stronger the photospheric fields at the
base of the filament, the stronger the field in the filament itself.
Similarly, the separation of the two opposite polarity regions
scales inversely with the field strength of the filament, the latter being larger
whenever the two polarities are closer together. Dense (highly packed)
fields in the photosphere correspond indeed to
the photospheric configuration found below AR filaments as observed
recently by \cite{lites05} and \cite{okamoto08}. In both
cases, they find opposite polarity ``abutted'' plage fields at the
NL, with sheared  horizontal fields in the hecto-Gauss range
and relatively high filling factors. These abutted field configurations seem to
also correspond with low-lying filaments structures. \cite{lites05}
comments that the height of the filaments on top of the abutted plage
fields is no more than 2.5\,Mm. Thus, high density (i.e., large
filling factor) horizontal plage fields near AR NLs, together with the
inferences from theoretical modeling (\cite{aulanier&demoulin98}, 2003),
 would indicate that fields of several hundred
Gauss can be expected in low-lying AR filaments.

Given the observation that energetic Coronal Mass Ejections are often associated with
AR filament eruptions (see, e.g., \cite{manchester08} and \cite{low01} for
a review), it is highly
desirable to develop diagnostic tools for direct measurements of
the AR filament magnetic field and its evolution, from its emergence
(see Okamoto et al. 2008) to the erupting phase. In this work, we present a
clear diagnostic tool of how this can be achieved using full Stokes polarimetry
of the \mbox{He {\sc i}} 10830\,\AA\ lines.
\begin{figure*}
\centering
\includegraphics[width=6cm,height=5.9cm]{cont_slit14_new.ps}
\includegraphics[width=6cm,height=5.9cm]{Si_stokesV14_blue_new.ps}
\\
\includegraphics[width=0.49\linewidth]{StokesI_new.ps}
\includegraphics[width=0.49\linewidth]{StokesQ_new.ps}
\includegraphics[width=0.49\linewidth]{StokesU_new.ps}
\includegraphics[width=0.49\linewidth]{StokesV_new.ps}
\caption{Top left: Slit reconstructed continuum intensity map centered at the NL
where the pores and penumbral-like formations can be identified. The
location of the slit for the time series map is displayed. Top right: \mbox{Si {\sc i}} Stokes V map normalized to the continuum intensity. Bottom: Stokes
profiles observed at the slit position indicated in the top images and averaged over the
time series (comprising 100 scans). The gray scale bar only applies to the polarization
signals. The Zeeman-like signature of Stokes $Q$ and $U$
are evident in this picture (\mbox{He {\sc i}} lines are centered at 0 and 1.2 \AA\
approximately. The zero in the wavelength scale corresponds to 10829.09 \AA). The arrows indicate the position of the presented Stokes profiles in Fig. \ref{fig:PBinversion1} and Fig. \ref{fig:PBinversion2}.
}
\label{fig:continuum_with_slit}
\end{figure*}

\section{Observations}

The observations described in this paper were carried out at the German Vacuum
Tower Telescope (VTT, Tenerife, Spain) on the 3rd and the 5th of July, 2005, using
the Tenerife Infrared Polarimeter (TIP-II, \cite{tip2}). TIP-II allows to measure
(almost) simultaneously the full Stokes vector for all the pixels along the
spectrograph slit. The slit (0\farcs5 wide and 35\arcsec\ long) was placed over
the target, a filament lying over the NL of active region NOAA 10781 situated
close to disk center (at coordinates N13-W4 around $\mu=0.97$ on July 3rd and
at N13-W29 or $\mu=0.92$ on July 5th) with the help of context H$\alpha$ and
continuum images. The H$\alpha$ frames showed the filament lying always on top
of bright plage regions visible immediately to either side of the absorption feature.
It is likely that the results presented in this paper particularly apply to
such an ``active'' filament configuration. SOHO/MDI (\cite{scherrer95})
magnetograms and continuum frames have been used to follow the evolution of the
active region as it crossed the visible surface.  The active region was
identified to be in its slow decay mode, encompassing a round symmetric leader
sunspot and a follower, spotless, plage region. At the location of the NL,
sporadic pore-like formations and penumbral-like configurations (with no evident
radially symmetric link to an umbral or pore structure as normal penumbra or partial
penumbrae do) were seen, especially on the 5th of July data-sets. Pores near the NL with both polarities
are identified and the penumbral-like region seems to correspond to field lines connecting
them. The NL was oriented in the N-E direction and showed a very compact
configuration, with the two polarities always remaining in close contact.  The
leader sunspot became invisible on the 6th of July after decreasing in size
while approaching the West solar limb.

Two observing strategies were used to map the NL region. Several spatial scans
that covered the whole AR were carried out in the course of the 3 days, and one
time series (with the slit fixed over the filament) was taken on July 5th. This
data-set was averaged in time in order to improve the signal-to-noise (S/N)
ratio of the spectral profiles.  The pixel size along the slit was
0\farcs17 and the scanning step (in the case of the rasters) was 0\farcs3
per step. The
exposure time per slit position was 8 seconds.  The adaptive optics system
(KAOS, \cite{kaos}) was locked on nearby pores during all the runs,
substantially improving the image quality of the observations carried out during
relatively poor seeing conditions, and providing a final estimated spatial
resolution of $\sim$1\arcsec.  Flat-field and dark current corrections were
performed for all the data-sets and, in order to compensate for the telescope
instrumental polarization, we also carried out the standard polarimetric
calibration (\cite{collados99}; \cite{collados2003}) for this instrument.

The observed spectral range spanned from 10825 to \mbox{10836 \AA}, with a spectral
sampling of $\sim$11.1 m\AA\,px$^{-1}$. However, a 3px binning in the spectral
domain was applied to all the data to increase the S/N while still preserving
subcritical sampling. The 10830 \AA\ spectral region is a powerful diagnostic
window of the solar atmospheric properties since it contains valuable
information coming simultaneously from the photosphere (carried by the
\mbox{Si {\sc i}} line at 10827 \AA ) and the chromosphere (encoded in the \mbox{He {\sc i}} 10830
\AA\, triplet). This He multiplet originates between a lower term $2^3S_1$ and
an upper term $2^3P_{2,1,0}$. Thus, it comprised three spectral lines,
namely a ``blue'' component at 10829.09 \AA\, ($J_L=1 \rightarrow J_U=0$) and a
``red'' component at $\sim 10830.3$ \AA\, that results from the remaining two
transitions ($J_L=1 \rightarrow J_U=1,2$), which appear completely blended at
typical solar atmospheric temperatures.
This multiplet is formed in the high chromosphere (\cite{avrett}) with no
contribution from photospheric levels, and serves as a unique diagnostic tool for
chromospheric magnetic fields. In filament structures, the height of formation of this
multiplet corresponds to the height of the opaque material inside them, which might
well correspond to typical coronal heights.

\subsection{Predominance of Zeeman-like signatures in AR filament He {\sc i} 10830\,\AA\
Stokes profiles}

The most striking finding encountered during the analysis of the data from this
campaign was the ubiquitous presence of Zeeman-like signatures in the Stokes
$Q$ and $U$ profiles of the \mbox{He {\sc i}} lines in the AR filament
(see Fig. \ref{fig:continuum_with_slit}). Recent
observations of this triplet carried out with the same instrument (Trujillo Bueno et al.
2002) have shown quiescent filament Stokes profiles dominated by atomic
level polarization (and its modification through the Hanle effect). These
forward-scattering signatures
(see Fig. 4 of Trujillo Bueno et al. 2002) correspond to one-lobe profiles that are positive
for the red component of the He triplet and negative for the blue one (in the positive Stokes $Q$
reference system) due, respectively, to selective emission and absorption
processes induced by the anisotropic illumination of the He atoms.
However, in the present case, the Stokes $Q$ and $U$ signals of the blue and red components
of the multiplet, exhibit the usual three-lobe profile that is expected from
the Zeeman effect (right panel of Fig. \ref{fig:continuum_with_slit}).
While some influence of atomic polarization in these profiles cannot be ruled out a priory,
it is clear that an explanation of their shapes should rely mainly on the Zeeman and
Paschen-Back effects. This fact, interesting in itself, will be briefly discussed below
(Sect. 3.3).

The continuum frame in the top panel of Figure \ref{fig:continuum_with_slit} was
reconstructed from the various slit scan positions (thus reflecting the real spatial
resolution of the data) of a map centered at the NL. This data-set was taken on the 5th of
July, 2005. The presence of pores and penumbral-like structures
is evident in this frame and their location corresponds to the NL. Right after this map,
a time series with the slit fixed over the NL
was performed in order to produce high S/N Stokes profiles.
The bottom four panels of Fig. \ref{fig:continuum_with_slit} correspond to the
Stokes parameters obtained after averaging over the full time series
(100 scans), which resulted
in a S/N of about 4000. The inversions performed later in this work were
carried
out on the spectral profiles extracted from this time-averaged data-set. The location of the NL
had to be guessed in real-time during the observations. The chosen slit position for the
time series is represented
by a vertical black line in the continuum frame of Fig. \ref{fig:continuum_with_slit}.
Posterior analysis has
shown that this location was a few arcsec off at one side of the actual NL, and that the
slit would have been more correctly placed at $x=18$\arcsec\
(referred to the abscissa coordinates in the continuum panel)
instead of at $x=20$\arcsec. Nevertheless, indications of polarity changes along the slit
in the \mbox{Si {\sc i}} line (evident around $y=13$\arcsec\ in the Stokes $V$
map of the same
figure) show that we were not very far away from it.

Inspection of the profiles at other slit locations, and of similar maps taken on the 3rd
of July, consistently exhibited Zeeman-dominated linear polarization signals in those regions
where the He line showed strong absorption features. However, in regions with weak He absorption
(typically near the boundaries of the filament), we often find linear profiles dominated by
scattering polarization signatures, similar to those observed by \cite{trujillo02}.
This indicates that, in general, the profiles obtained in AR filaments can
have significant contributions from various competing physical processes:
the atomic level polarization due to
the anisotropic illumination of the He atoms, the modification of these population imbalances
induced by the presence of a magnetic field inclined with respect to the axis of symmetry
(Hanle effect), and the Zeeman splitting characteristic of strong magnetic fields.
A detailed study with these physical ingredients for all the maps observed at this NL is
beyond the scope of this paper. In the present work, we concentrate on the implications of
the clear Zeeman-dominated signatures observed almost everywhere in this filament.

One point is worth mentioning after a simple visual inspection of Fig.
\ref{fig:continuum_with_slit}.  A comparison of the signs displayed by the
linear polarization profiles of the \mbox{He {\sc i}} and \mbox{Si {\sc i}} lines shows that
the orientation of the vector magnetic field is different in the filament than
in the underlying photosphere. In the Stokes $Q$ frame, the \mbox{He {\sc i}} 3-lobed
profiles have the same sign all along the slit whereas the \mbox{Si {\sc i}} line (and all the other smaller photospheric features at redder wavelengths)
exhibits a change in sign around $y=23$\arcsec. Conversely, in the Stokes $U$ frame, the
photospheric signals have a 3-lobed sign distribution that remains constant
almost all along the spatial domain while the \mbox{He {\sc i}} lines now show a sign reversal
near the middle of the slit. An analysis based on the formula for the azimuth angle by \cite{auer} gives an estimate of this difference between the azimuths inferred from the \mbox{Si {\sc i}} and \mbox{He {\sc i}} lines. We obtain values of the azimuth around $y = 5$\arcsec\ , which differ in $\sim 50^{\circ}$. This readily shows that the 10830 \AA\ spectral
region has a great potential to diagnose the orientation of the magnetic field
from the photosphere all the way up to the filament. We postpone this study to
a future paper.

\section{Vector magnetic field near the AR neutral line.}

Several analyses with various levels of complexity have been performed on the reduced data.
The first approach was a simple magnetograph-like analysis based on the
assumption of the weak-field approximation as formulated below.  This
method was applied to all the points in one of the maps obtained
during the campaign. We subsequently performed an analysis of the high S/N spectral profiles
obtained from the averaged time series with
a Milne-Eddington (ME) code and later with a more sophisticated inversion procedure
based on Principal Component Analysis (PCA) of a statistically generated database of
spectral profiles that account for the physics of atomic level polarization and the Hanle effect.
All these different methods consistently yield transverse field strengths in the filament
well above 500 G.

\subsection{Magnetograph analysis}

Typical Doppler widths for the red He {\sc i} line are in the range of 200-300 m\AA. The Land\' e
factors are 1.75 and 1.25 for the $J_L=1 \rightarrow J_U=1$ and the $J_L=1 \rightarrow J_U=2$
transitions, respectively. If the transitions are weighted with their line strengths,
an average Land\' e factor of $g_\mathrm{eff}=$1.42 for the red component of the Helium multiplet
is obtained. This Land\' e factor translates the above
Doppler widths to field strengths in the range of 2000-3000 G,
which is much stronger than the fields we expect for AR filaments. Together with the
assumption that the observed
signals are due to the Zeeman effect, these large Doppler widths justify the use of the
well-known weak field approximation (see, e.g., \cite{landi92}) as a first approach to
infer the magnetic field from the data.
In this approximation, the relation between the longitudinal field strength and Stokes
$V$ profile is given by:
\begin{equation}
V_o=f V_m= - f C B \cos\theta {{dI_m}\over{d\lambda}}.
\label{vmag}
\end{equation}
The subscript ``$o$'' stands for the observed profile, while ``$m$'' represents the profile generated
in the magnetic component that fills a fraction $f$ of the resolution element;
$B$ is the field strength, $\theta$
the angle between the line-of-sight (LOS) and the magnetic field direction and $C$ the
constant $4.67 \times 10^{-13} g_\mathrm{eff}\lambda_o^2$ (which forces the wavelength to be
expressed in \AA~and the field strength in G). The observed Stokes $I$ profile is given by:
\begin{equation}
I_o=f I_m +(1-f) I_{nm}
\end{equation}
with the subscript ``${nm}$'' referring to the non-magnetic component.
The last factor in the right-hand-side of Equation (\ref{vmag}) is the derivative of $I_m$
with respect to the wavelength. However, observations only provide the compound
profile $I_o$, product of the coexistence of magnetic and non-magnetic components in
the same resolution element. The present He {\sc i} analysis benefits from the
fact that the non-magnetic areas of the Sun display a very weak He absorption, so the
derivative of the Stokes $I_{nm}$ profile with respect to $\lambda$ can be neglected
(we note that, in non-magnetic regions, $dI_{nm}/d\lambda$ is found to peak at
one order of magnitude smaller values than
$dI_{m}/d\lambda$ as observed in the filament). This nicely eliminates
any dependence of the inferred longitudinal field on the unknown filling
factor. Following this argument of ignoring the derivative of the $I_{nm}$ profile
with respect to wavelength, the longitudinal field is directly inferred as:
\begin{equation}
B_{\parallel} =B \cos\theta= - {{1}\over{C}}
{{V_o(\lambda)}\over{{{dI_o}\over{d\lambda}}(\lambda)}}
\label{blos}
\end{equation}
The $\lambda$ dependence in the last ratio is written explicitly to emphasize that each point
within the profile provides an estimate of the longitudinal field. Thus, the way to
solve Equation (\ref{blos}) is through a least-square fit to all the observed
points within the profile.

\begin{figure*}
\centering
\includegraphics[width=6cm]{DoppHe14-2_new.ps}
\includegraphics[width=6cm]{BlosSi14_new.ps}
\includegraphics[width=6cm]{BlosHe14_new.ps} \\
\includegraphics[width=6cm]{coreHe14_new.ps}
\includegraphics[width=6cm]{BtransSi14_new.ps}
\includegraphics[width=6cm]{BtransHe14_new.ps}
\caption{From left to right and top to bottom: He {\sc i} Doppler width, \mbox{Si {\sc i}} LOS magnetic
field, He {\sc i} LOS magnetic field, He {\sc i} core intensity, \mbox{Si {\sc i}} transverse magnetic
field and He {\sc i} transverse magnetic field. The corresponding continuum intensity
frame can be seen in the top panel of Fig. \ref{fig:continuum_with_slit}.
}
\label{fig:magneto}
\end{figure*}
For transverse fields, a similar equation applies including
the quadratic dependence of the linear polarization signals on
the transverse component of the field
and the second derivative of the Stokes $I$ profile with respect to $\lambda$.
A similar argument gives:
\begin{equation}
B^2_{\perp} =B^2 \sin^2\theta=
{{4}\over{C^2}}{\sqrt{Q^2_o(\lambda_o)+U^2_o(\lambda_o)}
\over{{{d^2I_o}\over{d\lambda^2}}(\lambda_o)}}.
\label{btra}
\end{equation}
Like in the equation for the longitudinal field, the filling factor (of a
non-magnetic component) is absent.  On the other hand, and in contrast to Eq.
(\ref{blos}), the right hand side of Eq. (\ref{btra}) must be computed at the
central wavelength, $\lambda_o$, as this equation can only be formally derived
there (\cite{landi92}). Having only one wavelength point to compute the
transverse fields prevents the use of a least square approach, unlike in the
case of the longitudinal field equation. This, in turn, complicates the
determination of the transverse component because the Stokes $Q$ and $U$ signals are
often found to be at, or close to, the noise level. The quadratic average of
the linear polarization signals in \mbox{Eq. (\ref{btra})} hinders any possible
noise reduction and generates a veil of noise-induced transverse fields that
proves problematic in the use of this equation in many cases.
This is particularly true for the He {\sc i} line, whose second
derivative of the (broad) Stokes $I_o$ profile is very small. This
derivative in the denominator in \mbox{Equation (\ref{btra})}
results in values of the noise-induced transverse fields as large as
300-400 G, only slightly smaller than the actual values
measured in the filament region. Note that this 300-400 G noise level
corresponds to a S/N of 1000 in the map in \mbox{Fig. \ref{fig:magneto}}, whereas in
\mbox{Fig. \ref{fig:continuum_with_slit}} the S/N is four times larger.

Thus, and in order to obtain high S/N maps of the chromospheric
transverse fields, we averaged the unsigned $Q$ and $U$ signals over
three wavelengths, the central wavelength and two wavelengths at 278 m\AA\
on each side of line center, corresponding approximately to the locations of the peak
signals in the linear polarization profiles, so
\begin{equation}
\sqrt{Q^2_o(\lambda_o)+U^2_o(\lambda_o)}\approx\sqrt{{\bar Q}^2_o+{\bar U}^2_o},
\label{btra2}
\end{equation}
where ${\bar Q_o}$ and $\bar{U_o}$ are the averaged linear polarization signals over
these three wavelength points. We emphasize that this approach reduces
the noise level considerably
while at the same time preserves the distribution of transverse fields
already present when directly applying Eq. (\ref{btra}). The results
obtained from this approach have been validated by the ME inversions presented
in the next section.

A similar approach was followed with the \mbox{Si {\sc i}} line to obtain the
photospheric magnetic field. The Land\' e factor for this line is 1.5.
The spectral profiles arising from this transition are broad enough to partly justify
the application of the weak-field approximation. However, like for any other photospheric line,
the vector-magnetograph data obtained with this method have to be interpreted with caution.
In particular, the results for this line are affected by filling factor effects as is the
case for commonly used photospheric
magnetograph data. For the \mbox{Si {\sc i}} line, the exact formulation of Eq. (\ref{btra}) was
used. The intrinsically larger photospheric magnetic field strengths together with the
fact that the \mbox{Si {\sc i}} Stokes $I$ profile has a sizeable second derivative, makes
the determination of $B_{\perp}$ less affected by noise when using only the central wavelength
point of the linear polarization profiles for its determination.

Lastly, the Stokes $I$ profiles
of both \mbox{He {\sc i}} and \mbox{Si {\sc i}}  were fitted with
gaussian functions from which the line center, line width and strength were inferred.
In the case of the \mbox{He {\sc i}}, this approach proved very useful to identify the location of the
filament above the NL.

The results from this approach are presented in Figure \ref{fig:magneto} for
the map observed on the 5th of July. The He {\sc i} Doppler width and line core
frames show that the line becomes deeper and narrower inside the filament. This
could be an indication of the presence of a denser and cooler plasma than in
its immediate plage surroundings.  These frames also reveal that the filament
has a highly twisted topology, with filamentary threads running at almost 45
degrees from the direction defined by the NL. This becomes evident when
compared to the Si and He LOS magnetograms.  The twist in the filament is more
obvious in the bottom half of the frames than in the top part, where a more
diffuse linear topology is observed. These twisted signatures could be clearly
seen in the observations from the 5th of July, while they were hardly visible
in the maps obtained two days earlier. The \mbox{Si {\sc i}} LOS magnetogram
shows that the separation between the two polarities is less than 5\arcsec at
photospheric levels when the magnetogram is scaled at $\pm$1000 G. A tighter
scaling with a smaller threshold would provide a much narrower NL channel.  In
the \mbox{He {\sc i}} LOS magnetograms, with a scaling of $\pm$400 G, the
channel is practically absent, showing how intense the plage surrounding this
AR filament was. A thread-like structure in the longitudinal component running
parallel to the direction of the twist (i.e. at 45 degrees from the NL) in the
mid range of the frame is observed. Note that the AR is 22 degrees off disk
center so horizontal threads in the filament can give rise to sizeable Stokes
$V$ signals.  The intrinsic LOS magnetic field strength at the NL (we stress
that no filling factor effect is included in these estimates) is typically 100
G for each of the two polarities (in agreement with older measurements of AR
filaments that did not include transverse field measurements).  The neighboring
plage displays an intrinsic longitudinal field strength of 200-400 G. If we
assume an intrinsic photospheric field strength in the plage of \mbox{1400 G} (see
\cite{martinezpillet97}), this would indicate a filling factor of about 20\% in
this layer.

In both the photospheric and chromospheric magnetograms, the region where the
longitudinal field becomes weaker corresponds to the region with very strong
transverse fields. The somewhat wider photospheric NL channel corresponds to a
wider photospheric area occupied with strong transverse fields, while the
transverse fields in the chromosphere fill a narrower region. These transverse
fields seem to follow the structures observed in the He line core frame: a more
linear topology in the top part and a twisted formation in the bottom half.
What becomes readily surprising is the magnitude of the transverse fields
observed in the chromosphere (and to the extent that the weak field
approximation is valid, in the photosphere as well). The photosphere shows
fields in the range of up to kG in the region between abscissae $x=15$\arcsec\
to $x=20$\arcsec\ in the frames of Fig. \ref{fig:magneto}.  This area shows
penumbral-like structures and pores in the white light image.  It is clear
that these orphan penumbral-like regions are magnetically linked to the main body of
the AR filament and correspond to horizontal field lines at the photospheric
surface. The transverse fields derived from the \mbox{Si {\sc i}} line
undoubtedly show the twist configuration similar to that observed in the
\mbox{He {\sc i}} lines.  In the chromosphere, the transverse fields are also
strong. The present analysis yields transverse fields with strengths in the
range of 500-600 G, strongly concentrated in the narrow NL channel. Note that
this region, at best, a few arcsec wide, is below the low resolutions
of old observations. These large transverse magnetic fields are well above
most, if not all, of past field strength measurements in AR filaments.

The various maps for the observations on the 3rd of July support the hypothesis
that an AR filament corresponds to a narrow channel characterized with strong
($\sim 500$ G) transverse fields.  There is a clear spatial correlation between
strong absorption signatures in \mbox{He {\sc i}} (measured by the line core
intensities), and narrow (or ``abutted'') plage regions with strong transverse
fields in the chromosphere.  A complete description of all these maps is
postponed to a future paper. But it is interesting to point out the
following main differences with the map of Fig. \ref{fig:magneto}: the absence
of evident twisted threads, the absence of penumbral-like regions in
white-light maps, weaker and more diffuse transverse fields at the photosphere
(not reaching kG levels) and the presence of weaker \mbox{He {\sc i}}
absorption signatures in regions outside the intense NL channel that display
atomic polarization signals.

\begin{figure*}
\centering
\includegraphics[width=14cm]{invers13_015_p49_oplot_BW_robertosample_14ENE09.ps}
\caption{Stokes profiles of the \mbox{He {\sc i}} 10830 \AA\ triplet at a height
of 8 arcsec in Fig. \ref{fig:continuum_with_slit}
in the mid portion of the twisted region of the
filament. The dots represent the observed profiles obtained after averaging
the time series in order to achieve a larger S/N. The solid line corresponds
to the best fit achieved with a Milne-Eddington inversion code that takes
into account the IPB effect. The magnetic field obtained from
this particular fit is $B = 626\pm 10$ G,
$\theta = 80.2^{\circ}\pm 0.1^{\circ}$ and
$\chi =174.9^{\circ}\pm 0.1^{\circ}$ (the origin of azimuths being the local
solar radial direction).
The dot-dashed line shows the best fit obtained from an inversion in a slab
model of constant properties that accounts for the joint action of scattering
polarization
and the Hanle Zeeman effects.
This fit provided for the magnetic parameters
$B = 667\pm 16$ G, $\theta = 81.1^{\circ}\pm 0.3^{\circ}$ and $\chi = 174.8^{\circ}\pm 0.5^{\circ}$.
}
\label{fig:PBinversion1}
\end{figure*}

\begin{figure*}
\centering
\includegraphics[width=14cm]{invers13_015_p27_oplot_BW_robertosample_14ENE09.ps}
\caption{
Same as Fig. \ref{fig:PBinversion1} but for a profile near the bottom portion of the
twisted region of the filament (at a height of 5 arcsec).
The fitted parameters for the ME case are
$B = 627\pm 10$ G,
$\theta = 73.7^{\circ}\pm 0.3^{\circ}$ and
$\chi =5.7^{\circ}\pm 0.1^{\circ}$. The PCA-based algorithm including Hanle effect provided:
\mbox{$B = 664 \pm 15$ G}, $\theta = 74.9^{\circ}\pm 0.6^{\circ}$ and $\chi = 5.1^{\circ}\pm 0.8^{\circ}$.
}
\label{fig:PBinversion2}
\end{figure*}

\subsection{Milne-Eddington inversions}
The interest of the magnetograph analysis is that it is easy
computed over a complete map with a very low computational effort. However, in
order to validate the weak-field approximation inferences of the magnetic
field strength in the filament, we carried out a Milne-Eddington (ME) inversion
of the \mbox{He {\sc i}} Stokes profiles of the averaged time series. These
high S/N data correspond to the location of the vertical line at $x=20$\arcsec\
in the map of Fig. \ref{fig:continuum_with_slit}.  The inversion code we used
(MELANIE; \cite{MELANIE}) computes the Zeeman-induced Stokes spectra---in the
Incomplete Paschen-Back (IPB) effect regime---that emerge from a model
atmosphere described by the Milne-Eddington approximation. This assumes a
semi-infinite constant-property atmosphere whose source function varies
linearly with  optical depth, $\tau$.  This approach does not account for the
atomic-level polarization induced by the anisotropic radiation pumping of the
He atoms from the underlying photospheric continuum.

For magnetic fields between 400 and 1500 G the Zeeman splitting of the upper
J-levels of the \mbox{He {\sc i}} triplet is comparable to their energy separation. Thus, it
is crucial to compute the energy levels in the IPB regime to avoid an
under-estimation of the magnetic field strength when carrying out the inversion
to interpret observations (\cite{sn-pb}; \cite{sasso}).

The Milne-Eddington model atmosphere uses a set of eleven free parameters that
the inversion code modifies in an iterative manner in order to obtain the best fits to
the observed Stokes profiles. These free parameters are: the magnetic field strength ($B$),
its inclination ($\theta$) and azimuth ($\chi$) angles in the
reference system of the observer (the azimuth is measured
with respect to the local solar radial direction), the line strength ($\eta_0$),
the Doppler width ($\Delta\lambda_D$), a damping parameter, the LOS velocity, the source
function at $\tau=0$ and its gradient, a macroturbulence factor and a stray light
fraction $f$.
An initial guess model must be provided to the code. In order to prevent the inversion scheme
from getting locked into local minima, we carried out the inversions using several different
initializations, in which $B$, $\theta$ and $\chi$ were taken from
the above-described magnetograph analysis and the remaining parameters were obtained
from random perturbations of reasonable pre-set values.
All parameters were free except for the stray light fraction, which was set to zero
(i.e., assuming a filling factor of 1; we refer to the explanation in section 3.1).
We stress that, even if the data is contaminated with stray-light from the surroundings,
using a filling factor of unity allows a direct comparison with the results
given in the previous section.

Figures \ref{fig:PBinversion1} and \ref{fig:PBinversion2}
show the results of the ME inversion of the full Stokes
vector for two positions along the filament (only two
fits are shown but the complete set of profiles from
Fig.~\ref{fig:continuum_with_slit}
were inverted). The excellent performance of the ME inversion
tells us that there are barely any scattering polarization signatures in these
profiles, and that the formation physics of the multiplet is adequately described
in the IPB regime. As expected, the magnetic field strengths inferred from all the
inversions are of the order of 600--700\,G with very high
inclinations ($\sim 80^{\circ}$)
with respect to the line-of-sight, thus confirming the existence of
strong transverse
fields in this AR filament. The exact values obtained for
Fig.~\ref{fig:PBinversion1}
are $B_{\parallel}=109$\,G and $B_{\perp}=616$\,G while for
Fig.~\ref{fig:PBinversion2} we
obtain $B_{\parallel}=176$\,G and $B_{\perp}=602$\,G.
In the solar reference frame, the magnetic field
vector turns out to be close to horizontal, 59$^\circ$ and 99$^\circ$ of
inclination with respect to the local vertical, respectively.

The magnetic field strengths derived from the weak field
regime analysis are systematically under-estimated by \mbox{$\sim 100 - 150$ G}, but
the retrieved inclinations and azimuths are in very good agreement with those
obtained from the ME-IPB inversions. This means that the weak field approximation
yields a reliable magnetic field topology and a low-end value for the field strength.

\subsection{PCA-based atomic polarization inversion}

The observations described above were also independently inverted
using the Principal Component Analysis (PCA) approach described by
L{\'o}pez Ariste \& Casini (2002). In fact, pattern recognition techniques are particularly
well suited to attack ill-defined inversion problems characterized by
computationally intensive forward problems. This is exactly the case
for spectropolarimetric inversions in prominences and filaments,
where the Stokes profiles are often formed by the scattering
of resonant radiation. The computation of the emergent polarization in
such case requires the preliminary solution of the non-LTE problem of
atomic-level excitation by anisotropic illumination of the plasma
from the underlying photosphere. The presence of a magnetic field
further modifies the ensuing atomic polarization through the Hanle
effect (see, e.g. \cite{LL04}, for a review of atomic polarization effects).
This problem, which is totally by-passed in the ME approach
to spectropolarimetric inversion, completely dominates the numerical
computation of line polarization in radiation scattering. For this
reason, pattern recognition techniques provide a very attractive
strategy to Stokes inversion for radiation scattering, since the
numerically intensive forward problem is solved once and for all for
a comprehensive set of illumination, thermodynamic, and magnetic
conditions in the plasma for the problem at hand. The goal of these
techniques is thus to build universal databases of profiles that can
be searched for the solution to any given Stokes inversion problem.
Principal Component Analysis additionally provides a way of
compressing database information, by reducing it to a few principal
component profiles that contain all the fundamental physics of the
formation of the emergent Stokes profiles (\cite{LC02,CBL05}).

For the inversion of the observations illustrated in this paper, we created a
database of 250,000 profiles, spanning all possible orientations of the
magnetic field, with strengths between 0 and 2000\,G. The illumination
conditions were set by the radiation temperature and center-to-limb variation
(CLV) profile of the photospheric radiation at 1$\mu$m (\cite{cox}), and by
assuming a range of scattering heights between 0 and 0.06\,$R_\odot$ (42 Mm).
The LOS inclination in the database used for the inversion of the July 5, 2005,
observations spanned between 20$^\circ$ and 30$^\circ$. The thermal Doppler
width and micro-turbulent velocity, responsible for the overall profile
broadening, were accounted for by introducing an equivalent temperature
spanning between $10^4$ and $2\times 10^4$\,K.  Finally, the database Stokes
profiles were calculated by integrating the polarized radiation emerging from a
homogeneous slab with optical depth at line center varying between 0.2 and 1.4.
Very soon it was realized that such conventional scattering scenario could not
fit satisfactorily the observations, unless some additional depolarizing
mechanism could be accounted for to explain the surprisingly low level of
atomic polarization revealed by the Zeeman-like shape of the Stokes $Q$ and $U$
profiles. For this reason, in the creation of the inversion database we
introduced an ad-hoc weight factor for the anisotropy of the photosphere
radiation, ranging between 0 and 1. Extensive inversion tests that were run
over the entire filament map consistently gave anisotropy weight factors
significantly smaller than unity, with a predominance of values around 0.2,
thus confirming the presence of some unidentified depolarizing mechanism in the
formation of the observed profiles. The possible physical origin of such
depolarization has been extensively discussed by \cite{CML09}, but see also
Trujillo Bueno \& Asensio Ramos (2007).

The results of the PCA inversion are overplotted to the ME inversions as dot-dashed lines in Figs. \ref{fig:PBinversion1} and \ref{fig:PBinversion2}. The PCA inversion provides a magnetic field of $B=667$\,G in the first case and \mbox{$B = 664$\,G} in the latter. The inclination angles are almost the same as the ones obtained by the ME inversions, \mbox{$\sim 59^{\circ}$} and $\sim 99^{\circ}$ with respect to the local vertical, respectively. The inferred magnetic field azimuths are also in very good agreement.

\section{Conclusions}

The He {\sc i} 10830 \AA~lines observed in a filament at the NL of Active Region
NOAA 10781 exhibit linear polarization profiles dominated by the
Zeeman effect. Three different independent analyses (magnetograph approximation,
ME inversion and PCA inversion including atomic polarization) of the four Stokes
profiles consistently support inferred field strengths in the range 600-700 G at the formation
height of the Helium triplet. These fields are 3 to 7 times larger than those measured heretofore
in AR filaments. The field strengths found at the NL are largely horizontal with 500-600 G
transverse fields. The longitudinal component is typically measured to be in the range
of 100-200 G, in agreement with past measurements. It must be stressed that these previous observations
of AR filaments did not include transverse field estimates.
Thus, it is clear that the inclusion of the Stokes $Q$ and $U$ profiles in our analysis is what
has allowed the detection of such strong magnetic fields in AR filaments. While
the role played by the linear polarization signals was evident already from the magnetograph
analysis, a further test was made with the ME inversion to prove this point.
Using as input the observed $I$ and $V$ profiles but replacing the Stokes $Q$ and $U$ profiles
with noisy data, the ME inversions resulted in fields strengths in the range of 100-200 G.
We thus propose that the lack of full Stokes
polarimetry was the main reason why past measurements did not find the high field
strengths reported in this work. Note also that the spatial extent that displays such strong fields
is not more than a few arcsec wide, which would be hardly visible in low resolution
measurements. It is also clear that it becomes crucial to search for the
signatures of these large field strengths in different spectral windows such as H$\alpha$ and
the Ca {\sc ii} triplet.

A similar trend to infer higher field strengths (albeit in a lower range) is
presented in modern measurements of quiescent filaments that is also ascribed
to the use of full Stokes polarimetry (Casini et al. 2003). In the case
of quiescent filaments, observations are dominated by atomic polarization and its
modification through the Hanle effect. In the study presented
here, however,
the linear polarization profiles are dominated by the Zeeman effect as long as the profiles
are coming from the regions of strong Helium absorption. The reason why atomic polarization
signatures are almost absent from these profiles is not yet well understood. On the one hand,
the observed high field strengths are approaching the values where the Zeeman effect dominates
over atomic level polarization (\cite{tb07}), but even at this high field strengths,
one would have expected a clearer atomic polarization signature. On the other hand, some mechanism
to reduce the anisotropy of the radiation field could be present. For example,
Trujillo Bueno \& Asensio Ramos (2007) suggested that, if the
radiation comes from a high opacity region, the isotropy of the radiation field will
be such that a much reduced atomic polarization would be induced. A more recent proposal by \cite{CML09}
suggests that the presence of a randomly oriented field entangled with the main
filament field, and of a similar magnitude, could also explain the absence of atomic polarization
signatures in these profiles. It is interesting to note that the ad-hoc anisotropy
weight factor introduced in the PCA inversions was found to be in the range 0.1-0.5
over the filament region, with a predominant value of 0.2. Although the various atomic processes
generating polarization signals cannot be cleanly separated, this persisting value
of 0.2 is evidence of their presence in our profiles
(note that a weight factor of 1 corresponds to an anisotropic illumination described by the
standard CLV of the photospheric radiation field). We also
stress that the reason for the absence of atomic level polarization signatures is not simply due to
what is commonly referred to as a Van-Vleck configuration of the vector field. For example,
whereas the profiles inverted in Fig. \ref{fig:PBinversion1} give an inclination
(with respect to the local vertical) close to
the Van-Vleck angle (59$^\circ$; Van-Vleck angle corresponds to $\sim$ 55$^\circ$), those in
Fig. \ref{fig:PBinversion2} yield an inclination very far from it (99$^\circ$).

What is the origin of these strong transverse fields? This question relates
directly to the problem of filament formation and mass loading. Two basic
scenarios are commonly used to explain how these structures are formed:
photospheric (shearing) motions and flux emergence (see the recent review by
\cite{lites08}).  The first one uses photospheric plasma flows that move,
tangle and reconnect the  field lines of an already emerged active region to
form the filament directly in the corona. These processes include in one way or
another some form of magnetic cancellation and reconnection that provides a
source for mass upload of the filament. In this scenario, the presence of such
strong magnetic fields must be related to the existence of a dense plage
configuration at the NL and the low gradients inferred by
\cite{aulanier&demoulin03} for their model of AR filaments.  Our observations
pose the question of how filament field strengths in the range of 600-700 G may
be generated in this scenario from a surrounding ``abutted'' plage that has a
longitudinal field of no more than 400 G at the height of formation of the
Helium lines. The emergence of a flux rope from below the photosphere scenario has recently received strong
support from  the observations of \cite{okamoto08}. If the flux ropes are formed below the surface,
the answer to the observed field strengths could be related to the balance between
buoyancy forces and gravity acting over the flux system. Although this
balance is not yet fully understood, it is clear that the stronger the fields the easier
is for the flux system to emerge into the Corona and carry a significant amount of trapped
photospheric mass (\cite{Archontis04}).
Note that the photospheric transverse fields observed at the NL
are also very strong (including pore-like structures) and could represent the bottom part of
the flux rope system once emerged into the atmosphere.

It remains to be studied whether the observed strong transverse field strengths
presented in this paper are common to all AR filaments or only to those
surrounded by exceptionally dense plages. An extension of the present study to
other ARs with different degrees of activity is mandatory.

\begin{acknowledgements}
This work has been partly funded through project ESP2006-13030-C06-01.  The
National Center for Atmospheric Research (NCAR) is sponsored by the National
Science Foundation.  Help received by C. Kuckein during his stay at HAO/NCAR is
gratefully acknowledged. R. Manso Sainz has been partially supported by Spanish
Ministerio de Ciencia e Innovaci{\' o}n through project AYA2007-63881. H.
Socas-Navarro helped with the implementation of the MELANIE code and with the
interpretation of the obtained results. Comments on the manuscript
by B.C. Low and J. Trujillo Bueno are gratefully acknowledged.
\end{acknowledgements}


Local Group(s)

Abstract: The properties of the galaxies of the Local Group are reviewed, followed by a brief discussion of nearby groups. The galaxy groups in our vicinity - the M81 group, the Cen A group, and the IC 342/Maffei group - are in many respects Local Group analogs: Their luminosity functions, galaxy content, fractional galaxy type distribution, crossing times, masses, and zero-velocity surface radii are similar to those of the Local Group. Also, the nearby groups usually consist of two subgroups, some of which approach each other and may ultimately merge to form a fossil group. These poor groups contrast with the less evolved, loose and extended galaxy ``clouds'' such as the Scl group and the CVn I cloud. These are characterized by long crossing times, are dominated by gas-rich, late-type galaxies, and lack gas-deficient, low luminosity early-type dwarfs. These clouds may be groups still in formation. The local Hubble flow derived from the clouds and groups is very cold.

\section{Why are galaxy groups interesting?}
\label{sec:1}

The most conspicuous gatherings of mass and luminous matter in the
Universe are galaxy clusters.  However, the majority of nearby
galaxies -- about 85\% (\cite{Tully87, Kara05}) -- is observed to be
located outside of clusters and can be found mainly in galaxy groups.
While galaxy clusters reveal the location of the highest concentration
of visible and dark matter, these less massive galaxy agglomerations
trace the distribution of the filaments of the cosmic web and hence
the extended distribution of dark matter in less dense regions.  In
terms of luminous matter, the fraction of mass locked up in stars
increases with decreasing group size (\cite{Eke05}).  Most of the
stellar mass is in groups of the size of our Local Group, whereas
massive clusters only contain 2\% of the stellar mass in the Universe
(\cite{Eke05}).

Galaxy groups come in a variety of different morphologies, shapes, and
sizes, including, for instance, seemingly unbound ``clouds'' and
``spurs'' (\cite{Tully87}), loose groups, poor groups, compact groups,
rich groups, etc.  Precisely what constitutes a group is a question of
definition and depends primarily on how many galaxies down to a
certain limiting magnitude are found within a certain volume.  Common
usage often considers an agglomeration of some 30 galaxies within a
radius of one to two Mpc a galaxy group, but this broad definition is
usually adjusted to the practical needs of the actual application; for
instance, when searching for groups in redshift surveys (e.g.,
\cite{Lee04, Eke04}).

The different types of richness and compactness of groups permit us to
study galaxy evolution and environmental effects in lower-density
regions, and to contrast the results with the properties and the
evolutionary state of galaxies observed in the field or in rich
clusters.  The degree of clustering of groups increases with
increasing group luminosity (\cite{Padilla04}).  Within groups, the
earlier-type galaxies are more strongly clustered and tend to lie
closer to the centers of the groups (\cite{Girardi03}), similar to the
well-known morphology-density relation in galaxy clusters
(\cite{Dressler80, Binggeli87}). The final stage of evolution within
groups may be represented by fossil groups, which are dominated in
light and mass by a large central elliptical galaxy, presumably the
end product of major mergers (\cite{Jones03, D'Onghia05}).  The
different types of groups appear to epitomize different stages from
early to advanced structure formation.

Members of galaxy groups and ensembles of groups are kinematic tracers
that can be useful for mass determinations and that help to uncover
the local properties of large-scale galaxy flows, such as the local
Hubble flow.  It has even be suggested that they may reveal the
effects of dark energy on small scales ($>7$ Mpc; e.g., \cite{Kara02a,
Kara03a, Maccio05}).

The title of this contribution as assigned to me can be interpreted in
two ways:  Apart from the Local Group, one may understand it to cover
galaxy groups that are local, i.e., nearby, or one can interpret it as
referring to galaxy groups that are counterparts of the Local Group.
In fact, the nearby galaxy groups often combine these two properties.
Thus in this review I will concentrate mainly on the Local Group and
nearby Local Group analogs.

Our own Milky Way is a member of a poor groups of galaxies, the Local
Group.  The nearest neighboring group, the Sextans-Antlia group, is a
an even poorer group consisting only of dwarf galaxies -- sometimes
also referred to as a dwarf group \cite{vdB99a}, \cite{Tully02}.
Additional poor groups are located in our vicinity.  In a number of
ways we live in a fairly average extragalactic neighborhood.

\section{The Local Group}

The Local Group is the home to two large spiral galaxies, our Milky
Way and M31.  The Local Group is the only place where we can study the
ages, chemistry, star formation history, and kinematics of a range of
different galaxies in exceptional detail based on their resolved
stellar populations.  Here we can truly connect stellar and
extragalactic astrophysics.  We can compare our observational findings
with the predictions of cosmological models or, vice versa, test those
predictions through targeted observations.  Owing to our ability to
resolve and study even the oldest populations and very faint stars in
the closest galaxies and in our own Milky Way, we can conduct
near-field cosmology (\cite{Freeman02}), uncovering their detailed
evolutionary histories based on their resolved fossil stellar record.

\subsection{The Local Group census and galaxy distribution}

In spite of their proximity and in spite of the large efforts invested
into studying them over the past decades, the galaxies of the Local
Group continue to provide surprises.  As has also been found in other
groups, the number of faint galaxies in the Local Group lies below
theoretical expectations by about two orders of magnitude, a
deficiency also known as the missing satellite problem or the
cosmological substructure problem (\cite{Klypin99, Moore99}).  A
number of solutions to solve this problem have been proposed, but none
has proved fully satisfactory until now.  This problem remains one of
the key questions in structure formation in standard cold dark matter
(CDM) models.  Nonetheless, considerable progress has been made in
recent years in identifying new Local Group members.  Within a
zero-velocity radius of approximately 1 Mpc (\cite{Kara02a}) the
current Local Group census comprises at least 42 galaxies (including
the tidal streams of Sagittarius around the Milky Way and the great
stellar stream around M31, \cite{Ibata94, Ibata01}).

The three most luminous Local Group members are its spiral galaxies,
which form two subgroups:  The M31 plus M33 subgroup and the Milky Way
subgroup.  In terms of mass and luminosity, the Milky Way and M31 are
the two dominant galaxies and may be similarly massive
(\cite{Evans00}).  Each of these two galaxies is surrounded by an
entourage of mainly low-mass, gas-deficient galaxies.  In close
proximity to M31, one compact elliptical and three dwarf elliptical
(dE) galaxies are found.  M31 is the only Local Group galaxy with dE
companions.  The remaining low-mass early-type dwarfs in the Local
Group are all dwarf spheroidal (dSph) galaxies, the least massive
(total masses estimated to be a few times $10^7$ M$_{\odot}$), least
luminous ($M_V > -14$) type of galaxy known.  (For a more detailed
description of the different galaxy types of the Local Group, see
\cite{vdB00, Grebel01, GGH03}.) The dSphs and the dEs are almost all
found within a 300 kpc radius around the two dominant spirals.  This
radius also roughly corresponds to the size of the dark matter density
profile of $\sim$ L$_{\star}$ galaxies (e.g., \cite{McKay02,
Prada03}). The gas-rich, late-type irregular and dwarf irregular
(dIrr) galaxies, in contrast, show a much less concentrated
distribution and and are the most frequent galaxy type at larger
distances from the spirals.

This biased distribution is also apparent when considering the H\,{\sc
i} mass of dwarf and satellite galaxies as a function of distance to
the nearest principal galaxy:  The upper limits for the H\,{\sc i}
masses of dSphs are typically of the order of $10^5$ M$_{\odot}$ or
less, whereas dIrrs usually have H\,{\sc i} masses of at least $10^7$
M$_{\odot}$ (\cite{GGH03}).  In between, the so-called dIrr/dSph
transition-type galaxies are located.  These low-mass dwarf galaxies
share some of the properties of dIrrs (such as a measurable gas
content and recent star formation) and of dSphs (such as prominent old
populations and low luminosities).  Morphological segregation akin to
that within the Local Group is also observed in other groups and
indicates that environment plays an important role in shaping galaxy
properties (e.g., \cite{Einasto74}).  The distribution of the galaxies
in the Local Group is illustrated in, e.g., \cite{Grebel97, Grebel99}.

More than half of the Local Group's galaxies are dSph galaxies.  At
present, we know of 22 dSphs in the Local Group.   Three of the M31
dSph companions were only detected and confirmed approximately seven
years ago (\cite{Armand98, Armand99, KK99, GreGuh99}).  Two more M31
dSph companions were discovered during the past two years in the Sloan
Digital Sky Survey (SDSS) and confirmed through follow-up observations
by \cite{Zucker04a, Harbeck05, Zucker06a}. They are among the least
luminous, lowest surface brightness dwarfs known.  Altogether eight
confirmed or likely dSph companions of M31 are now known.  Additional
faint dSphs may yet to be found, whereas other features (e.g.,
\cite{Morrison03, Zucker04b}) are likely part of the giant stellar
stream around M31 (\cite{Ibata04, Fardal06}).  Recent additions to the
Milky Way dSph census include three extremely low surface brightness
dwarfs also found the SDSS (\cite{Willman05, Belokurov06a,
Zucker06b}), which increases the number of Galactic early-type
companions to 12.  The new dSph Boo is the faintest galaxy known so
far with $M_V \sim -5.7$ (\cite{Belokurov06a}), and the new CVn dSph
seems to have the lowest surface brightness of any of these dwarfs
($\mu_V \sim 28$ mag arcsec$^{-2}$; \cite{Zucker06b}).  The recent
detections beg the question whether there is a lower-mass cut-off for
dwarf galaxies with luminous baryonic matter.  Finally, some seven
years ago a new {\em isolated} Local Group dSph was discovered by
\cite{Whiting99}, raising the number of these rare objects to two.
Isolated dSphs are of particular interest since they seem to defy
commonly held ideas about the creation of dSphs through, e.g., tidal
or ram pressure stripping by massive galaxies.  On the other hand,
since we do not know the orbits of these more distant dSphs, we cannot
exclude that they may once have had close encounters with the massive
galaxies.

All recent galaxy discoveries added more objects to the faint end of
the Local Group's luminosity function.  It is generally believed that
the census is fairly complete for brighter galaxies (with the possible
exception of the zone of avoidance.)  Additional very low surface
brightness objects may be uncovered in the coming years using large
imaging sky surveys.  We may expect that nearby groups also host a
number of comparatively faint objects that have escaped detection so
far.  In spite of the recent impressive improvements in the census of
nearby galaxies (see, e.g., \cite{KKSG00, Kniazev04b} and references
therein) the numbers are still far too small to resolve the
substructure crisis.  It seems doubtful at present that new searches
will add the hundreds of objects required to solve this problem -- if
the missing dark matter halos contain luminous matter as well.

Another promising avenue is the search for stellar streams in the
halos of massive galaxies (e.g., \cite{Bullock01, Newberg02,
Newberg03, Yanny03, Martin04, Rocha04, Juric05, Pena05, Wyse06,
Belokurov06b}), which may ultimately permit us to constrain the number
of past accretion events, esp.\ once full phase-space information is
added -- one of the objectives of ESA's Gaia satellite mission
(\cite{Perryman01}).

Interestingly, the companions of the Milky Way do not seem to show a
random distribution. Instead, they appear to be arranged along one or
two polar great planes (e.g., \cite{Kunkel79, Lynden82, Majewski94,
Fusi95, Kroupa05}).  An investigation of the satellite distribution
around M31 revealed that most of its early-type companions lie along
one single, almost polar great plane with high statistical
significance (\cite{GreKolBra99, KG06}, see also \cite{McConn06}).  As
shown by \cite{KG06}, this plane does not coincide with the great
stellar stream around M31.  It is unclear whether these apparent
satellite anisotropies have a deeper physical significance.  If the
satellites orbit within the planes -- a requirement if the planes have
any physical reality -- then the satellites may have formed following
interactions or the break-up of a more massive progenitor (see, e.g.,
\cite{Kunkel79, Lynden82, Sawa05, Bournaud06}).  Another possibility
is that the planes are indicators of a prolate dark halo of the Milky
Way and M31 (e.g., \cite{Hartwick00, Navarro04, Zentner05}).  Or we
may be seeing the left-over effects of matter accretion along
large-scale dark matter filaments (see \cite{Knebe04, Zentner05,
Libeskind05} for a more detailed discussion).  \cite{KG06} found that
the M31 polar plane includes also M33 and points in the direction of
the M81 group.  -- In the absence of actual orbits it remains
difficult to assess the meaning of the observed polar alignments, if
any.  For the M31 subsystem the planar alignment comprises only a
subsample of the total number of its satellites (\cite{KG06}), which
may lend support to the suggestion by \cite{Zentner05} that such
distributions may also originate from random samples.  Unfortunately,
the distance uncertainties for galaxies in nearby groups are still too
large, preventing us from conducting a similar study there.

\subsection{A few remarks about star formation histories}

The star formation histories of the galaxies of the Local Group have
been reviewed fairly frequently (e.g., \cite{vdB00, Grebel00,
Grebel05}).  Here we will only summarize some interesting
characteristics that have emerged in recent years.

Old populations are ubiquitous in the galaxies of the Local Group, but
their fractions vary (\cite{Grebel00, GreGal04}).  For instances,
horizontal branches are an unmistakable tracer of old populations and
have been identified in all Local Group galaxies with sufficiently
deep photometry.  Here ``old'' refers to old Population II stars with
ages $> 10$ Gyr.  We have no evidence for the existence of possible
Population III stars in external galaxies.  For galaxies with deep
main-sequence photometry (mainly Galactic populations and satellites
of the Milky Way) there is evidence for a common epoch of early
(Population II) star formation.  Within the currently available
accuracy of the age dating techniques ($<$ 1 Gyr), the oldest
detectable populations in the Milky Way and its surroundings are
coeval (\cite{GreGal04}).  Moreover, there is no evidence for the
suppression of star formation in the low-mass galaxies of the Local
Group after the end of reionization (\cite{GreGal04}), contrasting
predictions of certain cosmological models that suggest that low-mass
halos experience complete photo-evaporation (\cite{Barkana99,
Ferrara00, Susa04}).

All Local Group galaxies show evidence for extended star formation
histories, but no two galaxies share the same star formation history,
not even within the same morphological type (e.g., \cite{Grebel97}).
The spirals, irregulars, and dIrrs all show evidence for active star
formation over a Hubble time.  In dwarf galaxies star formation
proceeded largely continuously with amplitude variations.  Strongly
episodic star formation with long pauses in between is only seen in
the Carina dSph (e.g., \cite{Monelli03}).  Even dSph galaxies with
seemingly entirely old populations show large abundance spreads of
more than 1 dex in [Fe/H] (e.g., \cite{Shetrone01}), which require
star formation and enrichment over several Gyr (\cite{Ikuta02,
Marcolini06, Fenner06}).  In contrast to many dIrrs, which are able to
continue to form stars for another Hubble time (\cite{Hunter97}),
dSphs do not show any ongoing star formation activity and appear to be
devoid of gas (\cite{Gallagher03} and references therein).  This is
unexpected, since even simple mass loss from red giants should lead to
detectable amounts of gas (e.g., \cite{Young00}).  The Milky Way is
surrounded by several dSphs with substantial intermediate-age (2 to 10
Gyr) star formation activity.  One of the Galactic companions, the
Fornax dSph, even ended its star formation activity only as recently
as about 200 Myr ago (\cite{GreSte99}).  For a more detailed
discussion of the gas loss problem see \cite{GGH03}.

Population gradients are obvious in spiral galaxies, but also many
dwarf galaxies show spatial variations in their star formation
histories (e.g., \cite{Grebel99, Grebel00}).  Essentially all galaxies
have extended ``halos'' of old Population II stars (e.g.,
\cite{Minniti99}).  DIrrs tend to also show extended intermediate-age
populations as traced by, e.g., carbon stars (e.g., \cite{Letarte02}).
Generally, the density distributions of different populations in
irregular galaxies become increasingly more regular and extended with
increasing age (e.g., \cite{Cioni00, Zaritsky00}), whereas the many
scattered young star-forming regions are responsible for the irregular
appearance.  Large star-forming regions may remain active for several
10 to 100 Myr, revealing a complex substructure in ages (e.g.,
\cite{Dohm97, GreBra98, Walborn99, GreChu00, Dieball00, Dieball02}).

DE and dSph galaxies have usually experienced continuous star
formation with decreasing intensity.  Star formation tends to be
longest-lasting in the centers of these galaxies, and in a number of
cases radial age (and possibly metallicity) gradients are observed
(see \cite{Grebel97, Stetson98, Hurley99, DaCosta00, Sarajedini02,
Tolstoy04, Koch06} for individual cases and \cite{Harbeck01} for a
comprehensive study).  Subpopulations may also be asymmetrically
distributed (e.g., \cite{Stetson98}).  Overall, the Galactic dSph
companions show a trend for increased intermediate-age population
fractions with increasing Galactocentric distance, possibly indicating
that star-forming material might have been removed earlier from closer
companions (\cite{vdB94, Grebel97}).  However, neither the fairly
isolated Local Group dSphs Cetus and Tucana fit this pattern, nor do
the M31 low-surface-brightness companions (\cite{GGH06}).

In the past star formation histories of galaxies with resolved stellar
populations were primarily based on photometry, and on the modelling
of the observed color-magnitude diagrams (as pioneered by
\cite{Tosi91, Greggio93}).  Limitations of modelling techniques are
discussed in \cite{Gallart05}.  Additional complications may arise
from rotation, multiplicity, unrecognized extinction effects, and
transformation problems (e.g., \cite{GR95, Roberts95, GRB96,
Zaritsky02, Koch04, Girardi04, Zaritsky04}).  While in many cases no
other information but photometric data are available, information on
special types of stars can provide important additional clues (e.g.,
\cite{Grebel97, Grebel99, Borissova00, Battinelli00, Dolphin02,
Harbeck04, Harbeck05}).  Special types of stars can be identified
photometrically if, e.g., they exhibit unique spectral features that
can be traced by special filter combinations (e.g., \cite{Cook86,
Cook89, GRdB92, Grebel97b, Majewski00, Nikolaev00, Keller01, Cioni03,
Helmi03, Harbeck04, Kerschbaum04, Kniazev06}) or if they are variable
(e.g., \cite{Mateo95, Kovacs95, Mateo98, Kovacs01, Dolphin02,
Pritzl05b}).

Photometric systems with improved metallicity sensitivity have been
employed to break the photometric age-metallicity degeneracy (e.g.,
Washington and Str\"omgren photometry, \cite{Geisler86, Geisler91,
GR92, GRvdR94, Hilker00, Cole00}).  However, this degeneracy can best
be addressed spectroscopically. This is now a realistic prospect for
nearby galaxies (e.g., \cite{Cote99, Cole00, Guha00, Shetrone01,
Pont04, Tolstoy04, Koch06, Koch06b}), largely thanks to the routine
availability of optical 8 to 10-m telescopes.

\subsection{A few remarks on abundances}

While the absorption-line measurements for stellar abundances tend to
require large telescopes and long integration times, emission lines
can relatively easily be measured with medium-sized telescopes.  In
galaxies with active star formation the metallicity of young
populations can be determined by measuring the emission lines of
H\,{\sc ii} regions.  This has not only been done for gas-rich
galaxies such as dIrrs in the Local Group and its surroundings (e.g.,
\cite{Skillman89, Saviane02, Lee03, Kniazev05}, but is routinely
carried out also for much more distant galaxies with sufficiently
strong emission lines (e.g., \cite{Kniazev04a}).  These measurements
yield the present-day abundances, and individual element abundance
ratios give us information about the modes of star formation.

Also nebular abundances derived from planetary nebulae can be obtained
with only a modest investment of telescope time, yielding metallicity
estimates for mainly intermediate-age populations (e.g.,
\cite{Richer98, Jacoby99}), which can be age-dated within certain
limits (\cite{Kniazev05, Buzzoni06}).  Combining ages with the results
from nebular abundances as well as stellar abundance measurements
supports the existence of an age-metallicity relation, i.e.,
increasing enrichment with decreasing age (e.g., \cite{DaCosta98,
Tolstoy03, Grebel04, Kniazev05, Koch06}).  The comparison with
chemical evolution models can then provide fairly detailed information
about the enrichment history as well as the relative importance of
closed-box evolution vs.\ outflows and galactic winds (e.g.,
\cite{Mighell98, DaCosta98, Lanfranchi03, Lanfranchi04, Hensler04,
Koch06}).

For a given age, a given dIrr galaxy is usually assumed to be
chemically well-mixed and homogeneous.  However, there are indications
in a few dIrrs for a metallicity spread in populations of similar ages
based on age-datable star clusters (\cite{DaCosta02}) or on differing
nebular abundances in H\,{\sc ii} regions (\cite{Kniazev05}).  For
dSphs, we do not yet have data that would permit us to quantify
possible abundance spreads within populations of the same age, but as
noted earlier dSphs exhibit large metallicity spreads overall, even in
galaxies dominated entirely by old populations (e.g.,
\cite{Shetrone01}).

Interestingly, high-resolution studies of individual red giants in
nearby dwarf galaxies indicate that at a given [Fe/H] their
[$\alpha$/Fe] ratios are on average lower by $\sim 0.2$ dex than those
of equally metal-poor stars in the Galactic halo (e.g., \cite{Hill00,
Shetrone01, Shetrone03, Boni04, Geisler05, Pritzl05}).  This makes
these {\em present-day} dwarf galaxies -- dIrrs and dSphs alike --
unlikely major contributors of the build-up of the Galactic halo and
hence unlikely survivors of a once much more numerous population of
primary Galactic building blocks.  However, at low [Fe/H] there is
consistency with the abundance ratios observed in the Galactic halo,
leaving very early accretion as a plausible option (\cite{GGH03,
Font06}).

On global scales, galaxies tend to follow a luminosity--metallicity
relation in the sense that more luminous galaxies are more metal-rich.
This trend is seen for all Local Group galaxies, but the
metallicity--luminosity relations of dSphs and dIrrs are offset from
each other (\cite{Richer98}), an offset that persists even when
comparing the same metallicity tracers and the same populations
(\cite{GGH03}): DSphs have higher mean stellar metallicities at a
given optical luminosity, which may imply more rapid star formation
and enrichment at early times as compared to dIrrs (\cite{GGH03}).
This suggests that the old populations in dSphs and dIrrs are
intrinsically different.  Transitions from gas-rich to gas-poor dwarfs
seem plausible only for present-day dIrr/dSphs to present-day dSphs.

\subsection{A few remarks on kinematics and dark matter}

The internal kinematics of galaxies are not only a valuable tool to
differentiate their components, populations, and evolutionary
histories, but also provide information about galaxy masses.  For the
large spiral galaxies, a variety of approaches have been used to
constrain their masses, including stellar and gas kinematics, the
kinematics of star clusters, and of satellites.  These data seem to
suggest that the Milky Way and M31 may be of similar mass
(\cite{Kochanek96, Evans00, Carignan06}).  This is somewhat
unexpected considering the larger luminosity, larger bulge, larger
number of globular clusters, and larger physical size of M31.

While the disks of spiral galaxies exhibit differential rotation, the
more massive dIrrs show solid-body rotation.  Low-mass dIrrs,
dIrr/dSph and dSph galaxies are dominated by random motions and do not
appear to be supported by rotation.  DSphs may contain large amounts
of dark matter (e.g., \cite{Gallagher94}).  This is inferred from the
high velocity dispersion and the resulting high mass-to-light ratios
derived under the assumption of virial equilibrium (\cite{Mateo97}).
Indirectly, a high dark matter content is also supported by the
smooth, symmetrical morphology of some nearby dSphs (\cite{Odenk01})
and by the observed lack of a significant depth extent
(\cite{Klessen03}).  The radial velocity dispersion profiles of dSphs
tend to be flat and fall off large radii (\cite{Wilkinson04,
Wilkinson06}).  While detailled modelling is still in progress, the
current data favor that dSphs share a common halo mass scale of about
$4\cdot 10^7$ M$_{\odot}$, are dark-matter dominated, and have cored
mass distributions (\cite{Wilkinson06b}).

Evidence for ongoing accretion events in the Local Group have been
mentioned already in Section 2.1.  Stellar streams may also be
contributed by disrupted globular clusters (e.g., \cite{Odenk01b,
Odenk03, Belokurov06c, Grillmair06}) and can provide valuable
information about the shape of the massive galaxy's halo (e.g.,
\cite{Murali99, Bullock01, Dehnen04, Pena05, Johnston05, Majewski06}).
The results can be compared to the halo shape derived from other
indicators such as the halo globular cluster mass density profile,
which seems to be partly of primordial origin (\cite{Parmentier05}).
--- A well-known example of ongoing interactions is the triple system
of the Magellanic Clouds and the Milky Way with the gaseous Magellanic
Bridge and Magellanic Stream as interaction signatures.  Tidal
interactions are also apparent in the S-shaped surface density profile
of the Galactic dSph Ursa Minor (\cite{Palma03}) and in the twisted
isophotes of the M31 dE companions M32 and NGC\,205 (\cite{Choi02}).
These galaxies and other dSphs are likely to be accreted eventually.
Moreover, dwarf-dwarf interactions and interactions with gas clouds
may play an important role (e.g., \cite{Wilcots98, deBlok00,
Coleman05}). The Local Group's wealth in low-mass, gas-deficient
early-type dwarfs as well as the radius-morphology relation may both
be indicative of the environmental impact on galaxy evolution.

\section{Other nearby groups and Local Group analogs}
\label{sec:2}

In our immediate cosmic neighborhood we find poor, loose groups and
``clouds'' or filaments.  The spatial and density distribution of
nearby galaxies is beautifully illustrated in \cite{TFatlas87}.  A few
years ago we initiated a project with the Hubble Space Telescope and
ground-based telescopes to study the properties of the groups and
clouds within a $\sim 5$ Mpc around the Local Group, i.e., in the
Local Volume (\cite{Gsnap00}).  This project led to an improved galaxy
census and resulted in knowledge of the approximate distances,
luminosities, luminous stellar content, and approximate metallicities
of a large number of nearby galaxies (\cite{Ksnap99, KKSG00, Ksnap00,
KM8100, Dolphin01, KM8101, KKR2501, Kara02b, Maka02, Sarajedini02,
Kara02c, MakaM8102, Kara03b, Kara03c, Kara03d, Lee03, Kniazev03,
Maka05, Kniazev05}; see also \cite{Cote97, Caldwell98, Jerjen00,
Barazza01, Parodi02, Rejkuba06}).  Moreover, kinematic properties and
global masses were derived, improving the characterization of these
nearby groups (see also \cite{Cote00, Kara02a, Kara03a, Tully05,
Kara05}).  Many of the results presented in the following stem from
this continuing project.

\cite{Plionis04} find that galaxy groups are considerably more
elongated than galaxy clusters, a trend that becomes most pronounced
in very poor groups of galaxies.  They suggest that the poorest groups
of galaxies are still in the process of being assembled through galaxy
infall, a scenario that is supported by the observed properties of
nearby, extended galaxy ``clouds'' like Sculptor and CVn I
(\cite{Jerjen98, Kara03b, Kara03c}):  These elongated clouds show
several subclumps and are not yet in dynamical equilibrium, since
their crossing times are of the order of half a Hubble time.  They
contain mainly early-type dwarfs, and their luminosity functions show
a lack of very low-luminosity dSph galaxies.  CVn I may be even less
evolved than Scl.  CVn I, Scl, and the Local Group form a 10 Mpc
filament that appears to be driven by the free Hubble flow
(\cite{Kara03c}).

In contrast to the clouds, the nearby, poor groups are Local Group
analogs in many respects.  They tend to be dominated by two massive
galaxies: In the M81 group, the dominant galaxies are the spirals M81
and NGC\,2403 (\cite{Kara02b}).  A prominent interaction between M81,
and the irregular galaxies M82 and NGC 3077 is currently in progress
(e.g., \cite{Yun94}), giving rise to impressive starburst phenomena in
M82 (e.g., \cite{OConnell95}) and possibly to the formation of tidal
dwarfs from material torn out during past close encounters (e.g.,
\cite{MakaM8102}).  In the Centaurus A group the spiral M83 and the
peculiar elliptical Cen A (NGC 5128) dominate the mass distribution
(\cite{Kara02c}).  NGC\,5128 seems to have experienced various
accretion events in the past (e.g., \cite{Malin83, Mirabel99}).  In
the IC\,342/Maffei group the main subgroups are centered around the
spiral IC\,342 and the elliptical Maffei I (\cite{Kara03d}).  In this
group the faint galaxy census is still highly incomplete due to the
high foreground extinction (e.g., \cite{Buta99}).  In any case, all of
these nearby groups reveal a ``binary'' substructure.

Just like the Local Group, the nearby groups also exhibit an increased
frequency of early-type dwarfs, particularly of dSphs, and a
comparable degree of morphological segregation (\cite{Grebel04,
Grebel05}).  The luminosity functions show the familiar rise at the
faint end (\cite{Kara02c, Kara03b}).  The differences in the
fractional dwarf galaxy type distribution between clouds and Local
Group analogs supports the idea of morphological transformations
induced by denser environments.  In accordance with this picture,
transition-type dIrr/dSph galaxies and dIrrs tend to be found at
larger distances from massive galaxies (e.g., \cite{GGH03}), and dwarf
S0 galaxies are located in the outskirts of groups (and clusters;
e.g., \cite{Beaulieu06, Lisker06}).  \cite{Mahdavi99} note that
star-forming galaxies are more likely on radial orbits and may be
falling in for the first time, whereas the orbits of galaxies
dominated by older populations are more consistent with isotropy.

In the M81 group, the M81 subgroup and the NGC\,2403 subgroup approach
each other, similar to what is seen in the Local Group between M31 and
the Milky Way.  A third, smaller subgroup around the less luminous
spiral NGC 4236 in the M81 group appears to be currently receding from
the other two subgroups (\cite{Kara02b}) and may actually even
constitute a small group outside of the M81 group.  In the Cen\,A
group the two dominant subgroups also appear to be moving away from
each other (\cite{Kara02c}).  We have insufficient data on the other
groups. -- Where subgroups approach each other the final result is
expected to be  one single large elliptical galaxy (see also
\cite{Forbes00}) and hence a fossil group (\cite{Jones03}).

The estimated crossing times in the nearby, poor groups range from 1.8
to 5.9 Gyr with a median around $\sim 2.3$ Gyr (\cite{Kara05}),
suggesting that they are closer to reaching dynamical equilibrium than
the unevolved clouds, although they are far from ``fossilization''.
The radii of the zero-velocity surfaces (beyond which galaxies are no
longer bound) of the nearby groups and of the Local Group are of the
order of 1 Mpc, and the total masses of the groups are approximately a
few $10^{12}$ M$_{\odot}$ (\cite{Kara02b, Kara02c, Kara03d}).  Within
groups, the mass is closely correlated with the luminous matter and
concentrated at the location of the massive galaxies (in agreement
also with the findings by \cite{McKay02, Prada03}).

On larger scales, the nearby groups and clouds nicely delineate
the local large-scale structure (\cite{TFatlas87}).  The centers of
the groups lie within a narrow layer with a thickness of only $\pm
0.33$ Mpc (\cite{Kara05}).  The centroids of the groups show a small
velocity dispersion with respect to the Hubble flow.  Overall the
local Hubble flow is remarkably cold (\cite{Kara02a, Kara03a}),
indicating a low local matter density.  \cite{Maccio05} suggest that
the local Hubble flow is best fit by $\Lambda$CDM models with
signatures of the impact of dark energy already becoming noticeable at
distances $> 7$ Mpc.

\section{Summary}

In the Local Group, old populations are ubiquitous in all galaxies,
but their fractions vary.  There appears to have been a common epoch
of early (Population II) star formation.  All galaxies show evidence
for extended star formation episodes, but no two galaxies share the
same detailed star formation history.  There is no obvious cessation
of star-formation activity after re-ionization in low-mass galaxies.
The apparent correlation between increasing intermediate-age
population fraction in dSph galaxies with increasing distance from the
Milky Way is not seen in the M31 dSph companions.  The morphological
segregation and H\,{\sc i} mass -- distance correlation in the Local
Group and other nearby groups hint at the importance of environment
and interactions.  This is also indicated by the observed ongoing
interactions and by the increased fraction of low-luminosity,
early-type dwarfs in groups as compared to loose, unvirialized
``clouds'' of galaxies.  However, there is an offset in the
metallicity-luminosity relation even for old populations that shows
dSphs to be too metal-rich for their luminosity in comparison to
dIrrs, making a simple transformation from dIrrs to dSphs unlikely.
The origin and nature of dSphs remains a puzzle.  Also, the meaning
(if any) of  the seemingly anisotropic distribution of the early-type
companions of the Milky Way and of M31 along a polar great plane
remains unclear as long as accurate orbits are lacking.

Nearby groups in the Local Volume are Local Group analogs in many
respects.  Their luminosity functions, galaxy content, and fractional
type distribution are reminiscent of the Local Group. They show
morphological segregation.  Just like the Local Group, they are
typically dominated by two luminous galaxies.  The crossing times,
group masses and radii resemble those of the Local Group.  These poor
groups may ultimately evolve into fossil groups once their subgroups
merge.  This stands in contrast to nearby galaxy ``clouds'', which
have long crossing times, are extended, have few early-type dwarfs and
are dominated by gas-rich late-type galaxies.  These clouds do not
show a turn-up at the faint end of the galaxy luminosity function.
Perhaps they are groups in formation. --- The local Hubble flow is
very quiet.

\input{grebelref}

\printindex
\end{document}

\begin{table}
\centering
\caption{Please write your table caption here}
\label{tab:1}
\begin{tabular}{lll}
\hline\noalign{\smallskip}
first & second & third  \\
\noalign{\smallskip}\hline\noalign{\smallskip}
number & number & number \\
number & number & number \\
\noalign{\smallskip}\hline
\end{tabular}
\end{table}
\begin{figure}
\centering
\includegraphics[height=4cm]{figure.eps}
\caption{Please write your figure caption here}
\label{fig:1}
\end{figure}

Local Group(s)


First results from 18-22cm VLBA polarisation observations of the MOJAVE-I AGNs

Abstract: We are in the process of obtaining VLBA polarisation data for the 135 MOJAVE-I Active Galactic Nuclei at four wavelengths in the 18-22cm band. These observations will enable studies of the evolution of the intensity and magnetic-field structures of these AGN jets as they propagate from parsec to kiloparsec scales, as well as studies of the thermal plasma present in the vicinity of the jets on these scales, manifest via Faraday rotation. A wide range of other multi-wavelength studies can also be carried out using these data. Preliminary results for selected sources from the first 3 of 9 observing sessions will be presented. We aim to have 18-cm intensity and polarisation images available via the MOJAVE website within 18 months after the last observing session.

\section{Studying the Jets of AGN on scales of 10-100 pc}

Active Galactic Nuclei (AGN) are compact regions at the centres of galaxies where there is believed to be a supermassive black hole which
attracts mass to it via its gravitational pull. The mass attracted to the black hole forms an accretion disk, and as mass spirals
from the accretion disk to the black hole, it heats up and begins to emit large amounts of high-energy radiation.
It is in part this intense emission from the accretion disk that makes the AGN considered to be `active'.\\

Much of the gas from the accretion disk spirals in towards the black hole until it passes the event horizon, after which no further emission can be
 detected. However, for reasons not fully understood, but which are thought to be connected with strong magnetic fields generated by the relativistic
 plasma in the accretion disk, in about 10\% of AGN some of the gas falling into the black hole is instead accelerated to relativistic speeds and sent
 flying away from the AGN along collimated relativistic `jets', presumably along the rotational axis of the system. These jets can be detected
 because of the large amounts of radio synchrotron radiation they emit. If the jet is orientated
towards the Earth then relativistic beaming boosts the electromagnetic flux from the jet received on Earth, making the jet more prominent.\\

The polarisation of the radiation detected from the jet is of particular interest because it provides information about the underlying synchrotron magnetic field: in optically thick/thin regions the emitted polarisation is parallel/perpendicular to this magnetic field. Changes in the polarisation angle of the radiation with
frequency due to Faraday rotation can give information about both the magnetic field and thermal electron density in the immediate vicinity of the jet. Multiwavelength VLBA observations at 18-22cm probe the jet on scales from parsecs to 100s of parsecs, making them a useful tool
to study Faraday rotation. Other jet physics that can be studied using such observations include jet collimation, bending and interaction with the surrounding medium and the evolution and morphology of the magnetic field along the jet. Specific proposed magnetic field structures, such as helical jet magnetic fields (expected due to the rotation of central accretion disk, jet outflow and shocks), can also be investigated. \\

\section{The Experiment : BG 196}

We are obtaining simultaneous VLBA observations at 1358, 1430, 1493 and 1665 MHz (22.1, 21.0, 20.1 and 18.0 cm) of the 135 AGN making up the MOJAVE-1 sample~\cite{lister}. The first observations were taken in early 2010 and the observations and processing are expected to be finished in 2011-2012. These observations are designed for intensity and polarisation studies of AGN jets on scales out to 100s of parsecs. The spectral index and Faraday rotation distributions can also be studied. The first 3 epochs (A,B and C) have been calibrated, and some preliminary results are presented here.

\section{Calibration and Imaging}

\begin{figure}
\label{cal_compare}
\begin{center}
\caption{J0841+7053 (Epoch C) with EVPA corrections applied at all four frequencies.}
\subfigure[1358 MHz]
{
  \includegraphics[scale=0.3]{J0841+7053_1_POL_UNCORR_EVPACORR_MAN.PS}
}
\subfigure[1430 MHz]
{
  \includegraphics[scale=0.3]{J0841+7053_2_POL_UNCORR_EVPACORR_MAN.PS}
}
\subfigure[1493 MHz]
{
  \includegraphics[scale=0.3]{J0841+7053_3_POL_UNCORR_EVPACORR_MAN.PS}
}
\subfigure[1665 MHz]
{
  \includegraphics[scale=0.3]{J0841+7053_4_POL_UNCORR_EVPACORR_MAN.PS}
}
\end{center}
\end{figure}

The preliminary calibration and polarisation (D-term) calibration were done in the Astronomical Image Processing System (AIPS) using standard techniques. In all cases, Los Alamos was used as the reference antenna. The polarisation angles for epochs A and B were calibrated using integrated (VLA) observations of compact polarised sources. However as no electric vector position angle (EVPA) calibrator was included in epoch C, this epoch was instead calibrated using a previous VLBI observation. Table 1 gives the resulting EVPA corrections for each frequency for each epoch. For comparison we also present the EVPA correction for the same frequencies for the observations of Gabuzda and Hallahan \cite{2004}, obtained in January 2004, which were reduced using the same reference antenna. The EVPA corrections are sometimes, but not always, remarkably stable over as long as 6 years. These EVPA corrections were tested to ensure that they resulted in self-consistant images for a range of sources. One example is shown in Figure 1, which shows the calibrated polarisation images for J0841+7053 observed at Epoch C: the polarisation structures are clearly consistant across the four frequencies. Semi-automated scripts were written in AIPS to speed up the creation of intensity, polarisation, spectral index and Faraday rotation maps.

\begin{table}

\begin{center}
  \begin{tabular}{| c | c | c | c | c | }
    \hline Frequency & Jan 2004 & Feb 2010 & Mar 2010 & May 2010 \\
     MHz & BG139 & Epoch A & Epoch B & Epoch C \\
      & $\Delta \chi^{\circ}$ & $\Delta \chi^{\circ}$ & $\Delta \chi^{\circ}$ & $\Delta \chi^{\circ}$ \\
    \hline 1358 & 128 & 130 & 91 & 132 \\
    \hline 1430 & 113 & 113 & 112 & 113 \\
    \hline 1493 & 100 & 147 & 84 & 96 \\
    \hline 1665 & 82 & 91 & 49 & 37 \\
    \hline
  \end{tabular}
\caption{EVPA corrections for epochs A, B and C and the earlier experiment BG 139.}
\end{center}
\label{cal_table}

\end{table}

\section{Results}

\begin{wrapfigure}{r}{8cm}
\caption{3C 120.}
\subfigure[Faraday Corrected Polarisation Map at 1358 MHz]
{
  \includegraphics[scale=0.3,angle=-90]{3C120_POL_PAPER.PS}
  \label{t3c_fcpol}
}
\subfigure[Rotation Measure Map]
{
  \includegraphics[scale=0.3,angle=-90]{3C120_RM_PAPER.PS}
  \label{t3c_rm}
}
\subfigure[RM of 3C 120 across the jet. About 40 mas from the core.]
{
  \includegraphics[scale=0.3]{rm_all.ps}
  \label{t3c_chi2}
}
\end{wrapfigure}

Intensity, polarisation, spectral index and Faraday rotation maps have been obtained for the first three epochs. The expected wavelength dependence for Faraday rotation is

\[ \chi_{obs} = \chi_{o} + RM \lambda^{2}\]
\[ RM \propto \int n_{e} \vec{B} \cdot d \vec{l}\]
\noindent
where $n_{e}$ is the electron density, $\vec{B}$ is the magnetic field and the integral is taken along the line of sight. The amount and direction of rotation are thus determined by the electron density and the line-of-sight magnetic field. Results for two individual sources are presented below.

\subsection{3C 120}

The well known AGN 3C 120 has an extended jet with complex polarisation structure. The polarisation map of 3C 120 at 1358 Mhz presented in Figure 2(a) has been corrected for Faraday rotation. The rotation measure distribution is shown in Figure 2(b). A tendancy for the RM to be higher along the upper edge of the jet of 3C 120 is visible; one possible origin for this transverse gradient in Faraday rotation is a helical magnetic field carried by the jet: the Faraday rotation gradient is due to a systematic change in the line-of-sight magnetic field across the jet. Figure 2(c) shows a rotation measure slice across the jet about 40 mas from the core, which displays a clear systematic gradient. Previous evidence that the jet of 3C 120 carries a helical field has been presented by G\'omez et al \cite{gomez}.

\subsection{J0121+1149}

\begin{figure}
\centering
\caption{Unusual behaviour in J0121+1149. Figures 3(a) and 3(c) show the faraday corrected polarisation and RM maps during May 2010, Figures 3(b) and 3(d) show the same maps during January 2004. Figure 3(e) shows $\chi$ vs. $\lambda^{2}$ plots for the three regions during May 2010. The $\chi$ errors are comparable to the size of the symbols.}
\subfigure[]{
\includegraphics[trim = 10mm 30mm 1mm 5mm, clip,scale=0.35]{J0121+1149_POL_FINAL_MAN.PS}
}
\subfigure[]{
\includegraphics[trim = 10mm 30mm 1mm 5mm, clip,scale=0.35]{J0121+1149_POL_FIN_OLD_MAN.PS}
}
\subfigure[]{
\includegraphics[trim = 0mm 10mm 0mm 13mm, clip,scale=0.3]{J0121+1149_RM2_MAN_no_rmr.PS}
}
\subfigure[]{
\includegraphics[trim = 0mm 20mm 0mm 10mm, clip,scale=0.3]{J0121+1149_RM_OLD_MAN.PS}
}
\subfigure[]{
\includegraphics[trim = 0mm 5mm 0mm 5mm, clip,scale=0.37]{man_j0121_3g.ps}
}
\label{j0121}
\end{figure}

The Faraday rotation map of this source, shown in Figure~\ref{j0121}, displays 3 distinct regions with different characteristic RM values: the core, inner jet and outer jet. After correcting for this Faraday rotation, the magnetic field structure is complex, but appears to become longitudinal in the outer jet. Figure~\ref{j0121} also shows a comparison of the polarisation images obtained at epoch C and in January 2004. The three distinct regions of the Faraday rotation are present at both epochs. The polarisation structure in the core and inner jet changed between the two epochs, while that in the outer jet remained constant. The $\chi$ vs. $\lambda^{2}$ plots in Figure~\ref{j0121} clearly show the linear dependence expected for Faraday rotation, with the RM decreasing markedly with distance from the core.

\section{Conclusions}

The images presented here show just two examples of the suitability of 18-22 cm VLBA observations for studying AGN jets on scales out to 10s or 100s of parsecs. We have worked out a scheme for calibrating and imaging the data for BG196 in a timely fashion, including the development of semi-automated scripts. We plan to provide calibrated data and images to the community within 18 months after the last observations, via the MOJAVE website \cite{lister}.

\section*{Acknowlegements}

This work was supported by the Irish Research Council for Science, Engineering and Technology (IRCSET).

Cosmology and Nongalactic Astrophysics (astro-ph.CO)