The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. We measured the tilt angle of the primary bipoles in 153 emerging active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 hours to define the polarities and measure the tilt angle up to four days after emergence. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that ARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Our finding that active regions emerge with an east-west alignment is surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation.

Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker et al. 2016). We first measured the position of each of the polarities eight hours after emergence and tracked their location forwards and backwards in time. We find that, on average, the polarities emerge with an east-west orientation and the separation speed between the polarities increases. At about 0.1 days after emergence, the average separation speed reaches a peak value of 229 +/- 11 m/s, and then starts to decrease, and about 2.5 days after emergence the polarities stop separating. We also find that the separation and the separation speed in the east-west direction are systematically larger for active regions with higher flux. Our results reveal two phases of the emergence process defined by the rate of change of the separation speed as the polarities move apart. Phase 1 begins when the opposite polarity pairs first appear at the surface, with an east-west alignment and an increasing separation speed. We define Phase 2 to begin when the separation speed starts to decrease, and ends when the polarities have stopped separating. This is consistent with the picture of Chen, Rempel, & Fan (2017): the peak of a flux tube breaks through the surface during Phase 1. During Phase 2 the magnetic field lines are straightened by magnetic tension, so that the polarities continue to move apart, until they eventually lie directly above their anchored subsurface footpoints.

We explore the general properties of near-surface flows around solar active regions. Helioseismic holography is applied to HMI Dopplergrams yielding nearly 5000 flow measurements of 336 unique active regions observed by the Solar Dynamics Observatory between 2010 and 2014. Ensemble averages of the flows, over subsets of regions sorted on the basis of magnetic flux, are performed. These averages show that converging flows, with speeds about 10 m/s and extending up to 10 degrees from the active region centers, are prevalent and have similar properties for all regions with magnetic flux above 10^21 Mx. Retrograde flows are also detected, with amplitudes around 10 m/s, which predominantly, but not exclusively, flank the polar side of the active regions. We estimate the expected contribution of these active-region flows to longitudinal averages of zonal and meridional flows and demonstrate the plausibility that they are responsible for at least some component of the time-varying global-scale flows. The reliability of our flow determination is tested using publicly available MHD simulations of both quiet-Sun convection and of a sunspot. While validating the overall methodology in general, the sunspot simulation demonstrates the presence of artifacts which may compromise quantitative flow inferences from some helioseismic measurements.

Here we use synthetic data to explore the performance of forward models and inverse methods for helioseismic holography. Specifically, this work presents the first comprehensive test of inverse modeling for flows using lateral-vantage (deep-focus) holography. We derive sensitivity functions in the Born approximation. We then use these sensitivity functions in a series of forward models and inversions of flows from a publicly available magnetohydrodynamic quiet-Sun simulation. The forward travel times computed using the kernels generally compare favorably with measurements obtained by applying holography, in a lateral-vantage configuration, on a 15-hour time series of artificial Dopplergrams extracted from the simulation. Inversions for the horizontal flow components are able to reproduce the flows in the upper 3Mm of the domain, but are compromised by noise at greater depths.

Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the propagation of a stochastic wave field through sunspot models with different properties, accounting for changes in the Wilson depression between 250 and 550 km and in the photospheric umbral magnetic field between 1500 and 3500 G. The results show that travel-time shifts at frequencies above approximately 3.50 mHz (depending on the phase-speed filter) are insensitive to the magnetic field strength. The travel time of these waves is determined exclusively by the Wilson depression and sound-speed perturbation. The travel time of waves with lower frequencies is affected by the direct effect of the magnetic field, although photospheric field strengths below 1500 G do not leave a significant trace on the travel-time measurements. These results could potentially be used to develop simplified travel-time inversion methods.

