Source author record

Hongchi Wang

Hongchi Wang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

CO Emission Delineating the Interface between the Milky Way Nuclear Wind Cavity and the Gaseous Disk

Based on the MWISP survey, we study high-z CO emission toward the tangent points, in which the distances of the molecular clouds (MCs) are well determined. In the region of l=12-26 deg and |b|<5.1 deg, a total of 321 MCs with |z|> 110 pc are identified, of which nearly 30 extreme high-z MCs (EHMCs at |z|> 260 pc) are concentrated in a narrow region of R_GC=2.6-3.1 kpc. The EHMC concentrations, together with other high-z MCs at R_GC=2.3-2.6 kpc, constitute molecular crater-wall structures surrounding the edges of the HI voids that are physically associated with the Fermi bubbles. Intriguingly, some large high-z MCs, which lie in the crater walls above and below the Galactic plane, show cometary structures with the head toward the plane, favouring the scenario that the entrained molecular gas moves with the multi-phase flows from the plane to the high-z regions. We suggest that the Milky Way nuclear wind has a significant impact on the Galactic gaseous disk. The powerful nuclear wind at ~3-6 Myr ago is likely responsible for the observational features, (1) the enhanced CO gas lying in the edges of the HI voids, (2) the deficiency of atomic and molecular gas within R_GC<3 kpc, (3) the possible connection between the EHMC concentrations and the 3-kpc arm, and (4) the elongated high-z MCs with the tail pointing away from the Galactic plane.

preprint2022arXiv

Dependence of Molecular Cloud Samples on Angular Resolution, Sensitivity, and Algorithms

In this work, we investigate the observational and algorithmic effects on molecular cloud samples identified from position-position-velocity (PPV) space. By smoothing and cutting off the high quality data of the Milky Way Imaging Scroll Painting (MWISP) survey, we extract various molecular cloud samples from those altered data with the DBSCAN (density-based spatial clustering of applications with noise) algorithm. Those molecular cloud samples are subsequently used to gauge the significance of sensitivity, angular/velocity resolution, and DBSCAN parameters. Two additional surveys, the FCRAO Outer Galaxy Survey (OGS) and the CfA-Chile 1.2 m complete CO (CfA-Chile) survey, are used to verify the MWISP results. We found that molecular cloud catalogs are not unique and the boundary and therefore the number shows strong variation with angular resolution and sensitivity. At low angular resolution (large beam sizes), molecular clouds merge together in PPV space, while low sensitivity (high cutoffs) misses small faint molecular clouds and takes bright parts of large molecular clouds as single ones. At high angular resolution and sensitivity, giant molecular clouds (GMCs) are resolved into individual clouds, and their diffuse components are also revealed. Consequently, GMCs are more appropriately interpreted as clusters or aggregates of molecular clouds, i.e., GMCs represent molecular cloud samples themselves.

preprint2022arXiv

Gas Column Density Distribution of Molecular Clouds in the Third Quadrant of the Milky Way

We have obtained column density maps for an unbiased sample of 120 molecular clouds in the third quadrant of the Milky Way mid-plane (b$\le |5|^{\circ}$) within the galactic longitude range from 195$^{\circ}$ to 225$^{\circ}$, using the high sensitivity $^{12}$CO and $^{13}$CO ($J=1-0$) data from the Milky Way Imaging Scroll Painting (MWISP) project. The probability density functions of the molecular hydrogen column density of the clouds, N-PDFs, are fitted with both log-normal (LN) function and log-normal plus power-law (LN+PL) function. The molecular clouds are classified into three categories according to their shapes of N-PDFs, i.e., LN, LN+PL, and UN (unclear), respectively. About 72\% of the molecular clouds fall into the LN category, while 18\% and 10\% into the LN+PL and UN categories, respectively. A power-law scaling relation, $σ_s\propto N_{H_2}^{0.44}$, exists between the width of the N-PDF, $σ_s$, and the average column density, $N_{H_2}$, of the molecular clouds. However, $σ_s$ shows no correlation with the mass of the clouds. A correlation is found between the dispersion of normalized column density, $σ_{N/\rm <N>}$, and the sonic Mach number, $\mathcal{M}$, of molecular clouds. Overall, as predicted by numerical simulations, the N-PDFs of the molecular clouds with active star formation activity tend to have N-PDFs with power-law high-density tails.

preprint2020arXiv

Distances to molecular clouds in the second Galactic quadrant

