In the outer heliosphere, beyond the solar wind termination shock, it is expected that the warped heliospheric current sheet forms a region of closely-packed, multiple, thin current sheets. Such a system may be subject to the ion-kinetic tearing instability, and hence generate magnetic islands and hot populations of ions associated with magnetic reconnection. Reconnection processes in this environment have important implications for local particle transport, and for particle acceleration at reconnection sites and in turbulence. We study this complex environment by means of three-dimensional hybrid simulations over long time scales, in order to capture the evolution from linear growth of the tearing instability to a fully developed turbulent state at late times. The final state develops from the highly ordered initial state via both forward and inverse cascades. Component and spectral anisotropy in the magnetic fluctuations is present when a guide field is included. The inclusion of a population of new-born interstellar pickup protons does not strongly affect these results. Finally, we conclude that reconnection between multiple current sheets can act as an important source of turbulence in the outer heliosphere, with implications for energetic particle acceleration and propagation.

Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfv\'en speed of the ramp, and motion in the sense of gyration of the reflected ions.