In this paper, we consider a four-point coupled boundary value problem for systems of the nonlinear semipositone fractional differential equation
\begin{gather*}\left\{ \begin{array}{ll}
\mathbf{D}_{0+}^\alpha u+\lambda f(t,u,v)=0,\quad 0<t<1, \lambda >0,\\
\mathbf{D}_{0+}^\alpha v+\lambda g(t,u,v)=0,\\
u^{(i)}(0)=v^{(i)}(0)=0, 0\leq i\leq n-2,\\
u(1)=av(\xi), v(1)=bu(\eta), \xi,\eta\in(0,1)
\end{array}\right.\end{gather*}
where λ is a parameter, a,b,ξ,η satisfy ξ,η∈(0,1), 0<abξη<1, α∈(n−1,n] is a real number and n≥3, and D0+α is the Riemann-Liouville's fractional derivative, and f,g are continuous and semipositone. We derive an interval on λ such that for any λ lying in this interval, the semipositone boundary value problem has multiple positive solutions