We present distances to 76 medium-sized molecular clouds and an extra large-scale one in the second Galactic quadrant ($104.75^\circ <l<150.25^\circ $ and $|b|<5.25^\circ$), 73 of which are accurately measured for the first time. Molecular cloud samples are drawn from $l$-$b$-$V$ space ($-95 < V_{\rm LSR}< 25$ \kms) with the density-based spatial clustering of applications with noise (DBSCAN) algorithm, and distances are measured with the background-eliminated extinction-parallax (BEEP) method using extinctions and Gaia DR2 parallaxes. The range of measured distances to 76 molecular clouds is from 211 to 2631 pc, and the extra large-scale molecular cloud appears to be a coherent structure at about 1 kpc, across about 40° ($\sim$700 pc) in the Galactic longitude.

preprint2020arXiv

Local Molecular Gas toward the Aquila Rift Region

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50" and a grid spacing of 30". The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

preprint2020arXiv

The Molecular Clouds associated with the H II Regions/Candidates between $l=207.7^\circ$ and $l=211.7^\circ$

Using the PMO-13.7 m millimeter telescope at Delingha in China, we have conducted a large-scale simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the sky region centered at $l$=$209.7^\circ$, $b$=$-$2.25$^\circ$ with a coverage of $4.0^\circ \times 4.5^\circ$. The majority of the emission in the region comes from the clouds with velocities lying in the range from $-$3 km s$^{-1}$ to 55 km s$^{-1}$, at kinematic distances from 0.5 kpc to 7.0 kpc. The molecular clouds in the region are concentrated into three velocity ranges. The molecular clouds associated with the ten H II regions/candidates are identified and their physical properties are presented. Massive stars are found within Sh2-280, Sh2-282, Sh2-283, and BFS54, and we suggest them to be the candidate excitation sources of the H II regions. The distributions of excitation temperature and line width with the projected distance from the center of H II region/candidate suggest that the majority of the ten H II regions/candidates and their associated molecular gas are three-dimensional structures, rather than two-dimensional structures.

preprint2019arXiv

The Column Density Structure of Orion A Depicted by N-PDF

We have conducted a large-field simultaneous survey of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ emission toward the Orion A giant molecular cloud (GMC) with a sky coverage of $\sim$ 4.4 deg$^2$ using the PMO-13.7 m millimeter-wavelength telescope. We use the probability distribution function of the column density (N-PDF) to investigate the distribution of molecular hydrogen in the Orion A GMC. The H$_2$ column density, derived from the $^{13}$CO emission, of the GMC is dominated by log-normal distribution in the range from $\sim$4$\times10^{21}$ to $\sim$1.5$\times10^{23}$ cm$^{-2}$ with excesses both at the low-density and high-density ends. The excess of the low-density end is possibly caused by an extended and low-temperature ($\sim$10 K) component with velocities in the range of 5$-$8 km s$^{-1}$. Compared with the northern sub-regions, the southern sub-regions of the Orion A GMC contain less gas with column density in $N_{H_2} > 1.25\times 10^{22}\ \rm{cm}^{-2}$. The dispersions of the N-PDFs of the sub-regions are found to correlate with the evolutionary stages of the clouds across the Orion A GMC. The structure hierarchy of Orion A GMC is explored with the DENDROGRAM algorithm, and it is found that the GMC is composed of two branches. All structures except one in the tree have virial parameters less than 2, indicating self-gravity is important on the spatial scales from $\sim$0.3 to $\sim$4 pc. Although power-laws and departures from log-normal distributions are found at the high-density end of N-PDFs of active star-forming regions, the N-PDFs of structures in the Orion A GMC are predominantly log-normal on scales from R$\sim$0.4 to 4 pc.

preprint2015arXiv

A deep near-infrared survey toward the Aquila molecular cloud - I. Molecular hydrogen outflows

We have performed an unbiased deep near-infrared survey toward the Aquila molecular cloud with a sky coverage of ~1 deg2. We identified 45 molecular hydrogen emission-line objects(MHOs), of which only 11 were previously known. Using the Spitzer archival data we also identified 802 young stellar objects (YSOs) in this region. Based on the morphology and the location of MHOs and YSO candidates, we associate 43 MHOs with 40 YSO candidates. The distribution of jet length shows an exponential decrease in the number of outflows with increasing length and the molecular hydrogen outflows seem to be oriented randomly. Moreover, there is no obvious correlation between jet lengths, jet opening angles, or jet H2 1-0 S(1) luminosities and spectral indices of the possible driving sources in this region. We also suggest that molecular hydrogen outflows in the Aquila molecular cloud are rather weak sources of turbulence, unlikely to generate the observed velocity dispersion in the region of survey.

preprint2014arXiv

Broad-band transmission spectrum and K-band thermal emission of WASP-43b as observed from the ground

