Magnetoelectric properties are studied by a combined experimental and theoretical study of a quasi-two-dimensional material composed of square cupolas, Ba(TiO)Cu4(PO4)4. The magnetization is measured up to above the saturation field, and several anomalies are observed depending on the field directions. We propose a S=1/2 spin model with Dzyaloshinskii-Moriya interactions, which well reproduces the full magnetization curves. Elaborating the phase diagram of the model, we show that the anomalies are explained by magnetoelectric phase transitions. Our theory also accounts for the scaling of the dielectric anomaly observed in experiments. The results elucidate the crucial role of the in-plane component of Dzyaloshinskii-Moriya interactions, which is induced by the noncoplanar buckling of square cupola. We also predict a 'hidden' phase and another magnetoelectric response both of which appear in nonzero magnetic field.PACS numbers: 77.80. Fm,75.85.+t,75.30.Kz Spatial asymmetry is a source of interesting phenomena in a broad range of condensed matter physics. A well-known example is the molecular asymmetry of water H 2 O, which leads to an electric polarization in each molecule. The asymmetry is at play also in magnets: the loss of inversion symmetry activates the asymmetric interactions through the relativistic spin-orbit coupling, such as the Dzyaloshinskii-Moriya (DM) interaction [1,2]. The asymmetric interactions lead to intriguing magnetism, e.g., weak ferromagnetism in antiferromagnets and spin-spiral ordering in helimagnets. They have also attracted growing interest as an origin of the magnetoelectric (ME) effect, that is, cross correlations between dielectricity and magnetism [3,4].Recently, an interesting series of chiral antiferromagnets, A(TiO)Cu 4 (PO 4 ) 4 (A = Ba, Sr) with space group P 42 1 2, was newly synthesized [5]. The materials have a quasi-two-dimensional structure, composed of an alternating array of Cu 4 O 12 clusters, as shown in Fig. 1(a). Each Cu 4 O 12 cluster consists of four corner-sharing CuO 4 plaquettes, forming a noncoplanar buckled structure termed (irregular) square cupola. The asymmetric unit can carry ME-active magnetic multipoles [6] associated with Cu spins. Indeed, a divergent anomaly of the dielectric constant was observed at the Néel temperature (T N =9.5 K) in magnetic fields applied along the [100] and [110] directions for A = Ba [7]. Although the ME response was argued by the magnetic quadrupole associated with noncoplanar antiferromagnetic ordering, the microscopic understanding is not fully obtained. It is highly desired to clarify how the unique asymmetry arising from the square cupolas affects the magnetic and dielectric properties in this series of compounds. [010][001] In this Letter, combining experimental and theoretical studies, we clarify the microscopic mechanism of ME behavior in A(TiO)Cu 4 (PO 4 ) 4 . First, from the magnetization measurement for the compound with A = Ba up to above the saturation field, we find several anomalies depending on the field dire...