X-ray polarimetry has seen a growing interest in recent years. Improvements in detector technology and focusing X-ray optics now enable sensitive astrophysical X-ray polarization measurements. These measurements will provide new insights into the processes at work in accreting black holes, the emission of X-rays from neutron stars and magnetars, and the structure of AGN jets. X-Calibur is a balloon-borne hard X-ray scattering polarimeter. An X-ray mirror with a focal length of 8 m focuses X-rays onto the detector, which consists of a plastic scattering element surrounded by Cadmium-Zinc-Telluride detectors, which absorb and record the scattered X-rays. Since X-rays preferentially scatter perpendicular to their polarization direction, the polarization properties of an X-ray beam can be inferred from the azimuthal distribution of scattered X-rays. A close alignment of the X-ray focal spot with the center of the detector is required in order to reduce systematic uncertainties and to maintain a high photon detection efficiency. This places stringent requirements on the mechanical and thermal stability of the telescope structure. During the flight on a stratospheric balloon, X-Calibur makes use of the Wallops Arc-Second Pointer (WASP) to point the telescope at astrophysical sources. In this paper, we describe the design, construction, and test of the telescope structure, as well as its performance during a 25-hour flight from Ft. Sumner, New Mexico. The carbon fiber-aluminum composite structure met the requirements set by X-Calibur and its design can easily be adapted for other types of experiments, such as X-ray imaging or spectroscopic telescopes.

Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.

Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted the sample to active regions that were visible in the continuum and emerged into quiet Sun largely avoiding pre-existing magnetic regions. As a reference data set we paired a control region (CR), with the same latitude and distance from central meridian, with each emerging active region (EAR). We call this data set, which is currently comprised of 105 emerging active regions observed between May 2010 and November 2012, the SDO Helioseismic Emerging Active Region (SDO/HEAR) survey. To demonstrate the utility of a data set of a large number of emerging active regions, we measure the relative east-west velocity of the leading and trailing polarities from the line-of-sight magnetogram maps during the first day after emergence. The latitudinally averaged line-of-sight magnetic field of all the EARs shows that, on average, the leading (trailing) polarity moves in a prograde (retrograde) direction with a speed of 121 +/- 22 m/s (-70 +/- 13 m/s) relative to the Carrington rotation rate in the first day. However, relative to the differential rotation of the surface plasma, the east-west velocity is symmetric, with a mean of 95 +/- 13 m/s. The SDO/HEAR data set will not only be useful for helioseismic studies, but will also be useful to study other features such as the surface magnetic field evolution of a large sample of EARs.

Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave propagation through sunspots. Helioseismic holography measurements made from the resulting simulated wavefields show qualitative agreement with observations of real sunspots. We use additional numerical experiments to determine, separately, the influence of the thermal structure of the sunspot and the direct effect of the sunspot magnetic field. We use the ray approximation to show that the travel-time shifts in the thermal (non-magnetic) sunspot model are primarily produced by changes in the wave path due to the Wilson depression rather than variations in the wave speed. This shows that inversions for the subsurface structure of sunspots must account for local changes in the density. In some ranges of horizontal phase speed and frequency there is agreement (within the noise level in the simulations) between the travel times measured in the full magnetic sunspot model and the thermal model. If this conclusion proves to be robust for a wide range of models, it would suggest a path towards inversions for sunspot structure.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

The Prime Focus Spectrograph (PFS) is a new optical/near-infrared multi-fiber spectrograph design for the prime focus of the 8.2m Subaru telescope. PFS will cover 1.3 degree diameter field with 2394 fibers to complement the imaging capability of Hyper SuprimeCam (HSC). The prime focus unit of PFS called Prime Focus Instrument (PFI) provides the interface with the top structure of Subaru telescope and also accommodates the optical bench in which Cobra fiber positioners are located. In addition, the acquisition and guiding cameras (AGCs), the optical fiber positioner system, the cable wrapper, the fiducial fibers, illuminator, and viewer, the field element, and the telemetry system are located inside the PFI. The mechanical structure of the PFI was designed with special care such that its deflections sufficiently match those of the HSC's Wide Field Corrector (WFC) so the fibers will stay on targets over the course of the observations within the required accuracy. In this report, the latest status of PFI development will be given including the performance of PFI components, the setup and performance of the integration and testing equipment.

Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence.

We use helioseismic holography to study the association of shallow flows with solar flare activity in about 250 large sunspot groups observed between 2010 and 2014 with the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory. Four basic flow parameters: horizontal speed, horizontal component of divergence, vertical component of vorticity, and a vertical kinetic helicity proxy, are mapped for each active region during its passage across the solar disk. Flow indices are derived representing the mean and standard deviation of these parameters over magnetic masks and compared with contemporary measures of flare X-ray flux. A correlation exists for several of the flow indices, especially those based on the speed and the standard deviation of all flow parameters. However, their correlation with X-ray flux is similar to that observed with the mean unsigned magnetic flux density over the same masks. The temporal variation of the flow indices are studied, and a superposed epoch analysis with respect to the occurrence to 70 M and X-class flares is made. While flows evolve with the passage of the active regions across the disk, no discernible precursors or other temporal changes specifically associated with flares are detected.