(Abridged) We observed one transit and one occultation of the hot Jupiter WASP-43b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2-meter telescope. From the transit event, we have independently derived WASP-43's system parameters with high precision, and improved the period to be 0.81347437(13) days. No significant variation in transit depths is detected, with the largest deviations coming from the i', H, and K bands. Given the observational uncertainties, the broad-band transmission spectrum can be explained by either a flat featureless straight line that indicates thick clouds, synthetic spectra with absorption signatures of atomic Na/K or molecular TiO/VO that indicate cloud-free atmosphere, or a Rayleigh scattering profile that indicates high-altitude hazes. From the occultation event, we have detected planetary dayside thermal emission in the K-band with a flux ratio of 0.197 +/- 0.042%, which confirms previous detections obtained in the 2.09 micron narrow band and Ks-band. The K-band brightness temperature 1878 +108/-116 K favors an atmosphere with poor day- to night-side heat redistribution. We also have a marginal detection in the i'-band (0.037 +0.023/-0.021%), which is either a false positive, a signature of non-blackbody radiation at this wavelength, or an indication of reflective hazes at high altitude.

preprint2014arXiv

Ground-based detection of the near-infrared emission from the dayside of WASP-5b

(Abridged) WASP-5b is a highly irradiated dense hot Jupiter orbiting a G4V star every 1.6 days. We observed two secondary eclipses of WASP-5b in the J, H and K bands simultaneously. Thermal emission of WASP-5b is detected in the J and K bands. The retrieved planet-to-star flux ratios in the J and K bands are 0.168 +0.050/-0.052% and 0.269+/-0.062%, corresponding to brightness temperatures of 2996 +212/-261K and 2890 +246/-269K, respectively. No thermal emission is detected in the H band, with a 3-sigma upper limit of 0.166%, corresponding to a maximum temperature of 2779K. On the whole, our J, H, K results can be explained by a roughly isothermal temperature profile of ~2700K in the deep layers of the planetary dayside atmosphere that are probed at these wavelengths. Together with Spitzer observations, which probe higher layers that are found to be at ~1900K, a temperature inversion is ruled out in the range of pressures probed by the combined data set. While an oxygen-rich model is unable to explain all the data, a carbon-rich model provides a reasonable fit but violates energy balance.

preprint2014arXiv

Herbig-Haro objects and mid-infrared outflows in the Vela C molecular cloud

We have performed a deep [SII]6717/6731 wide field Herbig-Haro (HH) object survey toward the Vela C molecular cloud with a sky coverage of about 2 deg2. In total, 18 new HH objects, HH 1090-1107, are discovered and the two previously known HH objects, HH 73-74, are also detected in our [SII] images. We also present an investigation of mid-infrared outflows in the Vela C molecular cloud using the Wide-field Infrared Survey Explorer images taken from AllWISE data release. Using the method suggested by Zhang & Wang, eleven extended green objects (EGOs) are identified to be the mid-infrared outflows, including 6 new mid-infrared outflows that have not been detected previously at other wavelengths and 5 mid-infrared counterparts of the HH objects detected in this work. Using the AllWISE Source Catalog and the source classification scheme suggested by Koenig et al., we have identified 56 young stellar object (YSO) candidates in the Vela C molecular cloud. The possible driving sources of the HH objects and EGOs are discussed based on the morphology of HH objects and EGOs and the locations of HH objects, EGOs and YSO candidates. Finally we associate 12 HH objects and five EGOs with 10 YSOs and YSO candidates. The median length of the outflows in Vela C is 0.35pc and the outflows seem to be oriented randomly.

preprint2014arXiv

Observed spectral energy distribution of the thermal emission from the dayside of WASP-46b

We aim to construct a spectral energy distribution (SED) for the emission from the dayside atmosphere of the hot Jupiter WASP-46b and to investigate its energy budget. We observed a secondary eclipse of WASP-46b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2m telescope. Eclipse depths of the acquired light curves were derived to infer the brightness temperatures at multibands that cover the SED peak. We report the first detection of the thermal emission from the dayside of WASP-46b in the K band at 4.2-sigma level and tentative detections in the H (2.5-sigma) and J (2.3-sigma) bands, with flux ratios of 0.253 +0.063/-0.060%, 0.194 +/- 0.078%, and 0.129 +/- 0.055%, respectively. The derived brightness temperatures (2306 +177/-187K, 2462 +245/-302K, and 2453 +198/-258K, respectively) are consistent with an isothermal temperature profile of 2386K, which is significantly higher than the dayside-averaged equilibrium temperature, indicative of very poor heat redistribution efficiency. We also investigate the tentative detections in the g'r'i' bands and the 3-sigma upper limit in the z' band, which might indicate the existence of reflective clouds if these tentative detections do not arise from systematics.

preprint2012arXiv

A Comparison of Approaches in Fitting Continuum SEDs

