  Motivated by the significance of modular symmetry in generating neutrino
masses and flavor mixings, we apply the modular $A_4$ symmetry in a 3+1 scheme
of active-sterile neutrino mixings. Neutrino oscillation observables in the
3$\sigma$ range are successfully reproduced through the vacuum expectation
value of the modulus $\tau$ in the Inverted Hierarchy(IH) only, whereas Normal
Hierarchy (NH) is ruled out. We also study phenomenologies related to the
effective neutrino masses $m_{\beta}$ in tritium beta decay and
$m_{\beta\beta}$ in neutrinoless double beta decay. Mixings between active
neutrinos and eV scale sterile neutrino are analyzed in detail. The model also
predicts the Dirac CP-violating phase $\delta_{CP}$ and Majonara phases
$\alpha$ and $\beta$. The best-fit values of the neutrino mixing angles and
ratios of two mass-squared differences are determined using minimum $\chi^2$
analysis. The best-fit values of the neutrino oscillation observables are
predicted as $\sin^2\theta_{23}=0.566,$ $\sin^2\theta_{12}=0.307,$ $
\sin^2\theta_{13}=0.023$ and $r = 0.172$. The Dirac and Majorana phases are
observed at $\delta_{CP}=350.19^o,\ \alpha=355.75^o$ and $\beta=333.10^o$. We
also observe that the predictions of effective neutrino mass parameters in the
3+1 scheme are significantly different from the three neutrino paradigm.