Installation of the ARIANNA Hexagonal Radio Array (HRA) on the Ross Ice Shelf of Antarctica has been completed. This detector serves as a pilot program to the ARIANNA neutrino telescope, which aims to measure the diffuse flux of very high energy neutrinos by observing the radio pulse generated by neutrino-induced charged particle showers in the ice. All HRA stations ran reliably and took data during the entire 2014-2015 austral summer season. A new radio signal direction reconstruction procedure is described, and is observed to have a resolution better than a degree. The reconstruction is used in a preliminary search for potential neutrino candidate events in the data from one of the newly installed detector stations. Three cuts are used to separate radio backgrounds from neutrino signals. The cuts are found to filter out all data recorded by the station during the season while preserving 85.4% of simulated neutrino events that trigger the station. This efficiency is similar to that found in analyses of previous HRA data taking seasons.

The ARIANNA collaboration completed the installation of the hexagonal radio array (HRA) in December 2014, serving as a pilot program for a planned high energy neutrino telescope located about 110 km south of McMurdo Station on the Ross Ice Shelf near the coast of Antarctica. The goal of ARIANNA is to measure both diffuse and point fluxes of astrophysical neutrinos at energies in excess of 1016 eV. Upgraded hardware has been installed during the 2014 deployment season and stations show a livetime of better than 90% between commissioning and austral sunset. Though designed to observe radio pulses from neutrino interactions originating within the ice below each detector, one station was modified to study the low-frequency environment and signals from above. We provide evidence that the HRA observed both continuous emission from the Galaxy and a transient solar burst. Preliminary work on modeling the (weak) Galactic signal confirm the absolute sensitivity of the HRA detector system.

The ARIANNA experiment seeks to observe the diffuse flux of neutrinos in the 10^8 - 10^10 GeV energy range using a grid of radio detectors at the surface of the Ross Ice Shelf of Antarctica. The detector measures the coherent Cherenkov radiation produced at radio frequencies, from about 100 MHz to 1 GHz, by charged particle showers generated by neutrino interactions in the ice. The ARIANNA Hexagonal Radio Array (HRA) is being constructed as a prototype for the full array. During the 2013-14 austral summer, three HRA stations collected radio data which was wirelessly transmitted off site in nearly real-time. The performance of these stations is described and a simple analysis to search for neutrino signals is presented. The analysis employs a set of three cuts that reject background triggers while preserving 90% of simulated cosmogenic neutrino triggers. No neutrino candidates are found in the data and a model-independent 90% confidence level Neyman upper limit is placed on the all flavor neutrino+antineutrino flux in a sliding decade-wide energy bin. The limit reaches a minimum of 1.9x10^-23 GeV^-1 cm^-2 s^-1 sr^-1 in the 10^8.5 - 10^9.5 GeV energy bin. Simulations of the performance of the full detector are also described. The sensitivity of the full ARIANNA experiment is presented and compared with current neutrino flux models.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 {\mu}m to 1.26 {\mu}m, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 {\mu}m to 0.89 {\mu}m will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase.

FOCCoS, Fiber Optical Cable and Connector System, has the main function of capturing the direct light from the focal plane of Subaru Telescope using optical fibers, each one with a microlens in its tip, and conducting this light through a route containing connectors to a set of four spectrographs. The optical fiber cable is divided in 3 different segments called Cable A, Cable B and Cable C. Multi-fibers connectors assure precise connection among all optical fibers of the segments, providing flexibility for instrument changes. To assure strong and accurate connection, these sets are arranged inside two types of assemblies: the Tower Connector, for connection between Cable C and Cable B; and the Gang Connector, for connection between Cable B and Cable A. Throughput tests were made to evaluate the efficiency of the connections. A lifetime test connection is in progress. Cable C is installed inside the PFI, Prime Focus Instrument, where each fiber tip with a microlens is bonded to the end of the shaft of a 2-stage piezo-electric rotatory motor positioner; this assembly allows each fiber to be placed anywhere within its patrol region, which is 9.5mm diameter.. Each positioner uses a fiber arm to support the ferrule, the microlens, and the optical fiber. 2400 of these assemblies are arranged on a motor bench plate in a hexagonal-closed-packed disposition.

