  Planar and linear arrays of SQUIDs (superconducting quantum interference
devices), operate as nonlinear magnetic metamaterials in microwaves. Such {\em
SQUID metamaterials} are paradigmatic systems that serve as a test-bed for
simulating several nonlinear dynamics phenomena. SQUIDs are highly nonlinear
oscillators which are coupled together through magnetic dipole-dipole forces
due to their mutual inductance; that coupling falls-off approximately as the
inverse cube of their distance, i.~e., it is non-local. However, it can be
approximated by a local (nearest-neighbor) coupling which in many cases
suffices for capturing the essentials of the dynamics of SQUID metamaterials.
For either type of coupling, it is numerically demonstrated that chimera states
as well as other spatially non-uniform states can be generated in SQUID
metamaterials under time-dependent applied magnetic flux for appropriately
chosen initial conditions. The mechanism for the emergence of these states is
discussed in terms of the multistability property of the individual SQUIDs
around their resonance frequency and the attractor crowding effect in systems
of coupled nonlinear oscillators. Interestingly, generation and control of
chimera states in SQUID metamaterials can be achieved in the presence of a
constant (dc) flux gradient with the SQUID metamaterial initially at rest.
