The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical MHD, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust \it{quantitative} measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to HI observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent MHD simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.

North-south (N-S) asymmetry of solar activity is a known statistical phenomenon but its significance is difficult to prove or theoretically explain. Here we consider each solar hemisphere as a separate dynamical system connected with the other hemisphere via an unknown coupling parameter. We use a non-linear dynamics approach to calculate the scale-dependent conditional dispersion (CD) of sunspots between hemispheres. Using daily Greenwich sunspot areas, we calculated the Neumann and Pearson chi-squared distances between CDs as indices showing the direction of coupling. We introduce an additional index of synchronization which shows the strength of coupling and allows us to discriminate between complete synchronization and independency of hemispheres. All indices are evaluated in a four-year moving window showing the evolution of coupling between hemispheres. We find that the driver-response interrelation changes between hemispheres have a few pulses during 130 years of Greenwich data with an at least 40 years-long period of unidirectional coupling. These sharp nearly simultaneous pulses of all causality indices are found at the decay of some 11-year cycles. The pulse rate of this new phenomenon of dynamic coupling is irregular: although the first two pulses repeat after 22-year Hale cycles, the last two pulses repeat after three and four 11-year cycles respectively. The last pulse occurs at the decay phase of Cycle 23 so the next pulse will likely appear during the decay of future Cycle 25 or later. This new phenomenon of dynamic coupling reveals additional constraints for understanding and modeling of the long-term solar activity cycles.