The Prime Focus Spectrograph (PFS) is a new optical/near-infrared multi-fiber spectrograph design for the prime focus of the 8.2m Subaru telescope. PFS will cover 1.3 degree diameter field with 2394 fibers to complement the imaging capability of Hyper SuprimeCam (HSC). The prime focus unit of PFS called Prime Focus Instrument (PFI) provides the interface with the top structure of Subaru telescope and also accommodates the optical bench in which Cobra fiber positioners are located. In addition, the acquisition and guiding (A&G) cameras, the optical fiber positioner system, the cable wrapper, the fiducial fibers, illuminator, and viewer, the field element, and the telemetry system are located inside the PFI. The mechanical structure of the PFI was designed with special care such that its deflections sufficiently match those of the HSC Wide Field Corrector (WFC) so the fibers will stay on targets over the course of the observations within the required accuracy.

The Cobra fiber positioner is being developed by the California Institute of Technology (CIT) and the Jet Propulsion Laboratory (JPL) for the Prime Focus Spectrograph (PFS) instrument that will be installed at the Subaru Telescope on Mauna Kea, Hawaii. PFS is a fiber fed multi-object spectrometer that uses an array of Cobra fiber positioners to rapidly reconfigure 2394 optical fibers at the prime focus of the Subaru Telescope that are capable of positioning a fiber to within 5um of a specified target location. A single Cobra fiber positioner measures 7.7mm in diameter and is 115mm tall. The Cobra fiber positioner uses two piezo-electric rotary motors to move a fiber optic anywhere in a 9.5mm diameter patrol area. In preparation for full-scale production of 2550 Cobra positioners an Engineering Model (EM) version was developed, built and tested to validate the design, reduce manufacturing costs, and improve system reliability. The EM leveraged the previously developed prototype versions of the Cobra fiber positioner. The requirements, design, assembly techniques, development testing, design qualification and performance evaluation of EM Cobra fiber positioners are described here. Also discussed is the use of the EM build and test campaign to validate the plans for full-scale production of 2550 Cobra fiber positioners scheduled to begin in late-2014.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms.

The subsurface properties of active regions prior to their appearance at the solar surface may shed light on the process of active region formation. Helioseismic holography has been applied to samples taken from two populations of regions on the Sun (pre-emergence and without emergence), each sample having over 100 members, that were selected to minimize systematic bias, as described in Paper I (Leka et al., 2012). Paper II (Birch et al., 2012) showed that there are statistically significant signatures in the average helioseismic properties that precede the formation of an active region. This paper describes a more detailed analysis of the samples of pre-emergence regions and regions without emergence, based on discriminant analysis. The property that is best able to distinguish the populations is found to be the surface magnetic field, even a day before the emergence time. However, after accounting for the correlations between the surface field and the quantities derived from helioseismology, there is still evidence of a helioseismic precursor to active region emergence that is present for at least a day prior to emergence.

We report on average subsurface properties of pre-emerging active regions as compared to areas where no active region emergence was detected. Helioseismic holography is applied to samples of the two populations (pre-emergence and without emergence), each sample having over 100 members, which were selected to minimize systematic bias, as described in Leka et al. We find that there are statistically significant signatures (i.e., difference in the means of more than a few standard errors) in the average subsurface flows and the apparent wave speed that precede the formation of an active region. The measurements here rule out spatially extended flows of more than about 15 m/s in the top 20 Mm below the photosphere over the course of the day preceding the start of visible emergence. These measurements place strong constraints on models of active region formation.

This first paper in a series describes the design of a study testing whether pre-appearance signatures of solar magnetic active regions were detectable using various tools of local helioseismology. The ultimate goal is to understand flux-emergence mechanisms by setting observational constraints on pre-appearance subsurface changes, for comparison with results from simulation efforts. This first paper provides details of the data selection and preparation of the samples, each containing over 100 members, of two populations: regions on the Sun that produced a numbered NOAA active region, and a "control" sample of areas that did not. The seismology is performed on data from the GONG network; accompanying magnetic data from SOHO/MDI are used for co-temporal analysis of the surface magnetic field. Samples are drawn from 2001 -- 2007, and each target is analyzed for 27.7 hr prior to an objectively determined time of emergence. The results of two analysis approaches are published separately: one based on averages of the seismology- and magnetic-derived signals over the samples, another based on Discriminant Analysis of these signals, for a statistical test of detectable differences between the two populations. We include here descriptions of a new potential-field calculation approach and the algorithm for matching sample distributions over multiple variables. We describe known sources of bias and the approaches used to mitigate them. We also describe unexpected bias sources uncovered during the course of the study and include a discussion of refinements that should be included in future work on this topic.

