We study inhomogeneous states in the t-J model using an unrestricted Gutzwiller approximation. We find that pa × pa checkerboard order, where p is a doping dependent number, emerges from Fermi surface instabilities of both the staggered flux phase and the Fermi liquid state with realistic band parameters. In both cases, the checkerboard order develops at wave vectors (±2π/pa, 0), (0, ±2π/pa) that are tied to the peaks of the wave-vector dependent susceptibility, and is of the Lomer-Rice-Scott type. The properties of such periodic, inhomogeneous states are discussed in connection to the checkerboard patterns observed by STM in underdoped cuprates.Spatially inhomogeneous, periodically modulated local density of states (LDOS) has been observed recently by STM in high-T c cuprates under a variety of conditions where the superconductivity is weakened 1,2,3,4,5 . The hope for a certain type of ordered state competing with superconductivity has been raised in connection to the pseudogap phenomena. The tunneling conductance map exhibits pa × pa checkerboard (CB) patterns corresponding to wave vectors (±2π/p, 0), (0, ±2π/p), where p takes on values between 4 and 5 (hereafter, we set the lattice constant a = 1). In and around a vortex core in Bi 2 Sr 2 CaCu 2 O 8+δ (BSCCO), p ∼ 4.3 1 ; p ∼ 4 in optimally doped BSCCO in zero magnetic field 2 ; p ≃ 4.7 in underdoped BSCCO in the pseudogap phase above T c 3 ; p ≃ 4 in lightly doped oxychlorides Ca 2−x Na x CuO 2 Cl 2 (Na-CCOC) at very low temperatures 4 ; and in substantially underdoped BSCCO, p ∼ 4.5 were observed in the dark regions of the conductance map 5 . The nature and the origin of the LDOS modulations, in particular the short coherence length and the role of dopants 6,7 , are unclear at present. Several inhomogeneous electronic states have been proposed theoretically, including pair density waves 8,9 ; hole Wigner crystal 10 ; a Wigner crystal of hole pairs embedded in a d-wave resonating valence bond (RVB) state 11 ; and valence-bond solid with or without charge order 12 . How certain inhomogeneous electronic states with CB order arise from the microscopic t-J model of doped Mott insulators is the focus of this work. Recent density matrix renormalization group calculations performed on small t-J clusters found approximate CB like patterns with strong one-dimensional stripe characters 13 . In this paper, we discuss a spatially unrestricted Gutzwiller approximation, which allows us to study large systems in the thermodynamic limit. We focus primarily on inhomogeneous solutions in the nonsuperconducting phase at moderate doping where the underlying homogeneous phase exhibits a Fermi surface (FS). Specifically, we consider the two situations displayed in Fig. 1 for the spectral intensity of the low energy single-particle excitations as measured by angleresolved photoemission (ARPES), in the staggered flux phase (SFP) 14 and the uniform short-range RVB state 15 with realistic band parameters. The quasiparticle scattering in the SFP (Case A in Fig.1a) is dominated by the...