We describe the self-assembly of multilayer hexagonal DNA arrays through highly regular interlayer packing. Slow cooling of a mixture of three single-stranded DNA sequences with various Mg2+ concentrations leads to the self-assembly of diverse multilayer architectures. The self-assembled aggregates were deposited onto mica surfaces and examined with atomic force microscopy. The size of the two-dimensional arrays and subsequent stacking to form multilayer structures are highly dependent on Mg2+ concentration. DNA bilayers and multilayers of defined shape are favored in 2-5 mM Mg2+ with an average lateral size of 700 nm. Arrays are much larger (up to 20 microm across) in 10-15 mM Mg2+, although multiple layers still make up 20-60% of the observed structures. Domains within single layer architectures were identified using Moiré pattern analysis. Distinct structural phases within the multilayer assemblies include two layers translated by 17.5 nm and interlayer rotations of 20 degrees and 30 degrees. Three layer assemblies have cubic close packing and taller multilayer architectures of 2D DNA sheets were also identified.