Supergranulation is one of the most visible length scales of solar convection and has been studied extensively by local helioseismology. We use synthetic data computed with the Seismic Propagation through Active Regions and Convection (SPARC) code to test regularized-least squares (RLS) inversions of helioseismic holography measurements for a supergranulation-like flow. The code simulates the acoustic wavefield by solving the linearized three-dimensional Euler equations in Cartesian geometry. We model a single supergranulation cell with a simple, axisymmetric, mass-conserving flow. The use of simulated data provides an opportunity for direct evaluation of the accuracy of measurement and inversion techniques. The RLS technique applied to helioseismic-holography measurements is generally successful in reproducing the structure of the horizontal flow field of the model supergranule cell. The errors are significant in horizontal-flow inversions near the top and bottom of the computational domain as well as in vertical-flow inversions throughout the domain. We show that the errors in the vertical velocity are due largely to cross talk from the horizontal velocity.

We use a publicly available numerical wave-propagation simulation of Hartlep et al. 2011 to test the ability of helioseismic holography to detect signatures of a compact, fully submerged, 5% sound-speed perturbation placed at a depth of 50 Mm within a solar model. We find that helioseismic holography as employed in a nominal "lateral-vantage" or "deep-focus" geometry employing quadrants of an annular pupil is capable of detecting and characterizing the perturbation. A number of tests of the methodology, including the use of a plane-parallel approximation, the definition of travel-time shifts, the use of different phase-speed filters, and changes to the pupils, are also performed. It is found that travel-time shifts made using Gabor-wavelet fitting are essentially identical to those derived from the phase of the Fourier transform of the cross-covariance functions. The errors in travel-time shifts caused by the plane-parallel approximation can be minimized to less than a second for the depths and fields of view considered here. Based on the measured strength of the mean travel-time signal of the perturbation, no substantial improvement in sensitivity is produced by varying the analysis procedure from the nominal methodology in conformance with expectations. The measured travel-time shifts are essentially unchanged by varying the profile of the phase-speed filter or omitting the filter entirely. The method remains maximally sensitive when applied with pupils that are wide quadrants, as opposed to narrower quadrants or with pupils composed of smaller arcs. We discuss the significance of these results for the recent controversy regarding suspected pre-emergence signatures of active regions.

The Prime Focus Spectrograph (PFS) of the Subaru Measurement of Images and Redshifts (SuMIRe) project has been endorsed by Japanese community as one of the main future instruments of the Subaru 8.2-meter telescope at Mauna Kea, Hawaii. This optical/near-infrared multi-fiber spectrograph targets cosmology with galaxy surveys, Galactic archaeology, and studies of galaxy/AGN evolution. Taking advantage of Subaru's wide field of view, which is further extended with the recently completed Wide Field Corrector, PFS will enable us to carry out multi-fiber spectroscopy of 2400 targets within 1.3 degree diameter. A microlens is attached at each fiber entrance for F-ratio transformation into a larger one so that difficulties of spectrograph design are eased. Fibers are accurately placed onto target positions by positioners, each of which consists of two stages of piezo-electric rotary motors, through iterations by using back-illuminated fiber position measurements with a wide-field metrology camera. Fibers then carry light to a set of four identical fast-Schmidt spectrographs with three color arms each: the wavelength ranges from 0.38 {\mu}m to 1.3 {\mu}m will be simultaneously observed with an average resolving power of 3000. Before and during the era of extremely large telescopes, PFS will provide the unique capability of obtaining spectra of 2400 cosmological/astrophysical targets simultaneously with an 8-10 meter class telescope. The PFS collaboration, led by IPMU, consists of USP/LNA in Brazil, Caltech/JPL, Princeton, & JHU in USA, LAM in France, ASIAA in Taiwan, and NAOJ/Subaru.

The Prime Focus Spectrograph (PFS) is a fiber fed multi-object spectrometer for the Subaru Telescope that will conduct a variety of targeted surveys for studies of dark energy, galaxy evolution, and galactic archaeology. The key to the instrument is a high density array of fiber positioners placed at the prime focus of the Subaru Telescope. The system, nicknamed "Cobra", will be capable of rapidly reconfiguring the array of 2394 optical fibers to the image positions of astronomical targets in the focal plane with high accuracy. The system uses 2394 individual "SCARA robot" mechanisms that are 7.7mm in diameter and use 2 piezo-electric rotary motors to individually position each of the optical fibers within its patrol region. Testing demonstrates that the Cobra positioner can be moved to within 5{\mu}m of an astronomical target in 6 move iterations with a success rate of 95 per cent. The Cobra system is a key aspect of PFS that will enable its unprecedented combination of high-multiplex factor and observing efficiency on the Subaru telescope. The requirements, design, and prototyping efforts for the fiber positioner system for the PFS are described here as are the plans for modular construction, assembly, integration, functional testing, and performance validation.

