We present a publicly-available toolkit of flight-proven hardware and software to retrieve 5 TB of data or small physical samples from a stratospheric balloon platform. Before launch, a capsule is attached to the balloon, and rises with it. Upon remote command, the capsule is released and descends via parachute, continuously transmitting its location. Software to predict the trajectory can be used to select a safe but accessible landing site. We dropped two such capsules from the SuperBIT telescope, in September 2019. The capsules took ~37 minutes to descend from ~30 km altitude. They drifted 32 km and 19 km horizontally, but landed within 300 m and 600 m of their predicted landing sites. We found them easily, and successfully recovered the data. We welcome interest from other balloon teams for whom the technology would be useful.

Future mm-wave and sub-mm space missions will employ large arrays of multiplexed Transition Edge Sensor (TES) bolometers. Such instruments must contend with the high flux of cosmic rays beyond our atmosphere that induce "glitches" in bolometer data, which posed a challenge to data analysis from the Planck bolometers. Future instruments will face the additional challenges of shared substrate wafers and multiplexed readout wiring. In this work we explore the susceptibility of modern TES arrays to the cosmic ray environment of space using two data sets: the 2015 long-duration balloon flight of the SPIDER cosmic microwave background polarimeter, and a laboratory exposure of SPIDER flight hardware to radioactive sources. We find manageable glitch rates and short glitch durations, leading to minimal effect on SPIDER analysis. We constrain energy propagation within the substrate through a study of multi-detector coincidences, and give a preliminary look at pulse shapes in laboratory data.

CMB full-sky temperature data show a hemispherical asymmetry in power nearly aligned with the Ecliptic. In real space, this anomaly can be quantified by the temperature variance in the northern and southern Ecliptic hemispheres, with the north displaying an anomalously low variance while the south appears consistent with expectations from the best-fitting theory, LCDM. While this is a well-established result in temperature, the low signal-to-noise ratio in current polarization data prevents a similar comparison. Even though temperature and polarization are correlated, polarization realizations constrained by temperature data show that the lack of variance is not expected to be present in polarization data. Therefore, a natural way of testing whether the temperature result is a fluke is to measure the variance of CMB polarization components. In anticipation of future CMB experiments that will allow for high-precision large-scale polarization measurements, we study how variance of polarization depends on LCDM parameters' uncertainties by forecasting polarization maps with Planck's MCMC chains. We find that, unlike temperature variance, polarization variance is noticeably sensitive to present uncertainties in cosmological parameters. This comes mainly from the current poor constraints on the reionization optical depth, tau, and the fact that tau drives variance at low multipoles. In this work we show how the variance of polarization maps generically depends on the cosmological parameters. We demonstrate how the improvement in the tau measurement seen between Planck's two latest data releases results in a tighter constraint on polarization variance expectations. Finally, we consider even smaller uncertainties on tau and how more precise measurements of tau can drive the expectation for polarization variance in a hemisphere close to that of the cosmic-variance-limited distribution.

At a fraction the total cost of an equivalent orbital mission, scientific balloon-borne platforms, operating above 99.7% of the Earth's atmosphere, offer attractive, competitive, and effective observational capabilities -- namely space-like resolution, transmission, and backgrounds -- that are well suited for modern astronomy and cosmology. SuperBIT is a diffraction-limited, wide-field, 0.5 m telescope capable of exploiting these observing conditions in order to provide exquisite imaging throughout the near-IR to near-UV. It utilizes a robust active stabilization system that has consistently demonstrated a 1 sigma sky-fixed pointing stability at 48 milliarcseconds over multiple 1 hour observations at float. This is achieved by actively tracking compound pendulations via a three-axis gimballed platform, which provides sky-fixed telescope stability at < 500 milliarcseconds and corrects for field rotation, while employing high-bandwidth tip/tilt optics to remove residual disturbances across the science imaging focal plane. SuperBIT's performance during the 2019 commissioning flight benefited from a customized high-fidelity science-capable telescope designed with exceptional thermo- and opto-mechanical stability as well as tightly constrained static and dynamic coupling between high-rate sensors and telescope optics. At the currently demonstrated level of flight performance, SuperBIT capabilities now surpass the science requirements for a wide variety of experiments in cosmology, astrophysics and stellar dynamics.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