We present a detailed comparison of two approaches, the use of a pre-calculated database and simulated annealing (SA), for fitting the continuum spectral energy distribution (SED) of astrophysical objects whose appearance is dominated by surrounding dust. While pre-calculated databases are commonly used to model SED data, only few studies to date employed SA due to its unclear accuracy and convergence time for this specific problem. From a methodological point of view, different approaches lead to different fitting quality, demand on computational resources and calculation time. We compare the fitting quality and computational costs of these two approaches for the task of SED fitting to provide a guide to the practitioner to find a compromise between desired accuracy and available resources. To reduce uncertainties inherent to real datasets, we introduce a reference model resembling a typical circumstellar system with 10 free parameters. We derive the SED of the reference model with our code MC3D at 78 logarithmically distributed wavelengths in the range [0.3um, 1.3mm] and use this setup to simulate SEDs for the database and SA. Our result shows directly the applicability of SA in the field of SED modeling, since the algorithm regularly finds better solutions to the optimization problem than a pre-calculated database. As both methods have advantages and shortcomings, a hybrid approach is preferable. While the database provides an approximate fit and overall probability distributions for all parameters deduced using Bayesian analysis, SA can be used to improve upon the results returned by the model grid.

preprint2012arXiv

An Algorithm for Preferential Selection of Spectroscopic Targets in LEGUE

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

preprint2012arXiv

LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) The survey science plan

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<δ<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

preprint2012arXiv

The LEGUE High Latitude Bright Survey Design for the LAMOST Pilot Survey

We describe the footprint and input catalog for bright nights in the LAMOST Pilot Survey, which began in October 2011. Targets are selected from two stripes in the north and south Galactic Cap regions, centered at $α$= 29$^\circ$, with 10$^\circ$ width in declination, covering right ascension of 135$^\circ-290^\circ$ and -30$^\circ$ to 30$^\circ$ respectively. We selected spectroscopic targets from a combination of the SDSS and 2MASS point source catalogs. The catalog of stars defining the field centers (as required by the Shack-Hartmann wavefront sensor at the center of the LAMOST field) consists of all V < 8m stars from the Hipparcos catalog. We employ a statistical selection algorithm that assigns priorities to targets based on their positions in multidimensional color/magnitude space. This scheme overemphasizes rare objects and de-emphasizes more populated regions of magnitude and color phase space, while ensuring a smooth, well-understood selection function. A demonstration of plate design is presented based on the Shack-Hartmann star catalog and an input catalog that was generated by our target selection routines.

preprint2012arXiv

The LEGUE Input Catalogue for Dark Night Observing in the LAMOST Pilot Survey

We outline the design of the dark nights portion of the LAMOST Pilot Survey, which began observations in October 2011. In particular, we focus on Milky Way stellar candidates that are targeted for the LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) survey. We discuss the regions of sky in which spectroscopic candidates were selected, and the motivations for selecting each of these sky areas. Some limitations due to the unique design of the telescope are discussed, including the requirement that a bright (V < 8) star be placed at the center of each plate for wavefront sensing and active optics corrections. The target selection categories and scientific goals motivating them are briefly discussed, followed by a detailed overview of how these selection functions were realized. We illustrate the difference between the overall input catalog - Sloan Digital Sky Survey (SDSS) photometry - and the final targets selected for LAMOST observation.

preprint2012arXiv

The selection of LEGUE disk targets for LAMOST's pilot survey

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<α<67^\circ$ and $42^\circ<δ<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

preprint2012arXiv

The site conditions of the Guo Shou Jing Telescope

The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.

preprint2009arXiv

Outflows in rho Ophiuchi as Seen with the Spitzer Infrared Array Camera

Using the IRAC images from the Spitzer c2d program, we have made a survey of mid-infrared outflows in the rho Ophiuchi molecular cloud. Extended objects that have prominent emission in IRAC channel 2 (4.5 micron) compared to IRAC channel 1 (3.6 micron) and stand out as green objects in the three-color images (3.6 micron in blue, 4.5 micron in green, 8.0 micron in red) are identified as mid-infrared outflows. As a result, we detected 13 new outflows in the rho Ophiuchi molecular cloud that have not been previously observed in optical or near-infrared. In addition, at the positions of previously observed HH objects or near-infrared emission, we detected 31 mid-infrared outflows, among which seven correspond to previously observed HH objects and 30 to near-infrared emission. Most of the mid-infrared outflows detected in the rho Ophiuchi cloud are concentrated in the L1688 dense core region. In combination with the survey results for Young Stellar Objects (YSOs) and millimeter and sub-millimeter sources, the distribution of mid-infrared outflows in the rho Ophiuchi molecular complex hints a propagation of star formation in the cloud in the direction from the northwest to the southeast as suggested by previous studies of the region.