There is a growing interest in developing dataâdriven reducedâorder models for atmospheric and oceanic flows that are trained on data obtained either from highâresolution simulations or satellite observations. The dataâdriven models are nonâintrusive in nature and offer significant computational savings compared to largeâscale numerical models. These lowâdimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equationâfree surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an endâtoâend nonâintrusive reducedâorder modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long shortâterm memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for longâterm forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.