BFORE is a high-altitude ultra-long-duration balloon mission to map the cosmic microwave background (CMB). During a 28-day mid-latitude flight launched from Wanaka, New Zealand, the instrument will map half the sky to improve measurements of the optical depth to reionization tau. This will break parameter degeneracies needed to detect neutrino mass. BFORE will also hunt for the gravitational wave B-mode signal, and map Galactic dust foregrounds. The mission will be the first near-space use of TES/mSQUID multichroic detectors (150/217 GHz and 280/353 GHz bands) with low-power readout electronics.

SPIDER is a balloon-borne instrument designed to map the polarization of the millimeter-wave sky at large angular scales. SPIDER targets the B-mode signature of primordial gravitational waves in the cosmic microwave background (CMB), with a focus on mapping a large sky area with high fidelity at multiple frequencies. SPIDER's first longduration balloon (LDB) flight in January 2015 deployed a total of 2400 antenna-coupled Transition Edge Sensors (TESs) at 90 GHz and 150 GHz. In this work we review the design and in-flight performance of the SPIDER instrument, with a particular focus on the measured performance of the detectors and instrument in a space-like loading and radiation environment. SPIDER's second flight in December 2018 will incorporate payload upgrades and new receivers to map the sky at 285 GHz, providing valuable information for cleaning polarized dust emission from CMB maps.

We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne SPIDER instrument. These FPUs are vital to SPIDER's primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a 16 x 16 grid of corrugated silicon feedhorns coupled to an array of aluminum-manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1,530 polarization sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be ~3 pW at 300 mK with a less than 6% variation across each array at one standard deviation. These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the SPIDER-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33<\ell<307$. No other limits exist over this full range of angular scales, and SPIDER improves upon the previous limit by several orders of magnitude, providing 95% C.L. constraints on $\ell (\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K ^2$ to 255 $\mu K ^2$ at 150 GHz for a thermal CMB spectrum. As linear CMB polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to obtain even stronger constraints on circular polarization.

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

BFORE is a NASA high-altitude ultra-long-duration balloon mission proposed to measure the cosmic microwave background (CMB) across half the sky during a 28-day mid-latitude flight launched from Wanaka, New Zealand. With the unique access to large angular scales and high frequencies provided by the balloon platform, BFORE will significantly improve measurements of the optical depth to reionization tau, breaking parameter degeneracies needed for a measurement of neutrino mass with the CMB. The large angular scale data will enable BFORE to hunt for the large-scale gravitational wave B-mode signal, as well as the degree-scale signal, each at the r~0.01 level. The balloon platform allows BFORE to map Galactic dust foregrounds at frequencies where they dominate, in order to robustly separate them from CMB signals measured by BFORE, in addition to complementing data from ground-based telescopes. The combination of frequencies will also lead to velocity measurements for thousands of galaxy clusters, as well as probing how star-forming galaxies populate dark matter halos. The mission will be the first near-space use of TES multichroic detectors (150/217 GHz and 280/353 GHz bands) using highly-multiplexed mSQUID microwave readout, raising the technical readiness level of both technologies.

We describe 280 GHz bolometric detector arrays that instrument the balloon-borne polarimeter SPIDER. A primary science goal of SPIDER is to measure the large-scale B-mode polarization of the cosmic microwave background in search of the cosmic-inflation, gravitational-wave signature. 280 GHz channels aid this science goal by constraining the level of B-mode contamination from galactic dust emission. We present the focal plane unit design, which consists of a 16$\times$16 array of conical, corrugated feedhorns coupled to a monolithic detector array fabricated on a 150 mm diameter silicon wafer. Detector arrays are capable of polarimetric sensing via waveguide probe-coupling to a multiplexed array of transition-edge-sensor (TES) bolometers. The SPIDER receiver has three focal plane units at 280 GHz, which in total contains 765 spatial pixels and 1,530 polarization sensitive bolometers. By fabrication and measurement of single feedhorns, we demonstrate 14.7$^{\circ}$ FHWM Gaussian-shaped beams with $<$1% ellipticity in a 30% fractional bandwidth centered at 280 GHz. We present electromagnetic simulations of the detection circuit, which show 94% band-averaged, single-polarization coupling efficiency, 3% reflection and 3% radiative loss. Lastly, we demonstrate a low thermal conductance bolometer, which is well-described by a simple TES model and exhibits an electrical noise equivalent power (NEP) = 2.6 $\times$ 10$^{-17}$ W/$\sqrt{\mathrm{Hz}}$, consistent with the phonon noise prediction.

