Given a domain above a Lipschitz graph, we establish solvability results for strongly elliptic second-order systems in divergence-form, allowed to have lower-order (drift) terms, with
L
p
L^p
-boundary data for
p
p
near
2
2
(more precisely, in an interval of the form
(
2
−
ε
,
2
(
n
−
1
)
n
−
2
+
ε
)
\big (2-\varepsilon ,\frac {2(n-1)}{n-2}+\varepsilon \big )
for some small
ε
>
0
\varepsilon >0
). The main novel aspect of our result is that the coefficients of the operator do not have to be constant, or have very high regularity; instead they will satisfy a natural Carleson condition that has appeared first in the scalar case. A significant example of a system to which our result may be applied is the Lamé system for isotropic inhomogeneous materials. We show that our result applies to isotropic materials with Poisson ratio
ν
>
0.396
\nu >0.396
.
Dealing with genuine systems gives rise to substantial new challenges, absent in the scalar case. Among other things, there is no maximum principle for general elliptic systems, and the De Giorgi–Nash–Moser theory may also not apply. We are, nonetheless, successful in establishing estimates for the square-function and the nontangential maximal operator for the solutions of the elliptic system described earlier, and use these as alternative tools for proving
L
p
L^p
solvability results for
p
p
near
2
2
.