  Local helioseismology provides a unique opportunity to investigate the
subsurface structure and dynamics of active regions and their effect on the
large-scale flows and global circulation of the Sun. We use measurements of
plasma flows in the upper convection zone, provided by the Time-Distance
Helioseismology Pipeline developed for analysis of solar oscillation data
obtained by Helioseismic and Magnetic Imager (HMI) on Solar Dynamics
Observatory (SDO), to investigate the subsurface dynamics of emerging active
region NOAA 11726. The active region emergence was detected in deep layers of
the convection zone about 12 hours before the first bipolar magnetic structure
appeared on the surface, and 2 days before the emergence of most of the
magnetic flux. The speed of emergence determined by tracking the flow
divergence with depth is about 1.4 km/s, very close to the emergence speed in
the deep layers. As the emerging magnetic flux becomes concentrated in sunspots
local converging flows are observed beneath the forming sunspots. These flows
are most prominent in the depth range 1-3 Mm, and remain converging after the
formation process is completed. On the larger scale converging flows around
active region appear as a diversion of the zonal shearing flows towards the
active region, accompanied by formation of a large-scale vortex structure. This
process occurs when a substantial amount of the magnetic flux emerged on the
surface, and the converging flow pattern remains stable during the following
evolution of the active region. The Carrington synoptic flow maps show that the
large-scale subsurface inflows are typical for active regions. In the deeper
layers (10-13 Mm) the flows become diverging, and surprisingly strong beneath
some active regions. In addition, the synoptic maps reveal a complex evolving
pattern of large-scale flows on the scale much larger than supergranulation.