We describe the cryogenic half-wave plate rotation mechanisms built for and used in Spider, a polarization-sensitive balloon-borne telescope array that observed the Cosmic Microwave Background at 95 GHz and 150 GHz during a stratospheric balloon flight from Antarctica in January 2015. The mechanisms operate at liquid helium temperature in flight. A three-point contact design keeps the mechanical bearings relatively small but allows for a large (305 mm) diameter clear aperture. A worm gear driven by a cryogenic stepper motor allows for precise positioning and prevents undesired rotation when the motors are depowered. A custom-built optical encoder system monitors the bearing angle to an absolute accuracy of +/- 0.1 degrees. The system performed well in Spider during its successful 16 day flight.

We describe the SPIDER flight cryostat, which is designed to cool six millimeter-wavelength telescopes during an Antarctic long-duration balloon flight. The cryostat, one of the largest to have flown on a stratospheric payload, uses liquid helium-4 to deliver cooling power to stages at 4.2 and 1.6 K. Stainless steel capillaries facilitate a high flow impedance connection between the main liquid helium tank and a smaller superfluid tank, allowing the latter to operate at 1.6 K as long as there is liquid in the 4.2 K main tank. Each telescope houses a closed cycle helium-3 adsorption refrigerator that further cools the focal planes down to 300 mK. Liquid helium vapor from the main tank is routed through heat exchangers that cool radiation shields, providing negative thermal feedback. The system performed successfully during a 17 day flight in the 2014-2015 Antarctic summer. The cryostat had a total hold time of 16.8 days, with 15.9 days occurring during flight.

An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

We present the results of integration and characterization of the SPIDER instrument after the 2013 pre-flight campaign. SPIDER is a balloon-borne polarimeter designed to probe the primordial gravitational wave signal in the degree-scale $B$-mode polarization of the cosmic microwave background. With six independent telescopes housing over 2000 detectors in the 94 GHz and 150 GHz frequency bands, SPIDER will map 7.5% of the sky with a depth of 11 to 14 $\mu$K$\cdot$arcmin at each frequency, which is a factor of $\sim$5 improvement over Planck. We discuss the integration of the pointing, cryogenic, electronics, and power sub-systems, as well as pre-flight characterization of the detectors and optical systems. SPIDER is well prepared for a December 2014 flight from Antarctica, and is expected to be limited by astrophysical foreground emission, and not instrumental sensitivity, over the survey region.

We present the technology and control methods developed for the pointing system of the SPIDER experiment. SPIDER is a balloon-borne polarimeter designed to detect the imprint of primordial gravitational waves in the polarization of the Cosmic Microwave Background radiation. We describe the two main components of the telescope's azimuth drive: the reaction wheel and the motorized pivot. A 13 kHz PI control loop runs on a digital signal processor, with feedback from fibre optic rate gyroscopes. This system can control azimuthal speed with < 0.02 deg/s RMS error. To control elevation, SPIDER uses stepper-motor-driven linear actuators to rotate the cryostat, which houses the optical instruments, relative to the outer frame. With the velocity in each axis controlled in this way, higher-level control loops on the onboard flight computers can implement the pointing and scanning observation modes required for the experiment. We have accomplished the non-trivial task of scanning a 5000 lb payload sinusoidally in azimuth at a peak acceleration of 0.8 deg/s$^2$, and a peak speed of 6 deg/s. We can do so while reliably achieving sub-arcminute pointing control accuracy.