The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models which differ in radius and total magnetic flux. A quiet Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m=0 and m=+/-1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p_1-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the emsemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase-shifts with the azimuthal order $m$ appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

The cause of enhanced acoustic power surrounding active regions, the acoustic halo, is not as yet understood. We explore the properties of the enhanced acoustic power observed near disk center from 21 to 27 January 2002, including AR 9787. We find that (i) there exists a strong correlation of the enhanced high frequency power with magnetic-field inclination, with greater power in more horizontal fields, (ii) the frequency of the maximum enhancement increases along with magnetic field strength, and (iii) the oscillations contributing to the halos show modal ridges which are shifted to higher wavenumber at constant frequency in comparison to the ridges of modes in the quiet-Sun.

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic field and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phase-sensitive holography demonstrates that they are significantly affected in a strong, moderately inclined magnetic field. This research indicates that the effects of the surface magnetic field are potentially very significant for local helioseismic analysis of active regions.

We perform helioseismic holography to assess the noise in p-mode travel-time shifts which would form the basis of inferences of large-scale flows throughout the solar convection zone. We also derive the expected travel times from a parameterized return (equatorward) flow component of the meridional circulation at the base of the convection zone from forward models under the assumption of the ray and Born approximations. From estimates of the signal-to-noise ratio for measurements focused near the base of the convection zone, we conclude that the helioseismic detection of the deep meridional flow including the return component may not be possible using data spanning an interval less than a solar cycle.

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled assuming the surface motion is elliptical, and we explore the properties of the elliptical motion on the magnetic field inclination. We also demonstrate that the phase shift of the outward propagating waves is opposite to the phase shift of the inward propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more inclined fields.

We present a comprehensive set of observations of the interaction of p-mode oscillations with sunspots using surface-focused seismic holography. Maps of travel-time shifts, relative to quiet-Sun travel times, are shown for incoming and outgoing p modes as well as their mean and difference. We compare results using phase-speed filters with results obtained with filters that isolate single p-mode ridges, and further divide the data into multiple temporal frequency bandpasses. The f mode is removed from the data. The variations of the resulting travel-time shifts with magnetic-field strength and with the filter parameters are explored. We find that spatial averages of these shifts within sunspot umbrae, penumbrae, and surrounding plage often show strong frequency variations at fixed phase speed. In addition, we find that positive values of the mean and difference travel-time shifts appear exclusively in waves observed with phase-speed filters that are dominated by power in the low-frequency wing of the p1 ridge. We assess the ratio of incoming to outgoing p-mode power using the ridge filters and compare surface-focused holography measurements with the results of earlier published p-mode scattering measurements using Fourier-Hankel decomposition.

We perform helioseismic holography on realistic solar convection simulations and compare the observed travel-time perturbations with the expected travel times from the horizontal flows in the simulations computed from forward models under the assumption of the Born approximation. We demonstrate reasonable agreement between the observed and model travel times which reinforces the validity of helioseismic holography in the detection of subsurface horizontal flows. From the variation of the signal-to-noise ratio with depth, we conclude that the helioseismic detection of individual flow structures with spatial scales of supergranulation or smaller is not possible for depths below about 5 Mm below the surface over time scales less than a day. Approximately half of the observed signal originates within the first 2 Mm below the surface. A consequence of this is a rapid decrease (and reversal in some cases) of the travel-time perturbations with depth due to the contribution to the measurements of oppositely directed surface flows in neighboring convective cells. This confirms an earlier interpretation of similar effects reported from observations.

We measure the mean frequencies of acoustic waves propagating toward and away from a sunspot employing a spot-centered Fourier-Hankel decomposition of p-mode amplitudes as measured from a set of observations made at the South Pole in 1991. We demonstrate that there is a significant frequency shift between the inward and outward traveling waves which is consistent with the Doppler effect of a radial outflow from the sunspot. For p-modes of temporal frequencies of 3 mHz it is observed that the frequency shift decreases slightly with spatial frequency, for modes with degree l between 160 to 600. From the l dependence of the frequency shift, we infer that the mean radial outflow within the observed annular region (which extends between 30 and 137 Mm from the spot) increases nearly linearly with depth, reaching a magnitude of about 200 m/s at a depth of 20 Mm. This outflow exhibits properties similar to flows recently reported by Lindsey, et al. (1996) using spatially sensitive local helioseismic techniques.