Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2\sigma$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at $1.42\,$GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen (Fialkov & Barkana 2019). We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on $141\,$hours of observations with LOFAR (Mertens et al. 2020) to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 68% C.L.) that the cosmological radio background can be at most 0.7% of the CMB at $1.42\,$GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-0.7$% over the CMB (at $1.42\,$GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data disfavour scenarios with inefficient X-ray sources and cases where the Universe was ionized by massive halos only.

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

We analyse the accuracy of radio interferometric gridding of visibilities with the aim to quantify the Epoch of Reionization (EoR) 21-cm power spectrum bias caused by gridding, ultimately to determine the suitability of different imaging algorithms and gridding settings for 21-cm power spectrum analysis. We simulate realistic LOFAR data, and construct power spectra with convolutional gridding and w-stacking, w-projection, image domain gridding and without w-correction. These are compared against directly Fourier transformed data. The influence of oversampling, kernel size, w-quantization, kernel windowing function and image padding are quantified. The gridding excess power is measured with a foreground subtraction strategy, for which foregrounds have been subtracted using Gaussian progress regression, as well as with a foreground avoidance strategy. Constructing a power spectrum that has a bias significantly lower compared to the expected EoR signals is possible with the tested methods, but requires a kernel oversampling factor > 4000 and, when using w-correction, > 500 w-quantization levels. These values are higher than typical values used for imaging, but are computationally feasible. The kernel size and padding factor parameters are less crucial. Among the tested methods, image domain gridding shows the highest accuracy with the lowest imaging time. LOFAR 21-cm power spectrum results are not affected by gridding. Image domain gridding is overall the most suitable algorithm for 21-cm EoR experiments, including for future SKA EoR analyses. Nevertheless, convolutional gridding with tuned parameters results in sufficient accuracy. This holds also for w-stacking for wide-field imaging. The w-projection algorithm is less suitable because of the kernel oversampling requirements, and a faceting approach is unsuitable due to the resulting spatial discontinuities.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2\sigma$ upper limit on the 21-cm power spectrum of $\Delta_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $\Delta_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes.

In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE; > 100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this gamma-ray emitting jet. The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and cover a time period of ~2 years. Through a segmented wavelet decomposition method we estimate the apparent displacement of individual plasma features; we then perform a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the SED. Proper motion is mostly detected along the edges of the flow, which appears strongly limb-brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ~$11.5$ c. This constrains the jet viewing angle to assume relatively small values ($\theta\lesssim10^{\circ}$). In the acceleration region, extending up to a de-projected distance of ~$4.8\times10^4$ Schwarzschild radii (~$11$ pc), the jet is collimating ($r\propto z^{0.40\pm 0.04}$), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated ($U_B/U_e>1$), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling.

We investigate the implications of interference detection for experiments that are pursuing a detection of the redshifted 21-cm signals from the Epoch of Reionization. Interference detection causes samples to be sporadically flagged and rejected. As a necessity to reduce the data volume, flagged samples are typically (implicitly) interpolated during time or frequency averaging or uv-gridding. This so-far unexplored systematic biases the 21-cm power spectrum, and it is important to understand this bias for current 21-cm experiments as well as the upcoming SKA Epoch of Reionization experiment. We analyse simulated data using power spectrum analysis and Gaussian process regression. We find that the combination of flagging and averaging causes tiny spectral fluctuations, resulting in `flagging excess power'. This excess power does not substantially average down over time and, without extra mitigation techniques, can exceed the power of realistic models of the 21-cm reionization signals in LOFAR observations. We mitigate the bias by i) implementing a novel way to average data using a Gaussian-weighted interpolation scheme; ii) using unitary instead of inverse-variance weighting of visibilities; and iii) using low-resolution forward modelling of the data. After these modifications, which have been integrated in the LOFAR EoR processing pipeline, the excess power reduces by approximately three orders of magnitude, and is no longer preventing a detection of the 21-cm signals.

(abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $\lambda = 3\,$mm) reach a resolution of about 50 $\mu$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the Global Millimeter VLBI Array. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy/beam, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the sources. The Gaussian model-fit-based estimates of brightness temperature ($T_\mathrm{b}$) at the jet base (core) and in moving regions (jet components) downstream from the core were compared to the estimates of $T_\mathrm{b}$ limits made directly from the visibility data, demonstrating a good agreement between the two methods. The apparent brightness temperature estimates for the jet cores in our sample range from $2.5 \times 10^{9}\,$K to $ 1.3\times 10^{12}\,$K. A population model with a single intrinsic value of brightness temperature, $T_\mathrm{0}$, is applied to reproduce the observed $T_\mathrm{b}$ distribution. It yields $T_\mathrm{0} = (3.77^{+0.10}_{-0.14}) \times 10^{11}\,$K for the jet cores, implying that the inverse Compton losses dominate the emission. In the jet components, $T_\mathrm{0} =(1.42^{+0.16}_{-0.19})\times 10^{11}\,$K is found, slightly higher than the equipartition limit of $\sim5\times 10^{10}\,$K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses dominate the observed changes of $T_\mathrm{b}$ along the jet.

Detecting and characterizing the Epoch of Reionization and Cosmic Dawn via the redshifted 21-cm hyperfine line of neutral hydrogen will revolutionize the study of the formation of the first stars, galaxies, black holes and intergalactic gas in the infant Universe. The wealth of information encoded in this signal is, however, buried under foregrounds that are many orders of magnitude brighter. These must be removed accurately and precisely in order to reveal the feeble 21-cm signal. This requires not only the modeling of the Galactic and extra-galactic emission, but also of the often stochastic residuals due to imperfect calibration of the data caused by ionospheric and instrumental distortions. To stochastically model these effects, we introduce a new method based on `Gaussian Process Regression' (GPR) which is able to statistically separate the 21-cm signal from most of the foregrounds and other contaminants. Using simulated LOFAR-EoR data that include strong instrumental mode-mixing, we show that this method is capable of recovering the 21-cm signal power spectrum across the entire range $k = 0.07 - 0.3 \ \rm{h\, cMpc^{-1}}$. The GPR method is most optimal, having minimal and controllable impact on the 21-cm signal, when the foregrounds are correlated on frequency scales $\gtrsim 3$\,MHz and the rms of the signal has $\sigma_{\mathrm{21cm}} \gtrsim 0.1\,\sigma_{\mathrm{noise}}$. This signal separation improves the 21-cm power-spectrum sensitivity by a factor $\gtrsim 3$ compared to foreground avoidance strategies and enables the sensitivity of current and future 21-cm instruments such as the {\sl Square Kilometre Array} to be fully exploited.

Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum.

Direct detection of the Epoch of Reionization (EoR) via the red-shifted 21-cm line will have unprecedented implications on the study of structure formation in the infant Universe. To fulfill this promise, current and future 21-cm experiments need to detect this weak EoR signal in the presence of foregrounds that are several orders of magnitude larger. This requires extreme noise control and improved wide-field high dynamic-range imaging techniques. We propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere, or equivalently in the $uvw$-domain. The method uses the one-to-one relation between spherical waves and spherical harmonics (SpH). It consistently handles signals from the entire sky, and does not require a $w$-term correction. The spherical-harmonics coefficients represent the sky-brightness distribution and the visibilities in the $uvw$-domain, and provide a direct estimate of the spatial power spectrum. Using these spectrally-smooth SpH coefficients, bright foregrounds can be removed from the signal, including their side-lobe noise, which is one of the limiting factors in high dynamics range wide-field imaging. Chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum, compared to a power spectrum computed directly from the images of the foreground visibilities where the wedge is clearly present. We illustrate our method using simulated LOFAR observations, finding an excellent reconstruction of the input EoR signal with minimal bias.

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

Very long baseline interferometry (VLBI) imaging of radio emission from extragalactic jets provides a unique probe of physical mechanisms governing the launching, acceleration, and collimation of relativistic outflows. The two-dimensional structure and kinematics of the jet in M\,87 (NGC\,4486) have been studied by applying the Wavelet-based Image Segmentation and Evaluation (WISE) method to 11 images obtained from multi-epoch Very Long Baseline Array (VLBA) observations made in January-August 2007 at 43 GHz ($\lambda = 7$ mm). The WISE analysis recovers a detailed two-dimensional velocity field in the jet in M\,87 at sub-parsec scales. The observed evolution of the flow velocity with distance from the jet base can be explained in the framework of MHD jet acceleration and Poynting flux conversion. A linear acceleration regime is observed up to $z_{obs} \sim 2$\,mas. The acceleration is reduced at larger scales, which is consistent with saturation of Poynting flux conversion. Stacked cross correlation analysis of the images reveals a pronounced stratification of the flow. The flow consists of a slow, mildly relativistic layer (moving at $\beta \sim 0.5\,c$), associated either with instability pattern speed or an outer wind, and a fast, accelerating stream line (with $\beta \sim 0.92$, corresponding to a bulk Lorentz factor $\gamma \sim 2.5$). A systematic difference of the apparent speeds in the northern and southern limbs of the jet is detected, providing evidence for jet rotation. The angular velocity of the magnetic field line associated with this rotation suggests that the jet in M87 is launched in the inner part of the disk, at a distance $r_0 \sim 5\, R_\mathrm{s}$ from the central engine. The combined results of the analysis imply that MHD acceleration and conversion of Poynting flux to kinetic energy play the dominant roles in collimation and acceleration of the flow in M\,87.

Velocity measurements made from multiple-epoch astronomical images of evolving objects with optically thin continuum emission (e.g. as relativistic jets or expanding supernova shells) may be confused as a result of the overlap of semi-transparent features moving at different speeds. Multi-scale wavelet decomposition can be effectively applied to identify and track such overlapping features, provided that their respective structural responses can be separated over the spatial scales used for the decomposition. We developed a new method that combines the stacked cross-correlation with the wavelet-based image segmentation and evaluation (WISE) technique of decomposition of two-dimensional structures, to separate and track dominant spatial responses of overlapping evolving features. The method is tested on a set of simulated images of a stratified relativistic jet, demonstrating the robust detection of both the faster spine and the slower sheath speeds. The method is applied to mutliple-epoch images from the MOJAVE survey, revealing two different superluminal streams inside the jet in 3C 273 and the acceleration of the flow in 3C 120. The method can be applied to densely monitored objects with composite structural evolution such as the parsec-scale jet in M87 or heavily resolved expanding supernova shells.

High-resolution Very-Long-Baseline Interferometry observations of relativistic jets are essential to constrain fundamental parameters of jet formation models. At a distance of 249 Mpc, Cygnus A is a unique target for such studies, being the only Fanaroff-Riley Class II radio galaxy for which a detailed sub-parsec scale imaging of the base of both jet and counter-jet can be obtained. Observing at millimeter wavelengths unveils those regions which appear self-absorbed at longer wavelengths and enables an extremely sharp view towards the nucleus to be obtained. We performed 7 mm Global VLBI observations, achieving ultra-high resolution imaging on scales down to 90 $\mu$as. This resolution corresponds to a linear scale of only ${\sim}$400 Schwarzschild radii (for $M_{\mathrm{BH}}=2.5 \times 10^9 M_{\odot}$). We studied the kinematic properties of the main emission features of the two-sided flow and probed its transverse structure through a pixel-based analysis. We suggest that a fast and a slow layer, with different acceleration gradients, exist in the flow. The extension of the acceleration region is large (${\sim} 10^4 R_{\mathrm{S}}$), indicating that the jet is magnetically-driven. The limb brightening of both jet and counter-jet and their large opening angles ($\phi_\mathrm{J}{\sim} 10^{\circ}$) strongly favor a spine-sheath structure. In the acceleration zone, the flow has a parabolic shape ($r\propto z^{0.55\pm 0.07}$). The acceleration gradients and the collimation profile are consistent with the expectations for a jet in "equilibrium'' (Lyubarsky 2009), achieved in the presence of a mild gradient of the external pressure ($p\propto z^{-k}, k\leq2$).}

Radio-loud AGN typically manifest powerful relativistic jets extending up to millions of light years and often showing superluminal motions organised in a complex kinematic pattern. A number of physical models are still competing to explain the jet structure and kinematics revealed by radio images using the VLBI technique. Robust measurements of longitudinal and transverse velocity field in the jets would provide crucial information for these models. This is a difficult task, particularly for transversely resolved jets in objects like 3C 273 and M87. To address this task, we have developed a new technique for identifying significant structural patterns (SSP) of smooth, transversely resolved flows and obtaining a velocity field from cross-correlation of these regions in multi-epoch observations. Detection of individual SSP is performed using the wavelet decomposition and multiscale segmentation of the observed structure. The cross-correlation algorithm combines structural information on different scales of the wavelet decomposition, providing a robust and reliable identification of related SSP in multi -epoch images. The algorithm enables recovering structural evolution on scales down to 0.25 full width at half maximum (FWHM) of the image point spread function (PSF). We present here examples of applying this algorithm to obtain the first detailed transverse velocity fields and to study the kinematic evolution in the parsec-scale jets in 3C 273 and M87.

Context. Images of spatially resolved astrophysical objects contain a wealth of morphological and dynamical information, and effective extraction of this information is of paramount importance for understanding the physics and evolution of these objects. Algorithms and methods employed presently for this purpose (such as, for instance, Gaussian model fitting) often use simplified approaches for describing the structure of resolved objects. Aims. Automated (unsupervised) methods for structure decomposition and tracking of structural patterns are needed for this purpose, in order to be able to deal with the complexity of structure and large amount of data involved. Methods. A new Wavelet-based Image Segmentation and Evaluation (WISE) method is developed for multiscale decomposition, segmentation, and tracking of structural patterns in astronomical images. Results. The method is tested against simulated images of relativistic jets and applied to data from long-term monitoring of parsec- scale radio jets in 3C 273 and 3C 120. Working at its coarsest resolution, WISE reproduces exceptionally well the previous results of model fitting evaluation of the structure and kinematics in these jets. Extending the WISE structure analysis to the finer scales provides the first robust measurements of two-dimensional velocity fields in these jets and indicates that the velocity fields are likely to reflect the evolution of Kelvin-Hemlholtz instability developing the flow.