We introduce the light-weight carbon fiber and aluminum gondola designed for the SPIDER balloon-borne telescope. SPIDER is designed to measure the polarization of the Cosmic Microwave Background radiation with unprecedented sensitivity and control of systematics in search of the imprint of inflation: a period of exponential expansion in the early Universe. The requirements of this balloon-borne instrument put tight constrains on the mass budget of the payload. The SPIDER gondola is designed to house the experiment and guarantee its operational and structural integrity during its balloon-borne flight, while using less than 10% of the total mass of the payload. We present a construction method for the gondola based on carbon fiber reinforced polymer tubes with aluminum inserts and aluminum multi-tube joints. We describe the validation of the model through Finite Element Analysis and mechanical tests.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

In January 2012, the 10m South Pole Telescope (SPT) was equipped with a polarization-sensitive camera, SPTpol, in order to measure the polarization anisotropy of the cosmic microwave background (CMB). Measurements of the polarization of the CMB at small angular scales (~several arcminutes) can detect the gravitational lensing of the CMB by large scale structure and constrain the sum of the neutrino masses. At large angular scales (~few degrees) CMB measurements can constrain the energy scale of Inflation. SPTpol is a two-color mm-wave camera that consists of 180 polarimeters at 90 GHz and 588 polarimeters at 150 GHz, with each polarimeter consisting of a dual transition edge sensor (TES) bolometers. The full complement of 150 GHz detectors consists of 7 arrays of 84 ortho-mode transducers (OMTs) that are stripline coupled to two TES detectors per OMT, developed by the TRUCE collaboration and fabricated at NIST. Each 90 GHz pixel consists of two antenna-coupled absorbers coupled to two TES detectors, developed with Argonne National Labs. The 1536 total detectors are read out with digital frequency-domain multiplexing (DfMUX). The SPTpol deployment represents the first on-sky tests of both of these detector technologies, and is one of the first deployed instruments using DfMUX readout technology. We present the details of the design, commissioning, deployment, on-sky optical characterization and detector performance of the complete SPTpol focal plane.

The SPTpol camera is a dichroic polarimetric receiver at 90 and 150 GHz. Deployed in January 2012 on the South Pole Telescope (SPT), SPTpol is looking for faint polarization signals in the Cosmic Microwave Background (CMB). The camera consists of 180 individual Transition Edge Sensor (TES) polarimeters at 90 GHz and seven 84-polarimeter camera modules (a total of 588 polarimeters) at 150 GHz. We present the design, dark characterization, and in-lab optical properties of the 150 GHz camera modules. The modules consist of photolithographed arrays of TES polarimeters coupled to silicon platelet arrays of corrugated feedhorns, both of which are fabricated at NIST-Boulder. In addition to mounting hardware and RF shielding, each module also contains a set of passive readout electronics for digital frequency-domain multiplexing. A single module, therefore, is fully functional as a miniature focal plane and can be tested independently. Across the modules tested before deployment, the detectors average a critical temperature of 478 mK, normal resistance R_N of 1.2 Ohm, unloaded saturation power of 22.5 pW, (detector-only) optical efficiency of ~ 90%, and have electrothermal time constants < 1 ms in transition.

The SPTpol camera is a two-color, polarization-sensitive bolometer receiver, and was installed on the 10 meter South Pole Telescope in January 2012. SPTpol is designed to study the faint polarization signals in the Cosmic Microwave Background, with two primary scientific goals. One is to constrain the tensor-to-scalar ratio of perturbations in the primordial plasma, and thus constrain the space of permissible inflationary models. The other is to measure the weak lensing effect of large-scale structure on CMB polarization, which can be used to constrain the sum of neutrino masses as well as other growth-related parameters. The SPTpol focal plane consists of seven 84-element monolithic arrays of 150 GHz pixels (588 total) and 180 individual 90 GHz single-pixel modules. In this paper we present the design and characterization of the 90 GHz modules.

We present the software system used to control and operate the South Pole Telescope. The South Pole Telescope is a 10-meter millimeter-wavelength telescope designed to measure anisotropies in the cosmic microwave background (CMB) at arcminute angular resolution. In the austral summer of 2011/12, the SPT was equipped with a new polarization-sensitive camera, which consists of 1536 transition-edge sensor bolometers. The bolometers are read out using 36 independent digital frequency multiplexing (\dfmux) readout boards, each with its own embedded processors. These autonomous boards control and read out data from the focal plane with on-board software and firmware. An overall control software system running on a separate control computer controls the \dfmux boards, the cryostat and all other aspects of telescope operation. This control software collects and monitors data in real-time, and stores the data to disk for transfer to the United States for analysis.

We consider the estimation of the regularization parameter for the simultaneous deblurring of multiple noisy images via Tikhonov regularization. We approach the problem in three ways. We first reduce the problem to a single-image deblurring for which the regularization parameter can be estimated through a classic generalized cross-validation (GCV) method. A modification of this function is used for correcting the undersmoothing typical of the original technique. With a second method, we minimize an average least-squares fit to the images and define a new GCV function. In the last approach, we use the classical $GCV$ on a single higher-dimensional image obtained by concatanating all the images into a single vector. With a reliable estimator of the regularization parameter, one can fully exploit the excellent computational characteristics typical of direct deblurring methods, which, especially for large images, makes them competitive with the more flexible but much slower iterative algorithms. The performance of the techniques is analyzed through numerical experiments. We find that under the independent homoscedastic and Gaussian assumptions made on the noise, the three approaches provide almost identical results with the first single image providing the practical advantage that no new software is required and the same image can be used with other deblurring algorithms.

We consider the simultaneous deblurring of a set of noisy images whose point spread functions are different but known and spatially invariant, and the noise is Gaussian. Currently available iterative algorithms that are typically used for this type of problem are computationally expensive, which makes their application for very large images impractical. We present a simple extension of a classical least-squares (LS) method where the multi-image deblurring is efficiently reduced to a computationally efficient single-image deblurring. In particular, we show that it is possible to remarkably improve the ill-conditioning of the LS problem by means of stable operations on the corresponding normal equations, which in turn speed up the convergence rate of the iterative algorithms. The performance and limitations of the method are analyzed through numerical simulations. Its connection with a column weighted least-squares approach is also considered in an appendix.

In this second paper in a series dedicated to developing efficient numerical techniques for the deblurring Cosmic Microwave Background (CMB) maps, we consider the case of asymmetric point spread functions (PSF). Although conceptually this problem is not different from the symmetric case, there are important differences from the computational point of view because it is no longer possible to use some of the efficient numerical techniques that work with symmetric PSFs. We present procedures that permit the use of efficient techniques even when this condition is not met. In particular, two methods are considered: a procedure based on a Kronecker approximation technique that can be implemented with the numerical methods used with symmetric PSFs but that has the limitation of requiring only mildly asymmetric PSFs. The second is a variant of the classic Tikhonov technique that works even with very asymmetric PSFs but that requires discarding the edges of the maps. We provide details for efficient implementations of the algorithms. Their performance is tested on simulated CMB maps.

Digital deblurring of images is an important problem that arises in multifrequency observations of the Cosmic Microwave Background (CMB) where, because of the width of the point spread functions (PSF), maps at different frequencies suffer a different loss of spatial resolution. Deblurring is useful for various reasons: first, it helps to restore high frequency components lost through the smoothing effect of the instrument's PSF; second, emissions at various frequencies observed with different resolutions can be better studied on a comparable resolution; third, some map-based component separation algorithms require maps with similar level of degradation. Because of computational efficiency, deblurring is usually done in the frequency domain. But this approach has some limitations as it requires spatial invariance of the PSF, stationarity of the noise, and is not flexible in the selection of more appropriate boundary conditions. Deblurring in real space is more flexible but usually not used because of its high computational cost. In this paper (the first in a series on the subject) we present new algorithms that allow the use of real space deblurring techniques even for very large images. In particular, we consider the use of Tikhonov deblurring of noisy maps with applications to {\it PLANCK}. We provide details for efficient implementations of the algorithms. Their performance is tested on Gaussian and non-Gaussian simulated CMB maps, and PSFs with both circular and elliptical symmetry. Matlab code is made available.